• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 11
    • 下载费用:30 金币  

    重庆时时彩个位杀1码: 静止卫星成像导航与配准恒星南北参数提取方法及系统.pdf

    摘要
    申请专利号:

    重庆时时彩单双窍门 www.4mum.com.cn CN201610971441.5

    申请日:

    2016.11.03

    公开号:

    CN106568437A

    公开日:

    2017.04.19

    当前法律状态:

    实审

    有效性:

    审中

    法律详情: 实质审查的生效IPC(主分类):G01C 21/02申请日:20161103|||公开
    IPC分类号: G01C21/02 主分类号: G01C21/02
    申请人: 上海卫星工程研究所
    发明人: 吕旺; 王田野; 董瑶海; 王皓; 沈毅力
    地址: 200240 上海市闵行区华宁路251号
    优先权:
    专利代理机构: 上海汉声知识产权代理有限公司 31236 代理人: 郭国中
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201610971441.5

    授权公告号:

    |||

    法律状态公告日:

    2017.05.17|||2017.04.19

    法律状态类型:

    实质审查的生效|||公开

    摘要

    本发明提供了一种静止卫星成像导航与配准恒星南北参数提取方法及系统,包括如下步骤:步骤1:构建列像元的灰度值与恒星影像南北向坐标分量的关系;步骤2:提取恒星划过列像元中心时的探测器像元灰度值;步骤3:计算恒星影像南北向坐标分量;步骤4:计算恒星影像中心划过探测器东西方向中心线时的南北向坐标分量。本发明在计算恒星影像坐标时,考虑了恒星影像的形状和亮度分布特性,大幅提高了识别精度;通过多列像元计算结果进行最小二乘拟合,降低了由探测器光学成像、电路噪声、卫星高频抖动等因素引起的星点位置高频误差。

    权利要求书

    1.一种静止卫星成像导航与配准恒星南北参数提取方法,其特征在于,包括如下步骤:
    步骤1:构建列像元的灰度值与恒星影像南北向坐标分量的关系;
    步骤2:提取恒星划过列像元中心时的探测器像元灰度值;
    步骤3:计算恒星影像南北向坐标分量;
    步骤4:计算恒星影像中心划过探测器东西方向中心线时的南北向坐标分量。
    2.根据权利要求1所述的静止卫星成像导航与配准恒星南北参数提取方法,其特征在
    于,所述步骤1包括:
    当恒星影像划过列像元中心时,设恒星影像南北向坐标分量为l,恒星影像的中心点坐
    标写为(l,L/2),其中,L为单个像元的尺寸;设列像元的第k个像元的灰度为hk,则列像元与
    恒星影像南北向坐标分量的关系如下:
    <mrow> <msub> <mi>h</mi> <mi>k</mi> </msub> <mo>=</mo> <munderover> <mo>&Integral;</mo> <mn>0</mn> <mi>L</mi> </munderover> <munderover> <mo>&Integral;</mo> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> <mo>&CenterDot;</mo> <mi>L</mi> </mrow> <mrow> <mi>k</mi> <mo>&CenterDot;</mo> <mi>L</mi> </mrow> </munderover> <mi>f</mi> <mrow> <mo>(</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>-</mo> <mi>l</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>y</mi> <mo>-</mo> <mi>L</mi> <mo>/</mo> <mn>2</mn> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>)</mo> </mrow> <mi>d</mi> <mi>x</mi> <mi>d</mi> <mi>y</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
    式中:函数f(·)表示恒星影像的点扩散函数,变量x表示恒星影像东西向坐标,变量y
    表示恒星影像南北向坐标。
    3.根据权利要求2所述的静止卫星成像导航与配准恒星南北参数提取方法,其特征在
    于,所述步骤2包括:恒星划过探测器的某一列像元时,探测器多次曝光记录了恒星划过列
    像元的整个过程,通过提取恒星影像引起该列像元总灰度值随时间的变化关系和曲线,使
    用半功率法计算恒星划过该列像元中心的时刻,并记录恒星影像处于列像元中心时各像元
    的灰度值;
    具体地,选择灰度值最高的三个列像元,并记录每个像元的灰度值。
    4.根据权利要求3所述的静止卫星成像导航与配准恒星南北参数提取方法,其特征在
    于,所述步骤3包括:
    假设步骤2中选取的三个列像元的排列序号分别为k、m和n,实测的列像元灰度值分别
    为hk、hm和hn;
    由步骤1得到恒星点扩散函数推导出来的三个列像元的灰度表达式分别如下:
    <mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>h</mi> <mi>k</mi> </msub> <mo>=</mo> <munderover> <mo>&Integral;</mo> <mn>0</mn> <mi>L</mi> </munderover> <munderover> <mo>&Integral;</mo> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> <mo>&CenterDot;</mo> <mi>L</mi> </mrow> <mrow> <mi>k</mi> <mo>&CenterDot;</mo> <mi>L</mi> </mrow> </munderover> <mi>f</mi> <mrow> <mo>(</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <mrow> <mi>x</mi> <mo>-</mo> <mi>l</mi> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mrow> <mi>y</mi> <mo>-</mo> <mi>L</mi> <mo>/</mo> <mn>2</mn> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>)</mo> </mrow> <mi>d</mi> <mi>x</mi> <mi>d</mi> <mi>y</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>h</mi> <mi>m</mi> </msub> <mo>=</mo> <munderover> <mo>&Integral;</mo> <mn>0</mn> <mi>L</mi> </munderover> <munderover> <mo>&Integral;</mo> <mrow> <mo>(</mo> <mi>m</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> <mo>&CenterDot;</mo> <mi>L</mi> </mrow> <mrow> <mi>m</mi> <mo>&CenterDot;</mo> <mi>L</mi> </mrow> </munderover> <mi>f</mi> <mrow> <mo>(</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <mrow> <mi>x</mi> <mo>-</mo> <mi>l</mi> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mrow> <mi>y</mi> <mo>-</mo> <mi>L</mi> <mo>/</mo> <mn>2</mn> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>)</mo> </mrow> <mi>d</mi> <mi>x</mi> <mi>d</mi> <mi>y</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>h</mi> <mi>n</mi> </msub> <mo>=</mo> <munderover> <mo>&Integral;</mo> <mn>0</mn> <mi>L</mi> </munderover> <munderover> <mo>&Integral;</mo> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> <mo>&CenterDot;</mo> <mi>L</mi> </mrow> <mrow> <mi>n</mi> <mo>&CenterDot;</mo> <mi>L</mi> </mrow> </munderover> <mi>f</mi> <mrow> <mo>(</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <mrow> <mi>x</mi> <mo>-</mo> <mi>l</mi> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mrow> <mi>y</mi> <mo>-</mo> <mi>L</mi> <mo>/</mo> <mn>2</mn> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>)</mo> </mrow> <mi>d</mi> <mi>x</mi> <mi>d</mi> <mi>y</mi> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
    恒星在探测器里的纵向坐标计算能够转化为最优化问题,最优化的目标是选取合适的
    恒星影像南北向坐标分量l,使点扩散函数法计算出的像元灰度值与实测的像元灰度偏差
    最??;目标函数如下:
    <mrow> <mtable> <mtr> <mtd> <mrow> <mi>Q</mi> <mo>=</mo> <msub> <mi>c</mi> <mn>1</mn> </msub> <mo>&CenterDot;</mo> <msup> <mrow> <mo>&lsqb;</mo> <mrow> <munderover> <mo>&Integral;</mo> <mn>0</mn> <mi>L</mi> </munderover> <munderover> <mo>&Integral;</mo> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> <mo>&CenterDot;</mo> <mi>L</mi> </mrow> <mrow> <mi>k</mi> <mo>&CenterDot;</mo> <mi>L</mi> </mrow> </munderover> <mi>f</mi> <mrow> <mo>(</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <mrow> <mi>x</mi> <mo>-</mo> <mi>l</mi> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mrow> <mi>y</mi> <mo>-</mo> <mi>L</mi> <mo>/</mo> <mn>2</mn> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>)</mo> </mrow> <mi>d</mi> <mi>x</mi> <mi>d</mi> <mi>y</mi> <mo>-</mo> <msub> <mi>h</mi> <mi>k</mi> </msub> </mrow> <mo>&rsqb;</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>c</mi> <mn>2</mn> </msub> <mo>&CenterDot;</mo> <msup> <mrow> <mo>&lsqb;</mo> <mrow> <munderover> <mo>&Integral;</mo> <mn>0</mn> <mi>L</mi> </munderover> <munderover> <mo>&Integral;</mo> <mrow> <mo>(</mo> <mi>m</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> <mo>&CenterDot;</mo> <mi>L</mi> </mrow> <mrow> <mi>m</mi> <mo>&CenterDot;</mo> <mi>L</mi> </mrow> </munderover> <mi>f</mi> <mrow> <mo>(</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <mrow> <mi>x</mi> <mo>-</mo> <mi>l</mi> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mrow> <mi>y</mi> <mo>-</mo> <mi>L</mi> <mo>/</mo> <mn>2</mn> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>)</mo> </mrow> <mi>d</mi> <mi>x</mi> <mi>d</mi> <mi>y</mi> <mo>-</mo> <msub> <mi>h</mi> <mi>m</mi> </msub> </mrow> <mo>&rsqb;</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>c</mi> <mn>3</mn> </msub> <mo>&CenterDot;</mo> <msup> <mrow> <mo>&lsqb;</mo> <mrow> <munderover> <mo>&Integral;</mo> <mn>0</mn> <mi>L</mi> </munderover> <munderover> <mo>&Integral;</mo> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> <mo>&CenterDot;</mo> <mi>L</mi> </mrow> <mrow> <mi>n</mi> <mo>&CenterDot;</mo> <mi>L</mi> </mrow> </munderover> <mi>f</mi> <mrow> <mo>(</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <mrow> <mi>x</mi> <mo>-</mo> <mi>l</mi> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mrow> <mi>y</mi> <mo>-</mo> <mi>L</mi> <mo>/</mo> <mn>2</mn> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>)</mo> </mrow> <mi>d</mi> <mi>x</mi> <mi>d</mi> <mi>y</mi> <mo>-</mo> <msub> <mi>h</mi> <mi>n</mi> </msub> </mrow> <mo>&rsqb;</mo> </mrow> <mn>2</mn> </msup> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
    式中:c1、c2和c3为对应的三个列像元的权重,c1、c2和c3的选取与对应列像元的灰度值
    大小相关,当Q值最小时,得出的l即为恒星影像南北向坐标分量。
    5.根据权利要求1所述的静止卫星成像导航与配准恒星南北参数提取方法,其特征在
    于,所述步骤4包括:
    步骤4.1:拟合恒星影像运动轨迹;
    步骤4.2:计算恒星影像中心划过探测器东西方向中心线时的南北向坐标分量。
    6.根据权利要求5所述的静止卫星成像导航与配准恒星南北参数提取方法,其特征在
    于,所述步骤4.1包括:
    恒星影像划过探测器面阵时,会划过多列像元,则认为恒星影像中心在探测器面阵中
    的做匀速直线运动,将恒星在探测器面阵内的南北向坐标分量与东西向坐标分量描述为线
    性关系:
    l=Kx+B (4)
    式中:K表示恒星影像运动的斜率、B分别表示恒星影像南北向坐标偏量;变量x表示恒
    星影像东西向坐标;变量l表示恒星影像东西向坐标。
    7.根据权利要求6所述的静止卫星成像导航与配准恒星南北参数提取方法,其特征在
    于,所述步骤4.2包括:计算出恒星影像中心处于每一列像元中心线的南北向坐标后,再结
    合探测器上像元布局几何关系得到如下的计算公式:
    <mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>x</mi> <mn>1</mn> </msub> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>x</mi> <mn>2</mn> </msub> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> <mtr> <mtd> <mtable> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> </mtable> </mtd> <mtd> <mtable> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> </mtable> </mtd> </mtr> <mtr> <mtd> <msub> <mi>x</mi> <mi>n</mi> </msub> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mi>K</mi> </mtd> </mtr> <mtr> <mtd> <mi>B</mi> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>l</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>l</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mtable> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> </mtable> </mtd> </mtr> <mtr> <mtd> <msub> <mi>l</mi> <mi>n</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
    式中:x1、x2、xn分别表示多次采样的恒星影像东西向坐标,下标表示采样次数;l1、l2、ln
    分别表示多次采样的恒星影像南北向坐标,下标表示采样次数。
    通过最小二乘法对系数进行估计,获取K、B的估计值,从而拟合出恒星影像中心在探测
    器面阵上划过的南北向坐标分量与东西向坐标分量关系。
    8.根据权利要求7所述的静止卫星成像导航与配准恒星南北参数提取方法,其特征在
    于,还包括恒星敏感东西方向参数获取步骤:据探测器面阵的几何特性,获取东西方向对称
    中心的列坐标,代入公式(5)获得恒星影像中心划过探测器东西方向中心线对应的南北向
    坐标。
    9.一种静止卫星成像导航与配准恒星南北参数提取系统,其特征在于,包括:
    列像元与恒星影像关系构建??椋河糜诠菇邢裨幕叶戎涤牒阈怯跋衲媳毕蜃攴?br />量的关系;
    灰度提取??椋河糜谔崛『阈腔邢裨行氖钡奶讲馄飨裨叶戎?;
    南北向坐标分量计算??橐唬河糜诩扑愫阈怯跋衲媳毕蜃攴至?;
    南北向坐标分量计算??槎河糜诩扑愫阈怯跋裰行幕讲馄鞫鞣较蛑行南呤钡?br />南北向坐标分量。

    说明书

    静止卫星成像导航与配准恒星南北参数提取方法及系统

    技术领域

    本发明涉及数字信号与图像处理技术,具体地,涉及一种静止卫星成像导航与配
    准恒星南北参数提取方法及系统。

    背景技术

    随着遥感卫星应用技术水平的提高,通常要求高精度静止卫星对地成像的图像被
    精确定位,并配准到地球标称网格。但卫星平台和遥感仪器受到空间力学环境、热环境变化
    等因素的影响,卫星会存在轨道漂移和姿态指向偏差,仪器会产生几何变形,影响遥感图像
    的定位与配准精度。卫星成像导航与配准技术可解决此类问题:卫星对遥感仪器的二维扫
    描机构进行实时角度增量补偿,将仪器观测点在地表的移动轨迹引导至预设路径,从而实
    现成像导航,使遥感图像配准到地球标称网格。静止遥感卫星在轨期间的仪器几何变形不
    可直接测量,可利用仪器对恒星进行敏感成像,获取恒星在探测器面阵内划过的运动轨迹,
    通过提取恒星运动规律和特征参数,与恒星在惯性系的实际位置对比可计算遥感仪器的在
    轨变形量。

    在计算恒星影像的纵向坐标时,常用的方法是重心法。但由于恒星在静止卫星探
    测器内的影像尺寸较小,甚至小于一个像元,若仍采用重心法,仅以像元的灰度为权重,不
    考虑恒星影像的形状和亮度分布特性,会引起一定的误差。

    本发明方法面向工程实际,将点扩散函数用在静止卫星成像导航与配准恒星敏感
    南北参数提取上,通过多帧数据信息处理和曲线拟合消除由探测器光学成像、电路噪声、卫
    星高频抖动等因素引起的星点位置高频误差,从而提高了恒星位置参数的识别精度。本发
    明可应用于我国第二代静止气象卫星风云四号等静止遥感卫星的成像定位与配准系统研
    制研发过程。J.L Fiorello等人编写的NOAA的研究报告(1989年,编号N90-13422)对图像导
    航配准概念和原理进行了介绍,提出了通过仪器观测恒星获取仪器变形参数的方法,但未
    阐述对图像数据的处理方法。

    经检索:

    专利号为:US4801202,名称为:Method and apparatus for radiometer star
    sensing,描述了静止气象卫星的辐射计在轨观测恒星的方法,声明了对具体操作步骤等的
    权力要求,但未阐述对恒星敏感图像数据的处理方法。

    专利号为:US4746976,名称为:Sightings by satellite for image
    navigation,介绍了静止气象卫星仪器在轨观测恒星的位置和数量,声明了对观星时机和
    方向的选取等方法的权力要求,但未阐述对多帧恒星敏感图像数据提取南北方向参数的内
    容。

    李晓坤、王淦泉、陈桂林在科学技术与工程第7卷第6期(2007年3月)的文章《风云
    四号气象卫星扫描成像仪——可见光通道恒星敏感》中,介绍了风云四号卫星扫描成像仪
    通过恒星敏感解算仪器内部形变的方法,提出通过加权分析获得恒星划过探测器的位置,
    但没有融合多列探测器的数据。

    发明内容

    针对现有技术中的缺陷,本发明的目的是提供一种静止卫星成像导航与配准恒星
    敏感南北参数提取方法。

    根据本发明提供的静止卫星成像导航与配准恒星南北参数提取方法,包括如下步
    骤:

    步骤1:构建列像元的灰度值与恒星影像南北向坐标分量的关系;

    步骤2:提取恒星划过列像元中心时的探测器像元灰度值;

    步骤3:计算恒星影像南北向坐标分量;

    步骤4:计算恒星影像中心划过探测器东西方向中心线时的南北向坐标分量。

    优选地,所述步骤1包括:

    当恒星影像划过列像元中心时,设恒星影像南北向坐标分量为l,恒星影像的中心
    点坐标写为(l,L/2),其中,L为单个像元的尺寸;设列像元的第k个像元的灰度为hk,则列像
    元与恒星影像南北向坐标分量的关系如下:


    式中:函数f(·)表示恒星影像的点扩散函数,变量x表示恒星影像东西向坐标,变
    量y表示恒星影像南北向坐标。

    优选地,所述步骤2包括:恒星划过探测器的某一列像元时,探测器多次曝光记录
    了恒星划过列像元的整个过程,通过提取恒星影像引起该列像元总灰度值随时间的变化关
    系和曲线,使用半功率法计算恒星划过该列像元中心的时刻,并记录恒星影像处于列像元
    中心时各像元的灰度值;

    具体地,选择灰度值最高的三个列像元,并记录每个像元的灰度值。

    优选地,所述步骤3包括:

    假设步骤2中选取的三个列像元的排列序号分别为k、m和n,实测的列像元灰度值
    分别为hk、hm和hn;

    由步骤1得到恒星点扩散函数推导出来的三个列像元的灰度表达式分别如下:


    恒星在探测器里的纵向坐标计算能够转化为最优化问题,最优化的目标是选取合
    适的恒星影像南北向坐标分量l,使点扩散函数法计算出的像元灰度值与实测的像元灰度
    偏差最??;目标函数如下:


    式中:c1、c2和c3为对应的三个列像元的权重,c1、c2和c3的选取与对应列像元的灰
    度值大小相关,当Q值最小时,得出的l即为恒星影像南北向坐标分量。

    优选地,所述步骤4包括:

    步骤4.1:拟合恒星影像运动轨迹;

    步骤4.2:计算恒星影像中心划过探测器东西方向中心线时的南北向坐标分量。

    优选地,所述步骤4.1包括:

    恒星影像划过探测器面阵时,会划过多列像元,则认为恒星影像中心在探测器面
    阵中的做匀速直线运动,将恒星在探测器面阵内的南北向坐标分量与东西向坐标分量描述
    为线性关系:

    l=Kx+B (4)

    式中:K表示恒星影像运动的斜率、B分别表示恒星影像南北向坐标偏量;变量x表
    示恒星影像东西向坐标;变量l表示恒星影像东西向坐标。

    优选地,所述步骤4.2包括:计算出恒星影像中心处于每一列像元中心线的南北向
    坐标后,再结合探测器上像元布局几何关系得到如下的计算公式:


    式中:x1、x2、xn分别表示多次采样的恒星影像东西向坐标,下标表示采样次数;l1、
    l2、ln分别表示多次采样的恒星影像南北向坐标,下标表示采样次数。

    通过最小二乘法对系数进行估计,获取K、B的估计值,从而拟合出恒星影像中心在
    探测器面阵上划过的南北向坐标分量与东西向坐标分量关系。

    优选地,还包括恒星敏感东西方向参数获取步骤:据探测器面阵的几何特性,获取
    东西方向对称中心的列坐标,代入公式(5)获得恒星影像中心划过探测器东西方向中心线
    对应的南北向坐标。

    根据本发明提供的静止卫星成像导航与配准恒星南北参数提取系统,包括:

    列像元与恒星影像关系构建??椋河糜诠菇邢裨幕叶戎涤牒阈怯跋衲媳毕蜃?br />标分量的关系;

    灰度提取??椋河糜谔崛『阈腔邢裨行氖钡奶讲馄飨裨叶戎?;

    南北向坐标分量计算??橐唬河糜诩扑愫阈怯跋衲媳毕蜃攴至?;

    南北向坐标分量计算??槎河糜诩扑愫阈怯跋裰行幕讲馄鞫鞣较蛑行南?br />时的南北向坐标分量。

    与现有技术相比,本发明具有如下的有益效果:

    本发明提供的静止卫星成像导航与配准恒星南北参数提取方法在计算恒星影像
    坐标时,考虑了恒星影像的形状和亮度分布特性,大幅提高了识别精度;通过多列像元计算
    结果进行最小二乘拟合,降低了由探测器光学成像、电路噪声、卫星高频抖动等因素引起的
    星点位置高频误差。

    附图说明

    通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、
    目的和优点将会变得更明显:

    图1为本发明提供的静止卫星成像导航与配准恒星南北参数提取方法的流程示意
    图;

    图2为恒星划过遥感探测器像元坐标示意图;

    图3为恒星影像和列像元的关系示意。

    具体实施方式

    下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术
    人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术
    人员来说,在不脱离本发明构思的前提下,还可以做出若干变化和改进。这些都属于本发明
    的?;し段?。

    本发明根据恒星影像的点扩散函数计算恒星影像在单列像元内的南北向坐标分
    量;利用多列像元南北向坐标分量的计算结果拟合恒星影像南北向运动轨迹,提取恒星敏
    感南北参数。根据点扩散函数对恒星影像和列像元的重叠面积进行积分,得出像元灰度与
    恒星在探测器南北向坐标分量的关系。再提取探测器像元的灰度值,进行最优化求解,得出
    星影像在单列像元内的南北向坐标分量;并根据静止遥感卫星仪器探测器面阵上各列像元
    的几何排布关系,拟合出恒星影像在探测器南北向的运动路径,并计算出恒星划过探测器
    面阵中心时的南北向坐标。

    具体地,根据本发明提供的静止卫星成像导航与配准恒星南北参数提取方法,包
    括如下步骤:

    步骤S1:构建列像元的灰度值与恒星影像南北向坐标分量的关系;

    步骤S2:提取恒星划过列像元中心时的探测器像元灰度值;

    步骤S3:计算恒星影像南北向坐标分量;

    步骤S4:计算恒星影像中心划过探测器东西方向中心线时的南北向坐标分量。

    所述步骤S1包括:

    当恒星影像划过列像元中心时(恒星影像和列像元的关系示意如附图3所示),设
    恒星影像南北向坐标分量为l,恒星影像的中心点坐标可写为(l,L/2),其中,L为单个像元
    的尺寸;设列像元的第k个列像元的灰度为hk,则列像元与恒星影像南北向坐标分量的关系
    如下:


    式中:函数f(·)表示恒星影像的点扩散函数,变量x表示恒星影像东西向坐标,变
    量y表示恒星影像南北向坐标。从上式可以看出,恒星影像的点扩散函数已知后,像元的灰
    度值与恒星影像的中心点坐标相关。

    所述步骤S2包括:

    恒星划过探测器的某一列像元时,探测器多次曝光,记录了恒星划过列像元的整
    个过程,通过提取恒星影像引起该列像元总灰度值随时间的变化关系和曲线,使用半功率
    法计算恒星划过该列像元中心的时刻,并记录恒星影像处于列像元中心时各像元的灰度
    值。

    一般情况下,恒星影像的尺寸和单个像元的尺寸大小接近,恒星影像一般会引起2
    ~3个像元的灰度值变化,因此应选择灰度值最高的三个像元(这三个像元应该是依次排
    列),并记录每个像元的灰度值。

    所述步骤S3包括:

    假设步骤S2中选取的三个列像元的排列序号分别为k、m和n,实测的列像元灰度值
    分别为hk、hm和hn。

    由步骤S1得到恒星点扩散函数推导出来的三个列像元的灰度表达式分别如下:


    恒星在探测器里的纵向坐标计算能够转化为最优化问题,最优化的目标是选取合
    适的恒星影像南北向坐标分量l,使点扩散函数法计算出的像元灰度值与实测的像元灰度
    偏差最??;目标函数可写为


    式中:c1、c2和c3为对应的三个列像元的权重,c1、c2和c3的选取与对应列像元的灰
    度值大小相关,当Q值最小时,得出的l即为恒星影像南北向坐标分量。

    所述步骤S4包括:

    恒星影像划过探测器面阵时,会划过多列像元,则认为恒星影像中心在探测器面
    阵中的做匀速直线运动,将恒星在探测器面阵内的南北向坐标分量与东西向坐标分量描述
    为线性关系:

    l=Kx+B (4)

    式中:K表示恒星影像运动的斜率、B分别表示恒星影像南北向坐标偏量;变量x表
    示恒星影像东西向坐标;变量l表示恒星影像东西向坐标。

    计算出恒星影像中心处于每一列像元中心线的南北向坐标后,再结合探测器上像
    元布局几何关系(附图2),得到如下的计算公式:


    式中:x1、x2、xn分别表示多次采样的恒星影像东西向坐标,下标表示采样次数;l1、
    l2、ln分别表示多次采样的恒星影像南北向坐标,下标表示采样次数。

    通过最小二乘法对系数进行估计,获取K、B的估计值,从而拟合出恒星影像中心在
    探测器面阵上划过的南北向坐标分量与东西向坐标分量关系。

    据探测器面阵的几何特性,获取东西方向对称中心的列坐标,代入公式(5)获得恒
    星影像中心划过探测器东西方向中心线对应的南北向坐标。

    以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述
    特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变化或修改,这并不影
    响本发明的实质内容。在不冲突的情况下,本申请的实施例和实施例中的特征可以任意相
    互组合。

    关 键 词:
    静止 卫星 成像 导航 恒星 南北 参数 提取 方法 系统
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    关于本文
    本文标题:静止卫星成像导航与配准恒星南北参数提取方法及系统.pdf
    链接地址://www.4mum.com.cn/p-6092978.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03