IRAF处理AGN数据步骤简介-测光与数据导出
已有 3816 次阅读
|个人分类:|系统分类:|关键词:IRAF
AGN 数据处理
图片测量appot 的测光程序为phot，phot为带参数的IRAF测光程序，因此运行之前，必须先修改相关参数。这里涉及到的要修改参数主要为 datapar,& photpar，fitskypar 和 phot，可用 epa 命令修改编辑。修改datapar参数ap & epa& datapar&&&&&&&&&&&&&&&&&&&&&&&&&&&& ［注：曾见在孔径测光中先进行了imexa来查看FWHM和天光的标准偏差］&&&&&&&&&&&&&&&&&&&&&&&&&scale&& ------&图象测量单位（每象素）fwhmpsf&----&&半宽& 值（由前面步骤得出）sigma&----&量背景计数标准偏差(取前面步骤所得的stddev值)datamin&-----&测量计数最小值(一般为0)datamax&-----&测量计数最大值（对一米镜，为65000）readnoi&------&CCD读出噪声（对一米镜，为3.9）epadu&-----&CCD增益（对一米镜，为4.0）obstime&----&观测时间（记录于文件头中的），各个天文天的记录名不同，有些为BJT，有些为DATE-OBS或者为DATE-STA的，要查看imhead中的信息后才能确定。itime&& ----&曝光时间（记录于文件头中的）&&&&&&&&&&&&&&&&&&&&&&&&&退出至 ap& 提示符下。修改photpar参数ap & epa& photpar&&&&&&&&&&&&&&&&&&&&&&&&&&& I R A FImage Reduction and Analysis FacilityPACKAGE = apphotTASK = photpars(weighti=& constant)Photometric weighting Scheme for wphot(apertur&=& 30.)& &list of aperture radii in scale unit(zmag&=& 25.)&&&& zero point of magnitude scale(mkapert=& no) &Draw apertures on the display(mode&=& ql)&&&&&&&&&&&&&&&&&&&&&&&&&apertur&------&孔径（取4～5倍的&半宽&值）&&&&&&&&&&&&&&&&&&&&&&&&&参照上页表修改完参数后用&shift& + & : &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& 切换至命令提示符&：&，用:& q&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& 命令退出至 ap& 提示符下。修改fitskypars参数ap & epa& fitskypar&&&&&&&&&&&&&&&&&&&&&&&&&& &&&&&&&&&&&&&&&&&&&&&&&&&annulus&------&观测内径（取5～6倍的&半宽&值）dannulu&----&观测内外径差sloreje&------&&剪切因子低限shiclip&----&&剪切因子低限注: 测量时不仅要采集测光对象的信息，同时还要采集对象周边的天空信息，为此这里要给出观测内径及观测内外径的差值。&&&&&&&&&&&&&&&&&&&&&&&&&退出至 ap& 提示符下。修改phot参数phot 为IRAF的测光程序，可完成手动测光和自动测光两种测光工作。自动测光的前提是提供 phot 测光目标的坐标值，为此手动测光是测光过程中必不可少的步骤，通过手动测光可获取测光目标的星等及坐标值等结果，为自动测光作出准备。Phot 为带参数的IRAF程序，使用前必须先设置相关参数，编辑命令及参数参考值如下所示：ap & epa& phot&&&&&&&&&&&&&&&&&&&&&&&&&&&&& &&&&&&&&&&&&&&&&&&&&&&&&&coords&------&观测目标坐标列表interac&------&测光方式注: 手动测光时坐标列表为空值，测光方式为交互式。&&&&&&&&&&&&&&&&&&&&&&&&&退出至 ap& 提示符下。手动测光ap & disp&& &文件名&.imh& 1&&&&&&&&&&&&&&&&&&& 使用手动测光时，首先应将待测光图片显示于ximtool命令打开的图片窗口中，然后再进行测光。输入命令ap & phot& &文件名&（无 .imh扩展名）&&&&&&&&&&& 后，就可开始测光了。测量中的一般顺序为：&对象&、 &比较星1（或A）&、&比较星2（B）&&&，方法是以测量光圈套住测量目标，按&空格键&获取测量数据，在所有的测量对象都获取数据后，在测量区域中以&q键&退出测量后，将光标置于命令窗口，这时会有几个选项可以选择：[Hit return to continue, n next image, q quit, w quit and save parameters] &enter键&继续此次测光，&q键&退出此次测量，&n键&写入测量结果并继续测光，&w键&存写测量结果并退出测光。测量所得结果将以 *.mag.1 的文件名保存于磁盘中。如此反复操作，直至测量完所有图片。读取待测目标坐标值为精确测量，需用系统的自动测光程序（即phot）再次测量相关的目标。因此必须提供测光程序以对象的坐标参数coords。通过手动测光后，相关的观测目标的坐标（xcenter, ycenter）实际上已记录于 *.mag.1 的文件中，只需将其读出并置于相应的坐标文件（如 *.coo.1 ）中即可。调用命令ap & epa& txdump&&&&&&&&&&&&&&&&&&&&&&&&&&&&& txdump就是读取文件的字段差数值。修改读取测量结果的IRAF程序 txdump 的参数，以获取测光目标的坐标值。txdump 的参数及参考值如下所示。&&&&&&&&&&&&&&&&&&&&&&&&&textfile&------&输入观测结果文件名fields&------&输出字段注: 坐标值获取时，需对每个手动测光结果（&文件名&.mag.1）分别取坐标（xcenter, ycenter）, 故此时的textfile 参数值为空。&&&&&&&&&&&&&&&&&&&&&&&&&退出至 ap& 提示符下，反复调用命令ap & txdump& &文件名&.imh.mag.1& && &文件名&.coo.1&& 将各观测图片的观测目标坐标值( xceter, ycenter )记录于文件 *.coo.1 中。(其中两个文件名中间的&&&为重定向符号)&指定待测自动测光的文件名列表调用系统的自动测光程序（即phot），除提供测光程序的坐标参数coords外，还要提供用于测光的图片名列表。这个文本文件可用前面生成t2 文件修改而得。调用命令ap & vi&&& t2&&&&&&&&&&&&&&&&&&&&&&&&&&&&& 修改 t2 内容，删除不用测量的flat 和bias 文件名，方法是将光标移至欲删除行，按&D键&一次或两次即可删除。修改完后用&shift& + & : &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& 切换至命令提示符&：&，用:& x&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& 退出至 ap& 提示符下。自动测光ap & epa& phot&&&&&&&&&&&&&&&&&&&&&&&&&&&& 再次改变IRAF自动测光程序 phot 的参数，进行自动测光。 Phot的参数值列表如下所示：&&&&&&&&&&&&&&&&&&&&&&&&&coords& ------&输入坐标文件，如果defaults不其作用，可以使用*.coo.1来包含你所要的坐标文件，也可以先将这些坐标文件名存储在一个文件用使用@t3来代替，只要image和coords中的文件对应一致即可。&&&&&&&&&&&&&&&&&&&&&&&&&参照上面表修改完参数后用&shift& + & : &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& 切换至命令提示符&：&，用:& g&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& 命令即开始自动测光，测量结果将存于 *.imh.mag.2 文件之中（go为直接运行任务）。至此，整个测光过程即告结束。最后，只需将测量结果存至相应的文本文件并读出就可以了。
数据读出利用IRAF的 txdump 程序可将测量结果 *.imh.mag.2 读出并存于文本文件中。数据读出和存放ap & epa& txdump&&&&&&&&&&&&&&&&&&&&&&&&&&&&& 再次修改txdump 程序的参数，参考值如下所示。&I R A FImage Reduction and Analysis FacilityPACKAGE = apphotTASK = txdumptextfile&=&&&&& Input apphot／daophot text datebase(s)fields&=& image, id,otime, mag, xcenter, ycenter,merr&&& Sky fileexpr&& =& yes)&&&& Boolean expression for record selection(headers=& no) & print the field headers(paremet=& yes)& print the parameters if headers if yes?(mode&=& ql)&&&&&&&&&&&&&&&&&&&&&&&&&textfile&------&输入观测结果文件名fields&------&输出字段&&&&&&&&&&&&&&&&&&&&&&&&&退出至 ap& 提示符下，再执行命令ap & txdump& *.mag.2& && &文件名&.txt&&&&&&&&&&&&& 完成测光数据的读出和存放。&(29) 退出IRAFap & bye&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& di & bye&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& cc & bye&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& cl & logout&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& $ & exit&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& (30) 退出Linux，键入命令$ & halt&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& 退至power down 提示下，关机。
参考文献：&向博主致谢！
转载本文请联系原作者获取授权，同时请注明本文来自任安炳科学网博客。链接地址：
上一篇：下一篇：
当前推荐数：0
评论 ( 个评论)
扫一扫，分享此博文
作者的其他最新博文
热门博文导读
Powered by
Copyright &internal error - server connection terminated
internal error - server connection terminated访问统计：
博客大巴模板：大战僵尸 | 作者：