pro prep

; Conversion du format .bin en format PGM des 15 images des spicules
; Preparation des donnees pour les programmes en C: spin.c etc

   im_init = bytarr(640,512)
   im      = bytarr(512,512)
   imdebruit=bytarr(512,512)  

   ;nb  = 14
   print,'number of files to process?'
   read, nb
   nb = nb - 1

   name1 = 'spi242'
   name2 = 'spi'

   for i = 0, nb do begin

   number = string( format='(2I)', i)
   no     = strcompress( number,/remove_all)
   s =  no
   if (i lt 10) then s = '0' + no

   file = name1 + s + '.bin'

   openr, 1, file
   readu, 1, im_init

   print, 'Initial image (640,512) from file ', file

   im = im_init(0:511,*)

   wdelete,0
   window, 0, xsize=512, ysize=512
   tvscl, im

   file = name1 + s +'.pgm'

   write_ppm, file, im
   print,'Initial image (512,512) format PGM is on ', file
  
   ; denoising 

   spawn,'spinbis'

   file = name2 + s +'.pgm'
   
   read_ppm, file, imdebruit
   print,'Processed image (512,512) format PGM is from file ',file 

   wdelete,1
   window,1, xsize=512, ysize=512
   tvscl, imdebruit

   close,1,2,3

   endfor

end