Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1376 1773 (N36W57 L=147)  1438    139        4 (1307/1729-1442/1801) -15.9 W49-W65 N34-N39
      2  694 1726 (N60W01 L= 91)   325     84        0 ( 646/1712- 728/1741) -11.3 E05-W05 N58-N61
      3  251 1111 (N29E52 L= 38)  1626    226        3 ( 155/1048- 350/1146) -14.0 E66-E42 N24-N32
      4 1891 1065 (S18W62 L=153)  1106    125        1 (1857/1011-1936/1111)  11.3 W56-W71 S20-S16
      5  219  884 (N16E59 L= 31)   550     97        3 ( 177/ 852- 264/ 919) -21.8 E64-E54 N14-N18
      6 1749  754 (S32W42 L=132)  2254    296        4 (1706/ 638-1780/ 901)  76.0 W36-W45 S39-S22
      7  365  826 (N10E47 L= 44)   616    158        1 ( 339/ 762- 399/ 894)  66.8 E48-E42 N06-N13
      8  628  814 (N04E28 L= 62)   243     73        0 ( 587/ 800- 665/ 829) -36.9 E30-E25 N02-N05
      9  235  776 (N10E61 L= 30)  1033    116        2 ( 191/ 721- 261/ 818) -50.2 E65-E56 N06-N12
     10 1527  486 (S42W18 L=109)   442     65        1 (1493/ 464-1569/ 509)  -8.1 W15-W22 S43-S40
     11  319  473 (S10E69 L= 22)   802     54        1 ( 288/ 449- 343/ 502)  16.7 E75-E64 S11-S08
     12  697  218 (S35E56 L= 35)  1699    169        1 ( 657/ 141- 726/ 282)  29.5 E74-E48 S42-S29
  

    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