CN109493382A - 一种基于像元内响应的恒星高精度位置提取方法 - Google Patents

一种基于像元内响应的恒星高精度位置提取方法 Download PDF

Info

Publication number
CN109493382A
CN109493382A CN201811176455.3A CN201811176455A CN109493382A CN 109493382 A CN109493382 A CN 109493382A CN 201811176455 A CN201811176455 A CN 201811176455A CN 109493382 A CN109493382 A CN 109493382A
Authority
CN
China
Prior art keywords
pixel
coordinate
fixed star
gray value
energy distribution
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.)
Granted
Application number
CN201811176455.3A
Other languages
English (en)
Other versions
CN109493382B (zh
Inventor
陈凡胜
孙胜利
于清华
林长青
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Institute of Technical Physics of CAS
Original Assignee
Shanghai Institute of Technical Physics of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai Institute of Technical Physics of CAS filed Critical Shanghai Institute of Technical Physics of CAS
Priority to CN201811176455.3A priority Critical patent/CN109493382B/zh
Publication of CN109493382A publication Critical patent/CN109493382A/zh
Application granted granted Critical
Publication of CN109493382B publication Critical patent/CN109493382B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/70—Determining position or orientation of objects or cameras
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/60—Analysis of geometric attributes
    • G06T7/66—Analysis of geometric attributes of image moments or centre of gravity
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10032—Satellite or aerial image; Remote sensing

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • Image Analysis (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明公开了一种基于像元内响应的恒星高精度位置提取方法，根据像元内响应分布，建立点目标位置与像元灰度之间的关系，实现恒星的高精度位置提取。其位置提取步骤是：1.对获得的图像进行去除背景和非均匀校正等处理。2.选取像元灰度值最大的十帧图像，根据像元的填充因子，求出该恒星目标经过光学系统后的能量分布函数。3.根据能量分布函数和像元填充因子得到一系列恒星目标在不同位置时的像元灰度值，并通过内插值的方法建立恒星位置与像元灰度的查找表。4.根据像元的灰度值，查找对应的像元位置，由于能量分布和像元响应的对称性，得到几组像元位置。5.使用质心法从步骤4中得到的几组像元位置中选出恒星位置。

Description

一种基于像元内响应的恒星高精度位置提取方法
技术领域
本发明属于空间光电遥感探测领域运动恒星位置提取方向，结合像元填充因子对恒星进行高精度位置提取，在空间遥感探测等领域有应用前景。
背景技术
对远距离目标的探测时，目标通常在像面上以点目标的形式出现，成像小于单个像元，由于焦平面成像阵列像素分辨率的限制，定位精度为单个像元。为了提高目标的定位精度，还需要进一步的计算。
由于像元的填充因子小于1，像元内响应的分布不均匀，目标成像在像元不同位置时，像元接收到的能量不同。接收到的能量差异会对使用传统的质心法进行点目标位置提取过程引入误差。在对恒星进行观测时，由于恒星相对于相机的距离不变，可认为像元灰度的变化是由恒星运动产生的。
发明内容
为了克服现有技术的不足，本发明提供了一种基于像元内响应的高精度恒星位置提取方法。
本发明所采用的技术方案是：根据高斯函数和像元填充因子建立点目标位置和像元灰度对应表，通过查找表和质心法得到恒星在像元的位置(x,y)。
与现有技术相比，本发明的特点是在利用像元灰度值进行恒星位置提取时，考虑了像元内响应不均匀对像元灰度值的影响。
一种基于像元内响应的恒星高精度位置提取方法，步骤如下：
步骤1：对图像进行预处理，包括去除背景和非均匀校正两个步骤。根据已有的黑体定标数据计算各个像元的非均匀校正系数，记为K(i,j)，i，j分别是该像元在探测器上所处的行列坐标，像元的背景灰度为Backgoung(i，j)。DN(i，j)是经过图像预处理后坐标为(i，j)像元的灰度数据，DN’(i,j)为原始像元灰度数据。
DN(i,j)＝((DN’(i,j)–Background(i，j))·K(i,j)
步骤2：恒星经光学系统成像后，可用高斯函数描述其能量分布。由于像元内响应的差异，当目标在单个像元上移动时，像元的灰度值发生变化，在多帧图像中找出单个像元灰度值最大的十帧，根据像元内感光区域的长和宽，通过反卷积的方法得到高斯函数的参数δ。在单个像元中，感光区域对光的响应系数为1，其他区域响应系数为0。单个像元的响应表示为pixresponse，其中，a和b分别是感光区域的长和宽，D为待探测的像元区域。g是恒星目标经过光学系统后的能量分布，A是点目标成像在像元中心时的能量，δ为标准差，x，y分别是探测器上某点的横纵坐标，xo，yo光斑中心在探测器上的横纵坐标，r是像元上的某点到光斑中心的距离。B是探测器的平均响应，I是该像元的灰度值。
I＝B∫∫D g·pixresponse dxdy
步骤3：根据步骤2得到的高斯函数以及像元的填充因子得到像元灰度分布随点目标位置移动的关系。
步骤4：通过步骤3中的关系，根据像元的灰度值，通过查表法得到恒星在像元内的位置。由于像元内响应和能量分布的对称性，得到多组位置。
步骤5：根据质心法的公式(如下)，求出一组(x′,y′)，以像元中心为原点，建立平面直角坐标系，判断(x′,y′)在单个像元所处的象限。从步骤4中所得的坐标中处于该象限的坐标，即为恒星的位置(x，y)。
附图说明
图1为基于像元内响应的恒星高精度位置提取方法的具体步骤。
图2为根据步骤3得到的像元随点目标位置移动的灰度分布。
具体实施方式
下面结合附图对本发明进一步说明。
步骤1：对图像进行预处理，去除背景和非均匀校正。根据已有的黑体定标数据计算各个像元的非均匀校正系数，记为K(i,j)，i，j分别是该像元在探测器上的横纵坐标。背景为Backgoung(i，j)。DN(i，j)是经过处理后坐标为(i，j)的像元灰度值，DN’(i,j)为原始像元灰度数据。
DN(i,j)＝((DN’(i,j)–Background(i，j))·K(i,j)
步骤2：点目标经光学系统成像后，可用高斯函数描述其能量分布。由于像元内响应的差异，当目标在单个像元上移动时，像元的灰度值也有一定的变化，在多帧图像中找出单个像元灰度值最大的十帧，根据像元内感光区域的长和宽，通过反卷积的方法得到高斯函数的参数δ。在单个像元中，感光区域对光的响应系数为1，其他区域响应系数为0。单个像元的响应表示为pixresponse，，其中，a和b分别是感光区域的长和宽，D为待探测的像元区域。g是恒星目标经过光学系统后的能量分布，A是点目标成像在像元中心时的能量，δ为标准差，x，y分别是探测器上某点的横纵坐标，xo，yo分别是光斑中心在探测器上的横纵坐标，r是像元上的某点到光斑中心的距离。B是探测器的平均响应，I是该像元的灰度值。
I＝B∫∫D g·pixresponse dxdy
步骤3：根据步骤2成像的高斯函数以及像元的填充因子得到像元灰度分布随点目标位置移动的关系。
步骤4：通过步骤3中的关系，根据像元的灰度值，通过查表法得到恒星在像元内的位置。由于像元响应和能量分布的对称性，得到的位置大于一组。
步骤5：根据质心法的公式(如下)，求出一组(x′,y′)，以灰度最大的像元中心为原点，建立平面直角坐标系，判断(x′,y′)在这个像元所处的象限。从步骤4中所得的坐标中处于该象限的坐标，即为点目标的位置。

Claims (1)

1.一种基于像元内响应的恒星高精度位置提取方法，其特征在于包括以下步骤：
步骤1：对图像进行预处理，包括去除背景和非均匀校正两个步骤。根据已有的黑体定标数据计算各个像元的非均匀校正系数，记为K(i,j)，i，j分别是该像元在探测器上的横纵坐标，像元的背景灰度为Backgoung(i，j)，DN(i，j)是经过图像预处理后坐标为(i，j)的像元灰度数据，DN’(i,j)为原始像元灰度数据，
DN(i,j)＝((DN’(i,j)–Background(i，j))·K(i,j)
步骤2：恒星经光学系统成像后，可用高斯函数描述其能量分布，由于像元内响应的差异，当目标在单个像元上移动时，像元的灰度值发生变化，在多帧图像中找出单个像元灰度值最大的十帧，根据像元内感光区域的长和宽，通过反卷积的方法得到高斯函数的参数δ；在单个像元中，感光区域对光的响应系数为1，其他区域响应系数为0；单个像元的响应表示为pixresponse，其中，a和b分别是感光区域的长和宽，D为待探测的像元区域；g是恒星目标经过光学系统后的能量分布，A是点目标成像在像元中心时的能量，δ为标准差，x，y分别是探测器上某点的横纵坐标，xo，yo光斑中心在探测器上的横纵坐标，r是像元上的某点到光斑中心的距离，B是探测器的平均响应，I是该像元的灰度值；
I＝B∫∫Dg·pixresponse dxdy
步骤3：根据步骤2得到的高斯函数以及像元的填充因子得到像元灰度分布随点目标位置移动的关系；
步骤4：通过步骤3中的关系，根据像元的灰度值，通过查表法得到恒星在像元内的位置，由于像元内响应和能量分布的对称性，得到的位置大于一组；
步骤5：根据如下的质心法公式，求出一组坐标(x′,y′)，以像元中心为原点，建立平面直角坐标系，判断恒星目标在单个像元内所处的象限。从步骤4中所得的坐标中选出处于该象限的坐标，即为恒星在像元中的位置；
CN201811176455.3A 2018-10-10 2018-10-10 一种基于像元内响应的恒星高精度位置提取方法 Active CN109493382B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811176455.3A CN109493382B (zh) 2018-10-10 2018-10-10 一种基于像元内响应的恒星高精度位置提取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811176455.3A CN109493382B (zh) 2018-10-10 2018-10-10 一种基于像元内响应的恒星高精度位置提取方法

Publications (2)

Publication Number Publication Date
CN109493382A true CN109493382A (zh) 2019-03-19
CN109493382B CN109493382B (zh) 2021-08-06

Family

ID=65689432

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811176455.3A Active CN109493382B (zh) 2018-10-10 2018-10-10 一种基于像元内响应的恒星高精度位置提取方法

Country Status (1)

Country Link
CN (1) CN109493382B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113916382A (zh) * 2021-09-14 2022-01-11 中国科学院上海技术物理研究所 一种基于像素内灵敏度模型的恒星能量提取方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101907463A (zh) * 2010-07-05 2010-12-08 中国人民解放军国防科学技术大学 一种星敏感器恒星像点位置提取方法
CN103913166A (zh) * 2014-03-12 2014-07-09 哈尔滨工程大学 一种基于能量分布的星点提取方法
US20180017394A1 (en) * 2016-07-12 2018-01-18 David Byron Douglas Planetary Surveillance System

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101907463A (zh) * 2010-07-05 2010-12-08 中国人民解放军国防科学技术大学 一种星敏感器恒星像点位置提取方法
CN103913166A (zh) * 2014-03-12 2014-07-09 哈尔滨工程大学 一种基于能量分布的星点提取方法
US20180017394A1 (en) * 2016-07-12 2018-01-18 David Byron Douglas Planetary Surveillance System

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
CARL CHRISTIAN LIEBE: "Accuracy Performance of Star Trackers", 《IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS》 *
连月勇等: "星点亚像元细分定位精度分析", 《测绘科学技术学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113916382A (zh) * 2021-09-14 2022-01-11 中国科学院上海技术物理研究所 一种基于像素内灵敏度模型的恒星能量提取方法
CN113916382B (zh) * 2021-09-14 2023-09-12 中国科学院上海技术物理研究所 一种基于像素内灵敏度模型的恒星能量提取方法

Also Published As

Publication number Publication date
CN109493382B (zh) 2021-08-06

Similar Documents

Publication Publication Date Title
CN109903241B (zh) 一种tof相机系统的深度图像校准方法及系统
CN108615244B (zh) 一种基于cnn和深度滤波器的图像深度估计方法及系统
CN107255521B (zh) 一种红外图像非均匀性校正方法及系统
CN108426585B (zh) 一种光场相机的几何标定方法
CN105716542B (zh) 一种基于柔性特征点的三维数据拼接方法
CN111288967B (zh) 一种基于机器视觉的远距离高精度位移检测方法
CN109272520B (zh) 一种联合运动指导与边缘检测的自适应红外焦平面非均匀校正方法
CN110211056B (zh) 基于局部中值直方图的自适应红外图像去条纹算法
US20110293146A1 (en) Methods for Estimating Peak Location on a Sampled Surface with Improved Accuracy and Applications to Image Correlation and Registration
CN108063932B (zh) 一种光度标定的方法及装置
CN109523585B (zh) 一种基于方向相位一致性的多源遥感影像特征匹配方法
CN107403410B (zh) 一种热红外图像的拼接方法
CN110223355B (zh) 一种基于双重极线约束的特征标志点匹配方法
CN110111292B (zh) 一种红外与可见光图像融合方法
CN102752504B (zh) 一种宽视场线阵ccd相机的相对辐射校正方法
CN103905746A (zh) 亚像素级图像偏移定位及叠加方法和装置以及视频设备
CN113012234A (zh) 基于平面变换的高精度相机标定方法
CN112489137A (zh) 一种rgbd相机标定方法及系统
CN107966137A (zh) 一种基于tdiccd拼接区图像的卫星平台颤振探测方法
WO2019048904A1 (en) STEREOSCOPIC DEPTH CARTOGRAPHY AND COMBINED PHASE DETECTION IN A DOUBLE-OPENING CAMERA
CN109084959B (zh) 一种基于双目测距算法的光轴平行度校正方法
CN109493382A (zh) 一种基于像元内响应的恒星高精度位置提取方法
CN108507564B (zh) 一种基于点扩散函数拟合的星敏感器质心定位方法
CN112419427A (zh) 用于提高飞行时间相机精度的方法
CN102885631B (zh) 一种运用于平板ccd探测器的畸变校正方法

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