Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1163 1823 (N49E18 L=154)   351     53        0 (1135/1803-1188/1838)   0.0 E20-E15 N46-N50
      2  638 1685 (N23E45 L=127)   808    202        0 ( 592/1589- 667/1766) -63.4 E49-E42 N15-N29
      3  943 1326 (N08E12 L=160)   210     27        0 ( 918/1311- 965/1342)  63.4 E13-E10 N06-N08
      4 1286 1204 (N09W10 L=182)   230    124        0 (1270/1156-1304/1261) -56.3 W08-W12 N06-N12
      5 1769  993 (N12W48 L=220)   337     68        1 (1747/ 959-1780/1031) -56.3 W46-W49 N09-N13
      6  169  924 (S30E65 L=107)   737     76        2 ( 132/ 901- 209/ 955)  13.0 E71-E58 S31-S28
      7  554  762 (S33E24 L=148)   331    136        0 ( 514/ 720- 591/ 821)  63.4 E26-E22 S36-S28
      8  255  792 (S38E54 L=118)   512     93        1 ( 210/ 770- 302/ 812)  13.0 E61-E47 S39-S35
      9  472  749 (S36E31 L=141)   273     93        2 ( 449/ 720- 511/ 785) -21.8 E32-E27 S37-S33
  

    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