pro makeqepics

	;indir='C:\figerdev\idtl\sipinrockwell\H1RG-18-SIPIN\cold2\dqe.15Jun04\temp\'
	;skip=1
	indir='\\rabbit\raid2\snaprockwell\H2RG-32-39\cold2\dqe.19May04\'
	skip=1
	outdir=indir

	x1=0
	x2=2047
	y1=0
	y2=2047

	cd,indir
	onlist=file_search('*_140K.fits')
	off='dqe_scan_bias_140K.fits'
	im1=readfits(off)*1.

	; this loop cycles over all images
	for i=0,n_elements(onlist)-1 do begin

		; this loop cycles over the number of types of pixels on a device, i.e. 1 type for H1RG-18-SIPIN and two
		; types for H1RG-021-SIPIN
		for j=0,skip-1 do begin

			; check to see if we are analyzing data from a detector that has separate types of pixels in alternating rows
			; if so, we should analyze the odd and even rows separately. on the second pass through this for loop, we should
			; start at row number 1.
			im2=readfits(onlist(i))*1.
			diff=im2(x1:x2,j:y2:skip,1)-im2(x1:x2,j:y2:skip,0)-(im1(x1:x2,j:y2:skip,1)-im1(x1:x2,j:y2:skip,0))
			meddiff=median(diff)
			min=meddiff*.8
			max=meddiff*1.2

			print,'processing '+onlist(i)

			; set text in JPG file name to indicate which pixels were processed
			if (skip eq 1) then ext=''
			if (skip eq 2 and j eq 0) then ext='evens'
			if (skip eq 2 and j eq 1) then ext='odds'

			jpgfile = outdir+strmid(onlist(i),0,strpos(onlist(i),'.fits'))+'_'+ext+'.jpg'
			; set up to plot to memory buffer.  This plot will later be
			; written from memory to a JPEG file.
			set_plot, 'z'
			sizeim=size(diff)

			device, set_resolution=[sizeim(1),sizeim(2)]
			device, set_font='Courier'
			fsub = jpgfile		; filename to be written at bottom of plot

			tv,bytscl(diff,min=min,max=max)

			jpgimg = tvrd()
			write_jpeg, jpgfile, $
			congrid(jpgimg, sizeim(1),sizeim(2), /interp, /center), quality=70

			device, /close

		endfor
	endfor

end