Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  933 1487 (N36W08 L=149)  1682    496        3 ( 755/1354-1144/1553)  -1.9 E07-W23 N32-N38
      2 1856 1258 (S05W65 L=206)   770    117        1 (1832/1189-1880/1310)  82.9 W63-W67 S09-S01
      3 1760  891 (S22W45 L=185)  2011    322        2 (1670/ 826-1861/ 993)  12.5 W36-W57 S24-S16
      4 1799  698 (S36W49 L=190)   484     68        0 (1775/ 673-1828/ 735)  26.6 W45-W52 S37-S34
      5  949  681 (S10E13 L=128)   224     64        0 ( 917/ 665- 981/ 699) -45.0 E14-E11 S12-S09
      6 1471  623 (S29W17 L=158)   426    100        1 (1439/ 602-1518/ 652)  -9.5 W14-W20 S30-S26
      7 1232  626 (S21W02 L=142)   360    113        0 (1177/ 613-1300/ 644) -33.7 E01-W05 S23-S18
      8 1348  500 (S33W06 L=147)   889    170        2 (1305/ 425-1388/ 575)  41.6 W01-W10 S37-S29
      9  869  442 (S23E26 L=115)   687    199        1 ( 809/ 374- 915/ 531)  31.6 E33-E20 S26-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