pro testfor0v


n1=640
n2=512
pwr=n2

   im_init   = bytarr(n1,n2)
   im        = bytarr(pwr,pwr)
   image_1   = bytarr(pwr)


   print,'number of files to process ?'
   read, nwindow

   nb = nwindow - 1


   name1 = 'spi_c242'

   for i = 0, nb do begin
   ;=====================

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

   file = name1 + ss + '.bin'

   openr, 1, file
   readu, 1, im_init
   close,1

   print, 'Read unformatted initial image (640,512) from file ', file


im=im_init(0:511,0:511)

set_plot,'x'
window,0,xsize=512,ysize=512
tv,im
;----

rep = 'yes'
print,'section in x?(yes/non)'
read,rep
print,'section number=?'
read,coupe

if rep eq 'yes' then image_1=im(coupe,*) else image_1 = im( *, coupe )

k=0
kk=0
for j=0,40 do begin
if image_1(j) eq 0B  then k=k+1
if image_1(511-j) eq 0B  then kk=kk+1
endfor

print,'i=',i,' k=',k,' kk=',kk

set_plot,'x'
window,1,xsize=pwr,ysize=pwr
plot,image_1,title='Original signal',yrange=[0,255]
;--------------------------------------------------

rep='yes'
print,'Print a copy?(yes/non)'
read,rep
if rep eq 'non' then goto,next

set_plot,'ps'
device,filename='dessin1.ps'
device,xsize=15.,ysize=15.
plot,image_1,title='Original signal',yrange=[0,255]
;--------------------------------------------------
device,/close

next:
answer='yes'
print,'Next file?(yes/non)'
read,answer
if answer eq 'non' then goto,endprog
; next file
endfor

endprog:
end