Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  432 1601 (N38E53 L=238)   586    125        1 ( 405/1576- 465/1628) -14.0 E56-E48 N36-N40
      2  201 1290 (N16E66 L=224)   925    206        1 ( 168/1259- 234/1332) -16.7 E72-E61 N14-N18
      3  672 1216 (N14E23 L=267)   728    233        1 ( 589/1146- 738/1303)  48.4 E28-E19 N09-N19
      4  538 1194 (N12E33 L=258)   381    103        3 ( 516/1160- 566/1240) -63.4 E34-E30 N09-N14
      5  161  914 (S08E68 L=222)  1413    132        3 ( 135/ 849- 178/ 972) -81.9 E72-E65 S11-S04
      6 1783  771 (S12W57 L=347)  3230    490        2 (1661/ 565-1890/ 960)  54.2 W48-W68 S25-N00
      7  569  886 (S07E29 L=261)   335    130        0 ( 527/ 860- 625/ 904)  -8.1 E32-E25 S08-S05
      8  883  645 (S21E08 L=282)  1260    318        2 ( 770/ 553- 966/ 776)  51.3 E16-E03 S27-S12
      9  382  616 (S26E47 L=243)   334    106        0 ( 370/ 571- 396/ 658)  71.6 E48-E44 S28-S22
     10  400  414 (S41E54 L=237)  1298    315        3 ( 353/ 384- 450/ 449)  20.6 E62-E46 S44-S37
  

    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