Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  308 1536 (N35E69 L=193)   829     78        2 ( 294/1492- 322/1586) -32.5 E76-E62 N31-N38
      2 1101 1445 (N30W06 L=267)   456    190        0 (1016/1414-1186/1483) -20.0 W00-W11 N28-N32
      3 1308 1301 (N21W19 L=280)   291    118        1 (1292/1248-1334/1349) -71.6 W18-W20 N17-N23
      4  198 1203 (N13E65 L=197)  1592    168        3 ( 155/1141- 244/1291)  68.2 E70-E58 N08-N18
      5 1886 1018 (N01W67 L=329)  3350    276       13 (1820/ 881-1935/1133) -65.0 W58-W77 S07-N08
      6  735  999 (N02E18 L=243)   216     70        1 ( 703/ 970- 760/1025) -36.9 E20-E16 N00-N03
      7  225  926 (S04E59 L=203)   689    115        0 ( 204/ 872- 254/ 974)  49.4 E62-E55 S07-N00
      8  351  798 (S11E47 L=214)   758    169        2 ( 302/ 736- 383/ 872)  48.8 E53-E44 S15-S06
      9  437  625 (S23E43 L=219)  1570    227        4 ( 383/ 543- 501/ 695)  47.7 E50-E36 S28-S17
     10 1716  565 (S28W56 L=318)   783    148        0 (1687/ 500-1736/ 637) -68.2 W52-W58 S32-S22
     11 1457  518 (S30W32 L=293)  2337    457        5 (1223/ 450-1607/ 604)  13.6 W14-W43 S34-S24
     12  424  400 (S40E57 L=204)  1032     68        2 ( 396/ 371- 456/ 429)  18.4 E65-E50 S42-S37
  

    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