CN106570902B - 基于探测器像素响应频谱获取的psf相对质心计算方法 - Google Patents
基于探测器像素响应频谱获取的psf相对质心计算方法 Download PDFInfo
- Publication number
- CN106570902B CN106570902B CN201610965410.9A CN201610965410A CN106570902B CN 106570902 B CN106570902 B CN 106570902B CN 201610965410 A CN201610965410 A CN 201610965410A CN 106570902 B CN106570902 B CN 106570902B
- Authority
- CN
- China
- Prior art keywords
- psf
- pixel response
- detector
- pixel
- target
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/22—Matching criteria, e.g. proximity measures
Landscapes
- Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明涉及一种基于探测器像素响应频谱获取的PSF相对质心计算方法,包括:利用四步相移法获得阵列探测器中每个像素响应函数的空间频谱,然后利用所得到的像素响应函数的空间频谱实现参考图PSF的重采样,进而通过拟合得到参考图PSF与目标图PSF的相对质心的最优估计。本发明的方法充分考虑到阵列探测器的像素响应函数对质心定位的影响,由于基于像素响应函数频谱对PSF进行精确重建,避免了PSF绝对质心的解算,能够大幅度提高质心定位精度。
Description
技术领域
本发明涉及天文学和空间技术领域,特别涉及一种基于探测器像素响应频谱获取的PSF(点扩散函数)相对质心计算方法。
背景技术
图像传感器的像素响应函数(PRF)一般通过将一个聚焦之后的光斑在图像传感器的像素感光面扫描得到。利用光斑聚焦到感光面上进行扫描,我们可以得到光斑位于像素内部不同位置的响应输出,从而得到像素响应的二维分布。对于高精度的质心定位来说,图像传感器的像素响应函数的二维分布是非常重要的。成像系统的点源响应被定义为点扩散函数(PSF)。如果图像传感器的像素响应函数已知,那么我们可以通过将图像传感器的像素响应函数与观测点扩散函数卷积获得精确的点源位置。这种方法被用于2009年3月份发射升空的Kepler空间望远镜。Kepler空间望远镜将标定的像素响应函数用于望远镜姿态的精确确定和焦平面目标的跟踪。在国际空间探测和天文技术领域,众多的空间天文卫星项目将图像传感器的像素响应函数的测量和标定作为卫星项目的关键技术,比如,HST,Kepler,JASMINE,Euclid等等。
随着探测技术和天文技术的发展,改进星点像质心定位方法以提高质心定位精度变得越来越重要。对于高精度的光度测量和天体测量,PSF的质心定位精度直接影响测光精度和天体测量精度。尤其是对于微角秒天体测量,PSF的相对质心定位精度需要达到微像素量级(10-6像素)。很多研究人员致力于对高精度质心定位的研究,并发展了许多经典的数值质心定位算法,如高斯拟合法、修正矩方法、中值法和寻导法等,并在此基础上后续延伸出一些经过改进的算法,如可调阈值的修正矩方法,递推的高斯曲面拟合技术等。这些算法都在不同程度上提高了测量精度。这些算法虽然不同程度地提高了测量精度,但由于仅关注于图像数据进行处理,在实际测量中精度仅能够达到亚像素量级。
传统的质心定位方法仅仅局限于算法的改进,而没有考虑到图像传感器像素物理特性的标定,没有从根本上解决问题。
发明内容
本发明的目的在于克服现有的质心定位方法定位精度有限的缺陷,从而提供一种能够有效提高质心定位精度的方法。
为了实现上述目的,本发明提供了一种基于探测器像素响应频谱获取的PSF相对质心计算方法,包括:
利用四步相移法获得阵列探测器中每个像素响应函数的空间频谱,然后利用所得到的像素响应函数的空间频谱实现参考图PSF的重采样,进而通过拟合得到参考图PSF与目标图PSF的相对质心的最优估计。
上述技术方案中,该方法具体包括:
步骤1)、利用阵列探测器采集一幅星图,所述星图中的星点像需满足奈奎斯特采样定理;
步骤2)、对步骤1)获得的星图进行预处理,所述预处理包括暗场改正、本底改正和平场改正,获得预处理之后的星图;
步骤3)、在步骤2)所得到的经过预处理的星图中,计算每颗星点像PSF灰度值峰值像素坐标,并以峰值所在像素为中心选择一个矩形区域,所述矩形区域包括星点像PSF绝大部分能量;
步骤4)、对阵列探测器进行四步相移正弦条纹照明,计算阵列探测器的像素响应函数的空间频谱;
步骤5)、从步骤3)得到的已为星点像选择矩形区域的星图中任意选取两颗星,将其中一个星点像PSF作为参考图,另一幅星点像PSF作为目标图;然后利用步骤4)得到的像素响应函数的空间频谱和参考图的强度分布获得参考图的原始图像频谱;接着利用傅里叶平移定理对参考图的原始图像频谱图在频域空间做相移,并在步骤3)得到的目标图的矩形区域内对经过相移的参考图进行重采样;
步骤6)、对相移后重采样的参考图PSF与目标图PSF进行匹配,计算其残差,最终获得参考图PSF与目标图PSF之间的质心间距。
上述技术方案中,所述矩形区域的大小由探测器的信噪比和质心定位精度的需求共同决定:在PSF信号强度一定的情况下,噪声越大取得范围越小,噪声越小,取得范围越大;质心定位精度要求越高,矩形区域越大。
上述技术方案中,在步骤4)中,所述四步相移正弦条纹采用下式表示:
其中,I为像素输出的ADU值,B为直流分量,A为交流分量,kx为x方向的空间频率,ky为y方向的空间频率,为相位初值,其取值分别为0,π/2,π,3π/2;当用四步相移正弦条纹照明时,探测器像素阵列的输出为像素响应函数和正弦条纹的卷积,即:
其中,Q(x,y)为像素响应函数在实域的表达式;
采用四步相移正弦条纹照明所得到的频域的像素响应函数的表达式为:
其中,表示傅里叶变换,PRF(x,y)表示像素响应函数,就是频域的像素响应函数,i为虚数单位。
上述技术方案中,在所述的步骤5)中,参考图重采样所得到结果的表达式为:
其中,和表示图像的傅里叶变换和逆变换,表示重采样之后的参考星,PSFref(x,y)由参考图像解出,即通过如下关系解出:
为参考星在步骤3)所选择矩形区域内的数据。
上述技术方案中,在步骤6)中,采用非线性最小二乘法或马尔科夫链蒙特卡洛算法对相移后重采样的参考图PSF与目标图PSF进行匹配。
上述技术方案中,采用非线性最小二乘法对相移后重采样的参考图PSF与目标图PSF进行匹配包括:
通过加权非线性最小二乘拟合得到参考星PSFref和目标星PSFtar之间的相对质心:
其中,(Δx,Δy)为参考星PSF和目标星PSF之间的相对质心。
本发明的优点在于:
本发明的方法充分考虑到阵列探测器的像素响应函数对质心定位的影响,由于基于像素响应函数频谱对PSF进行精确重建,避免了PSF绝对质心的解算,能够大幅度提高质心定位精度。
附图说明
图1是本发明的基于探测器像素响应频谱获取的PSF相对质心计算方法的流程图。
具体实施方式
现结合附图对本发明作进一步的描述。
本发明的方法利用四步相移法获得阵列探测器每个像素响应函数的频谱,并且利用所得像素响应函数频谱实现PSF的重采样,进而通过非线性最小二乘拟合得到PSF的相对质心的最优估计。
参考图1,本发明的基于探测器像素响应频谱获取的PSF相对质心计算方法包括以下步骤:
步骤1)、利用CCD或CMOS等阵列探测器采集一幅星图,星图可以是实验室模拟的星图或天文望远镜等设备采集的真实星图,且星点像即恒星的点扩散函数要满足奈奎斯特采样定理;
步骤2)、对步骤1)获得的星图进行预处理,所述预处理包括暗场改正、本底改正和平场改正,获得预处理之后的星图;
步骤3)、利用步骤2)所得到的经过预处理的星图,计算每颗星点像PSF灰度值峰值像素坐标,并以峰值所在像素为中心选择一个矩形区域,这个矩形区域包括星点像PSF绝大部分能量;其中,所述矩形区域的大小由探测器的信噪比和质心定位精度的需求共同决定。在PSF信号强度一定的情况下,噪声越大取得范围越小,噪声越小,取得范围越大;另外,质心定位精度要求越高,矩形区域越大。
对同一套光学系统和探测器,不同星点所形成的矩形区域的形状和大小是相同的;对不同的光学系统和探测器,星点矩形区域的形状和大小会改变。
步骤4)、对阵列探测器进行四步相移正弦条纹照明,计算阵列探测器的像素响应函数的空间频谱;
步骤5)、从步骤3)得到的已为星点像选择矩形区域的星图中任意选取两颗星,将其中一个星点像PSF作为参考图(即图1中的PSF1),另一幅星点像PSF作为目标图(即图1中的PSF2);然后利用步骤4)得到的像素响应函数的空间频谱和参考图的强度分布获得参考图的原始图像频谱;由于星点像满足采样定理,该参考图的原始图像频谱精确记录了星点像的完整信息,因此利用傅里叶平移定理对参考图的原始图像频谱图在频域空间做相移(实域空间即为任意平移),并在步骤3)得到的目标图的矩形区域对经过相移的参考图进行重采样;
步骤6)、对相移后重采样的参考图PSF与目标图PSF进行匹配,计算其残差,最终获得参考图PSF与目标图PSF之间的质心间距。
其中,在对相移后重采样的参考图PSF与目标图PSF进行匹配时可采用非线性最小二乘法、马尔科夫链蒙特卡洛(MCMC)算法等现有技术中已有的参数估计方法。
下面对本发明方法中的各个步骤做进一步的说明。
在步骤1)中探测器平面上的像素点记录的离散化的点扩散函数PSFdata可以由点扩散函数PSF和像素响应函数PRF的卷积表示如下:
恒星的点扩散函数要满足奈奎斯特采样定理,即:
2e≤λf/D (2)
其中,e表示像素尺寸,λ表示波长,f表示焦距,D表示孔径直径。
在步骤4)中,四步相移正弦条纹可以用下式表示:
其中,I为像素输出的ADU值,B为直流分量,A为交流分量,kx为x方向的空间频率,ky为y方向的空间频率,为相位初值,其取值分别为0,π/2,π,3π/2。当用四步相移正弦条纹照明时,探测器像素阵列的输出为像素响应函数和正弦条纹的卷积,即:
其中,Q(x,y)为像素响应函数在实域的表达式。
在步骤4)中,我们可以通过下式四步相移得到频域的像素响应函数:
其中,表示傅里叶变换,PRF(x,y)表示像素响应函数,就是频域的像素响应函数,i为虚数单位。
在步骤5)中选择两个PSF,在之前的步骤3)中,这两个PSF已经以峰值所在像素为中心扣取矩形区域,所述矩形区域包括星点像绝大部分能量,在步骤5)中将其中一个记为参考星PSFref,另一个图像记为目标星PSFtar。因为恒星点扩散函数是有限带宽的信号,离散化的点扩散函数PSFdata满足奈奎斯特采样定理且像素响应函数PRF经过标定,所以我们可以在探测器平面的任意位置对PSFdata进行重采样,精确地重建参考星PSF或目标星PSF。
其中,和表示图像的傅里叶变换和逆变换,表示重采样之后的参考星,PSFref(x,y)可由参考图像解出,即通过如下关系解出:
由上式可知,我们只要知道了像素响应函数的傅里叶变换(即频域的像素响应函数,也就是前述步骤5)描述中所提到的像素响应函数的空间频谱)、参考星在步骤3)所选择矩形区域内的数据(也就是前述步骤5)描述中所提到的参考星的强度分布)就能实现上述采样图像的重建,而可以通过步骤4)获得,而为已知值。
步骤6)中,通过加权非线性最小二乘拟合即可得到参考星PSFref和目标星PSFtar之间的相对质心:
其中,(Δx,Δy)为参考星PSF和目标星PSF之间的相对质心。
最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (6)
1.一种基于探测器像素响应频谱获取的PSF相对质心计算方法,包括:
利用四步相移法获得阵列探测器中每个像素响应函数的空间频谱,然后利用所得到的像素响应函数的空间频谱实现参考图PSF的重采样,进而通过拟合得到参考图PSF与目标图PSF的相对质心的最优估计;
该方法具体包括:
步骤1)、利用阵列探测器采集一幅星图,所述星图中的星点像需满足奈奎斯特采样定理;
步骤2)、对步骤1)获得的星图进行预处理,所述预处理包括暗场改正、本底改正和平场改正,获得预处理之后的星图;
步骤3)、在步骤2)所得到的经过预处理的星图中,计算每颗星点像PSF灰度值峰值像素坐标,并以峰值所在像素为中心选择一个矩形区域,所述矩形区域包括星点像PSF绝大部分能量;
步骤4)、对阵列探测器进行四步相移正弦条纹照明,计算阵列探测器的像素响应函数的空间频谱;
步骤5)、从步骤3)得到的已为星点像选择矩形区域的星图中任意选取两颗星,将其中一个星点像PSF作为参考图,另一幅星点像PSF作为目标图;然后利用步骤4)得到的像素响应函数的空间频谱和参考图的强度分布获得参考图的原始图像频谱;接着利用傅里叶平移定理对参考图的原始图像频谱图在频域空间做相移,并在步骤3)得到的目标图的矩形区域内对经过相移的参考图进行重采样;
步骤6)、对相移后重采样的参考图PSF与目标图PSF进行匹配,计算其残差,最终获得参考图PSF与目标图PSF之间的质心间距。
2.根据权利要求1所述的基于探测器像素响应频谱获取的PSF相对质心计算方法,其特征在于,所述矩形区域的大小由探测器的信噪比和质心定位精度的需求共同决定:在PSF信号强度一定的情况下,噪声越大取得范围越小,噪声越小,取得范围越大;质心定位精度要求越高,矩形区域越大。
3.根据权利要求1所述的基于探测器像素响应频谱获取的PSF相对质心计算方法,其特征在于,在步骤4)中,所述四步相移正弦条纹采用下式表示:
其中,I为像素输出的ADU值,B为直流分量,A为交流分量,kx为x方向的空间频率,ky为y方向的空间频率,为相位初值,其取值分别为0,π/2,π,3π/2;当用四步相移正弦条纹照明时,探测器像素阵列的输出为像素响应函数和正弦条纹的卷积,即:
其中,Q(x,y)为像素响应函数在实域的表达式;
采用四步相移正弦条纹照明所得到的频域的像素响应函数的表达式为:
其中,表示傅里叶变换,PRF(x,y)表示像素响应函数,就是频域的像素响应函数,i为虚数单位。
4.根据权利要求3所述的基于探测器像素响应频谱获取的PSF相对质心计算方法,其特征在于,在所述的步骤5)中,参考图重采样所得到结果的表达式为:
其中,和表示图像的傅里叶变换和逆变换,表示重采样之后的参考星,PSFref(x,y)由参考图像解出,即通过如下关系解出:
为参考星在步骤3)所选择矩形区域内的数据。
5.根据权利要求4所述的基于探测器像素响应频谱获取的PSF相对质心计算方法,其特征在于,在步骤6)中,采用非线性最小二乘法或马尔科夫链蒙特卡洛算法对相移后重采样的参考图PSF与目标图PSF进行匹配。
6.根据权利要求5所述的基于探测器像素响应频谱获取的PSF相对质心计算方法,其特征在于,采用非线性最小二乘法对相移后重采样的参考图PSF与目标图PSF进行匹配包括:
通过加权非线性最小二乘拟合得到参考星PSFref和目标星PSFtar之间的相对质心:
其中,(Δx,Δy)为参考星PSF和目标星PSF之间的相对质心。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610965410.9A CN106570902B (zh) | 2016-11-04 | 2016-11-04 | 基于探测器像素响应频谱获取的psf相对质心计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610965410.9A CN106570902B (zh) | 2016-11-04 | 2016-11-04 | 基于探测器像素响应频谱获取的psf相对质心计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106570902A CN106570902A (zh) | 2017-04-19 |
CN106570902B true CN106570902B (zh) | 2019-09-24 |
Family
ID=58536257
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610965410.9A Active CN106570902B (zh) | 2016-11-04 | 2016-11-04 | 基于探测器像素响应频谱获取的psf相对质心计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106570902B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108507564B (zh) * | 2018-06-29 | 2020-06-16 | 中国科学院国家空间科学中心 | 一种基于点扩散函数拟合的星敏感器质心定位方法 |
CN109738855B (zh) * | 2018-12-28 | 2020-09-18 | 中国科学院国家空间科学中心 | 一种卫星编队的星间基线方向的测量系统 |
CN113063371B (zh) * | 2021-03-23 | 2021-09-21 | 广东工业大学 | 面向正弦条纹的非线性自矫正结构光三维测量方法及系统 |
CN113552648A (zh) * | 2021-07-09 | 2021-10-26 | 中国测绘科学研究院 | 空间目标的天文定位及自动跟踪方法、系统及电子设备 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101013065A (zh) * | 2006-03-21 | 2007-08-08 | 北京航空航天大学 | 一种基于像素频率的星敏感器高精度校准方法 |
CN103913166A (zh) * | 2014-03-12 | 2014-07-09 | 哈尔滨工程大学 | 一种基于能量分布的星点提取方法 |
CN105023281A (zh) * | 2015-05-25 | 2015-11-04 | 中国科学院空间科学与应用研究中心 | 基于点扩散函数波前修正的星点像质心计算方法 |
CN105738073A (zh) * | 2016-02-03 | 2016-07-06 | 中国科学院国家空间科学中心 | 一种在空间频率域进行像素响应函数测量的方法 |
CN105841925A (zh) * | 2016-03-22 | 2016-08-10 | 中国科学院国家空间科学中心 | 一种基于探测器像素响应傅里叶谱获取的图像重建方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6516134B2 (ja) * | 2013-06-29 | 2019-05-22 | 堀 健治 | 位相変換作用を持つフィルター、レンズ、結像光学系及び撮像システム |
-
2016
- 2016-11-04 CN CN201610965410.9A patent/CN106570902B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101013065A (zh) * | 2006-03-21 | 2007-08-08 | 北京航空航天大学 | 一种基于像素频率的星敏感器高精度校准方法 |
CN103913166A (zh) * | 2014-03-12 | 2014-07-09 | 哈尔滨工程大学 | 一种基于能量分布的星点提取方法 |
CN105023281A (zh) * | 2015-05-25 | 2015-11-04 | 中国科学院空间科学与应用研究中心 | 基于点扩散函数波前修正的星点像质心计算方法 |
CN105738073A (zh) * | 2016-02-03 | 2016-07-06 | 中国科学院国家空间科学中心 | 一种在空间频率域进行像素响应函数测量的方法 |
CN105841925A (zh) * | 2016-03-22 | 2016-08-10 | 中国科学院国家空间科学中心 | 一种基于探测器像素响应傅里叶谱获取的图像重建方法 |
Non-Patent Citations (2)
Title |
---|
A focal plane metrology;Haitao Li, Baoquan Li, Yang Cao, Ligang Li;《Proc. SPIE 10155, Optical Measurement Technology and Instrumentation》;20161019;全文 * |
基于PSF相关的星点定位算法;郝云彩等;《空间控制技术与应用》;20151231;第41卷(第6期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN106570902A (zh) | 2017-04-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106570902B (zh) | 基于探测器像素响应频谱获取的psf相对质心计算方法 | |
Intema et al. | The GMRT 150 MHz all-sky radio survey-First alternative data release TGSS ADR1 | |
CN109632104B (zh) | 一种红外阵列辐射源校准装置及校准方法 | |
West et al. | A survey and examination of subpixel measurement techniques | |
CN106803892B (zh) | 一种基于光场测量的光场高清晰成像方法 | |
Varenius et al. | Subarcsecond international LOFAR radio images of Arp 220 at 150 MHz-A kpc-scale star forming disk surrounding nuclei with shocked outflows | |
McCormac et al. | DONUTS: A science frame autoguiding algorithm with sub-pixel precision, capable of guiding on defocused stars | |
CN111551129B (zh) | 大口径平面镜的中、低阶面形检测装置、系统及存储介质 | |
CN105716725B (zh) | 一种基于叠层扫描的相位差波前探测和图像复原方法 | |
Moore et al. | Development of a new underwater bathymetric laser imaging system: L-bath | |
US20110032586A1 (en) | Light microscope with novel digital method to achieve super-resolution | |
NL2010457C2 (en) | Hartmann wavefront measuring instrument adapted for non-uniform light illumination. | |
CN112034485A (zh) | 利用飞行时间相机的反射率感测 | |
Li et al. | Self-measurements of point-spread function for remote sensing optical imaging instruments | |
Lal et al. | Milliarcsec-scale radio structure of a matched sample of Seyfert 1 and Seyfert 2 galaxies | |
JP2014531160A (ja) | 電気光学センサ用のぼけキャリブレーションシステム及び移動マルチ焦点マルチターゲットコンステレーションを用いる方法 | |
Baldwin et al. | The point spread function in Lucky Imaging and variations in seeing on short timescales | |
Haist et al. | Towards one trillion positions | |
Renhorn et al. | High spatial resolution hyperspectral camera based on exponentially variable filter | |
JP2017208016A (ja) | 移動体検出装置、観測システム及び移動体検出方法 | |
Xie et al. | Frequency Estimation Method for Wideband Microwave Camera | |
CN106530351B (zh) | 一种基于图像传感器像素内量子效率获取的质心定位方法 | |
Stafford et al. | Experimental demonstration of a stripmap holographic aperture ladar system | |
Ketchazo et al. | A new technique of characterization of the intrapixel response of astronomical detectors | |
He et al. | On-line beam diagnostics based on single-shot beam splitting phase retrieval |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |