Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1704 1550 (N22W68 L=272)  1503    237        2 (1588/1518-1799/1582) -21.0 W56-W82 N16-N27
      2 1498 1437 (N21W41 L=245)   279     52        1 (1474/1414-1518/1459) -90.0 W38-W42 N19-N22
      3  856 1264 (N24E06 L=198)  1117    263        4 ( 777/1171- 922/1382) -68.2 E10-E03 N17-N32
      4 1631 1292 (N10W46 L=250)   326    107        0 (1590/1261-1680/1310)   0.0 W41-W50 N08-N10
      5 1870 1003 (S13W62 L=266)   479     68        0 (1851/ 961-1886/1034)  45.0 W59-W65 S15-S11
      6 1354  805 (S12W16 L=220)   471    128        1 (1329/ 753-1390/ 854)  51.3 W14-W19 S15-S09
      7 1007  795 (S06E05 L=199)   363    146        0 ( 978/ 739-1065/ 835)  33.7 E07-E01 S09-S04
      8  779  736 (S06E20 L=184)   661    195        3 ( 698/ 689- 866/ 781) -37.9 E24-E15 S10-S02
      9 1625  729 (S24W35 L=239)   318     98        1 (1581/ 703-1653/ 767)  18.4 W31-W37 S25-S21
     10 1760  713 (S29W49 L=253)   493     89        0 (1726/ 678-1796/ 748)  23.2 W45-W52 S30-S26
     11 1014  477 (S27E12 L=192)   636    293        0 ( 938/ 379-1124/ 555)  27.9 E19-E02 S33-S23
     12 1210  455 (S33W01 L=205)   389     99        0 (1174/ 421-1253/ 499)  38.7 E01-W04 S35-S30
  

    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