pro stdnoiseaiv,a, hdr
;
;	generates standard plots for noise analysis from aiv
;       images that are 64 col x 1024 rows
;
; SCCS variables for IDL use
; 
; %W% %H% :NRL Solar Physics
;
;
s=size(a)
nx=s(1)
ny=s(2)
!p.multi=[0,3,2,0,1]
print, !P.multi
!y.title='DN'

s=stdev(a(*,0))
plot,a(*,0),/ynozero,title='Line 1, Std Dev = '+strtrim(string(s),2), xtitle='Column Number'
print,'Standard deviation, line 0   = ',s

s=stdev(a(*,511))
plot,a(*,511),/ynozero,title='Line 512, Std Dev = '+strtrim(string(s),2), xtitle='Column Number'
print,'Standard deviation, line 512 = ',s

col = a(14,*)
s=stdev(col)
plot,col,/ynozero,title='Column 15, Std Dev = '+strtrim(string(s),2), xtitle='Line Number'
print,'Standard deviation, Column 15 = ',s

col = a(19,*)
s=stdev(col)
plot,col,/ynozero,title='Column 20, Std Dev = '+strtrim(string(s),2), xtitle='Line Number'
print,'Standard deviation, Column 20 = ',s

!P.multi(0)=1
col = a(49,*)
s=stdev(col)
plot,col,/ynozero,title='Column 50, Std Dev = '+strtrim(string(s),2), xtitle='Line Number'
print,'Standard deviation, Column 50 = ',s

sidebarleb,hdr
return


window, 1
!p.multi=[0,3,2,0,1]
sd = fltarr(ny)
for i=0,ny-1 do sd(i)=average(a(*,i))
!x.title='Line Number'
!y.title='DN'
plot,sd,title='Line Average',/ynozero
empty
for i=0,ny-1 do sd(i)=stdev(a(*,i))
plot,sd,title='Std Dev of Line',/ynozero
empty
sd = fltarr(nx)
for i=0,nx-1 do sd(i)=average(a(i,*))
!x.title='Column Number'
plot,sd,title='Column Average',/ynozero
empty
for i=0,nx-1 do sd(i)=stdev(a(i,*))
plot,sd,title='Std Dev of Column',/ynozero
empty
sidebarleb,hdr
window, 2
!p.multi = [0,3,2,0,1]
!p.title='Histogram of Line 512'
!x.title='Channel No'
!y.title='Frequency'
!p.psym=10
plot_io,histogram(a(*,511),min=0)>1
!p.title='Histogram of Entire Image'
plot_io,histogram(a,min=0)>1
dif=a(15,*)-a(16,*)
s=stdev(dif)/1.414
!x.title='Row Number'
!p.title='Column 15 - 16, Std Dev = '+strtrim(string(s,format='(f5.2)'),2)
plot,dif,/ynozero
nc = nx/2
dif=a(*,nc)-a(*,nc+1)
s=stdev(dif)/1.414
!x.title='Row Number'
!p.title='Row '+strtrim(string(nc),2)+' - '+strtrim(string(nc+1),2)
!p.title=!p.title+', Std Dev = '+strtrim(string(s,format='(f5.2)'),2)
plot,dif,/ynozero
sidebarleb,hdr
!p.multi=0
!p.title=''
!x.title=''
!y.title=''
empty
return
end