Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  992 1790 (N52W20 L=331)   456    107        0 ( 936/1767-1042/1815)   4.8 W13-W25 N50-N52
      2 1307 1677 (N35W37 L=347)   360     53        1 (1279/1650-1335/1697)   9.5 W33-W39 N33-N36
      3 1553 1620 (N26W53 L=  3)  1711    221        3 (1455/1564-1636/1654) -19.7 W43-W60 N21-N29
      4 1732 1507 (N16W63 L= 13)   452     60        0 (1696/1494-1763/1519)  -8.1 W58-W67 N14-N16
      5 1433 1462 (N19W35 L=346)   268     62        0 (1397/1443-1468/1477) -21.8 W32-W37 N17-N20
      6  541 1277 (N25E26 L=285)  2135    555        2 ( 320/1212- 806/1335)  -9.0 E45-E07 N22-N29
      7 1842 1263 (S01W63 L= 13)   634     89        0 (1801/1232-1870/1303)   0.0 W57-W68 S01-N01
      8 1641 1102 (S06W40 L=351)   237     32        1 (1615/1081-1660/1120) -26.6 W37-W41 S06-S04
      9 1880  903 (S22W61 L= 12)   661     71        1 (1852/ 871-1902/ 942) -78.7 W57-W64 S24-S19
     10 1282  780 (S17W11 L=322)   361    178        0 (1246/ 712-1300/ 868)  58.0 W07-W12 S21-S12
     11 1109  784 (S14W00 L=311)   196     80        0 (1092/ 740-1130/ 819)  53.1 E01-W02 S16-S12
     12  456  427 (S24E55 L=256)   584     96        1 ( 441/ 376- 472/ 484)  35.0 E61-E51 S26-S19
  

    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