JSOC/proj/flatfield/pzt_flat_IDL/dark_rem_inorbit.pro - view

Return to dark_rem_inorbit.pro CVS log

Up to [Development] / JSOC / proj / flatfield / pzt_flat_IDL

File: [Development] / JSOC / proj / flatfield / pzt_flat_IDL / dark_rem_inorbit.pro (download)
Revision: 1.1, Fri Feb 18 00:21:17 2011 UTC (12 years, 7 months ago) by richard
Branch: MAIN
CVS Tags: Ver_LATEST, Ver_9-5, Ver_9-41, Ver_9-4, Ver_9-3, Ver_9-2, Ver_9-1, Ver_9-0, Ver_8-8, Ver_8-7, Ver_8-6, Ver_8-5, Ver_8-4, Ver_8-3, Ver_8-2, Ver_8-12, Ver_8-11, Ver_8-10, Ver_8-1, Ver_8-0, Ver_7-1, Ver_7-0, Ver_6-4, Ver_6-3, Ver_6-2, Ver_6-1, Ver_6-0, Ver_5-14, Ver_5-13, HEAD

IDL package for calculating pzt flatfields
2011.02.17

pro dark_rem_inorbit, imx, headers, time=time, imr, headim, filename=filename, xlow=xlow, xhigh=xhigh, ylow=ylow, yhigh=yhigh, drk=drk



ss=size(imx)

nbin=ss[1]

ns=ss[3]

if not keyword_set(xlow) then xlow=0
if not keyword_set(xhigh) then xhigh=nbin-1

if not keyword_set(ylow) then ylow=0
if not keyword_set(yhigh) then yhigh=nbin-1

nbinx=xhigh-xlow+1
nbiny=yhigh-ylow+1

; Identify dark images
exposure=getpar(headers, 'HSHIEXP')


; give exposure by hand


idx_image=where(exposure ne 0)
idx_dark=where(exposure eq 0, n_dark)





headim=headers[idx_image]


;;;;
n_image=n_elements(idx_image)


imr=fltarr(nbinx, nbiny, n_image)
dark=fltarr(nbinx,nbiny)
if keyword_set(drk) then dark=drk

dark_e=fltarr(nbin/16,nbin/16, n_image)



time=getshs(headers)
;;;







  for i=xlow, xhigh do begin
     
 
   for j=ylow, yhigh do begin

    if not keyword_set(drk) then dark[i-xlow,j-ylow]=median(imx[i,j,idx_dark])

  endfor

endfor





  for i=0, n_image-1 do begin

   imr[*,*,i]=imx[xlow:xhigh,ylow:yhigh,idx_image[i]] - dark

  endfor





end

Karen Tian