Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1652 1554 (N23W56 L=221)   925    163        0 (1589/1495-1709/1605) -53.1 W52-W59 N18-N26
      2  629 1510 (N36E21 L=145)   829    233        1 ( 545/1451- 745/1546)   3.6 E27-E11 N33-N38
      3 1274 1106 (N01W16 L=181)   198     76        0 (1252/1073-1302/1145)  36.9 W14-W18 N00-N03
      4 1569  957 (S11W33 L=198)   327    110        0 (1545/ 900-1591/1007)  56.3 W30-W35 S14-S08
      5 1211  731 (S20W07 L=172)   506    168        0 (1144/ 677-1261/ 798)  33.7 W02-W11 S22-S16
      6 1495  750 (S23W26 L=191)   234     53        0 (1483/ 722-1509/ 789) -76.0 W25-W27 S24-S20
      7  291  628 (S13E59 L=106)  1475    207        5 ( 254/ 558- 323/ 709)  63.4 E64-E56 S17-S07
      8 1011  660 (S21E06 L=159)   267    108        0 ( 963/ 642-1070/ 683) -26.6 E08-E02 S22-S19
      9  352  609 (S15E54 L=111)   551    105        0 ( 335/ 561- 374/ 665)  71.6 E56-E51 S17-S11
     10  691  597 (S20E28 L=137)   287     76        0 ( 667/ 562- 718/ 633)  38.7 E30-E25 S22-S18
     11  401  383 (S30E65 L=100)  1419    181        1 ( 346/ 322- 460/ 453) -84.3 E68-E61 S34-S24
     12  498  363 (S33E56 L=109)   573    100        0 ( 476/ 321- 518/ 416)  28.6 E62-E51 S35-S29
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y