[发明专利]一种CCD像元内部响应差异的矫正方法

权利要求书

1.一种CCD像元内部响应差异的矫正方法，其特征在于步骤如下：

(1)根据观测星的测量流量误差与观测星归一化后的中心坐标的相关性，构建像元内部响应差异二维矫正图；

(2)测量观测目标星的流量和中心位置，归一化观测目标星的中心位置；

(3)根据观测目标星归一化后的中心位置通过内插像元内部响应差异二维矫正图，计算出流量矫正量；

(4)根据流量矫正量对观测目标星测量流量进行误差矫正。

2.根据权利要求1所述的一种CCD像元内部响应差异的矫正方法，其特征在于：所述构建像元内部响应差异二维矫正图步骤如下：

(1.1)对选择的固定天区进行连续的时序测光观测，得到一组时间序列的图像；

(1.2)对步骤(1.1)中观测得到的所有图像进行CCD本底去除、平场矫正和暗电流去除处理；

(1.3)对步骤(1.2)处理后的每幅图像j中的所有星进行孔径测光，获取每颗观测星k的流量flux_kj和信噪比(S/N)_kj；同时测出每颗观测星的中心位置(xc_kj，yc_kj)，筛选出信噪比大于50的所有观测星；

(1.4)根据步骤(1.3)筛选后的观测星的流量flux_kj和中心位置xc_kj，yck_j，得到每颗观测星在整个测光时间序列上的归一化后的流量fluxnorm_kj，进而得到这一组测光时间序列的流量中值fluxMedian_k，根据二者的比值计算出流量比率RatioFlux_kj；另外，还需要完成归一化中心位置坐标(xcPixIn_kj，ycPixIn_kj)的处理，计算公式为

xcPixInkj=xckj-int(xckj)ycPixInkj=yckj-int(yckj);]]>

(1.5)将实际CCD靶面进行分组，根据观测星的中心位置(xc_kj，yc_kj)，获得每颗观测星在实际CCD靶面上的相应分组；

(1.6)对步骤(1.5)中的每个组分别进行下述处理：

首先，将1个像元至少划分为一个10×10的网格，根据步骤(1.4)得到的归一化中心位置坐标值(xcPixIn_kj，ycPixIn_kj)，得出归一化中心位置坐标值落入每个子网格中相应图像的观测星；

然后，根据上述确定的每个子网格中相应图像的观测星，确定落入该子网格中所有观测星的流量比率；并计算出其中值median(RatioFlux)，用该中值代表所在子网格的值，即得到上述划分网格的像元内响应差异二维矫正图。

3.根据权利要求1所述的一种CCD像元内部响应差异的矫正方法，其特征在于：所述步骤(3)中的内插采用三次样条函数内插以提高计算的流量矫正量的精度。