noao
nproto
#!rm tmp0.tbl
#!ls FMSA*a.fits | tail -n 25 > tmp0.tbl
#imdel tmp1
#imcomb @tmp0.tbl tmp1 reject=sigclip lsigma=2 hsigma=2
#imfunc tmp1 tmp1 abs
#fit1d tmp1 tmp1 fit axis=2 inter- function=legendre order=1 high_reject=2 low_reject=3 niterate=5 grow=0
#imrepl tmp1 0. upper=0.
#imrepl tmp1 65000. lower=65000.
#!rm tmp.tbl tmp1.tbl
#imhist tmp1 z1=-5 z2=65005 nbins=6501 hist_type=cumulative listout+ > tmp.tbl
### どの程度夜光輝線が残っているかでマスクする個数を調整して下さい。
#!../thresh tmp.tbl 50000 tmp1 > tmp1.tbl
#cl < tmp1.tbl
#imstati tmp1 upper=0.5
#displ tmp1 1

!rm tmp1.tbl
!ls FMSA*x.fits | tail -n 25 > tmp1.tbl
#imarith @tmp1.tbl * tmp1 @tmp1.tbl
!rm tmp2.tbl tmp3.tbl tmp4.tbl
!sed -e 's/x\.fits/b\.fits/g' tmp1.tbl > tmp2.tbl
!sed -e 's/fits/fits[100:1950,*]/g' tmp1.tbl > tmp3.tbl
#!sed -e 's/fits/fits[41:480,*]/g' tmp1.tbl > tmp3.tbl
#!sed -e 's/fits/fits[481:920,*]/g' tmp1.tbl > tmp3.tbl
#!sed -e 's/fits/fits[1021:1500,*]/g' tmp1.tbl > tmp3.tbl
#!sed -e 's/fits/fits[1501:1980,*]/g' tmp1.tbl > tmp3.tbl
imdel @tmp2.tbl,raster,tmp1
blkav @tmp3.tbl @tmp2.tbl 2000 1
fit1d @tmp2.tbl @tmp2.tbl difference axis=2 inter- function=legendre order=1 high_reject=1 niterate=20 grow=0
imcomb @tmp2.tbl tmp1 reject=minmax nlow=24 nhigh=0
irmosaic @tmp2.tbl,tmp1 raster.fits tmp4.tbl 26 1 nxoverlap=0
displ raster 1 zs+
!rm tmp1.tbl tmp2.tbl tmp3.tbl
listpix raster[1:1,*] > tmp1.tbl
!cut -d" " -f3 tmp1.tbl > tmp2.tbl
!cut -f1 ../irs1_slit_spine.tbl > tmp3.tbl
!paste tmp3.tbl tmp2.tbl > raster11.tbl
!rm tmp1.tbl tmp2.tbl tmp3.tbl
listpix raster[2:2,*] > tmp1.tbl
!cut -d" " -f3 tmp1.tbl > tmp2.tbl
!cut -f1 ../irs1_slit_spine.tbl > tmp3.tbl
!paste tmp3.tbl tmp2.tbl > raster12.tbl
!rm tmp1.tbl tmp2.tbl tmp3.tbl
listpix raster[3:3,*] > tmp1.tbl
!cut -d" " -f3 tmp1.tbl > tmp2.tbl
!cut -f1 ../irs1_slit_spine.tbl > tmp3.tbl
!paste tmp3.tbl tmp2.tbl > raster13.tbl
!rm tmp1.tbl tmp2.tbl tmp3.tbl
listpix raster[4:4,*] > tmp1.tbl
!cut -d" " -f3 tmp1.tbl > tmp2.tbl
!cut -f1 ../irs1_slit_spine.tbl > tmp3.tbl
!paste tmp3.tbl tmp2.tbl > raster14.tbl
!rm tmp1.tbl tmp2.tbl tmp3.tbl
listpix raster[5:5,*] > tmp1.tbl
!cut -d" " -f3 tmp1.tbl > tmp2.tbl
!cut -f1 ../irs1_slit_spine.tbl > tmp3.tbl
!paste tmp3.tbl tmp2.tbl > raster15.tbl
!rm tmp1.tbl tmp2.tbl tmp3.tbl
listpix raster[6:6,*] > tmp1.tbl
!cut -d" " -f3 tmp1.tbl > tmp2.tbl
!cut -f1 ../irs1_slit_spine.tbl > tmp3.tbl
!paste tmp3.tbl tmp2.tbl > raster21.tbl
!rm tmp1.tbl tmp2.tbl tmp3.tbl
listpix raster[7:7,*] > tmp1.tbl
!cut -d" " -f3 tmp1.tbl > tmp2.tbl
!cut -f1 ../irs1_slit_spine.tbl > tmp3.tbl
!paste tmp3.tbl tmp2.tbl > raster22.tbl
!rm tmp1.tbl tmp2.tbl tmp3.tbl
listpix raster[8:8,*] > tmp1.tbl
!cut -d" " -f3 tmp1.tbl > tmp2.tbl
!cut -f1 ../irs1_slit_spine.tbl > tmp3.tbl
!paste tmp3.tbl tmp2.tbl > raster23.tbl
!rm tmp1.tbl tmp2.tbl tmp3.tbl
listpix raster[9:9,*] > tmp1.tbl
!cut -d" " -f3 tmp1.tbl > tmp2.tbl
!cut -f1 ../irs1_slit_spine.tbl > tmp3.tbl
!paste tmp3.tbl tmp2.tbl > raster24.tbl
!rm tmp1.tbl tmp2.tbl tmp3.tbl
listpix raster[10:10,*] > tmp1.tbl
!cut -d" " -f3 tmp1.tbl > tmp2.tbl
!cut -f1 ../irs1_slit_spine.tbl > tmp3.tbl
!paste tmp3.tbl tmp2.tbl > raster25.tbl
!rm tmp1.tbl tmp2.tbl tmp3.tbl
listpix raster[11:11,*] > tmp1.tbl
!cut -d" " -f3 tmp1.tbl > tmp2.tbl
!cut -f1 ../irs1_slit_spine.tbl > tmp3.tbl
!paste tmp3.tbl tmp2.tbl > raster31.tbl
!rm tmp1.tbl tmp2.tbl tmp3.tbl
listpix raster[12:12,*] > tmp1.tbl
!cut -d" " -f3 tmp1.tbl > tmp2.tbl
!cut -f1 ../irs1_slit_spine.tbl > tmp3.tbl
!paste tmp3.tbl tmp2.tbl > raster32.tbl
!rm tmp1.tbl tmp2.tbl tmp3.tbl
listpix raster[13:13,*] > tmp1.tbl
!cut -d" " -f3 tmp1.tbl > tmp2.tbl
!cut -f1 ../irs1_slit_spine.tbl > tmp3.tbl
!paste tmp3.tbl tmp2.tbl > raster33.tbl
!rm tmp1.tbl tmp2.tbl tmp3.tbl
listpix raster[14:14,*] > tmp1.tbl
!cut -d" " -f3 tmp1.tbl > tmp2.tbl
!cut -f1 ../irs1_slit_spine.tbl > tmp3.tbl
!paste tmp3.tbl tmp2.tbl > raster34.tbl
!rm tmp1.tbl tmp2.tbl tmp3.tbl
listpix raster[15:15,*] > tmp1.tbl
!cut -d" " -f3 tmp1.tbl > tmp2.tbl
!cut -f1 ../irs1_slit_spine.tbl > tmp3.tbl
!paste tmp3.tbl tmp2.tbl > raster35.tbl
!rm tmp1.tbl tmp2.tbl tmp3.tbl
listpix raster[16:16,*] > tmp1.tbl
!cut -d" " -f3 tmp1.tbl > tmp2.tbl
!cut -f1 ../irs1_slit_spine.tbl > tmp3.tbl
!paste tmp3.tbl tmp2.tbl > raster41.tbl
!rm tmp1.tbl tmp2.tbl tmp3.tbl
listpix raster[17:17,*] > tmp1.tbl
!cut -d" " -f3 tmp1.tbl > tmp2.tbl
!cut -f1 ../irs1_slit_spine.tbl > tmp3.tbl
!paste tmp3.tbl tmp2.tbl > raster42.tbl
!rm tmp1.tbl tmp2.tbl tmp3.tbl
listpix raster[18:18,*] > tmp1.tbl
!cut -d" " -f3 tmp1.tbl > tmp2.tbl
!cut -f1 ../irs1_slit_spine.tbl > tmp3.tbl
!paste tmp3.tbl tmp2.tbl > raster43.tbl
!rm tmp1.tbl tmp2.tbl tmp3.tbl
listpix raster[19:19,*] > tmp1.tbl
!cut -d" " -f3 tmp1.tbl > tmp2.tbl
!cut -f1 ../irs1_slit_spine.tbl > tmp3.tbl
!paste tmp3.tbl tmp2.tbl > raster44.tbl
!rm tmp1.tbl tmp2.tbl tmp3.tbl
listpix raster[20:20,*] > tmp1.tbl
!cut -d" " -f3 tmp1.tbl > tmp2.tbl
!cut -f1 ../irs1_slit_spine.tbl > tmp3.tbl
!paste tmp3.tbl tmp2.tbl > raster45.tbl
!rm tmp1.tbl tmp2.tbl tmp3.tbl
listpix raster[21:21,*] > tmp1.tbl
!cut -d" " -f3 tmp1.tbl > tmp2.tbl
!cut -f1 ../irs1_slit_spine.tbl > tmp3.tbl
!paste tmp3.tbl tmp2.tbl > raster51.tbl
!rm tmp1.tbl tmp2.tbl tmp3.tbl
listpix raster[22:22,*] > tmp1.tbl
!cut -d" " -f3 tmp1.tbl > tmp2.tbl
!cut -f1 ../irs1_slit_spine.tbl > tmp3.tbl
!paste tmp3.tbl tmp2.tbl > raster52.tbl
!rm tmp1.tbl tmp2.tbl tmp3.tbl
listpix raster[23:23,*] > tmp1.tbl
!cut -d" " -f3 tmp1.tbl > tmp2.tbl
!cut -f1 ../irs1_slit_spine.tbl > tmp3.tbl
!paste tmp3.tbl tmp2.tbl > raster53.tbl
!rm tmp1.tbl tmp2.tbl tmp3.tbl
listpix raster[24:24,*] > tmp1.tbl
!cut -d" " -f3 tmp1.tbl > tmp2.tbl
!cut -f1 ../irs1_slit_spine.tbl > tmp3.tbl
!paste tmp3.tbl tmp2.tbl > raster54.tbl
!rm tmp1.tbl tmp2.tbl tmp3.tbl
listpix raster[25:25,*] > tmp1.tbl
!cut -d" " -f3 tmp1.tbl > tmp2.tbl
!cut -f1 ../irs1_slit_spine.tbl > tmp3.tbl
!paste tmp3.tbl tmp2.tbl > raster55.tbl
!rm tmp1.tbl tmp2.tbl tmp3.tbl
listpix raster[26:26,*] > tmp1.tbl
!cut -d" " -f3 tmp1.tbl > tmp2.tbl
!cut -f1 ../irs1_slit_spine.tbl > tmp3.tbl
!paste tmp3.tbl tmp2.tbl > rastermx.tbl