tmpscr(2/2): irs1_sample1a1_lrxp.fits

display range: -40,40 Subtract "bkg" image Fit and subtract the background using object mask Merge 9 1D images

imcomb tmp4a,tmp4b bkg.fits offset=tmp.offset
median bkg bkg 27 9
blkav bkg bkg 1 9
magnify bkg bkg 1. 25. y1=1.01 y2=9.0 interpo=linear boundar=nearest fluxcon-
#Correct the spectral distorsion
imdel tmps.fits
imarith tmpe.fits - bkg tmps.fits
!../align tmps.fits irs1_sample1a1_lrx1.fits irs1_cont_lra.dat
imdel tmps.fits
imarith tmpe_1.fits - bkg tmps.fits
!../align tmps.fits irs1_sample1a1_lrx2.fits irs1_cont_lra.dat
imdel tmps.fits
imarith tmpe_2.fits - bkg tmps.fits
!../align tmps.fits irs1_sample1a1_lrx3.fits irs1_cont_lra.dat
imdel tmps.fits
imarith tmpe_3.fits - bkg tmps.fits
!../align tmps.fits irs1_sample1a1_lrx4.fits irs1_cont_lra.dat
imdel tmps.fits
imarith tmpe_4.fits - bkg tmps.fits
!../align tmps.fits irs1_sample1a1_lrx5.fits irs1_cont_lra.dat
imdel tmps.fits
imarith tmpe_5.fits - bkg tmps.fits
!../align tmps.fits irs1_sample1a1_lrx6.fits irs1_cont_lra.dat
imdel tmps.fits
imarith tmpe_6.fits - bkg tmps.fits
!../align tmps.fits irs1_sample1a1_lrx7.fits irs1_cont_lra.dat
imdel tmps.fits
imarith tmpe_7.fits - bkg tmps.fits
!../align tmps.fits irs1_sample1a1_lrx8.fits irs1_cont_lra.dat
imdel tmps.fits
imarith tmpe_8.fits - bkg tmps.fits
!../align tmps.fits irs1_sample1a1_lrx9.fits irs1_cont_lra.dat
!echo "0.5 1 1 1 1 1 1 1 0.5" > extract.dat
#Prepare object mask
imdel tmp1,tmp2
imcomb irs1_sample1a1_lrx4,irs1_sample1a1_lrx5,irs1_sample1a1_lrx6 tmp1
imcomb irs1_sample1a1_lrx1,irs1_sample1a1_lrx2,irs1_sample1a1_lrx3,irs1_sample1a1_lrx7,irs1_sample1a1_lrx8,irs1_sample1a1_lrx9 tmp2
imarith tmp1 - tmp2 tmp1
imfunc tmp1 tmp1 abs
fit1d tmp1 tmp1 fit axis=1 inter- function=legendre order=1
imrepl tmp1 -1. upper=100
imrepl tmp1  0. lower=100
imarith tmp1 * -1 tmp1
imdel mask,mask2
imcopy tmp1 mask.fits
imarith mask.fits * irs1_sample1a1_lr_mask.fits mask.fits
#Fit and subtract the background using object mask
imdel tmp1
imarith irs1_sample1a1_lrx1 / irs1_cont_lra irs1_sample1a1_lrx1
!../fixbadpix irs1_sample1a1_lrx1.fits tmp1.fits y
fit1d tmp1 tmp1 fit axis=2 inter- function=spline3 order=10 high_r=3 low_r=3 nite=10 grow=0
imarith irs1_sample1a1_lrx1 - tmp1 irs1_sample1a1_lrx1
imarith irs1_sample1a1_lrx1 * irs1_cont_lra irs1_sample1a1_lrx1
imdel tmp1
imarith irs1_sample1a1_lrx2 / irs1_cont_lra irs1_sample1a1_lrx2
!../fixbadpix irs1_sample1a1_lrx2.fits tmp1.fits y
fit1d tmp1 tmp1 fit axis=2 inter- function=spline3 order=10 high_r=3 low_r=3 nite=10 grow=0
imarith irs1_sample1a1_lrx2 - tmp1 irs1_sample1a1_lrx2
imarith irs1_sample1a1_lrx2 * irs1_cont_lra irs1_sample1a1_lrx2
imdel tmp1
imarith irs1_sample1a1_lrx3 / irs1_cont_lra irs1_sample1a1_lrx3
!../fixbadpix irs1_sample1a1_lrx3.fits tmp1.fits y
fit1d tmp1 tmp1 fit axis=2 inter- function=spline3 order=10 high_r=3 low_r=3 nite=10 grow=0
imarith irs1_sample1a1_lrx3 - tmp1 irs1_sample1a1_lrx3
imarith irs1_sample1a1_lrx3 * irs1_cont_lra irs1_sample1a1_lrx3
imdel tmp1
imarith irs1_sample1a1_lrx4 / irs1_cont_lra irs1_sample1a1_lrx4
!../fixbadpix irs1_sample1a1_lrx4.fits tmp1.fits y
fit1d tmp1 tmp1 fit axis=2 inter- function=spline3 order=10 high_r=3 low_r=3 nite=10 grow=0
imarith irs1_sample1a1_lrx4 - tmp1 irs1_sample1a1_lrx4
imarith irs1_sample1a1_lrx4 * irs1_cont_lra irs1_sample1a1_lrx4
imdel tmp1
imarith irs1_sample1a1_lrx5 / irs1_cont_lra irs1_sample1a1_lrx5
!../fixbadpix irs1_sample1a1_lrx5.fits tmp1.fits y
fit1d tmp1 tmp1 fit axis=2 inter- function=spline3 order=10 high_r=3 low_r=3 nite=10 grow=0
imarith irs1_sample1a1_lrx5 - tmp1 irs1_sample1a1_lrx5
imarith irs1_sample1a1_lrx5 * irs1_cont_lra irs1_sample1a1_lrx5
imdel tmp1
imarith irs1_sample1a1_lrx6 / irs1_cont_lra irs1_sample1a1_lrx6
!../fixbadpix irs1_sample1a1_lrx6.fits tmp1.fits y
fit1d tmp1 tmp1 fit axis=2 inter- function=spline3 order=10 high_r=3 low_r=3 nite=10 grow=0
imarith irs1_sample1a1_lrx6 - tmp1 irs1_sample1a1_lrx6
imarith irs1_sample1a1_lrx6 * irs1_cont_lra irs1_sample1a1_lrx6
imdel tmp1
imarith irs1_sample1a1_lrx7 / irs1_cont_lra irs1_sample1a1_lrx7
!../fixbadpix irs1_sample1a1_lrx7.fits tmp1.fits y
fit1d tmp1 tmp1 fit axis=2 inter- function=spline3 order=10 high_r=3 low_r=3 nite=10 grow=0
imarith irs1_sample1a1_lrx7 - tmp1 irs1_sample1a1_lrx7
imarith irs1_sample1a1_lrx7 * irs1_cont_lra irs1_sample1a1_lrx7
imdel tmp1
imarith irs1_sample1a1_lrx8 / irs1_cont_lra irs1_sample1a1_lrx8
!../fixbadpix irs1_sample1a1_lrx8.fits tmp1.fits y
fit1d tmp1 tmp1 fit axis=2 inter- function=spline3 order=10 high_r=3 low_r=3 nite=10 grow=0
imarith irs1_sample1a1_lrx8 - tmp1 irs1_sample1a1_lrx8
imarith irs1_sample1a1_lrx8 * irs1_cont_lra irs1_sample1a1_lrx8
imdel tmp1
imarith irs1_sample1a1_lrx9 / irs1_cont_lra irs1_sample1a1_lrx9
!../fixbadpix irs1_sample1a1_lrx9.fits tmp1.fits y
fit1d tmp1 tmp1 fit axis=2 inter- function=spline3 order=10 high_r=3 low_r=3 nite=10 grow=0
imarith irs1_sample1a1_lrx9 - tmp1 irs1_sample1a1_lrx9
imarith irs1_sample1a1_lrx9 * irs1_cont_lra irs1_sample1a1_lrx9
#Merge 9 images
!echo "0.5 1 1 1 1 1 1 1 0.5" > extract.dat
imdel irs1_sample1a1_lrxp.fits,irs1_sample1a1_lr_mask.fits
!../pileup irs1_sample1a1_lrx
displ irs1_sample1a1_lrxp 1