Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1174 1886 (N72W22 L=142)   872     98        1 (1140/1853-1212/1924)  15.9 W12-W61 N66-N80
      2 1694 1480 (N34W56 L=177)  2162    390        3 (1550/1352-1822/1570) -33.1 W43-W68 N24-N40
      3 1371 1420 (N29W24 L=144)   339     97        1 (1346/1364-1392/1463) -74.1 W22-W24 N25-N32
      4  356 1406 (N24E54 L= 66)   456    108        3 ( 328/1363- 380/1456)  90.0 E56-E50 N20-N27
      5  798 1387 (N25E17 L=103)   362    109        1 ( 744/1359- 845/1412)  29.7 E20-E13 N22-N26
      6 1620 1364 (N26W43 L=164)   446     78        2 (1583/1336-1654/1391)   8.1 W39-W46 N23-N27
      7 1637 1058 (N06W41 L=162)   769    180        3 (1594/ 978-1666/1149)  73.3 W37-W44 N01-N11
      8  212  833 (S13E62 L= 58)  2916    285        6 ( 165/ 729- 283/ 963) -70.3 E70-E54 S19-S05
      9  822  842 (S09E12 L=108)   388    251        1 ( 728/ 762- 885/ 919)  45.0 E18-E08 S14-S04
     10 1748  854 (S06W52 L=173)   627     58        3 (1722/ 822-1776/ 882)  76.0 W49-W55 S08-S04
     11  321  757 (S17E51 L= 69)   718     94        4 ( 291/ 715- 359/ 796)  80.5 E54-E47 S20-S14
     12  388  636 (S25E47 L= 73)  3126    393        7 ( 242/ 532- 519/ 747)  21.6 E77-E33 S34-S16
     13  750  673 (S20E17 L=103)   635    139        3 ( 702/ 638- 811/ 725) -35.5 E20-E12 S22-S16
     14 1660  632 (S21W49 L=169)   560    121        1 (1614/ 596-1713/ 659)  26.6 W45-W54 S23-S18
     15 1018  558 (S27W01 L=121)   326    104        0 ( 990/ 514-1057/ 596)  63.4 W00-W03 S30-S24
     16 1103  468 (S34W08 L=128)   603    105        2 (1049/ 439-1148/ 497) -20.6 W03-W11 S35-S31
     17 1362  480 (S32W28 L=148)   260     95        0 (1322/ 469-1405/ 496)  -8.1 W24-W31 S32-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