UBAND 数据处理方法 - BATC home page.PDFVIP

UBAND 数据处理方法 蒋兆基 2011,6,11 。原始数据仪器星等+SDSS 星等 。测光方法： DAO_PHOTO+BERTIN_SEX 。批处理： 数据按天存放于 15 个目录，多窗口进行运算。（~48 小时） 。处理包： 0key， 1key_sky, 2key, 4key_all, 6key. （一键完成） Uband 批处理流程： 综合大家思想，用 C 编写“包包”，有的“死亡”，有的新生，有的发展。 同时发展了一些“工具”包。 拷贝实际巡天 fits 文件，在各天的数据目录内运行，批处理。 0key ： 统计出当月的 supbias, 拷贝到各天的数据目录内。 1key_sky: 分割图像，减 OVERSCAN，减SUPBIAS 图像。 用当天巡天图像生成 SUPERSKY，作为平场图像。 巡天图像除以平场，得到 p*_?.fit （* 代表天区编号，？代表 1,2,3,4 号 CCD ） 2key ： 一次拍摄的四幅 CCD 图像独立位置定标； （~0.5” ） 四幅图像利用统一的投影中心+BOK 焦面改正，位置定标。 （~0.2” ） 4key_all: 求出星表。 p*_?.fit ( DAO_PHOTO + SEX ) p*_?.cat 6key ： 仪器星等星表加上 SL （斯隆）信息。 加上 DAO_PHOTO 的变孔径测光星等。加上坏像元标记。 用已知的所有 USL 位置， 在我们的 fits 图上跑一下。 p*_?.cat + p*.usl r*_?.cat “工具”包： 1. ucac3 （全天星表）,usl （南银冠，斯隆五色星等） 模拟星空+星表输出 u3 ： 快一年以前的东西了，uband 位置定标参考系。 usl ： 斯隆 u 星等——模拟星空（目前观测天区）。键入 usl 有说明。 2 。uband 观测图像，实时求出其中心坐标， （不到 8 秒，观测时换天区用） getcoord *.fits (最原始图像） 其中，用到的命令：d4,e4,coord8, 是uband 批处理“标准”程序。 未减 Bias ，未除 Flat ， 初级位置定标， 精度足够，误差在2 角秒以内。 3. z_8 拼成（马赛克）大图， （0.2~1.0” ） 0key: 巡天 fits 文件减完 OVERSCAN 后，再减一个按月统计出的 SUPBIAS 图像。 该 BIAS 图像也是预先减过OVERSCAN 的。 做法： 把某月的 bias 文件拷到某一目录，运行 0key ! 用“中值”统计法，求出。 结果分别存放在 10bias，11bias，12bias 目录中。 可比较， 10 月的bias 形状 略不同于 11,12 月的。 bias-bias 目录，各月bias 相减图像，仅供观看。 1key_sky: （利用天光做平场） 包含单个程序： d4 转换原始图像--》 4 个 4032+20*2 =4076*4096 图像 e4b 减 OVERSCAN，统计天光，编辑FITS 头信息，图像旋转，减 SUPBIAS 图像。 f4 天光做平场。 p4 除以平场。 2key: BOK_uband CCD 特点： 1。视场大，中间有呈十字，CCD 拼接盲区 2. 非 Schmidt 系统 在 RMS0.5 的判据下：在半平方度的一幅图上，通常约有 130 颗星。 传统做法，只能匹配上不到 50 颗。外边角处差。 现在能对上 110 余颗。 误差分布均匀， RMS~0.2 X,Y 和 RA,DEC 之间的转换： 逐步进行： 天球投影（纯几何），矩阵转换（6 底片常数）。 或反之。 现在上述两者间增加了一步。 乘上 ~1 的系数，（和R 有关） 若近轴图像+Schmidt 系统，可省略此步。 例： 若为 Schmidt 系统， 公式经典, 来回转化无误差。 若为其他望远镜， 公式，系数不得而