Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  742 1533 (N43E13 L=156)   407    125        1 ( 712/1496- 788/1577) -49.4 E15-E09 N39-N46
      2 1661 1437 (N20W53 L=222)  3155    504       10 (1505/1367-1807/1505) -32.5 W41-W66 N12-N26
      3 1231 1291 (N19W17 L=187)   418    144        2 (1158/1263-1291/1310) -14.0 W12-W20 N16-N20
      4  269 1264 (N29E54 L=116)   365     73        0 ( 244/1250- 300/1281)  11.3 E56-E50 N28-N30
      5  609 1100 (N16E25 L=144)   357     69        2 ( 586/1063- 633/1131) -90.0 E27-E23 N13-N18
      6 1788 1003 (S08W53 L=222)  1697    438        8 (1736/ 848-1851/1106)  51.8 W46-W59 S16-S02
      7 1905  952 (S15W66 L=236)  1427    262        2 (1873/ 883-1931/1035)  45.0 W61-W72 S18-S10
      8 1099  865 (S04W02 L=172)   594    131        2 (1025/ 832-1159/ 893)  -8.1 E01-W06 S06-S02
      9 1143  621 (S20W01 L=171)  1836    527        7 (1017/ 432-1236/ 844)  53.1 E10-W09 S32-S07
     10 1563  732 (S21W31 L=201)   582    143        2 (1522/ 706-1603/ 762)   0.0 W28-W34 S22-S18
     11  786  505 (S23E24 L=145)   409    231        0 ( 768/ 423- 808/ 612)  74.1 E25-E21 S29-S15
     12  391  465 (S22E57 L=112)  1009    290        5 ( 349/ 418- 430/ 516)  26.6 E64-E50 S25-S18
  

    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