Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1443 1761 (N46W53 L=127)  2298    285        3 (1392/1675-1506/1837)  65.2 W49-W59 N38-N51
      2 1808 1328 (N11W64 L=138)  1678    296        1 (1740/1218-1887/1439) -58.0 W60-W71 N02-N18
      3 1441 1371 (N21W33 L=107)   309    106        3 (1400/1333-1484/1406) -63.4 W30-W35 N17-N23
      4 1731 1348 (N14W56 L=130)   439     91        1 (1707/1324-1763/1373)   0.0 W52-W59 N12-N15
      5  420 1336 (N32E40 L= 34)   328    130        0 ( 379/1319- 471/1360) -18.4 E44-E35 N30-N33
      6  240 1284 (N30E57 L= 17)  2295    496        5 ( 135/1253- 336/1323)  -2.6 E69-E47 N27-N32
      7 1459 1254 (N13W31 L=105)   263    100        1 (1411/1219-1494/1287)  15.9 W27-W34 N11-N14
      8  530 1067 (N14E31 L= 43)   334     86        1 ( 490/1038- 564/1102) -38.7 E34-E28 N12-N16
      9 1711  899 (S12W45 L=119)   375     25        0 (1693/ 877-1732/ 926) -71.6 W43-W46 S13-S10
     10 1776  669 (S29W52 L=126)  1531    233        5 (1703/ 599-1835/ 744)  26.6 W44-W57 S33-S24
     11 1620  568 (S33W37 L=112)   807    188        4 (1561/ 534-1686/ 598)   5.2 W32-W43 S34-S30
     12  921  499 (S26E14 L= 60)  1008    231        1 ( 819/ 466-1024/ 543) -17.1 E20-E07 S27-S21
  

    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