Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1193 1890 (N73W27 L=145)   697     68        0 (1174/1860-1217/1921)  16.6 W18-W55 N68-N79
      2 1694 1493 (N34W57 L=176)  2045    336        1 (1575/1387-1809/1572) -33.0 W46-W68 N26-N40
      3 1391 1418 (N29W25 L=144)   386    112        0 (1363/1364-1412/1462) -81.9 W23-W26 N25-N32
      4  820 1385 (N25E15 L=104)   431    147        0 ( 763/1347- 876/1419)  35.5 E18-E11 N22-N27
      5  379 1393 (N23E51 L= 68)   322     59        1 ( 356/1362- 405/1415)  53.1 E53-E48 N20-N24
      6 1635 1367 (N26W45 L=164)   423     98        1 (1598/1340-1669/1403)   8.1 W41-W48 N24-N28
      7 1659 1057 (N06W43 L=162)   530    175        3 (1634/ 983-1687/1139)  84.3 W41-W45 N01-N11
      8  208  851 (S12E63 L= 56)  2321    231        6 ( 176/ 754- 242/ 954)  77.0 E68-E58 S18-S05
      9 1750  858 (S06W52 L=171)   518     32        0 (1724/ 836-1776/ 877)  11.3 W49-W55 S07-S04
     10  337  757 (S17E49 L= 70)   680     86        1 ( 310/ 720- 377/ 788) -90.0 E51-E45 S19-S15
     11  487  688 (S21E37 L= 82)  1225    223        2 ( 412/ 617- 556/ 769)  35.5 E44-E30 S26-S15
     12  741  683 (S20E18 L=101)   469    110        0 ( 706/ 649- 782/ 734) -51.3 E19-E14 S21-S16
     13 1653  623 (S21W48 L=167)   333     62        0 (1627/ 596-1674/ 652)  90.0 W46-W50 S23-S19
     14  299  575 (S30E61 L= 58)  1105    128        1 ( 247/ 542- 337/ 627)  15.9 E73-E52 S32-S25
     15 1038  554 (S27W03 L=121)   364    101        0 (1012/ 509-1076/ 597)  63.4 W00-W05 S30-S24
     16 1745  569 (S25W62 L=181)   522     76        0 (1707/ 552-1782/ 583)   5.7 W57-W68 S26-S24
     17 1377  479 (S32W29 L=148)   266     82        0 (1344/ 469-1424/ 497)  -8.1 W25-W32 S32-S30
     18 1121  467 (S33W09 L=128)   548    104        1 (1065/ 437-1161/ 492) -15.9 W05-W12 S35-S31
  

    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