Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  494 1338 (N33E28 L=178)  1170    169        4 ( 414/1296- 561/1377) -26.6 E34-E22 N29-N36
      2  139  927 (N16E66 L=140)   559     97        0 ( 121/ 895- 163/ 962)  23.2 E69-E62 N13-N17
      3  592  769 (S02E31 L=175)  1964    391        3 ( 490/ 655- 744/ 845)   0.0 E41-E22 S06-N01
      4  737  593 (S15E27 L=178)  2006    433        6 ( 619/ 458- 808/ 742)  33.0 E40-E19 S21-S07
      5  518  586 (S11E42 L=164)   286     89        1 ( 495/ 549- 541/ 620)  26.6 E45-E39 S12-S09
      6  876  492 (S24E22 L=183)   571    160        1 ( 827/ 445- 959/ 528)  11.3 E26-E16 S26-S22
      7  491  450 (S19E50 L=156)   985    244        1 ( 429/ 381- 551/ 524) -79.7 E53-E47 S23-S12
  

    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