Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1114 1888 (N57E35 L=  1)   866     58        1 (1086/1861-1142/1914) -10.3 E41-E28 N53-N60
      2 1480 1715 (N56W13 L= 49)   679     90        1 (1453/1668-1515/1738)  73.3 W09-W17 N51-N61
      3 1689 1576 (N53W40 L= 76)   704     73        1 (1665/1538-1712/1613)  81.9 W36-W44 N49-N56
      4 1051 1429 (N19E10 L= 26)   312    118        1 (1000/1402-1102/1452)   7.1 E13-E05 N18-N20
      5 1622 1256 (N26W31 L= 67)   429    135        0 (1567/1220-1649/1306)  54.5 W27-W33 N21-N29
      6 1633 1105 (N17W35 L= 71)   588    159        0 (1609/1019-1666/1170)  90.0 W33-W36 N11-N21
      7 1767 1043 (N18W48 L= 84)   558     50        2 (1738/1018-1793/1066) -11.3 W45-W50 N16-N19
      8 1892  967 (N18W65 L=101)  1618    117        7 (1856/ 921-1935/1019)  66.8 W59-W74 N14-N21
      9  244  928 (S29E54 L=342)  1050    208        2 ( 209/ 811- 280/1005) -61.4 E58-E49 S35-S24
     10 1713  960 (N11W45 L= 81)   481     84        1 (1665/ 926-1752/ 978)  14.0 W40-W48 N09-N13
     11 1826  898 (N12W58 L= 94)   924    125        1 (1786/ 850-1852/ 969)  69.4 W54-W61 N08-N16
     12  562  893 (S24E25 L= 11)   238     40        2 ( 536/ 875- 589/ 910)  33.7 E26-E22 S24-S22
     13  351  858 (S31E42 L=354)   392    111        0 ( 304/ 829- 402/ 888)  29.1 E46-E37 S33-S28
     14  188  834 (S36E63 L=333)   728     72        2 ( 166/ 797- 221/ 869)  76.0 E67-E57 S38-S34
     15 1197  411 (S35W33 L= 69)   640     80        1 (1154/ 392-1247/ 440)  38.7 W29-W36 S36-S32
     16  985  416 (S41W19 L= 54)   353     89        0 ( 932/ 401-1032/ 430)  23.2 W14-W21 S42-S39
     17  502  394 (S62E27 L=  9)   450     52        0 ( 471/ 368- 531/ 414)  -6.3 E31-E22 S63-S59
  

    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