CN114593736A - 一种摆扫式卫星的地理定位方法、定位误差分析方法及其系统 - Google Patents

一种摆扫式卫星的地理定位方法、定位误差分析方法及其系统 Download PDF

Info

Publication number
CN114593736A
CN114593736A CN202210105023.3A CN202210105023A CN114593736A CN 114593736 A CN114593736 A CN 114593736A CN 202210105023 A CN202210105023 A CN 202210105023A CN 114593736 A CN114593736 A CN 114593736A
Authority
CN
China
Prior art keywords
coordinates
coordinate
focal plane
satellite
earth
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
CN202210105023.3A
Other languages
English (en)
Other versions
CN114593736B (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.)
Central South University
Original Assignee
Central South University
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 Central South University filed Critical Central South University
Priority to CN202410371124.4A priority Critical patent/CN118274841A/zh
Priority to CN202210105023.3A priority patent/CN114593736B/zh
Publication of CN114593736A publication Critical patent/CN114593736A/zh
Application granted granted Critical
Publication of CN114593736B publication Critical patent/CN114593736B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/20—Instruments for performing navigational calculations
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C11/00—Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
    • G01C11/04—Interpretation of pictures
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Multimedia (AREA)
  • Computing Systems (AREA)
  • Algebra (AREA)
  • Manufacturing & Machinery (AREA)
  • Automation & Control Theory (AREA)
  • Image Processing (AREA)

Abstract

本发明公开了一种摆扫式卫星的地理定位方法、定位误差分析方法及其系统，本发明的定位方法是利用摆扫式卫星的成像几何模型计算出每个目标像点的光轴指向，进而基于光轴指向与地球椭圆相交以获取目标像点的对应地面点坐标，实现定位目的。本发明的定位误差分析方法是利用检查点的像方坐标恢复成像光轴指向，并利用检查点的物方坐标及摆扫式卫星的成像几何模型逆推出相机坐标系下光轴指向，最后将两个光轴指向转换为焦平面上齐次坐标差值以分析定位误差。本发明提供的定位技术，其基于严密成像几何模型提高了定位精度以及可靠性；定位误差分析技术是以一种全新的手段来分析摆扫式卫星的定位误差。

Description

