Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1577 1558 (N44W22 L=138)  1069    127        2 (1514/1524-1639/1581)  13.0 W15-W28 N40-N47
      2 1466 1438 (N31W15 L=131)   368     83        2 (1437/1403-1498/1471) -31.0 W12-W17 N29-N32
      3  312 1411 (S02E59 L= 57)   789     92        2 ( 290/1368- 336/1461) -45.0 E63-E56 S04-N01
      4 1480 1209 (N17W21 L=137)   304    135        0 (1443/1171-1537/1256) -14.0 W17-W25 N15-N18
      5 1750 1191 (N25W42 L=158)   325     48        0 (1728/1167-1775/1213)   0.0 W39-W45 N24-N26
      6  208 1051 (S25E60 L= 56)  1297    204        2 ( 123/ 988- 273/1117)  25.3 E73-E52 S29-S19
      7 1785 1045 (N18W48 L=165)   422     95        1 (1755/1010-1817/1090) -15.9 W44-W51 N16-N20
      8 1751  987 (N14W46 L=163)   945    243        1 (1696/ 903-1807/1089) -21.8 W39-W54 N11-N17
      9  333  983 (S26E45 L= 71)   459    134        2 ( 303/ 925- 353/1045) -50.2 E48-E42 S29-S23
     10  370  808 (S36E40 L= 77)   339     72        0 ( 340/ 781- 392/ 844)  59.0 E42-E37 S38-S33
     11 1613  722 (S05W45 L=161)  1066    286        1 (1497/ 641-1683/ 808)  65.2 W40-W47 S12-N01
     12  976  697 (S26W07 L=123)  1194    234        4 ( 886/ 631-1084/ 737)  56.3 W02-W12 S32-S20
     13  721  555 (S43E05 L=111)   369     68        2 ( 685/ 528- 751/ 578)  36.9 E07-E02 S44-S40
     14 1747  541 (S09W68 L=184)   485     78        0 (1713/ 518-1781/ 566)  51.3 W64-W70 S11-S06
     15  901  535 (S38W09 L=125)   268     63        0 ( 868/ 511- 931/ 557)   0.0 W05-W11 S38-S36
     16  731  469 (S48E01 L=115)   944    120        2 ( 677/ 443- 795/ 517)  29.1 E05-W04 S51-S45
     17 1603  473 (S18W58 L=175)   390     59        1 (1578/ 451-1628/ 496)  63.4 W55-W60 S19-S15
  

    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