Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1396 1616 (N32W43 L= 94)   272     76        0 (1357/1593-1434/1632) -45.0 W40-W44 N29-N33
      2  530 1478 (N46E27 L= 24)   955    248        1 ( 462/1427- 605/1541) -38.7 E34-E19 N40-N52
      3 1610 1475 (N17W51 L=102)   562    126        1 (1581/1423-1630/1531)  60.3 W47-W53 N14-N21
      4 1647  953 (S12W37 L= 88)   268     96        0 (1621/ 919-1663/1002)  45.0 W34-W39 S14-S09
      5  194  965 (N18E61 L=350)   400    160        0 ( 157/ 936- 235/ 992)  10.3 E66-E55 N16-N19
      6 1354  805 (S13W14 L= 65)   724    162        6 (1309/ 748-1388/ 868)  35.5 W10-W17 S16-S09
      7  463  641 (S04E45 L=  6)   426     95        1 ( 446/ 597- 486/ 695)  81.9 E46-E43 S07-N00
      8 1355  640 (S23W11 L= 62)   445    161        0 (1285/ 613-1422/ 679) -41.2 W07-W15 S26-S19
      9  631  625 (S08E33 L= 18)   999    215        0 ( 532/ 580- 729/ 651) -26.6 E40-E27 S12-S05
     10  821  458 (S22E27 L= 24)   478    144        1 ( 804/ 398- 843/ 522)  63.4 E29-E24 S26-S18
     11  584  384 (S23E47 L=  4)  1723    375        4 ( 537/ 282- 619/ 506)  46.8 E57-E41 S30-S14
     12  372  436 (S17E64 L=347)   793    230        2 ( 341/ 402- 413/ 462)  22.6 E70-E58 S19-S14
  

    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