Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1109 1909 (N60W49 L=229)  1109    108        1 (1075/1867-1150/1948)   9.5 W37-W67 N56-N62
      2 1423 1759 (N38W54 L=235)   384     53        0 (1391/1748-1456/1770) -15.9 W50-W58 N36-N38
      3  368 1170 (N23E41 L=140)   259     69        0 ( 345/1136- 389/1205) -59.0 E42-E39 N20-N25
      4 1898 1141 (S10W67 L=247)   683     86        0 (1879/1096-1916/1178)  26.6 W62-W71 S11-S07
      5 1778 1112 (S08W52 L=232)   335     64        0 (1765/1079-1795/1154) -90.0 W50-W53 S10-S05
      6  391 1096 (N18E39 L=141)   264     62        1 ( 365/1079- 414/1117)  26.6 E41-E37 N17-N19
      7  824  739 (S10E17 L=163)   323    114        0 ( 808/ 676- 843/ 788)  81.9 E18-E15 S14-S07
      8  375  641 (S08E51 L=129)  1028    165        0 ( 343/ 559- 409/ 725)  79.7 E53-E49 S13-S02
      9 1233  681 (S22W06 L=186)   370    141        0 (1172/ 643-1290/ 722)  14.0 W01-W09 S23-S20
     10 1712  649 (S36W43 L=224)   324     40        1 (1687/ 631-1731/ 669)  26.6 W40-W45 S37-S34
     11  958  398 (S35E20 L=161)   918    252        1 ( 861/ 371-1077/ 434) -31.6 E25-E11 S39-S30
  

    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