pro image_plot,a,hdr
;+
;Name:
;       IMAGE_PLOT
;
;Purpose:
;    generates standard plots for image analysis
;
;Usage:
;       IMAGE_PLOT,A,Hdr
;
;Inputs:
;       A   = image
;       Hdr = FITS header
;
;-
naxis1 = fxpar(hdr,'NAXIS1')
naxis2 = fxpar(hdr,'NAXIS2')
!p.multi=[0,3,2,1,1]
!x.title='Column Number'
!y.title='DN'
for i=0,2 do begin
    case I of
        0:  row = naxis1 /2
        1:  row = naxis1 /4
        2:  row = 3 * naxis1 /4
    endcase
    rs = strtrim(string(row),2)
    s = stdev(a(*,row))
    plot,a(*,row-1),/ynozero,title='Line  +rs+', $
         StdDev =  '+strtrim(string(s),2)
endfor
!x.title='Row Number'
col = naxis2 /2
cs = strtrim(string(col),2)
s=stdev(a(col-1,*))
plot,a(col-1,*),/ynozero,title='Column  +cs+', $
     Std Dev = '+strtrim(string(s),2)
lasco_sidebar,hdr
row = naxis1 /2
!p.multi=[0,3,2,1,1]
!p.title='Histogram of Line  +strtrim(string(row),2)
!y.title='Frequency'
!x.title='Channel No'
!p.psym=0
plot_io,histogram(a(*,row),min=0)>.1
!p.title='Histogram of Entire Image'
plot_io,histogram(a,min=0)>.1
sd = fltarr(1024)
for i=0,1023 do sd(i)=average(a(*,i))
!x.title='Row Number'
!y.title='DN'
plot,sd,title='Line Average',/ynozero
empty
for i=0,1023 do sd(i)=stdev(a(*,i))
plot,sd,title='Std Dev',/ynozero
lasco_sidebar,hdr
!p.multi=0
empty
!p.title=''
!x.title=''
!y.title=''
!p.psym=0
return
end