Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  795 1679 (N26E34 L=287)   668     92        1 ( 768/1639- 837/1729) -20.6 E38-E30 N23-N29
      2  531 1675 (N20E53 L=268)  1087    199        2 ( 428/1628- 594/1724)  45.0 E59-E49 N15-N24
      3 1213 1375 (N17W02 L=322)   275     69        1 (1194/1333-1234/1404)  90.0 W00-W03 N14-N19
      4 1807 1237 (N28W48 L=  9)  1301    250        1 (1765/1122-1842/1347) -59.7 W43-W53 N22-N34
      5 1097  878 (S13W08 L=329)   208     55        0 (1077/ 853-1127/ 897)  63.4 W06-W09 S15-S11
      6  892  645 (S32W03 L=323)   268     71        0 ( 854/ 621- 925/ 667)  59.0 W00-W04 S34-S29
      7 1120  539 (S32W21 L=341)   498    120        0 (1077/ 507-1170/ 595) -12.5 W15-W24 S33-S29
      8  596  558 (S46E17 L=303)   275     63        0 ( 567/ 542- 633/ 579)   0.0 E20-E13 S47-S45
      9  631  355 (S61E10 L=310)   397     42        1 ( 605/ 333- 656/ 379)   8.1 E13-E06 S63-S58
     10  793  290 (S61W11 L=332)   638    114        1 ( 755/ 250- 839/ 334) -26.6 W06-W19 S62-S56
  

    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