Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  526 1340 (N27E34 L=156)   276     63        1 ( 492/1325- 556/1352) -11.3 E36-E31 N26-N28
      2  323 1085 (N12E49 L=141)   747    167        2 ( 305/1009- 356/1168) -73.3 E50-E45 N06-N16
      3 1614  938 (S05W38 L=228)   838    188        2 (1538/ 897-1693/ 971) -13.0 W32-W45 S07-S02
      4 1835  881 (S12W60 L=250)   710     88        1 (1805/ 854-1873/ 918) -24.0 W56-W65 S13-S08
      5  452  746 (S10E41 L=149)   345     68        1 ( 428/ 721- 481/ 774) -53.1 E42-E38 S11-S07
      6 1628  630 (S26W42 L=231)   584    182        0 (1599/ 545-1663/ 716)  79.7 W40-W43 S31-S19
      7 1034  659 (S18E02 L=188)   590    158        1 ( 980/ 608-1100/ 700)  29.7 E05-W02 S21-S15
      8 1474  607 (S25W29 L=218)   656    212        2 (1396/ 529-1535/ 669)  45.0 W23-W33 S30-S21
      9  421  616 (S18E46 L=143)   446    119        1 ( 380/ 591- 470/ 656) -15.9 E49-E42 S19-S15
     10  900  528 (S26E12 L=178)   527     76        1 ( 858/ 503- 942/ 559)  26.6 E15-E08 S27-S24
     11  827  393 (S36E20 L=169)   537    120        0 ( 800/ 331- 858/ 437)  48.4 E24-E16 S41-S32
  

    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