Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  741 1612 (N20E35 L=135)   886    252        4 ( 679/1550- 822/1677) -14.9 E42-E27 N18-N22
      2  314 1567 (N08E69 L=101)   662     65        1 ( 289/1543- 336/1592)  90.0 E72-E66 N06-N10
      3  719 1571 (N17E35 L=135)   239     51        1 ( 697/1555- 743/1592) -26.6 E36-E32 N15-N18
      4 1835 1127 (N25W51 L=221)   529    106        2 (1807/1092-1863/1167) -15.9 W47-W54 N23-N26
      5 1039  808 (S17W07 L=177)   599    188        1 ( 955/ 761-1135/ 865)  55.0 W03-W10 S21-S11
      6 1721  678 (S03W56 L=225)  1573    294        4 (1667/ 557-1758/ 822) -65.0 W50-W60 S10-N05
      7 1239  625 (S21W25 L=195)   691    230        2 (1142/ 567-1345/ 720)  -3.2 W15-W33 S22-S18
      8  731  509 (S45E03 L=167)  1037    290        0 ( 603/ 476- 869/ 540)  24.2 E13-W07 S49-S40
      9 1519  325 (S28W60 L=230)  1655    344        2 (1419/ 295-1638/ 353)  27.9 W52-W70 S31-S22
  

    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