Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1288 1733 (N39W43 L=203)  1330    284        1 (1172/1678-1407/1782)   0.0 W30-W57 N36-N40
      2 1444 1538 (N22W41 L=201)   330    102        0 (1406/1505-1500/1564) -51.3 W38-W43 N19-N24
      3  503 1341 (N33E28 L=132)   743    196        2 ( 412/1318- 601/1358) -10.6 E35-E19 N31-N35
      4 1121 1273 (N15W11 L=171)   305     90        0 (1072/1254-1164/1287) -31.0 W08-W13 N13-N16
      5 1637 1266 (N02W44 L=204)   249     25        0 (1613/1252-1658/1282) -14.0 W41-W45 N01-N03
      6 1812 1225 (S04W58 L=218)   454     63        0 (1785/1190-1836/1258) -53.1 W55-W59 S06-S02
      7  703 1237 (N22E15 L=145)   291     68        1 ( 667/1217- 738/1254)   0.0 E17-E12 N21-N23
      8 1745 1047 (S13W47 L=207)  1736    271        4 (1680/ 969-1806/1138) -61.7 W42-W53 S18-S05
      9  277 1006 (N18E51 L=110)  1711    307        2 ( 143/ 945- 373/1095)   4.2 E67-E40 N14-N21
     10 1845 1029 (S16W57 L=217)   532     63        2 (1826/1000-1867/1063)  31.0 W54-W60 S17-S14
     11 1445  953 (S10W23 L=183)   227     72        0 (1429/ 912-1461/ 985)  53.1 W21-W24 S12-S08
     12 1574  802 (S23W29 L=189)   334     56        0 (1554/ 773-1590/ 840)  71.6 W27-W30 S24-S20
     13 1399  728 (S22W16 L=176)   773    182        0 (1333/ 644-1452/ 804)  35.0 W10-W20 S26-S18
     14 1185  633 (S23W00 L=161)  1051    268        2 (1145/ 524-1247/ 755)  50.2 E04-W06 S29-S17
     15 1426  604 (S31W16 L=176)   507    130        0 (1369/ 575-1477/ 647)   6.3 W11-W20 S32-S29
     16  840  455 (S26E26 L=134)   536    145        1 ( 806/ 387- 862/ 526)  41.6 E31-E22 S30-S21
     17 1636  452 (S50W36 L=196)  1103    112        1 (1569/ 429-1695/ 484) -10.0 W27-W44 S52-S47
  

    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