pro STDIMGPLOT2,a,hdr
;+
;Name:
;       STDIMGPLOT2
;
;Purpose:
;    generates standard plots for analysis of 64x1024 images 
;
;Usage:
;       STDIMGPLOT2,A,Hdr
;
;Inputs:
;       A   = image
;       Hdr = FITS header
;
; 
; %W% %H% :LASCO IDL LIBRARY
;
;-
naxis1 = fxpar(hdr,'NAXIS1')
naxis2 = fxpar(hdr,'NAXIS2')
IF ((naxis1 NE 64) OR (naxis2 NE 1024))  THEN BEGIN
   PRINT,'ERROR, STDIMGPLOT2:  Wrong image size.  Should be 64x1024'
   RETURN
ENDIF
!p.multi=[0,3,2,1,1]
!x.title='Column Number'
!y.title='DN'
row=512
    rs = strtrim(string(row),2)
    arow = a(*,row)
    good = WHERE(arow NE 0)
    s = stdev(arow(good))
    av = median(arow(good))
IF (av GT 500) THEN BEGIN
   ymin = (av-150)>0
   ymax = av+150
ENDIF ELSE BEGIN
   ymin = (av-100)>0
   ymax = av+100
ENDELSE
    plot,a(*,row),/ynozero,title='Line  '+rs+' StdDev =  '+strtrim(string(s),2)+'  Avg = '+strtrim(av,2), yrange=[ymin,ymax]

filename = FXPAR(hdr, 'FILENAME')
date_obs = FXPAR(hdr, 'DATE-OBS')
time_obs = STRMID(FXPAR(hdr, 'TIME-OBS'),0,8)
camera   = STRTRIM(FXPAR(hdr, 'DETECTOR'),2)
port     = STRTRIM(FXPAR(hdr, 'READPORT'),2)

!x.title='Row Number'
for i=0,2 do begin
    case I of
        0:  col = 5
        1:  col = 15
        2:  col = 30
    endcase
    cs = strtrim(string(col),2)
    acol = a(col,*)
    good = WHERE(acol NE 0)
    md = MEDIAN(acol(good))
    av = average(acol(good))
    good = WHERE((acol NE 0) AND ((acol-md) LT 10))
    s = stdev(acol(good))
    av = average(acol(good))
    plot,a(col,*),/ynozero,title='Column  '+cs+' StdDev = '+strtrim(string(s),2)+'  Avg = '+strtrim(av,2), yrange=[ymin,ymax]
    PRINT,filename+' '+date_obs+' '+time_obs+' '+camera+' '+port,col,av,s
endfor

XYOUTS, 0.7, 0.9,  /NORMAL, 'FILENAME = '+filename
XYOUTS, 0.7, 0.85, /NORMAL, 'DATE     = '+date_obs
XYOUTS, 0.7, 0.8,  /NORMAL, 'TIME     = '+time_obs
XYOUTS, 0.7, 0.75, /NORMAL, 'CAMERA   = '+camera
XYOUTS, 0.7, 0.7,  /NORMAL, 'PORT     = '+port

av=median(a)
TV, /normal, BYTSCL(REBIN(a,naxis1/2,naxis2/2), av-10, av+10), 0.93, 0.0
!p.multi=0
empty
!p.title=''
!x.title=''
!y.title=''
!p.psym=0
return
end