Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  443 1489 (N31E48 L=179)   340     54        1 ( 422/1465- 472/1515) -15.9 E51-E44 N29-N32
      2 1720 1371 (N25W54 L=281)   358     52        0 (1691/1353-1750/1385)   7.1 W50-W58 N24-N26
      3  443  703 (S19E40 L=187)   346    128        0 ( 395/ 685- 499/ 720)   0.0 E44-E35 S20-S17
      4  205  666 (S23E71 L=156)  1219     66        4 ( 177/ 627- 234/ 701)  21.8 E80-E64 S26-S20
      5 1710  649 (S20W53 L=280)   489     75        0 (1679/ 625-1732/ 679)  21.8 W49-W55 S21-S18
      6  452  608 (S25E42 L=185)   446    144        0 ( 434/ 546- 490/ 660)  53.1 E44-E37 S29-S21
      7 1164  367 (S41W13 L=240)   457     59        2 (1139/ 339-1192/ 394)  78.7 W10-W15 S43-S38
      8  912  323 (S46E08 L=219)   857    170        2 ( 827/ 281- 976/ 354) -17.1 E15-E02 S49-S43
  

    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