基于矢量累积的恒星拖尾星点提取方法
专利汇可以提供基于矢量累积的恒星拖尾星点提取方法专利检索，专利查询，专利分析的服务。并且本 发明 涉及基于矢量累积的 恒星 拖尾星点提取方法。首先，从原始图像中提取特征图像；然后分别计算拖尾星点的旋转 角 和长度；之后，将整幅星图沿着拖尾建造矢量以拖尾步长进行累积，提高拖尾星点的 信噪比 ，得到 能量 增强的目标，获得结果图像；最后对恒星目标通过传统质心法对目标进行检测并进行 位置 估计。本发明能够有效提高对拖尾星点 定位 精度 ；可以增强拖尾星点的信噪比和极限星等；在 飞行器 快速转动条件下可有效工作；在星敏感器非线性转动的条件下可稳定工作。，下面是基于矢量累积的恒星拖尾星点提取方法专利的具体信息内容。
1.基于矢量累积的恒星拖尾星点提取方法，针对星敏感器转动条件下恒星目标信噪比下降导致的定姿精度下降问题，以及拖尾恒星加速转动条件下提点精度误差增加的问题，其特征在于，首先，从原始图像中提取特征图像；然后分别计算拖尾星点的旋转角和长度；
之后，通过矢量累积算法，将整幅星图沿着拖尾建造矢量以拖尾步长进行累积，提高拖尾星点的信噪比，得到能量增强的目标，获得结果图像；最后对恒星目标通过传统质心法对目标进行检测并进行位置估计。
2.根据权利要求1所述的基于矢量累积的恒星拖尾星点提取方法，其特征在于，所述为计算拖尾星点的旋转角和长度，建立基于线性拖尾模型的星点数学模型：
从姿态方程W＝AV出发分析拖尾星点的形状，提出一种分析拖尾星点的线性模型，建立拖尾星点的数学模型，其中A是星敏感器输出的飞行器姿态矩阵，W和V分别表示体坐标系和惯性坐标系下的单位矢量，W和V的值可由以下方程算得：
其中(α,δ)表示恒星在天球上的赤经和赤纬，(x,y)表示星点在焦平面上位置，f表示相机的焦距；
设(u,v)表示星点在焦平面上的速度， t为时间，而相机的焦距与时
间无关，因此， 在大多数情况下，(x2+y2)/f小于0.01，于是，以时间为自变量对姿态方程求导可得近似方程：
在视场较小和飞行器姿态变化较小条件下，星点在焦平面上每个位置运动速度认为近似相等，即拖尾星点形状在焦平面上每个位置都是近似相同的，因此，提出一种基于线性拖尾模型的星点数学模型如下：
其中I0表示星敏感器上的星点总能量，(-x0,0)为曝光开始时焦平面上星点位置，(x0,
0)为曝光结束时焦平面上星点位置，σ为星点的高斯半径，拖尾星点的中心位于坐标原点，在曝光时间沿着x轴方向移动，这个模型为标准星点数学模型。
3.根据权利要求1所述的基于矢量累积的恒星拖尾星点提取方法，其特征在于，所述矢量累积算法，包括五个步骤：
1)提取星图Is中高对比度的局部图像If；
2)计算拖尾星点旋转角；
用Random变换计算拖尾星点的旋转角，Radon变换将二维图像转化到直线参数空间，星图中每条直线在直线参数空间上都有一个峰值，这个特性用于图像处理中的直线检测，在二维欧几里德空间中的基本Radon变换为：
其中I为二维图像灰度分布函数，θ为投影角，δ为冲击函数，若进行基本Radon变换，当特征图像投影在拖尾方向时，能量最集中，沿着别的方向能量较分散，如果投影方向垂直于拖尾方向，投影能量最分散，
但是，由于图像的形状为矩形，在相同θ条件下积分路径随着ρ的不同而不同，即，投影值不仅仅与灰度值相关，而且与积分路径相关，基本Radon变换结果不能准确反映图像上的直线信息，提出一种归一化的Radon变换算法解决这个问题，如公式(3.2)所示，积分路径被归一化，于是，在Radon参数空间中的投影值在相同积分路径下进行比较；
在归一化Radon变换的基础上，定义目标函数f(θ)：
目标函数通过以下方式建立：(I)归一化投影曲线；(II)对该曲线进行积分，θ越接近拖尾旋转角时，函数值越小，通过归一化Radon变换计算实际拖尾模型和标准拖尾模型之间的旋转角，定义解算出的旋转角为θc，实际旋转角为θa，其中：
当θ＝θc时：
当θ≠θc时，
通过下式计算旋转角：
θc＝arg{min[f(θ)]} (3.7)，
3)计算拖尾星点长度；
待解算拖尾长度的图像为拖尾方向为水平的特征图像，用标准星点模型进行分析，积分长度用L表示，在水平方向上对星点数学模型公式(2.4)进行积分并对其归一化：
将公式(2.5)、(2.6)代入公式(3.8)，得到：
分析沿x轴方向的积分，x轴方向的子函数为：
不同L下条件下，g1(x)有不同性质，因此定义函数：
提出一个以L为自变量的基于残差平方和的目标函数e(L)来估计fv(L)以得到高精度的拖尾星点长度，当该函数值达到最小时，L值最接近实际拖尾星点长度：
其中，参数 和 通过下式计算：
其中，L∈[1,L0]，推导过程：设计分段函数fp(L)以估计函数fv(L)：
模型的参数使用最小二乘法来估计，
4)矢量累积；
整幅星图通过拖尾星点长度沿着旋转角方向累积，得到结果图像Ir，在解算的旋转角θc和长度Lc基础上，将原始图像上沿着图像中心旋转θc，并以长度Lc进行累积，表达式为：
Ir(x,y)＝∫ζIs(ξ,η)ds (3.17)
其中，ξ为一条直线，起点为(x,y)，终点为(x+cosθcLc,y+sinθcLc)，经过矢量累积后，获得结果图像Ir(x,y)，
5)使用质心法处理结果图像Ir求解星点坐标；
其中(x,y)为图像坐标点，uxy为此点的强度值，Ω为星点所在的区域，由矢量累积方法可知，理论上，质心位置为拖尾星点的起始点位置。
基于矢量累积的恒星拖尾星点提取方法
技术领域
背景技术
星敏感器的内部处理算法包括图像预处理、星点位置估计、星图识别、星图匹配和姿态估计。星点位置估计是星图识别和匹配的基础。只有在星点位置估计准确且可靠的条件下，星敏感器才能输出高精度的姿态信息。在动态条件下，星点会产生图像拖尾现象。图像拖尾现象对飞行器姿态估计有非常重要的影响。
发明内容
其中(α,δ)表示恒星在天球上的赤经和赤纬，(x,y)表示星点在焦平面上位置，f表示相机的焦距；
设(u,v)表示星点在焦平面上的速度， t为时间，而相机的焦距与
时间无关，因此， 在大多数情况下，(x2+y2)/f小于0.01。于是，以时间为自变量对姿态方程求导可得近似方程：
在视场较小和飞行器姿态变化较小条件下，星点在焦平面上每个位置运动速度可以认为近似相等，即拖尾星点形状在焦平面上每个位置都是近似相同的。这些发现证明了单帧星图上的所有星点拖尾方向和长度都是相同的。因此，提出一种分析拖尾星点数学模型。
0)为曝光结束时焦平面上星点位置，σ为星点的高斯半径，拖尾星点的中心位于坐标原点，在曝光时间沿着x轴方向移动。这个模型为标准星点数学模型；
上面公式可以写成两个独立函数的乘积形式：I(x,y)＝I1(x)I2(x)，其中：
星点数学模型的一个例子如附图2所示。在附图2(a)和(b)中分别显示了线性拖尾星点模型的三维曲面和二维投影图像。附图2(c)为I1(x)函数形状，附图2(d)为I2(y)函数形状。
由于实际拖尾星点存在平移和旋转，因此，我们可通过标准拖尾模型的几何变换(旋转和平移)来描述实际拖尾星点。我们的模型是一个在小视场和飞行器姿态变化较小条件下的线性模型，并且此线性模型有明确的表达式。
1)提取星图Is中高对比度的局部图像If；
2)计算拖尾星点旋转角；
本发明用Random变换计算拖尾星点的旋转角。Radon变换将二维图像转化到直线参数空间。星图中每条直线在直线参数空间上都有一个峰值。这个特性可以用于图像处理中的直线检测。在二维欧几里德空间中的基本Radon变换为：
其中I为二维图像灰度分布函数，θ为投影角，δ为冲击函数，ρ为(x,y)平面内直线到原点的距离。若进行基本Radon变换，当特征图像投影在拖尾方向时，能量最集中，如附图3(a)所示。沿着别的方向能量较分散，如附图3(b)所示。如果投影方向垂直于拖尾方向，投影能量最分散，如附图3(c)所示。
其中：
目标函数通过以下方式建立：(I)归一化投影曲线；(II)对该曲线进行积分。θ越接近拖尾旋转角时，函数值越小。
如附图4所示，阴影部分的面积为f(θ)函数的值。图中的三条曲线分别为rN0(ρ,θa)，rN0(ρ,θ1)，rN0(ρ,θ2)归一化函数。很容易得知，当θ＝θa时，r(ρ,θ)曲线的能量最集中。因此，阴影部分的面积最小。数学证明如下：
如附图5所示，模糊星点的图像在θ＝θc和θ≠θc的条件下被分别绘制在(a)和(b)图上。
然后，在正下方绘制出归一化投影曲线r(ρ,θc)和r(ρ,θ)。
当θ＝θc时：
当θ≠θc时，
最后，通过下式计算旋转角：
θc＝arg{min[f(θ)]} (3.7)
将公式(2.5)、(2.6)代入公式(3.8)，得到：
分析沿x轴方向的积分，x轴方向的子函数为：
不同积分步长L条件下，g1(x)有不同性质，因此定义函数：
提出一个以L为自变量的基于残差平方和的目标函数e(L)来估计fv(L)以得到高精度的拖尾星点长度。当该函数值达到最小时，L值最接近实际拖尾星点长度：
其中，参数 和 通过下式计算：
其中，L∈[1,L0]。推导过程如下：
设计分段函数fp(L)以估计函数fv(L)：
模型的参数使用最小二乘法来估计。
求解上式得到参数 和
附图6给出一个示例。拖尾的实际长度为20.0，使用总残差函数计算出的长度为21.0。
该结果在误差范围之内。
Ir(x,y)＝∫ζIs(ξ,η)ds (3.17)
其中，ζ为一条直线，起点为(x,y)，终点为(x+cosθcLc,y+sinθcLc)。经过矢量累积后，获得结果图像Ir(x,y)。此图像的特点为：该图像上的拖尾星点经过能量累积，成为一个单峰函数，峰值为拖尾恒星的起始点。附图7为矢量累积的示意图。左图圆圈中心为拖尾星点的端点，右图圆圈中心为对原始图像进行累积后的峰值位置。这两个圆圈在图像上的位置相同。
5)对原始图像进行方向累积后的结果图像上，使用质心法计算星点坐标；
其中(x,y)为图像坐标点，uxy为此点的强度值，Ω为星点所在的区域。由矢量累积方法可知，理论上，质心位置为拖尾星点的起始点位置。
下面从理论上推导此偏差。在解算得到拖尾旋转角和星点长度的条件下，星点的质心可以通过下式获得，可得公式：
其中，x2＝x1+cos△θc△Lc,y2＝y1+sin△θc△Lc,△θc＝θc-θr和△Lc＝Lc-Lr。θc和θr分别为拖尾星点的解算旋转角和实际旋转角。Lc和Lr分别为拖尾星点的长度解算值和实际值；
Ir(x,y)的解析表达式为：
当x1<-x0且x2<-x0,
Ir(x,y)＝G1(x1,x2) (3.21)
当x1<-x0且–x0≤x2
当x1<-x0且x2≥x0,
Ir(x,y)＝G1(x1,-x0)+G2(-x0,x0)+G3(x0,x2) (3.23)
当–x0≤x1
当–x0≤x1
当x1≥x0且x2≥x0,
Ir(x,y)＝G3(x1,x2) (3.26)为简化计算，分别在三个区间[-∞,-x0]，[-x0,x0]和[x0,∞]上定义三个函数G1(x1,x2)、G2(x1,x2)、G3(x1,x2)。下面推导这三个函数：
当△θc为零时,公式(3.20)可以被简化为：
当x1<-x0且x2<-x0时，得：
当-x0≤x1
当△θc不为零时，可得：
当x1<-x0且x2<-x0,
当–x0≤x1
在前文中证明了在星图中所有的条痕都可以使用相同的模型建模。然后分析了图像上所有星点的精度。在不同的星点旋转角和长度条件下进行仿真实验以验证Ir(x,y)的解析表达式。由于篇幅所限，只列出其中一颗星的精度结果作为示例，如附图9和附图10所示。需要注意的是，所有恒星的精度分析结果都是相似的；
仿真拖尾星点的基本设计参数如下：I0为50.0，x0为10.0，σ为5.0。附图9和附图10显示了δx和δy随着△Lc和△θc的变化。│δy│随着△θc增加而增加。同时│δy│随着△Lc增加而增加。
当△Lc小于等于零时，│δx│随着△θc值的增加而增加。但是当△Lc大于零时，│δx│随着△θc的增加先减少后增加。根据公式(3.20)至(3.24)解算出理论值用于对比。从附图9和附图10这两幅图像中可知，理论值与实验值符合的很好。
2、可以增强拖尾星点的信噪比和极限星等；
3、在飞行器快速转动条件下可有效工作；
4、在星敏感器非线性转动的条件下可稳定工作。
附图说明
图3特征图像在不同方向的投影示意图；
图4 f(θ)函数示意图；
图5 f(θ)函数性质的数学证明
图6总残差值e(L)曲线示例图；
图7矢量累积示意图；
图8拖尾星点位置估计精度的模型示意图；
图9以△Lc和△θc为变量的函数δx；
图10以△Lc和△θc为变量的函数δy。
具体实施方式
1)从姿态方程出发分析拖尾星点的形状，提出一种分析拖尾星点的线性模型，建立拖尾星点数学模型；
2)实行矢量累积算法估计星点位置，包括五个子步骤：提取星图中高对比度的局部图像；通过归一化Radon变换计算实际拖尾模型和标准拖尾模型之间的旋转角；通过这个旋转角旋转局部图像，计算拖尾星点的长度；整幅星图通过拖尾星点长度沿着旋转角方向进行累积得到结果图像；使用质心法处理结果图像求解星点坐标；
3)提出一个数学模型估计拖尾星点位置的精度，分析解算的旋转角偏差和拖尾星点长度偏差带来的拖尾星点位置偏差误差。
| 标题 | 发布/更新时间 | 阅读量 |
|---|---|---|
| 基于恒星时角的姿态控制方法和系统 | 2020-05-11 | 13 |
| 用于恒星模拟器成像系统的宽带通滤光片 | 2020-05-15 | 901 |
| 恒星挤出机 | 2020-05-11 | 620 |
| 一种恒星自行修正方法 | 2020-05-11 | 635 |
| 一种紫外恒星地球模拟器 | 2020-05-12 | 188 |
| 恒星及绕其运行的行星排列表 | 2020-05-14 | 604 |
| 静止气象卫星恒星敏感的成像试验方法 | 2020-05-15 | 957 |
| 一种自动恒星转角自锁结构 | 2020-05-12 | 620 |
| 恒星日与太阳日示教板 | 2020-05-14 | 628 |
| 一种恒星轴调差差速器 | 2020-05-14 | 197 |
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台，是提供专利分析，专利查询，专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据，免费查、免费专利分析。
专利汇分析报告产品可以对行业情报数据进行梳理分析，涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度，系统通过AI智能系统对图表进行解读，只需1分钟，一键生成行业专利分析报告。
