Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1108 1320 (N18E01 L=282)  3196    753        4 ( 842/1025-1334/1589)  17.0 E24-W18 N04-N30
      2 1705 1567 (N50W49 L=331)   573    134        0 (1669/1555-1735/1584)  12.5 W43-W52 N48-N52
      3 1439 1555 (N41W19 L=302)   340     83        0 (1403/1540-1476/1570)  26.6 W16-W22 N39-N42
      4 1873 1080 (N19W61 L=344)  1533    267        2 (1828/1011-1923/1152) -28.3 W55-W69 N15-N22
      5  145 1006 (S18E66 L=216)   791    160        0 ( 125/ 950- 163/1065)  74.1 E68-E64 S21-S14
      6 1372  939 (N00W22 L=305)   397    133        0 (1328/ 896-1426/ 992) -23.2 W18-W26 S01-N02
      7  408  830 (S24E38 L=245)   620    157        1 ( 343/ 797- 475/ 863)   5.2 E43-E32 S25-S22
      8  941  791 (S17E01 L=282)  1782    480        3 ( 723/ 745-1166/ 851)  22.9 E14-W12 S21-S09
      9 1389  373 (S34W43 L=325)  5713    748       10 (1083/ 325-1640/ 453)  26.6 W21-W60 S43-S23
     10 1022  367 (S43W17 L=299)   362     54        0 ( 988/ 351-1056/ 387)  26.6 W14-W19 S45-S41
  

    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