Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1261 1840 (N47W54 L=254)   442     31        1 (1243/1817-1282/1861)  12.5 W49-W58 N45-N48
      2 1048 1571 (N35W17 L=217)  1923    620        0 ( 787/1506-1332/1618) -16.6 E02-W36 N27-N39
      3 1604 1440 (N12W48 L=249)  1041    174        1 (1547/1371-1659/1513) -78.7 W46-W50 N07-N17
      4  162 1194 (N32E66 L=135)   681     89        0 ( 122/1172- 210/1219) -12.5 E75-E57 N29-N34
      5 1241 1012 (S02W12 L=212)   448    126        1 (1196/ 976-1304/1033) -15.9 W08-W15 S04-N00
      6 1547  815 (S21W27 L=228)   662    146        1 (1503/ 761-1572/ 883)  45.0 W23-W29 S23-S17
      7 1787  612 (S43W52 L=252)  1625    252        1 (1687/ 544-1862/ 707)  12.8 W39-W62 S45-S37
      8  834  499 (S23E25 L=175)   457    170        0 ( 807/ 418- 875/ 558)  38.7 E30-E20 S27-S19
      9  654  445 (S22E39 L=161)  1214    248        2 ( 618/ 338- 699/ 550)  60.3 E43-E33 S29-S15
     10  400  512 (S13E56 L=144)   329     57        0 ( 384/ 486- 421/ 538)  18.4 E58-E52 S14-S11
     11 1371  343 (S50W06 L=207)  1692    311        0 (1237/ 285-1509/ 408)  -4.4 E06-W20 S51-S47
  

    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