Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  866 1276 (N23E03 L= 82)   532    138        2 ( 791/1246- 922/1311) -15.9 E07-W00 N21-N24
      2 1748  989 (S17W44 L=129)   474    166        0 (1702/ 930-1805/1038)  15.3 W39-W50 S18-S15
      3  291  963 (N18E49 L= 36)  1731    328        4 ( 174/ 931- 466/1023) -21.0 E61-E34 N13-N23
      4  698  901 (N06E21 L= 64)   528    180        1 ( 658/ 847- 731/ 969)  50.2 E24-E19 N03-N09
      5  533  831 (N06E33 L= 52)   281     94        0 ( 485/ 812- 573/ 847) -38.7 E36-E31 N03-N07
      6  180  770 (N09E67 L= 18)   723     69        1 ( 156/ 735- 204/ 805) -59.0 E69-E64 N06-N11
      7 1709  713 (S32W38 L=123)   466    102        0 (1673/ 675-1742/ 764)  29.7 W34-W41 S34-S29
      8  512  652 (S04E40 L= 45)   938    315        1 ( 388/ 579- 602/ 761)  12.8 E52-E30 S06-N00
      9 1565  485 (S44W23 L=108)   496     92        0 (1520/ 469-1612/ 505) -26.6 W18-W28 S46-S40
     10 1231  334 (S43E08 L= 77)   908    297        2 (1115/ 265-1342/ 416)  10.3 E19-W03 S46-S40
     11  924  201 (S44E38 L= 47)   545    109        0 ( 879/ 184- 978/ 220)  -5.2 E43-E32 S45-S42
  

    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