Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  998 1657 (N35E09 L=177)   505    129        0 ( 928/1632-1063/1677)  -5.2 E14-E03 N33-N36
      2  533 1321 (N09E34 L=152)   352     86        1 ( 502/1284- 561/1351) -45.0 E36-E31 N06-N10
      3 1598 1226 (N13W35 L=221)   245     48        1 (1580/1206-1617/1256) -56.3 W33-W36 N11-N14
      4 1855  905 (S02W61 L=247)   414     62        0 (1837/ 870-1876/ 935)  90.0 W58-W63 S04-N00
      5  437  480 (S44E44 L=142)   590    160        0 ( 395/ 452- 480/ 526)  21.8 E48-E38 S46-S40
      6  996  369 (S48W06 L=192)   508     61        1 ( 967/ 331-1018/ 402) -63.4 W03-W08 S51-S45
      7 1172  309 (S50W24 L=210)   573    115        1 (1140/ 277-1210/ 348) -29.1 W19-W29 S52-S47
  

    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