Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  925 1819 (N44E38 L=322)  1249    143        4 ( 863/1769- 979/1857)  -4.1 E45-E29 N41-N47
      2  684 1722 (N30E47 L=312)   372     57        1 ( 661/1696- 712/1755) -80.5 E49-E45 N27-N33
      3 1075 1657 (N36E16 L=343)   440     88        2 (1048/1620-1114/1696) -20.6 E19-E11 N34-N38
      4  381 1582 (N14E64 L=296)  1400    267        1 ( 326/1485- 450/1682) -77.9 E67-E61 N07-N21
      5  501 1592 (N17E52 L=307)   456    103        1 ( 461/1575- 550/1618)   0.0 E56-E47 N16-N18
      6 1476 1135 (N14W24 L= 23)   256    104        0 (1431/1098-1510/1160)  -9.5 W20-W26 N12-N15
      7  603  894 (S21E22 L=337)   414    206        0 ( 562/ 819- 651/ 976) -41.2 E26-E17 S24-S16
      8 1053  623 (S26W13 L= 12)   493    141        3 ( 992/ 598-1129/ 649)  40.6 W09-W16 S28-S22
      9 1362  254 (S39W59 L= 58)   776    112        1 (1341/ 217-1385/ 306) -14.0 W50-W66 S41-S36
     10 1437  253 (S37W67 L= 66)   904     49        1 (1406/ 225-1461/ 278)  45.0 W62-W72 S38-S34
     11 1288  252 (S42W53 L= 52)   480     27        1 (1270/ 231-1311/ 276)   0.0 W48-W56 S43-S40
  

    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