- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
SouthGalacticCapu-bandSkySurvey程序说明范舟2010-8-2蒋兆基.doc
South Galactic Cap u-band Sky Survey程序说明
范舟 2010-8-2
蒋兆基老师对老CCD图像进行了分析,并基于BATC原有pipeline写了新的u波段巡天数据处理批作业程序。
所有的程序都放在了uband server 159.226.170.53:/vega2/bok/
100503.tar # 测光批作业程序
首先我们先将其解压缩,得到一个文件夹100503/
一、图像格式转化程序,将.fits格式转化为.jpg格式的文件:
edo edo.fits (.fits文件名)
屏幕提示:
!edo.gif
!convert edo.gif edo.jpg
!rm edo.gif
然后我们可得到edo.jpg文件。
二、测光程序
为了更精确的得到测光的结果,我们采用了IRAF/DAOPHOT和burtin SEXtractor结合的测光方法。
位置定标
由于观测得到的图像位置和方向都不对,我们必须对其进行过图像的位置定标,具体来说,图像要经过平移和旋转的矫正。
coord8 edo
屏幕显示:
! wait, auto match stars. Now, pick out CCD stars
!....
! Auto shift 0.7 arc_minutes
! By using 22 stars to calculate 8 ceof, RMS: 0.48 (arcsec)
! produce new_fits_file: uedo.fit
这样我们就会得到一副图像uedo.fit,此图的位置都进行过了改正,而且fits头也进行了标准化的调整,增加了必须的信息。
2.坏线修正
badline uedo
屏幕显示
!modified bad_line of: uedo.fit
这样uedo.fit就经过了坏线处理。
3.Burtin SEXtracter测光
# 该测光同时给出了固定孔径的孔径测光以及等亮度测光
sex uedo.fit
屏幕显示:
!Measuring from: edo_a1413 / 4064 x 4064 / 32 bits FLOATING POINT data
!(M+D) Background: 89.4817 RMS: 10.1477 / Threshold: 10.1477
!Objects: detected 2123 / sextracted 1538
! All done
该程序会产生两个文件test.cat和check.fits
# test.cat: 使用不同方法(固定孔径的孔径测光以及等亮度测光)得到的测光结果。
# check.fits: 对图像中每个源分别画出了SEXtractor等亮度星等的轮廓(椭圆)和固定孔径的测光半径(圆)。
4.删除图像边缘10行之内的源,得到可信度较高的.coo文件
chbercoo uedo
屏幕显示:
!uedo.coo produced!
!all: 1271
!cp uedo.coo test.coo
生成了test.coo文件,图像中比较可靠的源的位置
5.IRAF/DAOPHOT 测光。
此测光方法基于BATC的测光批处理程序,将图像分成3X3块,然后分别对每一块进行PSF拟合。
pip231 uedo
最终得到了如下的文件:
uedo.mes #各种测光参数,??以更改
uedo.lst #psf星的参数
uedo.als #所有源的x,y和测光星等及误差
uedo.fits #原始图像减去psf拟合轮廓后残留的图像
uedo.cat6 #psf星的测光结果
coo2als uedo
屏幕显示
produce coo2als.dat
得到coo2als.dat文件。
我们怀疑很多找出的源可能是宇宙线,所以我们需要找出可能的宇宙线
checkmis uedo
!prodced miss.dat 26
得到了miss.dat文件
我们将sharpness (uedo.als最后一列) -10的源认为是宇宙线,得到宇宙线的源的列表文件cosmos.dat
getcosmos uedo
屏幕显示
!prodced cosmos.dat 107个可能的宇宙线源
!prodced uedo.coo 1164 !剔除宇宙线后真实的源
得到了cosmos.dat(宇宙线列表)和uedo.coo(干净的源列表)
由于原来宇宙线的轮廓可能会影响我们测光时拟合的参数,因此我们需要进行迭代
pip231 uedo
得到了1163个源
coo2als uedo
文档评论(0)