Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  870 1760 (N33E35 L=264)   539     77        1 ( 847/1725- 896/1798) -63.4 E37-E32 N30-N36
      2  376 1685 (N18E75 L=225)  1010     90        1 ( 343/1654- 406/1720) -59.0 E79-E70 N15-N20
      3 1232 1596 (N32E04 L=295)  1257    350        0 (1075/1553-1335/1673)  54.8 E12-W00 N25-N42
      4 1210 1381 (N18W01 L=300)   338    114        0 (1176/1339-1247/1433)  82.9 W00-W01 N14-N22
      5 1566 1251 (N21W26 L=325)   918    240        0 (1488/1201-1647/1288)  22.6 W21-W33 N15-N23
      6 1897  938 (N15W64 L=  4)  1909    214        0 (1874/ 846-1937/1008) -84.8 W61-W71 N09-N20
      7 1084  911 (S12W06 L=306)   346    147        0 (1032/ 867-1154/ 945)  -6.3 W02-W11 S12-S10
      8 1010  721 (S24W08 L=307)   624    307        0 ( 913/ 669-1117/ 811)  58.4 W03-W11 S29-S16
      9 1223  683 (S20W21 L=321)   358    105        1 (1172/ 650-1270/ 709)   8.1 W17-W24 S21-S18
     10 1561  622 (S13W45 L=345)   473    133        0 (1503/ 585-1619/ 656)  54.5 W42-W48 S16-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