Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  824 1695 (N30E35 L=235)  1029    144        3 ( 779/1634- 873/1752) -32.5 E40-E29 N26-N33
      2  485 1628 (N17E56 L=214)   843    129        3 ( 424/1599- 536/1665)  32.0 E61-E50 N14-N19
      3  735 1493 (N15E31 L=240)   272     75        0 ( 710/1457- 760/1528)  90.0 E31-E29 N11-N17
      4 1687 1257 (N28W37 L=308)   349    125        0 (1643/1211-1732/1295) -12.5 W32-W41 N26-N29
      5 1853 1162 (N28W57 L=327)   787    100        2 (1823/1120-1879/1208) -38.7 W52-W60 N25-N30
      6  159 1147 (S18E69 L=201)   634     77        0 ( 123/1125- 194/1163)  -4.1 E77-E63 S18-S16
      7  426  946 (S24E36 L=234)   272     88        0 ( 405/ 912- 441/ 992)  78.7 E37-E34 S26-S21
      8  273  969 (S26E51 L=220)   307     47        0 ( 249/ 948- 302/ 987)  31.0 E53-E48 S28-S24
      9  581  747 (S32E21 L=250)   494    165        0 ( 517/ 720- 660/ 782)  35.0 E25-E15 S35-S27
     10  979  512 (S35W14 L=284)   234     58        1 ( 949/ 498-1011/ 530)  36.9 W11-W15 S36-S33
     11 1662  418 (S19W70 L=341)  2220    245        3 (1615/ 319-1749/ 484) -67.4 W63-W82 S24-S12
     12 1430  352 (S30W54 L=325)   528     95        0 (1384/ 334-1482/ 371)  51.3 W50-W57 S32-S27
     13 1314  247 (S41W57 L=327)   472     72        0 (1272/ 237-1354/ 260)  20.6 W52-W60 S42-S38
  

    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