一种摆扫式卫星的地理定位方法、定位误差分析方法及其 系统
技术领域
本发明属于卫星影像处理技术领域，具体涉及一种摆扫式卫星的地理定位方法、定位误差分析方法及其系统。
背景技术
摆扫式卫星影像的地理定位是卫星影像处理的基础环节，该环节将会为摆扫式卫星影像生成相应的地理定位查找表，即建立摆扫式卫星影像与地面坐标间的一一对应关系。由于摆扫式卫星通常采用多元并扫的成像机制，需要采用大量的活动部件，以获取大的视场角。45°镜是一种常用的摆扫式成像方案，但该机械结构将会使得多元并扫影像产生显著的像旋效应，因此需要进一步引入其他光学部件如K镜等实现对影像像旋的消除。然而，由于卫星发射时的巨大冲力及在轨后运行环境的差异，使得各个部件之间的相对关系发生变化。这些变化使得成像光线不再符合理想摆扫模型，制约了摆扫式卫星地理定位的精度。摆扫式卫星影像的地理定位查找表是通过构建成像几何模型，与地球椭球相交以获取像点的对应地面点坐标。不准确的成像几何模型将导致定位结果中存在定位误差，因此，构建精确的成像几何模型是摆扫式卫星的地理定位的关键性环节。
现有两种方法可以分析摆扫式相机的地理定位误差，分别为：在物方空间和像方空间。在第一种情况下，摆扫式影像将在特定投影中进行校正，从更高精度的参考影像中提取高精度地面控制点，在特定的投影中确定绝对地理定位精度。但摆扫式相机的地面采样间隔(GSD)随着卫星和物体之间的距离而变化。因此，像方空间中的地理定位误差取决于视角。此外，考虑到地图投影的失真，很难为超过2900公里的条带选择合适的地图投影方法。另一种是基于像方的定位误差分析方法，即将物方坐标重新投影到像方空间来确定像平面的误差。但地理定位查找表建立的像方到物方的映射过程，由于物方坐标反算像方坐标过程较为复杂，实现较为困难。
针对上述成像几何模型的精度将影响摆扫式卫星的地理定位精度的问题，本发明提供一种摆扫式卫星的地理定位方法，此外，针对上述摆扫式相机的地理定位误差分析问题，本发明提供一种定位误差分析方法予以解决。
发明内容
本发明的目的是针对上述现有技术中至少存在的部分技术问题，提供一种摆扫式卫星的地理定位方法、定位误差分析方法及其系统。本发明提供的地理定位方法，其基于严密的成像几何模型进行定位，提高了定位结果的精度以及可靠性；本发明提供的定位误差分析方法，其引入齐次坐标，提供了一种全新的误差分析技术，且分析过程也利用了严密的成像几何模型，进一步保证了误差分析结果的可靠性。
一方面，本发明提供的一种摆扫式卫星的地理定位方法，其包括以下步骤：
步骤1：构建摆扫式卫星的成像几何模型，所述成像几何模型表示为如下公式或其等价变形公式：
Figure BDA0003492104450000021
或
Figure BDA0003492104450000022
式中，坐标(x，y)表示像点在焦平面坐标系的坐标，f为主距，
Figure BDA0003492104450000023
为卫星本体坐标到轨道坐标的变换矩阵，
Figure BDA0003492104450000024
为轨道坐标到地心惯性系ECI坐标的变换矩阵，
Figure BDA0003492104450000025
为卫星本体坐标到地心惯性系ECI坐标的变换矩阵，
Figure BDA0003492104450000026
为地心惯性系ECI坐标到地心地固坐标系ECF的变换矩阵，Rmkt为扫描组件的变换矩阵，m为比例因子，[Xs Ys Zs]T表示卫星在地心地固坐标系下的位置，T为矩阵转置符号，
Figure BDA0003492104450000027
表示地心地固坐标系中物方坐标；
步骤2：利用所述成像几何模型计算出每个目标像点的光轴指向；
步骤3：基于光轴指向与地球椭球或数字高程模型相交以获取目标像点的对应地面点物方坐标。
进一步可选地，所述像点恢复的成像光线，表示为：
Figure BDA0003492104450000028
式中，
Figure BDA0003492104450000029
表示相机坐标系下像点的光轴指向，rN为像点所在帧的探元标号，μ表示探元的像元大小，f为主距大小，(x0，y0)是焦平面CCD(电荷耦合器件)的中心点坐标，M为探元数目。
进一步可选地，所述扫描组件包括K镜和扫描镜，所述扫描组件的变换矩阵Rmkt表示为如下公式或其等价变形公式：
Figure BDA00034921044500000210
式中，Rm(θ)表示扫描镜绕X轴旋转θ角度后的反射矩阵，Rk表示K镜的反射矩阵，
Figure BDA00034921044500000211
表示相机坐标系光轴指向变换到镜头坐标系的变换矩阵。
进一步可选地，若所述K镜具有三个反射镜，对应的安装误差用绕X轴，Y轴，Z轴三个角度φ-α-γ分别表示，所述K镜的反射矩阵表示为如下公式或其等价变形公式：
Figure BDA0003492104450000031
或
Figure BDA0003492104450000032
且满足：
Rki＝Rφ·Rα·Rγ
式中，Rθ/2+φ、Rθ/2均表示K镜的反射矩阵，Rki为旋转矩阵，Rφ，Rα，Rγ分别表示φ-α-γ方向对应的旋转矩阵，Rk0为扫描镜指向天底时K镜的反射矩阵。
进一步可选地，当像点扫描坐标c位于图像的中心时，视为扫描角θ等于0，当匀速扫描时，扫描角θ表示为：
Figure BDA0003492104450000033
当非匀速扫描时，扫描角θ表示为：
Figure BDA0003492104450000034
式中，W是图像的宽度，t0是采样间隔，
Figure BDA0003492104450000035
为扫描镜的角速度，t表示帧内的扫描时间，ai，wi，和φi(i＝1，2，...，n)分别是振幅、角频率和相位，i表示正弦曲线编号。
进一步可选地，所述方法还包括利用步骤1-步骤3的方式以预设间隔进行像点定位，基于像点及其对应的地面点坐标构建地理定位查找表。
第二方面，本发明提供一种摆扫式卫星的地理定位误差分析方法，其包括以下步骤：
利用检查点的像点坐标计算出对应的焦平面齐次坐标，所述焦平面齐次坐标与相机坐标系下的光轴指向的关系满足：
Figure BDA0003492104450000036
其中，
Figure BDA0003492104450000037
表示相机坐标系下的光轴指向，定义为：
Figure BDA0003492104450000041
其中，坐标(x，y)表示像点在焦平面坐标系的坐标，f为主距；
利用检查点的物方坐标以及摆扫式卫星的成像几何模型逆转换得到焦平面系统中的坐标，再计算出对应的焦平面齐次坐标，其中，按照如下组公式或者其等价公式进行转换得到焦平面系统中的坐标：
Figure BDA0003492104450000042
或者
Figure BDA0003492104450000043
式中，[Xs Ys Zs]T表示卫星在地心地固坐标系下的位置，T为矩阵转置符号，
Figure BDA0003492104450000044
表示地心地固坐标系中的物方坐标，
Figure BDA0003492104450000045
表示本体坐标系下的光轴指向，Rmkt为扫描组件的变换矩阵，
Figure BDA0003492104450000046
为地心惯性系ECI坐标到卫星本体坐标的变换矩阵，
Figure BDA0003492104450000047
为地心地固ECF坐标到地心惯性系ECI坐标的变换矩阵，
Figure BDA0003492104450000048
为轨道坐标到卫星本体坐标的变换矩阵，
Figure BDA0003492104450000049
为地心惯性系ECI坐标到轨道坐标的变换矩阵；
最后利用基于像点坐标对应的焦平面齐次坐标与基于物方坐标对应的焦平面齐次坐标的差值分析定位误差。
第三方面，本发明提供一种基于所述地理定位方法或基于所述地理定位误差分析方法的摆扫式卫星系统，其设有扫描镜、以及设置于扫描镜与望远镜之间的K镜，经过K镜反射作用，光线被望远镜聚焦并投射到焦平面上。
第四方面，本发明提供一种摆扫式卫星的地理定位，其包括：
成像几何模型构建模块，用于构建摆扫式卫星的成像几何模型，所述成像几何模型表示为如下公式或其等价变形公式：
Figure BDA00034921044500000410
或
Figure BDA00034921044500000411
式中，坐标(x，y)表示像点在焦平面坐标系的坐标，f为主距，
Figure BDA0003492104450000051
为卫星本体坐标到轨道坐标的变换矩阵，
Figure BDA0003492104450000052
为轨道坐标到地心惯性系ECI坐标的变换矩阵，
Figure BDA0003492104450000053
为卫星本体坐标到地心惯性系ECI坐标的变换矩阵，
Figure BDA0003492104450000054
为地心惯性系ECI坐标到地心地固坐标系ECF的变换矩阵，Rmkt为扫描组件的变换矩阵，m为比例因子，[Xs Ys Zs]T表示卫星在地心地固坐标系下的位置，T为矩阵转置符号，
Figure BDA0003492104450000055
表示地心地固坐标系中物方坐标；
目标像点的光轴指向计算模块，利用所述成像几何模型计算出每个目标像点的光轴指向；
定位模块，用于基于光轴指向与地球椭球或数字高程模型相交以获取目标像点的对应地面点物方坐标。
第五方面，本发明提供一种基于所述地理定位误差分析方法的系统，其包括：
基于像点坐标的齐次坐标计算模块，用于利用检查点的像点坐标计算出对应的焦平面齐次坐标；
基于物方坐标的齐次坐标计算模块，用于利用检查点的物方坐标以及摆扫式卫星的成像几何模型逆转换得到焦平面系统中的坐标，再计算出对应的焦平面齐次坐标；
误差分析模块，利用基于像点坐标对应的焦平面齐次坐标与基于物方坐标对应的焦平面齐次坐标的差值分析定位误差。
第六方面，本发明提供一种可读存储介质，其存储了计算机程序，所述计算机程序被处理器调用以实现：
一种摆扫式卫星的地理定位方法的步骤；
或者所述计算机程序被处理器调用以实现：
一种摆扫式卫星的地理定位误差分析方法的步骤。
有益效果
1.本发明提供的地理定位方法构建了严密的摆扫式卫星的成像几何模型，得以保证摆扫式卫星的地理定位的精度，提高定位结果的可靠性。
2.本发明提供的摆扫式卫星的地理定位误差分析方法，其引入了焦平面齐次坐标，并利用检查点的像方坐标计算得到的齐次坐标与物方坐标计算出的焦平面坐标进行误差评估，提供了一种全新的全局定位误差评估的思路以及技术手段，可以更直观的揭示出摆扫式相机误差来源。
附图说明
图1是本发明实施例提供的摆扫式卫星的成像示意图；
图2是本发明提供的传统定位模型的定位误差分布图；
图3是本发明提供的扫描方向的残差中值的示意图；
图4是采用本发明所述地理定位方法实现时的误差分布图。
具体实施方式
本发明提供的一种摆扫式卫星的地理定位方法，是针对扫描组件精确建模，进而构建出严密的摆扫式卫星的成像几何模型，最终基于摆扫式卫星的成像几何模型得到定位结果。下述实施例中将以设有三个反射镜的K镜和45°扫描镜为例，其他可行的实施例中，在不脱离本发明构思的基础上，对K镜和扫描镜进行结构调整也是可行的。本发明提供的一种摆扫式卫星的地理定位误差分析方法，是一种用于评估定位结果误差的技术手段，其引入了焦平面齐次坐标，并利用检查点的像方坐标计算得到的齐次坐标与物方坐标计算出的焦平面坐标进行误差评估，提供了一种全新的定位误差评估的思路以及技术手段。
实施例1：
本实施例以设有三个反射镜的K镜和45°扫描镜为例。如图1所示，MERSI2使用45°扫描镜来捕获垂轨向的地面图像。但是，对于具有多探元的传感器，反射镜会导致图像旋转错位。为了纠正这种图像旋转，采用了具有三个反射镜的K镜。经过多次反射后，光线被望远镜聚焦并投射到焦平面上。为了分析成像几何，引入了仪器框架坐标系，其X轴为45°镜扫描的旋转轴，Z轴指向地心，Y轴由右手定则确定。线性探元垂直于焦平面中的扫描方向。在焦平面坐标系中，x轴与飞行方向相反，y轴由右手定则确定。
基于上述结构，本实施例中的摆扫式卫星的成像几何模型的推理过程如下：
1.构建摆扫式相机镜头坐标系下光轴指向，具体参照公式(3)。
通过45°扫描镜来捕获垂轨向的地面图像影像，针对像点(r，c)，r表示像点所在图像坐标中的行数，c为像点扫描方向坐标。利用行数r和传感器探元数量M的商的整数部分表示像点(r，c)对应的帧N，行数r和传感器探元数量M的商的余数部分表示像点所在帧的探元标号rN。
进而根据像点所在帧的探元标号rN确定对应探元在焦平面上的位置，利用探元的安装和像元大小，恢复像素(r，c)在焦平面坐标系中的坐标(x，y)为：
Figure BDA0003492104450000061
式中，μ表示探元的像元大小，(x0，y0)是焦平面CCD(电荷耦合器件)的中心点坐标。
再利用二维像点和相机主距f恢复相机坐标系下的光轴指向
Figure BDA0003492104450000071
具体表示为：
Figure BDA0003492104450000072
关于相机坐标系的光轴指向与镜头坐标系之间的变换矩阵
Figure BDA0003492104450000073
其可以将相机坐标系的光轴指向变换到镜头坐标系中，变换矩阵
Figure BDA0003492104450000074
表示为：
Figure BDA0003492104450000075
2.构建摆扫式相机扫描组件模型的变换矩阵，具体参照公式(18)。
由于K镜和45°扫描镜的旋转，扫描角度取决于扫描时间t。1km分辨率的探元以t0的间隔进行采样，而250m分辨率的采样频率是1km分辨率的四倍。其中，扫描时间t的定义为：
Figure BDA0003492104450000076
式中，W是影像宽度，t0为采样间隔。
45°扫描镜以角速度
Figure BDA0003492104450000079
扫描地球，当像点扫描坐标c位于图像的中心时，定义扫描角θ等于0。本发明考虑匀速扫描和非匀速扫描，其中，匀速扫描时，扫描角定义为：
Figure BDA0003492104450000077
当扫描镜非匀速扫描时，扫描角定义为：
Figure BDA0003492104450000078
式中，ai，ωi，和φi(i＝1，2，...，n)分别是振幅，角频率和相位。应当说明的是，本发明分为匀速扫描和非匀速扫描，是对扫描过程的进一步优化，可以进一步保证构建的成像几何模型的精度，但是其并非本发明的唯一实现方式，在不脱离本发明构思的基础，即构建基于扫描组件的成像几何模型的基础上，是否对扫描过程进行优化会影响到模型精度，但是不影响本发明的技术构思的实现。
为了构建摆扫式相机扫描组件模型的变换矩阵，首先考虑构建45°扫描镜的转动模型，即反射矩阵Rm(θ)，过程为：
反射光线r′由法向量
Figure BDA0003492104450000081
和入射光线r决定。根据反射定律，反射光线r′存在：
r′＝r-2n·(nT·r)＝(I-2n·nT)·r＝Rm·r (7)
反射矩阵Rm表示为：
Figure BDA0003492104450000082
考虑到45°扫描镜绕X轴旋转θ，则旋转后的法向量为：
Figure BDA0003492104450000083
对应的反射矩阵Rm(θ)表示为：
Rm(θ)＝R(θ)·Rm·RT(θ) (10)
在无摆动成像时，45°扫描镜的法向量为：
Figure BDA0003492104450000084
因此，仪器坐标下45°扫描镜的反射矩阵为：
Figure BDA0003492104450000085
从公式(10)以及公式(11)来看，仪器坐标下45°扫描镜的反射矩阵存在两种形式。
关于K镜，K镜的旋转反射矩阵Rk(θ/2)由三个旋转反射矩阵组成，具体可以表示为：
Figure BDA0003492104450000086
当45°扫描镜指向天底时，三个反射面的法向量分别为[cosβ 0 sinβ]T、[0 0 -1]T和[-cosβ 0 sinβ]T，K镜的反射矩阵为：
Figure BDA0003492104450000087
根据式(13)，K镜的反射矩阵与角度β无关。由式(12)可知，K镜的旋转反射矩阵为：
Figure BDA0003492104450000091
K镜的安装误差可以用绕X轴，Y轴，Z轴三个角度φ-α-γ表示，旋转矩阵Rki定义如下：
Figure BDA0003492104450000092
三个反射镜可以作为一个单元进行变换。因此，K镜的反射矩阵变为：
Figure BDA0003492104450000093
或者将其简化为：
Figure BDA0003492104450000094
式中，Rθ/2+φ、Rθ/2均表示K镜的反射矩阵，其为本领域的常规设定。本实施例结合K镜和45°扫描镜，扫描组件的变换矩阵为：
Figure BDA0003492104450000095
3.基于摆扫式相机镜头坐标系下光轴指向、扫描组件模型的变换矩阵构建成像几何模型。
其中，由帧N内插帧的成像时间tN，结合帧内成像时间t，计算每个像点的时间t’。
为了实现严密的成像几何模型的构建，将相机坐标系转换到本体坐标系中的光轴指向表示为：
Figure BDA0003492104450000096
即从上述公式可知，利用扫描组件的变换矩阵，将像素在相机坐标系转换到本体坐标系中的光轴指向。此外，相机的姿态定义了从本体坐标系到轨道坐标系的变换。在Y-Z-X旋转下，用三个欧拉角
Figure BDA0003492104450000097
w，κ确定旋转矩阵，具体表示为：
Figure BDA0003492104450000101
利用地心惯性(ECI)坐标中的位置和速度，可以计算从轨道坐标到ECI坐标的变换矩阵
Figure BDA0003492104450000102
ECI中的光轴指向可以通过IERS约定转换为地心地固(ECF)坐标。通过上述变换，本发明构建的成像几何模型表示为如下公式或其等价变形公式：
Figure BDA0003492104450000103
或
Figure BDA0003492104450000104
式中，坐标(x，y)表示像点在焦平面坐标系的坐标，f为主距，
Figure BDA0003492104450000105
为卫星本体坐标到轨道坐标的变换矩阵，
Figure BDA0003492104450000106
为轨道坐标到ECI坐标的变换矩阵，
Figure BDA0003492104450000107
为卫星本体坐标到ECI坐标的变换矩阵，
Figure BDA0003492104450000108
为ECI坐标到地心地固ECF的变换矩阵，Rmkt为扫描组件的变换矩阵，m为比例因子，[Xs Ys Zs]T表示卫星的位置，T为矩阵转置符号，
Figure BDA0003492104450000109
表示地心地固坐标系中物方坐标。
应当理解，本发明涉及的公式是基于几何的坐标转换，因此，公式的表达并未唯一表达式，可以对其进行等价变形，在不脱离本发明构思的基础上，对其进行等价替换均属于本发明的保护范围。
基于上述推理的成像几何模型，本发明提供的一种摆扫式卫星的地理定位方法，其包括以下步骤：
步骤1：构建摆扫式卫星的成像几何模型。
步骤2：利用所述成像几何模型计算出每个目标像点的光轴指向。其中，将目标像点的参数代入所述成像几何模型得到目标像点的光轴指向。
步骤3：基于光轴指向与地球椭球或数字高程模型相交以获取目标像点的对应地面点坐标。
其中，恢复的光线与地球表面相交以确定比例因子m，并计算ECF中的笛卡尔坐标。由于视角不同，地形起伏会给图像带来差异。因此，采用7.5弧秒空间分辨率全球数字高程模型(DEM)，GMTED2010。为了避免发散和遮挡引起的问题，选择光线追踪法。笛卡尔坐标到大地坐标的变换表示为四次方程的解。为了获得更好的计算效率和准确性，应用了迭代方法。由于“基于光轴指向与地球椭球或数字高程模型相交以获取目标像点的对应地面点坐标”的实现过程是现有技术可以实现的，且本发明对其并未进行优化，因此，对其具体过程不进行详细陈述。
本发明为了星载摆扫式相机建立地理定位查找表，进而按照上述方法以一定间隔进行像点坐标计算。由于地理定位查找表的构建规则以及形式均是现有的，故对其不进行具体的陈述。
实施例2：
本实施例提供一种摆扫式卫星的地理定位误差分析方法，其包括以下步骤：
S1：利用检查点的像点坐标计算出对应的焦平面齐次坐标。其中，每个像素定义了光轴指向，它有一个自由度。为分析MERSI2的系统误差，在焦平面上引入齐次坐标(px py f/μ)，其定义为：
Figure BDA0003492104450000111
其中，(px，py)的单位是像素，1个像素等于1个IFOV。考虑到扫描方向的驻留角较小，py的空间分辨率在实像空间中较大。此外，px在飞行方向，py在扫描方向。即按照公式(2)以及(22)计算出像点对应的焦平面齐次坐标。
S2：利用检查点的物方坐标以及摆扫式卫星的成像几何模型逆转换得到焦平面系统中的坐标，再计算出对应的焦平面齐次坐标，其中，按照如下组公式或者其等价公式进行转换得到焦平面系统中的坐标，再利用公式(22)计算出物方坐标对应的焦平面齐次坐标：
Figure BDA0003492104450000112
Figure BDA0003492104450000121
或者
Figure BDA0003492104450000122
Figure BDA0003492104450000123
式中，[Xs Ys Zs]T表示卫星的位置，T为矩阵转置符号，
Figure BDA0003492104450000124
表示地心地固坐标系中像素的光轴指向，
Figure BDA0003492104450000125
表示本体坐标系下的光轴指向，
Figure BDA0003492104450000126
表示相机坐标系下的光轴指向，Rmkt为扫描组件的变换矩阵，
Figure BDA0003492104450000127
为地心惯性系ECI坐标到卫星本体坐标的变换矩阵，
Figure BDA0003492104450000128
为地心地固ECF到ECI坐标的变换矩阵，
Figure BDA0003492104450000129
为ECI坐标到轨道坐标的变换矩阵，
Figure BDA00034921044500001210
为轨道坐标到卫星本体坐标的变换矩阵。
最后利用基于像点坐标对应的焦平面齐次坐标与基于物方坐标对应的焦平面齐次坐标的差值分析定位误差。
Figure BDA00034921044500001211
公式(27)等价变形式的一种形式为：
Figure BDA00034921044500001212
其中，Px1、py1为基于像点坐标对应的焦平面齐次坐标；px2、py2为基于物方坐标对应的焦平面齐次坐标。Δpx，Δpy为两个焦平面齐次坐标的差值，应当理解，Δpx，Δpy的数值大小反应了误差大小，误差大小程度与Δpx，Δpy的数值大小的关系是可以根据实际需要以及实验精度为依据进行调整和设置。XC、YC、ZC可以视为是像点对应在相机坐标系下的光轴指向，x、y可以视为是物方坐标经过摆扫式卫星的成像几何模型逆转换得到焦平面作差。本实施例中，在焦平面的齐次坐标下，像点主点误差表现为扫描方向的平移量，主距误差表现为帧内的缩放，转角误差直接表现为扫描方向的非线性误差，如图3所示，姿态滚动角误差为平移误差，而俯仰角误差引起飞行方向类似余弦曲线误差，偏航角误差引起飞行方向正弦曲线误差。
验证：
精度验证采用的FY-3D MERSI2数据可从国家卫星气象中心免费获得。采用了数据集FY3D_MERSI_GBAL_L1_20200517_0505_0250M_MS。图2为传统定位模型的定位误差分布图，扫描方向的残差中值如图3所示，从图中可以得到看出原始成像几何模型中存在显著的姿态误差与扫描误差。采用本发明的定位方法，其误差分布图如图4，相较于图2，采用本发明所述方法的RMSE减少到0.32像素。
实施例3：
本实施例提供一种基于所述地理定位方法或者基于所述地理定位误差分析方法的摆扫式卫星系统，其内设有扫描镜、以及设置于扫描镜与望远镜之间的K镜，经过K镜反射作用，光线被望远镜聚焦并投射到焦平面上。具体的结构可以参照图1。
实施例4：
本实施例提供一种基于所述地理定位方法的系统，其包括：
成像几何模型构建模块，用于构建摆扫式卫星的成像几何模型；
目标像点的光轴指向计算模块，利用所述成像几何模型计算出每个目标像点的光轴指向；
定位模块，用于基于光轴指向与地球椭球或数字高程模型相交以获取目标像点的对应地面点物方坐标。
或者提供一种基于所述地理定位误差分析方法的系统，其包括：
基于像点坐标的齐次坐标计算模块，用于利用检查点的像点坐标计算出对应的焦平面齐次坐标；
基于物方坐标的齐次坐标计算模块，用于利用检查点的物方坐标以及摆扫式卫星的成像几何模型逆转换得到焦平面系统中的坐标，再计算出对应的焦平面齐次坐标；
误差分析模块，利用基于像点坐标对应的焦平面齐次坐标与基于物方坐标对应的焦平面齐次坐标的差值分析定位误差。
在一些实现方式中，基于所述地理定位方法的系统和基于所述地理定位误差分析方法的系统可以是同一个系统，即实现定位的基础上还实现了定位误差分析，实现不同功能的功率模块可以共用一个模块，即一个功能模块可以实现上述两种或更多的功能；也可以是每个功能对应一个功能模块，本发明对此不进行具体的限定。
其中，各个功能模块单元的实现过程可以参照对应方法的描述。
上述功能模块单元的划分仅仅是一种逻辑功能的划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。同时，上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。
实施例5：
本实施例提供一种可读存储介质，其存储了计算机程序，所述计算机程序被处理器调用以实现：一种摆扫式卫星的地理定位方法的步骤，具体是执行：
步骤1：构建摆扫式卫星的成像几何模型。
步骤2：利用所述成像几何模型计算出每个目标像点的光轴指向。
步骤3：基于光轴指向与地球椭球或数字高程模型相交以获取目标像点的对应地面点坐标。
在一些实现方式中，计算机程序被处理器调用以实现：利用步骤1-步骤3的方式以预设间隔进行像点定位，基于像点及其对应的地面点坐标构建地理定位查找表。
或者，所述计算机程序被处理器调用以实现：一种摆扫式卫星的地理定位误差分析方法的步骤。具体是执行：
利用检查点的像点坐标计算出对应的焦平面齐次坐标，所述焦平面齐次坐标与相机坐标系下的光轴指向的关系满足：
利用检查点的物方坐标以及摆扫式卫星的成像几何模型逆转换得到焦平面系统中的坐标，再计算出对应的焦平面齐次坐标；最后利用基于像点坐标对应的焦平面齐次坐标与基于物方坐标对应的焦平面齐次坐标的差值分析定位误差。
各个步骤的具体实现过程请参照前述方法的阐述。
所述可读存储介质为计算机可读存储介质，其可以是前述任一实施例所述的控制器的内部存储单元，例如控制器的硬盘或内存。所述可读存储介质也可以是所述控制器的外部存储设备，例如所述控制器上配备的插接式硬盘，智能存储卡(Smart Media Card，SMC)，安全数字(Secure Digital，SD)卡，闪存卡(Flash Card)等。进一步地，所述可读存储介质还可以既包括所述控制器的内部存储单元也包括外部存储设备。所述可读存储介质用于存储所述计算机程序以及所述控制器所需的其他程序和数据。所述可读存储介质还可以用于暂时地存储已经输出或者将要输出的数据。
基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分，或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的可读存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
实施例6：
本实施例提供一种电子终端，其包括：至少一个或多个处理器以及一个或多个存储器，所述存储器存储了计算机程序，计算机程序被处理器调用以实现：一种摆扫式卫星的地理定位方法的步骤，具体是执行：
步骤1：构建摆扫式卫星的成像几何模型。
步骤2：利用所述成像几何模型计算出每个目标像点的光轴指向。
步骤3：基于光轴指向与地球椭球或数字高程模型相交以获取目标像点的对应地面点坐标。
在一些实现方式中，计算机程序被处理器调用以实现：利用步骤1-步骤3的方式以预设间隔进行像点定位，基于像点及其对应的地面点坐标构建地理定位查找表。
或者计算机程序被处理器调用以实现：一种摆扫式卫星的地理定位误差分析方法的步骤。具体是执行：
利用检查点的像点坐标计算出对应的焦平面齐次坐标，所述焦平面齐次坐标与相机坐标系下的光轴指向的关系满足：
利用检查点的物方坐标以及摆扫式卫星的成像几何模型逆转换得到焦平面系统中的坐标，再计算出对应的焦平面齐次坐标；最后利用基于像点坐标对应的焦平面齐次坐标与基于物方坐标对应的焦平面齐次坐标的差值分析定位误差。
该电子终端还包括：通信接口，用于与外界设备进行通信，进行数据交互传输。
其中，存储器可能包含高速RAM存储器，也可能还包括非易失性除颤器，例如至少一个磁盘存储器。
如果存储器、处理器和通信接口独立实现，则存储器、处理器和通信接口可以通过总线相互连接并完成相互间的通信。所述总线可以是工业标准体系结构总线，外部设备互联总线或扩展工业标准体系结构总线等。所述总线可以分为地址总线、数据总线、控制总线等。
可选的，在具体实现上，如果存储器、处理器和通信接口集成在一块芯片上，则存储器、处理器即通信接口可以通过内部接口完成相互之间的通信。
各个步骤的具体实现过程请参照前述方法的阐述。
应当理解，在本发明实施例中，所称处理器可以是中央处理单元(CentralProcessing Unit，CPU)，该处理器还可以是其他通用处理器、数字信号处理器(DigitalSignal Processor，DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现成可编程门阵列(Field-Programmable GateArray，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。存储器可以包括只读存储器和随机存取存储器，并向处理器提供指令和数据。存储器的一部分还可以包括非易失性随机存取存储器。例如，存储器还可以存储设备类型的信息。
需要强调的是，本发明所述的实例是说明性的，而不是限定性的，因此本发明不限于具体实施方式中所述的实例，凡是由本领域技术人员根据本发明的技术方案得出的其他实施方式，不脱离本发明宗旨和范围的，不论是修改还是替换，同样属于本发明的保护范围。

Claims (10)

1.一种摆扫式卫星的地理定位方法，其特征在于：包括以下步骤：
步骤1：构建摆扫式卫星的成像几何模型，所述成像几何模型表示为如下公式或其等价变形公式：
Figure FDA0003492104440000011
或
Figure FDA0003492104440000012
式中，坐标(x，y)表示像点在焦平面坐标系的坐标，f为主距，
Figure FDA0003492104440000013
为卫星本体坐标到轨道坐标的变换矩阵，
Figure FDA0003492104440000014
为轨道坐标到地心惯性系ECI坐标的变换矩阵，
Figure FDA0003492104440000015
为卫星本体坐标到地心惯性系ECI坐标的变换矩阵，
Figure FDA0003492104440000016
为地心惯性系ECI坐标到地心地固坐标系ECF的变换矩阵，Rmkt为扫描组件的变换矩阵，m为比例因子，[Xs Ys Zs]T表示卫星在地心地固坐标系下的位置，T为矩阵转置符号，
Figure FDA0003492104440000017
表示地心地固坐标系中物方坐标；
步骤2：利用所述成像几何模型计算出每个目标像点的光轴指向；
步骤3：基于光轴指向与地球椭球或数字高程模型相交以获取目标像点的对应地面点物方坐标。
2.根据权利要求1所述的地理定位方法，其特征在于：像点恢复的成像光线，表示为：
Figure FDA0003492104440000018
式中，
Figure FDA0003492104440000019
表示相机坐标系下像点的光轴指向，rN为像点所在帧的探元标号，μ表示探元的像元大小，f为主距大小，(x0，y0)是焦平面CCD(电荷耦合器件)的中心点坐标，M为探元数目。
3.根据权利要求1所述的地理定位方法，其特征在于：所述扫描组件包括K镜和扫描镜，所述扫描组件的变换矩阵Rmkt表示为如下公式或其等价变形公式：
Figure FDA00034921044400000110
式中，Rm(θ)表示扫描镜绕X轴旋转θ角度后的反射矩阵，Rk表示K镜的反射矩阵，
Figure FDA00034921044400000111
表示相机坐标系光轴指向变换到镜头坐标系的变换矩阵。
4.根据权利要求3所述的地理定位方法，其特征在于：若所述K镜具有三个反射镜，对应的安装误差用绕X轴，Y轴，Z轴三个角度φ-α-γ分别表示，所述K镜的反射矩阵表示为如下公式或其等价变形公式：
Figure FDA0003492104440000021
或
Figure FDA0003492104440000022
且满足：
Rki＝Rφ·Rα·Rγ
式中，Rθ/2+φ、Rθ/2均表示K镜的反射矩阵，Rki为旋转矩阵，Rφ，Rα，Rγ分别表示φ-α-γ方向对应的旋转矩阵，Rk0为扫描镜指向天底时K镜的反射矩阵。
5.根据权利要求3所述的地理定位方法，其特征在于：当像点扫描坐标c位于图像的中心时，视为扫描角θ等于0，当匀速扫描时，扫描角θ表示为：
Figure FDA0003492104440000023
当非匀速扫描时，扫描角θ表示为：
Figure FDA0003492104440000024
式中，W是图像的宽度，t0是采样间隔，
Figure FDA0003492104440000025
为扫描镜的角速度，t表示帧内的扫描时间，ai，wi，和φi(i＝1，2，...，n)分别是振幅、角频率和相位，i表示正弦曲线编号。
6.根据权利要求1所述的方法，其特征在于：还包括利用步骤1-步骤3的方式以预设间隔进行像点定位，基于像点及其对应的地面点坐标构建地理定位查找表。
7.一种摆扫式卫星的地理定位误差分析方法，其特征在于：包括以下步骤：
利用检查点的像点坐标计算出对应的焦平面齐次坐标，所述焦平面齐次坐标与相机坐标系下的光轴指向的关系满足：
Figure FDA0003492104440000026
其中，
Figure FDA0003492104440000027
表示相机坐标系下的光轴指向，定义为：
Figure FDA0003492104440000031
其中，坐标(x，y)表示像点在焦平面坐标系的坐标，f为主距；
利用检查点的物方坐标以及摆扫式卫星的成像几何模型逆转换得到焦平面系统中的坐标，再计算出对应的焦平面齐次坐标，其中，按照如下组公式或者其等价公式进行转换得到焦平面系统中的坐标：
Figure FDA0003492104440000032
或者
Figure FDA0003492104440000033
式中，[Xs Ys Zs]T表示卫星在地心地固坐标系下的位置，T为矩阵转置符号，
Figure FDA0003492104440000034
表示地心地固坐标系中的物方坐标，
Figure FDA0003492104440000035
表示本体坐标系下的光轴指向，Rmkt为扫描组件的变换矩阵，
Figure FDA0003492104440000036
为地心惯性系ECI坐标到卫星本体坐标的变换矩阵，
Figure FDA0003492104440000037
为地心地固ECF坐标到地心惯性系ECI坐标的变换矩阵，
Figure FDA0003492104440000038
为轨道坐标到卫星本体坐标的变换矩阵，
Figure FDA0003492104440000039
为地心惯性系ECI坐标到轨道坐标的变换矩阵；
最后利用基于像点坐标对应的焦平面齐次坐标与基于物方坐标对应的焦平面齐次坐标的差值分析定位误差。
8.一种摆扫式卫星的地理定位系统，其特征在于：包括：
成像几何模型构建模块，用于构建摆扫式卫星的成像几何模型，所述成像几何模型表示为如下公式或其等价变形公式：
Figure FDA00034921044400000310
或
Figure FDA00034921044400000311
式中，坐标(x，y)表示像点在焦平面坐标系的坐标，f为主距，
Figure FDA00034921044400000312
为卫星本体坐标到轨道坐标的变换矩阵，
Figure FDA00034921044400000313
为轨道坐标到地心惯性系ECI坐标的变换矩阵，
Figure FDA00034921044400000314
为卫星本体坐标到地心惯性系ECI坐标的变换矩阵，
Figure FDA00034921044400000315
为地心惯性系ECI坐标到地心地固坐标系ECF的变换矩阵，Rmkt为扫描组件的变换矩阵，m为比例因子，[Xs Ys Zs]T表示卫星在地心地固坐标系下的位置，T为矩阵转置符号，
Figure FDA0003492104440000041
表示地心地固坐标系中物方坐标；
目标像点的光轴指向计算模块，利用所述成像几何模型计算出每个目标像点的光轴指向；
定位模块，用于基于光轴指向与地球椭球或数字高程模型相交以获取目标像点的对应地面点物方坐标。
9.一种基于权利要求7所述地理定位误差分析方法的系统，其特征在于：包括：
基于像点坐标的齐次坐标计算模块，用于利用检查点的像点坐标计算出对应的焦平面齐次坐标；
基于物方坐标的齐次坐标计算模块，用于利用检查点的物方坐标以及摆扫式卫星的成像几何模型逆转换得到焦平面系统中的坐标，再计算出对应的焦平面齐次坐标；
误差分析模块，利用基于像点坐标对应的焦平面齐次坐标与基于物方坐标对应的焦平面齐次坐标的差值分析定位误差。
10.一种可读存储介质，其特征在于：存储了计算机程序，所述计算机程序被处理器调用以实现：
权利要求1-6任一项所述地理定位方法的步骤；
或者所述计算机程序被处理器调用以实现：
权利要求7所述的地理定位误差分析方法的步骤。
CN202210105023.3A 2022-01-27 2022-01-27 一种摆扫式卫星的地理定位方法、定位误差分析方法及其系统 Active CN114593736B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202410371124.4A CN118274841A (zh) 2022-01-27 2022-01-27 一种摆扫式卫星的地理定位误差分析方法及其系统
CN202210105023.3A CN114593736B (zh) 2022-01-27 2022-01-27 一种摆扫式卫星的地理定位方法、定位误差分析方法及其系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210105023.3A CN114593736B (zh) 2022-01-27 2022-01-27 一种摆扫式卫星的地理定位方法、定位误差分析方法及其系统

Related Child Applications (1)

Application Number Title Priority Date Filing Date
CN202410371124.4A Division CN118274841A (zh) 2022-01-27 2022-01-27 一种摆扫式卫星的地理定位误差分析方法及其系统

Publications (2)

Publication Number Publication Date
CN114593736A true CN114593736A (zh) 2022-06-07
CN114593736B CN114593736B (zh) 2024-04-19

Family

ID=81806805

Family Applications (2)

Application Number Title Priority Date Filing Date
CN202210105023.3A Active CN114593736B (zh) 2022-01-27 2022-01-27 一种摆扫式卫星的地理定位方法、定位误差分析方法及其系统
CN202410371124.4A Pending CN118274841A (zh) 2022-01-27 2022-01-27 一种摆扫式卫星的地理定位误差分析方法及其系统

Family Applications After (1)

Application Number Title Priority Date Filing Date
CN202410371124.4A Pending CN118274841A (zh) 2022-01-27 2022-01-27 一种摆扫式卫星的地理定位误差分析方法及其系统

Country Status (1)

Country Link
CN (2) CN114593736B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115164938A (zh) * 2022-07-04 2022-10-11 湖北工业大学 一种线阵摆扫式光学遥感卫星在轨几何定标方法
CN115311363A (zh) * 2022-10-08 2022-11-08 中南大学 一种基于高分卫星影像辅助的足印相机定向方法及系统
CN115856891A (zh) * 2023-03-01 2023-03-28 中国科学院空天信息创新研究院 一种高分辨率星载sar非均匀波束扫描成像方法
CN116058774A (zh) * 2023-01-03 2023-05-05 中山千寻光学有限公司 一种内窥镜装置及成像控制方法

Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2653918A1 (fr) * 1989-10-31 1991-05-03 Europ Agence Spatiale Procede de rectification en temps reel des images des satellites meteorologiques geostationnaires.
US5032023A (en) * 1990-07-02 1991-07-16 General Electric Company Optical fiber based sensor for a variable depth range camera
DE4206304A1 (de) * 1991-02-28 1992-09-03 Ricoh Kk Bilderzeugender spiegel fuer eine optische abtastung mit konstanter geschwindigkeit und optischer scanner mit einem solchen spiegel
CN101692000A (zh) * 2009-09-18 2010-04-07 中国科学院长春光学精密机械与物理研究所 影响摆扫式航空相机两条图像间实际重叠率因素的处理方法
US20110317228A1 (en) * 2010-06-24 2011-12-29 Canon Kabushiki Kaisha Optical scanning apparatus provided with resonance-driven swing mirror and image forming apparatus
RU2457504C1 (ru) * 2011-04-15 2012-07-27 Открытое акционерное общество "Научно-исследовательский институт оптико-электронного приборостроения" (ОАО "НИИ ОЭП") Способ обзора пространства оптико-электронной системой
CN103063200A (zh) * 2012-11-28 2013-04-24 国家测绘地理信息局卫星测绘应用中心 高分辨率光学卫星正射纠正影像生成方法
DE102013203312A1 (de) * 2013-02-27 2014-08-28 Sirona Dental Systems Gmbh Extraoraler Dentalscanner
KR101690950B1 (ko) * 2015-11-19 2016-12-29 국방과학연구소 최소 기준점을 이용한 푸시브롬 스캔 방식 스트립 위성영상의 기하보정 방법 및 이의 시스템
CN107490391A (zh) * 2017-06-06 2017-12-19 上海卫星工程研究所 基于高带宽陀螺测量的天基探测基准标定方法
WO2020093577A1 (zh) * 2018-11-07 2020-05-14 山东理工大学 一种直升机载激光雷达激光扫描姿态角稳定方法与装置
CN111524196A (zh) * 2020-04-30 2020-08-11 中国科学院微小卫星创新研究院 一种摆扫大幅宽光学卫星在轨几何标定方法
CN111538051A (zh) * 2020-04-30 2020-08-14 中国科学院微小卫星创新研究院 一种摆扫大幅宽光学卫星精准处理方法
CN113610924A (zh) * 2021-10-09 2021-11-05 中南大学 顾及地形效应的线性逼近高精度坐标反算方法、设备和介质
CN115311365A (zh) * 2022-07-08 2022-11-08 中国科学院上海技术物理研究所 一种长线列摆扫相机高精度在轨几何定位方法

Patent Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2653918A1 (fr) * 1989-10-31 1991-05-03 Europ Agence Spatiale Procede de rectification en temps reel des images des satellites meteorologiques geostationnaires.
US5032023A (en) * 1990-07-02 1991-07-16 General Electric Company Optical fiber based sensor for a variable depth range camera
DE4206304A1 (de) * 1991-02-28 1992-09-03 Ricoh Kk Bilderzeugender spiegel fuer eine optische abtastung mit konstanter geschwindigkeit und optischer scanner mit einem solchen spiegel
CN101692000A (zh) * 2009-09-18 2010-04-07 中国科学院长春光学精密机械与物理研究所 影响摆扫式航空相机两条图像间实际重叠率因素的处理方法
US20110317228A1 (en) * 2010-06-24 2011-12-29 Canon Kabushiki Kaisha Optical scanning apparatus provided with resonance-driven swing mirror and image forming apparatus
RU2457504C1 (ru) * 2011-04-15 2012-07-27 Открытое акционерное общество "Научно-исследовательский институт оптико-электронного приборостроения" (ОАО "НИИ ОЭП") Способ обзора пространства оптико-электронной системой
CN103063200A (zh) * 2012-11-28 2013-04-24 国家测绘地理信息局卫星测绘应用中心 高分辨率光学卫星正射纠正影像生成方法
DE102013203312A1 (de) * 2013-02-27 2014-08-28 Sirona Dental Systems Gmbh Extraoraler Dentalscanner
KR101690950B1 (ko) * 2015-11-19 2016-12-29 국방과학연구소 최소 기준점을 이용한 푸시브롬 스캔 방식 스트립 위성영상의 기하보정 방법 및 이의 시스템
CN107490391A (zh) * 2017-06-06 2017-12-19 上海卫星工程研究所 基于高带宽陀螺测量的天基探测基准标定方法
WO2020093577A1 (zh) * 2018-11-07 2020-05-14 山东理工大学 一种直升机载激光雷达激光扫描姿态角稳定方法与装置
CN111524196A (zh) * 2020-04-30 2020-08-11 中国科学院微小卫星创新研究院 一种摆扫大幅宽光学卫星在轨几何标定方法
CN111538051A (zh) * 2020-04-30 2020-08-14 中国科学院微小卫星创新研究院 一种摆扫大幅宽光学卫星精准处理方法
CN113610924A (zh) * 2021-10-09 2021-11-05 中南大学 顾及地形效应的线性逼近高精度坐标反算方法、设备和介质
CN115311365A (zh) * 2022-07-08 2022-11-08 中国科学院上海技术物理研究所 一种长线列摆扫相机高精度在轨几何定位方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
PAN, H.等: "Systematic Geolocation Errors of FengYun-3D MERSI-II", IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING, 4 March 2022 (2022-03-04), pages 1, XP011905585, DOI: 10.1109/TGRS.2022.3156999 *
叶思熔;江万寿;李金龙;刘晓波;: "某大视场机载摆扫红外扫描仪几何成像仿真与误差分析", 红外与激光工程, no. 04, 25 April 2017 (2017-04-25), pages 134 - 139 *
夏中秋;黄巧林;何红艳;杨媛丽;孙世君;岳春宇;: "空间相机摆扫成像立体定位精度仿真分析", 中国空间科学技术, vol. 37, no. 03, 31 December 2017 (2017-12-31), pages 117 - 125 *
陈海燕;尹球;胡勇;: "行扫描光学遥感器扫描参数对遥感图像定位结果影响模型", 上海航天, no. 03, 25 June 2012 (2012-06-25), pages 23 - 27 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115164938A (zh) * 2022-07-04 2022-10-11 湖北工业大学 一种线阵摆扫式光学遥感卫星在轨几何定标方法
CN115164938B (zh) * 2022-07-04 2024-06-07 湖北工业大学 一种线阵摆扫式光学遥感卫星在轨几何定标方法
CN115311363A (zh) * 2022-10-08 2022-11-08 中南大学 一种基于高分卫星影像辅助的足印相机定向方法及系统
CN115311363B (zh) * 2022-10-08 2022-12-23 中南大学 一种基于高分卫星影像辅助的足印相机定向方法及系统
CN116058774A (zh) * 2023-01-03 2023-05-05 中山千寻光学有限公司 一种内窥镜装置及成像控制方法
CN116058774B (zh) * 2023-01-03 2023-10-03 中山千寻光学有限公司 一种内窥镜装置及成像控制方法
CN115856891A (zh) * 2023-03-01 2023-03-28 中国科学院空天信息创新研究院 一种高分辨率星载sar非均匀波束扫描成像方法

Also Published As

Publication number Publication date
CN118274841A (zh) 2024-07-02
CN114593736B (zh) 2024-04-19

Similar Documents

Publication Publication Date Title
CN114593736B (zh) 一种摆扫式卫星的地理定位方法、定位误差分析方法及其系统
CN110500995B (zh) 利用rpc参数建立高分辨率卫星影像等效几何成像模型的方法
CN106403902B (zh) 一种星地协同的光学卫星在轨实时几何定位方法及系统
US10803621B2 (en) Method and device for building camera imaging model, and automated driving system for vehicle
CN107202582A (zh) 静止卫星在轨实时成像导航与配准方法
CN111398980B (zh) 一种机载LiDAR数据处理的方法及装置
CN109696182A (zh) 一种星载推扫式光学传感器内方位元素定标方法
CN108828623B (zh) 静止气象卫星成像仪的地球固定网格映射方法
CN110006452B (zh) 高分六号宽视场相机相对几何定标方法及系统
CN106709944B (zh) 卫星遥感图像配准方法
Fraser et al. Sub-metre geopositioning with Ikonos GEO imagery
CN115311365B (zh) 一种长线列摆扫相机高精度在轨几何定位方法
CN107292831A (zh) 静止卫星成像导航与配准的姿态运动补偿方法
CN111912430B (zh) 高轨光学卫星的在轨几何定标方法、装置、设备及介质
CN109828248B (zh) 星载微波遥感仪器角度变形对视向量影响的分析方法
Poli General model for airborne and spaceborne linear array sensors
CN110779517A (zh) 激光雷达的数据处理方法及装置、存储介质和计算机终端
CN112816184A (zh) 一种光学遥感卫星的无控定标方法及装置
CN111598930A (zh) 彩色点云生成方法、装置及终端设备
CN112379344B (zh) 信号补偿方法及装置、设备、存储介质
CN113421300A (zh) 确定鱼眼相机图像中物体实际位置的方法及装置
CN111751007A (zh) 热成像测温方法、装置及存储介质
Tonolo et al. Georeferencing of EROS-A1 high resolution images with rigorous and rational function model
CN118129798B (zh) 天回归轨道遥感卫星的几何定标方法、装置、设备及介质
CN109918724A (zh) 星载遥感仪器光路组件位移变形对视向量影响的分析方法

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