Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1405 1545 (N39W17 L=149)   290     41        0 (1382/1525-1431/1560)  45.0 W14-W18 N37-N40
      2 1075 1434 (N24E05 L=127)   597    170        0 (1016/1363-1142/1495) -31.0 E10-W00 N21-N27
      3 1730 1378 (N35W46 L=178)  2889    440        4 (1601/1295-1872/1467)  -4.2 W33-W61 N30-N39
      4 1581 1195 (N20W33 L=165)   403     48        1 (1549/1177-1616/1211)  11.3 W30-W35 N18-N20
      5 1537 1073 (N11W31 L=163)   739    271        1 (1502/ 997-1598/1185) -55.0 W28-W36 N07-N18
      6 1332  918 (S02W20 L=152)   484    119        1 (1269/ 885-1377/ 963)  45.0 W17-W22 S04-N01
      7 1320  642 (S18W26 L=158)   312     36        0 (1296/ 622-1349/ 665) -14.0 W24-W28 S19-S16
      8 1729  626 (S11W59 L=191)   691    172        1 (1686/ 601-1777/ 656)   0.0 W53-W64 S12-S09
  

    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