Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  432 1638 (N14E60 L=281)   517     67        0 ( 417/1603- 448/1673) -53.1 E62-E57 N12-N16
      2 1122 1630 (N32E13 L=328)   244     51        2 (1088/1615-1150/1647)  14.0 E14-E10 N30-N33
      3 1596 1573 (N46W23 L=  4)   796    144        0 (1568/1525-1617/1620)  73.3 W21-W25 N40-N50
      4  531 1445 (N05E42 L=299)   328    114        0 ( 485/1411- 566/1466)   0.0 E44-E38 N03-N05
      5 1770 1423 (N42W42 L= 23)  1230    172        1 (1717/1384-1820/1451)  10.3 W36-W48 N37-N44
      6 1851 1259 (N34W51 L= 31)  2517    342        5 (1780/1133-1881/1358) -58.4 W42-W54 N25-N39
      7  203 1239 (S13E63 L=278)  2128    404        5 ( 107/1202- 292/1265)   0.0 E78-E52 S14-S11
      8 1703 1072 (N16W40 L= 21)   326     95        0 (1672/1022-1724/1110) -38.7 W36-W42 N13-N17
      9  647  970 (S19E21 L=320)   249    119        0 ( 619/ 911- 669/1035) -45.0 E23-E17 S22-S16
     10 1622  926 (N05W38 L= 19)   658    153        0 (1564/ 875-1676/ 975)  -5.2 W32-W43 N03-N07
     11 1175  854 (S12W13 L=354)   258    148        0 (1158/ 787-1194/ 928) -60.3 W11-W16 S14-S07
     12 1836  864 (N10W58 L= 39)  1582    261        3 (1764/ 829-1903/ 896)  17.4 W51-W67 N05-N13
     13 1739  611 (S06W61 L= 41)  2156    329        3 (1588/ 486-1864/ 768)  61.4 W54-W67 S17-N05
     14 1447  685 (S13W35 L= 16)   541    150        0 (1394/ 651-1503/ 752)  82.9 W31-W36 S15-S07
     15  836  653 (S33E00 L=341)   561    202        1 ( 747/ 627- 940/ 701)  48.0 E05-W04 S36-S26
     16 1230  565 (S26W27 L=  8)   860    194        1 (1140/ 524-1304/ 633)  -4.1 W18-W32 S27-S23
     17 1405  417 (S28W45 L= 26)   396    120        1 (1370/ 389-1446/ 450)   0.0 W41-W49 S28-S27
  

    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