CN105510913B - 基于类光学像方改正的异源光学和sar遥感影像联合定位方法 - Google Patents

基于类光学像方改正的异源光学和sar遥感影像联合定位方法 Download PDF

Info

Publication number
CN105510913B
CN105510913B CN201510767975.1A CN201510767975A CN105510913B CN 105510913 B CN105510913 B CN 105510913B CN 201510767975 A CN201510767975 A CN 201510767975A CN 105510913 B CN105510913 B CN 105510913B
Authority
CN
China
Prior art keywords
remote sensing
image
point
sar
optics
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.)
Expired - Fee Related
Application number
CN201510767975.1A
Other languages
English (en)
Other versions
CN105510913A (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.)
Hubei University of Technology
Original Assignee
Hubei University of Technology
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 Hubei University of Technology filed Critical Hubei University of Technology
Priority to CN201510767975.1A priority Critical patent/CN105510913B/zh
Publication of CN105510913A publication Critical patent/CN105510913A/zh
Application granted granted Critical
Publication of CN105510913B publication Critical patent/CN105510913B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9021—SAR image post-processing techniques
    • G01S13/9027—Pattern recognition for feature extraction
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/904—SAR modes

Abstract

一种基于类光学像方改正的异源光学和SAR遥感影像联合定位方法，包括：1）利用严格几何模型建立立体格网，计算每个格网点在SAR遥感影像上对应的像点坐标，建立虚拟控制点格网，解算SAR遥感影像的RPC参数；2）获取光学和SAR遥感影像间连接点和地面控制点对应像点的像素坐标；3）建立对应的像方误差改正模型，利用基于类光学像方改正的异源光学和SAR遥感影像RFM区域网平差方法，实现精确联合定位。本发明采用类光学像方改正模型对其系统误差进行改正，在近似统一的平差模型下实现异源光学和SAR遥感影像的联合定位，可以取得与严格几何模型相当甚至更优的定位精度，极大降低对地面控制点的要求。

Description

基于类光学像方改正的异源光学和SAR遥感影像联合定位 方法
技术领域
本发明属于遥感影像处理技术领域，涉及一种基于类光学像方改正的异源光学和SAR遥感影像联合定位方法。
背景技术
合成孔径雷达(Synthetic Aperture Radar，SAR)，是一种是用微波探测地表目标的主动式成像传感器。星载SAR，具有全天候、全天时、大范围、强穿透性等独特优点，其获取的SAR遥感影像，可作为光学遥感影像的有益补充，已广泛地应用于地形测绘、国土资源调查、农作物估产、水文监测和灾害评估等领域。异源光学和SAR遥感影像精确联合定位，是诸多应用的重要前提和关键环节，对其研究具有重要的现实意义。
就几何定位而言，构成广义立体像对实现对地目标定位的充要条件是“不同位置拍摄”和“影像存在一定重叠”，与影像是否属于同类传感器无关。利用异源光学和SAR遥感影像进行联合定位，不仅可以降低遥感影像立体成像的要求，而且可以充分利用单景异源遥感影像的互补信息，为此，已有较多学者对其进行深入研究。目前，研究方法主要可分为如下三类：
1)先空间后方交会-再前方交会的方法，即首先利用大量均匀分布的地面控制点分别对每景影像进行定向，然后利用前方交会技术进行定位，如2004年Toutin学者利用其提出的几何模型对Landsat-7ETM+、SPOT4HRV、ASTER VNIR、RADARSAT、ERS等多源影像进行了联合定位试验、2008年邢帅等人尝试将光学与SAR卫星遥感影像联立构成“立体像对”进行复合式定位等。
2)基于严格几何模型的区域网平差方法，如2006年Toutin学者对Landsat-7ETM+、SPOT4HRV、ASTER VNIR、RADARSAT、ERS等影像，采用25个控制点进行平差试验，2008年邢帅对多星多源传感器SPOT、IKONOS、ERS及其RADARSAT影像采用严格几何模型进行了联合区域网平差实验研究等。
3)基于有理函数模型(Rational Function Model，RFM)的区域网平差方法，如秦绪文、张过等人直接将RFM模型用于多源遥感影像的区域网平差，并取得较好的定位结果。
然而，上述方法在实际应用中主要存在如下几点问题：
1)先空间后方交会-再前方交会以及基于严格几何模型的区域网平差方法，对地面控制点要求较高。采用第一种方法时，对每景影像进行定向，均需要野外布设大量、均匀分布的地面控制点；而采用第二种方法时，可一定程度减少所需的控制点数量，但对控制点的分布要求较高，一般需尽量四角布设。
2)基于严格几何模型的区域网平差方法，求解比较复杂。采用第一、二种方法时，当涉及不同类型传感器影像的联合定位时，需要利用不同传感器对应的严格几何模型建立误差方程进行求解，涉及较多复杂的物理参数，求解比较复杂。
3)基于RFM的区域网平差模型，虽具有与传感器无关、形式简单、便于计算等优点，但其像方改正模型主要针对于光学遥感影像。而对于SAR遥感影像，则直接沿用光学遥感影像的仿射变换误差补偿模型。尽管SAR遥感影像亦具有视场角窄、飞行高度高的特点，使其定位误差均具有一定系统性，但由于SAR侧视斜距投影的成像方式，造成地形起伏引起的像点位移较光学影像会严重几倍甚至几十倍，这在某种意义上将破坏误差系统性，从而使得用于光学影像的误差补偿模型不再完全适用于SAR遥感影像。
发明内容
本发明针对目前异源光学和SAR遥感影像联合定位方法存在的缺陷，提出了一种基于类光学像方改正的异源光学和SAR遥感影像联合定位方法，可在有限地面控制点条件下，实现异源光学和SAR遥感影像的精确几何定位，大大降低了野外工作量，具有显著的经济和社会效益。
为了达到上述目的，本发明采用如下技术方案：
一种基于类光学像方改正的异源光学和SAR遥感影像联合定位方法，其特征在于：所述基于类光学像方改正的异源光学和SAR遥感影像联合定位方法包括以下步骤：
1)解算异源光学和SAR遥感影像的RPC参数；
2)利用影像匹配或者人工量测的方法，获取光学和SAR遥感影像间连接点和地面控制点对应像点的像素坐标；
3)待完成步骤2)后，基于步骤1)所解算得到的异源光学和SAR遥感影像的RPC参数以及基于类光学像方改正的异源光学和SAR遥感影像区域网平差，实现精确联合定位。
作为优选，本发明所采用的步骤1)的具体实现方式是：
1.1)构建SAR遥感影像在物方均匀分布的虚拟控制点格网；
1.2)根据步骤1.1)所建立的虚拟控制点格网求解SAR遥感影像有理函数多项式模型；
1.3)对步骤1.2)所建立得到的有理函数多项式模型进行线性化，同时采用谱修正迭代算法计算得到异源光学和SAR遥感影像的RPC参数。
作为优选，本发明所采用的步骤1.1)的具体实现方式是：
1.1.1)根据SAR遥感影像的四个角点的像素坐标求解对应物方点的大地经纬度坐标，确定SAR遥感影像的覆盖范围；
1.1.2)读取SAR遥感影像的覆盖范围内的免费全球DEM数据，确定最大高程值以及最小高程值，并根据给定的高程分层数建立三维立体格网；
1.1.3)对步骤1.1.2)中所建立得到的三维立体格网中各格网点，由其大地经纬度坐标计算各格网点在SAR遥感影像上对应的像点坐标，并进一步对像点坐标进行从侧视斜距投影到光学透视投影的转换，建立物方均匀分布的虚拟控制点格网。
作为优选，本发明所采用的步骤1.1.1)中求解对应物方点的大地经纬度坐标的具体实现方式是：
由某像点像素坐标(x,y)和人工设定的投影面高程h，利用SAR严格几何模型并引入地球椭球模型，计算该像点对应物方点的地心旋转坐标系坐标(X,Y,Z)，通过坐标转换可得到对应的大地经纬度坐标(lat,lon,h)；
所述SAR严格几何模型是根据距离条件和多普勒频率条件来表达雷达图像瞬时构像的几何关系：
(X-Xs)VSX+(Y-Ys)VSY+(Z-Zs)VSZ＝0
其中：
RSlant为卫星天线中心S到地面点P的斜距；
Rp以及Vp分别为地面点P的物方坐标(X,Y,Z)以及速度矢量(Vpx,Vpy,Vpz)；
Rs以及Vs分别为卫星天线中心S成像瞬间对应的物方空间坐标(XS,YS,ZS)以及速度矢量(VSX,VSY,VSZ)；
所述地球椭球模型的函数表达式是：
其中：
re,rp分别为地球椭球的长、短半轴；
h为人工设定的投影面高程；
在计算过程中，首先从SAR遥感影像头文件中读取所需的模型参数，然后以物方点的地心旋转坐标系坐标(X,Y,Z)为未知数，对上述公式线性化，建立误差方程的矩阵形式为：
V＝BX-L,P
式中：
为(X,Y,Z)的改正数对应的系数矩阵；
X为(X,Y,Z)的改正数向量；
L＝-[F1 F2 F3]T为常数向量；
P为权矩阵；
根据误差方程，按照最小二乘法迭代求解解算出(X,Y,Z)，经坐标转换即可得到大地经纬度坐标(lat,lon,h)。
作为优选，本发明所采用的步骤1.1.3)中由其大地经纬度坐标计算各格网点在SAR遥感影像上对应的像点坐标的具体实现方式是：
由某物方点的大地经纬度坐标(lat,lon,h)，转换为地心旋转坐标系坐标(X,Y,Z)，通过SAR严格几何模型的间接定位公式，计算该物方点在SAR遥感影像上对应像点的像素坐标(x,y)：
所述SAR严格几何模型是根据距离条件和多普勒频率条件来表达雷达图像瞬时构像的几何关系：
其中：
RSlant为卫星天线中心S到地面点P的斜距；
Rp以及Vp分别为地面点P的物方坐标(X,Y,Z)以及速度矢量(Vpx,Vpy,Vpz)；
Rs以及Vs分别为卫星天线中心S成像瞬间对应的物方空间坐标(XS,YS,ZS)以及速度矢量(VSX,VSY,VSZ)；
fD为多普勒频率；
λ为卫星天线所接收的雷达的波长；
其中，斜距RSlant为距离向坐标y的函数；
对于斜距影像，RSlant＝R0+yMy，R0为第一斜距，My为斜距影像的距离向分辨率，参数R0和My可以从SAR遥感影像头文件中直接读取；y为像点的像素坐标(x,y)在斜距向的值；
对于SAR遥感影像，方位向上每一行对应不同的卫星位置和速度矢量，表示为成像时刻t的函数：
t＝t0+(x-1)Δt
其中：
Δt为每行SAR遥感影像成像的时间间隔，Δt是根据SAR遥感影像的头文件提供的脉冲重复频率PRF得到，Δt＝1/PRF；
t为对应的成像时刻；
x为方位向坐标；
t0为参考时刻；
ai,bi,ci(i＝0,1,2)为多项式拟合系数，由SAR遥感影像的头文件中提供的若干个卫星状态向量拟合得到；
星载SAR遥感影像一般都进行了多普勒归零处理，且地物在地心旋转坐标系中速度为零，为此多普勒条件方程(2)可简化为：
(X-Xs)VSX+(Y-Ys)VSY+(Z-Zs)VSZ＝0 (4)
联立公式(1)和公式(4)即为距离多普勒数学模型；
对距离多普勒数学模型线性化，建立误差方程的矩阵形式为：
V＝AX-L,P (5)
其中：
为未知数t、y改正数的系数；
X为未知数t、y改正数的向量；
L＝-[F1 F2]T为常数向量，其中，T表示矩阵转置，F1以及F2分别代表多普勒对应两个方程计算得到的值；
P为权矩阵；
根据误差方程(5)，按照最小二乘法迭代求解可解算出未知数(t,y)，进而由公式t＝t0+(x-1)Δt可得到像点像素坐标(x,y)。
作为优选，本发明在步骤1.1.3)中对像点坐标进行从侧视斜距投影到光学透视投影的转换的具体实现方式是：
假设O点为卫星天线中心所在的位置，卫星天线中心距离整个测区平均高程面的高度均为H，O在地面上的投影为E点；A点为SAR遥感影像的近地点所对应的地面点，而B点为远地点所对应的地面点；OO'为SAR遥感影像成像面，DD'为等效焦距为f时的光学影像成像面；
则近地点A对应侧视角为：
θ0＝arccos(H/R0)
其中：
R0为近地点斜距OA的距离；
H为卫星天线航高；
远地点B的斜距为：
OB＝R0+MxW
其中：
Mx为SAR遥感影像的斜距分辨率；
W为SAR遥感影像的宽度；
根据公式(7)，则有近地点A与远地点B之间的距离AB是：
令近地点A与远地点B在等效焦距为f时的光学影像成像面DD'上的投影长度ab为W，则等效焦距f为：
设物方点P高程为hp，物方点P在SAR遥感影像对应像点坐标为(x,y)，像点坐标(x,y)转换为透视中心投影坐标为(x',y')，则有：
作为优选，本发明所采用的步骤1.2)中有理函数多项式模型的具体定义是：
其中：
NumL(P,L,H)、DenL(P,L,H)、Nums(P,L,H)以及Dens(P,L,H)为三次多项式，三次多项式的系数为RPC参数；
(P,L,H)为正则化的地面坐标，(X,Y)为正则化的影像坐标，计算公式如下：
其中：
LAT_OFF、LAT_SCALE、LONG_OFF、LONG_SCALE、HEIGHT_OFF以及HEIGHT_SCALE为地面点坐标的正则化参数；
SAMP_OFF、SAMP_SCALE、LINE_OFF以及LINE_SCALE为影像坐标的正则化参数；
Lat表示纬度，Lon表示经度，Height表示高程，此三项即代表地面点的空间坐标；
Sample代表影像列坐标；Line代表影像行坐标。
作为优选，本发明所采用的步骤1.3)的具体实现方式是：
以RPC参数为未知数，对步骤1.2)所建立得到的有理函数多项式模型进行线性化，列立误差方程：
V＝BX-L,P
其中：
x＝[ai bi ci di]T
P为权矩阵
采用谱修正迭代算法，计算得到SAR遥感影像的RPC参数。
作为优选，本发明所采用的步骤3)的具体实现方式是：
3.1)计算和给定各未知参数的初始值；
3.2)利用异源光学和SAR遥感影像的RPC参数，对步骤2)中所得到的异源光学和SAR遥感影像之间的连接点和地面控制点的每个连接点和地面控制点，以连接点和地面控制点对应的物方坐标，以及各图像对应的仿射变换系数为未知参数列立误差方程；
3.3)利用最小二乘原理，可解算出未知参数的改正数；
3.4)将步骤3.3)解算得到的未知参数的改正数与上次迭代的未知参数值相加，得到本次迭代更新后的未知参数值；当第一次迭代时，未知参数值为初始值；
3.5)利用更新后的未知参数值按步骤3.2)～3.4)重新进行计算，直至达到收敛条件，即完成异源光学和SAR遥感影像的联合定位。。
作为优选，本发明所采用的步骤3.2)的具体实现方式是：
对于光学遥感影像，直接采用像方仿射变换模型，即：
其中：
(x,y)为物方点在光学遥感影像上对应像点的像素坐标；
(Sample,Line)为物方点利用所在光学遥感影像的RPC参数投影到光学遥感影像上的像素坐标；
(pxi,pyi)(i＝0,1,2)是待求光学遥感影像对应的仿射变换系数；其中，px1、py2的初始值为1，其他系数初始值均为0；
对于SAR遥感影像，采用光学遥感影像的像方仿射变换模型，即：
其中：
(x′,y′)为对物方点在SAR遥感影像上对应像点的像素坐标进行透视中心投影转换后的像素坐标；
(sample,line)为利用所在SAR遥感影像的RPC参数计算的地面点投影到等效光学影像面上的像素坐标；
(px′i,py′i)(i＝0,1,2)是SAR遥感影像对应的仿射变换系数；
对于光学遥感影像上的每个连接点和地面控制点的像素坐标，按公式(6)进行线性化，列立误差方程：
V1＝A1t+B1X-L1,P1 (8)
对于SAR遥感影像上的每个连接点和地面控制点的像素坐标，按公式(7)进行线性化，列立误差方程：
V2＝A2t+B2X-L2,P2 (9)
对于地面控制点，将地面控制点的物方坐标视为带权观测值进行处理，还需列立如下误差方程：
V3＝EX-L3,P3 (10)
联立公式(8)～公式(10)，即为类光学像方改正的异源光学和SAR遥感影像区域网平差模型；
其中：
t为各影像对应的仿射变换系数改正数向量；
X为连接点以及地面控制点对应的物方坐标改正数向量；
Ai和Bi(i＝1,2)分别为未知数t和X对应的系数矩阵；
Li(i＝1,2,3)为对应误差方程的常数向量；
Pi为(i＝1,2,3)为相应的权矩阵。
本发明的优点在于：
本发明主要分为如下两个部分：一是异源遥感影像RPC参数计算，二是基于类光学像方改正的异源光学和SAR遥感影像区域网平差。本发明首先对于SAR遥感影像，采用与地形无关的方法，利用严格几何模型建立立体格网，计算每个立体格网点在SAR遥感影像上对应的像点坐标，并进行从侧视斜距投影到光学透视投影的转换，根据每个立体格网点在SAR遥感影像上对应的像点坐标建立虚拟控制点格网，在此基础上，采用谱修正迭代算法，解算SAR遥感影像的RPC参数。其次，利用影像匹配或者人工量测的方法，获取光学以及SAR遥感影像间连接点和地面控制点对应像点的像素坐标。最后，针对不同的光学、SAR遥感影像，建立对应的像方误差改正模型，利用基于类光学像方改正的异源光学和SAR遥感影像RFM区域网平差方法，实现其精确联合定位。本发明针对SAR遥感图像，建立可顾及其系统误差特性的类光学像方改正模型，可以获得更优的定位精度。充分利用RFM模型传感器无关性，形式简单、解算方便等优点，建立了近似统一的异源光学和SAR遥感影像区域网平差模型。本方法具有RFM模型区域网平差方法受地面控制点分布影响较小的优点，对于地形复杂的困难地区，可大大降低地面控制点测量的难度。本发明为异源光学和SAR遥感影像的联合定位提供了一种新方法，创新性地提出了一种基于类光学像方改正的RFM模型区域平差方法，可充分顾及SAR遥感影像的误差特性，非常适用于缺少地面控制点的实际情况。实践表明，当地面控制点数量较少、分布不佳时，可以获得相比严格几何模型更优的定位精度。为此，本方法在借鉴RFM区域网平差模型方法的基础上，针对SAR遥感影像误差特性，提出了一种基于类光学像方改正的异源光学和SAR遥感影像联合定位方法，具有非常重要的现实意义。
附图说明
图1为SAR遥感影像像点从侧视斜距投影转换为透视中心投影的示意图。
具体实施方式
下面结合附图对本发明做进一步详细描述。
实施例1：
以下结合图1对本发明做进一步详细描述。一种基于类光学像方改正的异源光学和SAR遥感影像联合定位方法，各步骤详细阐述如下：
步骤1)异源光学和SAR遥感影像的RPC参数计算
目前，光学遥感影像均提供相应的RPC参数，可直接利用。下面将重点说明本方法中SAR遥感影像对应的RPC参数求解方法：
方式一)由某物方点的大地经纬度坐标(lat,lon,h)，转换为地心旋转坐标系坐标(X,Y,Z)，通过SAR严格几何模型的间接定位公式，计算该物方点在SAR遥感影像上对应像点的像素坐标(x,y)。
SAR严格几何模型，是根据距离条件和多普勒频率条件来表达雷达图像瞬时构像的几何关系，即：
式中：
RSlant为卫星天线中心S到地面点P的斜距；
Rp以及Vp分别为地面点P的物方坐标(X,Y,Z)以及速度矢量(Vpx,Vpy,Vpz)；
Rs以及Vs分别为卫星天线中心S成像瞬间对应的物方空间坐标(XS,YS,ZS)以及速度矢量(VSX,VSY,VSZ)；
fD为多普勒频率；
λ为卫星天线所接收的雷达的波长；
其中，斜距RSlant为距离向坐标y的函数；
对于斜距影像，RSlant＝R0+yMy，R0为第一斜距，My为斜距影像的距离向分辨率，参数R0和My可以从SAR遥感影像头文件中直接读取；y为像点的像素坐标(x,y)在斜距向的值；
对于SAR遥感影像，方位向上每一行对应不同的卫星位置和速度矢量，可表示为成像时刻t的函数：
t＝t0+(x-1)Δt
其中：
Δt为每行SAR遥感影像成像的时间间隔，Δt是根据SAR遥感影像的头文件提供的脉冲重复频率PRF得到Δt＝1/PRF；
t为对应的成像时刻；
x为方位向坐标；
t0为参考时刻；
ai,bi,ci(i＝0,1,2)为多项式拟合系数，由SAR遥感影像的头文件中提供的若干个卫星状态向量拟合得到。
目前，星载SAR遥感影像一般都进行了多普勒归零处理，且地物在地心旋转坐标系中速度为零，为此多普勒条件方程(2)可简化为：
(X-Xs)VSX+(Y-Ys)VSY+(Z-Zs)VSZ＝0 (4)
联立公式(1)和公式(4)即为距离多普勒数学模型。
本发明实施例中，首先从SAR遥感影像头文件中读取所需的模型参数，然后将时间t和列号y作为未知数，对距离多普勒数学模型线性化，建立误差方程的矩阵形式为：
V＝AX-L,P (5)
式中：
为未知数t、y改正数的系数；
X为未知数t、y改正数的向量；
L＝-[F1 F2]T为常数向量，其中，T表示矩阵转置，F1以及F2分别代表多普勒对应两个方程计算得到的值；
P为权矩阵。
根据误差方程(5)，按照最小二乘法迭代求解可解算出未知数(t,y)，进而由公式t＝t0+(x-1)Δt可得到像点像素坐标(x,y)。
方式二：由某像点像素坐标(x,y)和人工设定的投影面高程h，利用SAR严格几何模型并引入地球椭球模型，计算该像点对应物方点的地心旋转坐标系坐标(X,Y,Z)，通过坐标转换可得到对应的大地经纬度坐标(lat,lon,h)。
所引入的地球椭球模型表示如下：
式中：
re,rp分别为地球椭球的长、短半轴；
h为人工设定的投影面高程。
实施例的计算过程中，首先从SAR遥感影像头文件中读取所需的模型参数，然后以物方点的地心旋转坐标系坐标(X,Y,Z)为未知数，对公式(1)、式(4)和式(6)线性化，建立误差方程的矩阵形式为：
V＝BX-L,P (7)
式中：
为(X,Y,Z)的改正数对应的系数矩阵；
X为(X,Y,Z)的改正数向量；
L＝-[F1 F2 F3]T为常数向量；
P为权矩阵。
根据误差方程(7)，按照最小二乘法迭代求解可解算出(X,Y,Z)，经坐标转换即可得到大地经纬度坐标(lat,lon,h)。
1.2)根据方式一以及方式二所得结果，对物方点在SAR遥感影像上对应像点的像素坐标进行从侧视斜距投影到光学透视投影的转换，构建SAR遥感影像在物方均匀分布的虚拟控制点格网。
具体包括以下步骤：
影像覆盖范围的确定：根据影像四个角点的像素坐标，按方式二所述方法，求出对应物方点的大地经纬度坐标，从而确定影像覆盖范围。
立体格网的建立：读取影像覆盖范围内的免费全球DEM数据，确定出的最大和最小高程值，并根据给定的高程分层数建立三维立体格网。
虚拟控制点格网的建立：对三维立体格网中各格网点，由其大地经纬度坐标，根据方式一中所述方法计算各格网点在SAR遥感影像上对应的像点坐标，并进一步对像点坐标进行从侧视斜距投影到光学透视投影的转换，从而建立物方均匀分布的虚拟控制点格网。
下面将以图1为例对本发明的步骤1)进行详细说明：
图1给出了SAR遥感影像像点从侧视斜距投影转换为透视中心投影的示意图。O为SAR天线所在的位置，并假设天线距离整个测区平均高程面的高度为H，其在地面上的投影为E点。A点为SAR遥感影像的近地点所对应的地面点，而B点为远地点所对应的地面点。OO'为SAR遥感影像成像面，DD'为等效焦距为f时的光学影像成像面。设近地点斜距OA的距离为R0，天线航高为H，故可以计算近地点对应侧视角：
θ0＝arccos(H/R0) (8)
远地点的斜距OB＝R0+MxW，其中Mx为SAR遥感影像的斜距分辨率，W为SAR遥感影像的宽度，则有：
令近地点与远地点在光学影像成像面成像面上的投影长度ab为W，则等效焦距f为：
设物方点P高程为hp，其在SAR遥感影像对应像点坐标为(x,y)，像点坐标(x,y)转换为透视中心投影坐标为(x',y')，则有：
1.3)根据步骤1.2)建立的虚拟控制点格网求解SAR遥感影像有理函数多项式模型参数。
RFM模型的具体定义关系式如下：
式中：
NumL(P,L,H)、DenL(P,L,H)、Nums(P,L,H)以及Dens(P,L,H)为三次多项式，三次多项式的系数即称为RPC参数；
(P,L,H)为正则化的地面坐标，(X,Y)为正则化的影像坐标，计算公式如下：
式中：
LAT_OFF、LAT_SCALE、LONG_OFF、LONG_SCALE、HEIGHT_OFF和HEIGHT_SCALE为地面点坐标的正则化参数；
SAMP_OFF、SAMP_SCALE、LINE_OFF和LINE_SCALE为影像坐标的正则化参数。Lat表示纬度，Lon表示经度，Height表示高程，此三项即可代表地面点的空间坐标。
Sample代表影像列坐标，Line代表影像行坐标。
本实施例采用最小二乘方法进行求解，以RPC参数为未知数，对有理函数多项式模型进行线性化，列立误差方程：
V＝BX-L,P (14)
式中
x＝[ai bi ci di]T
P为权矩阵 (15)
采用谱修正迭代算法，可计算出所述的RPC参数。
步骤2)基于类光学像方改正的异源光学和SAR遥感影像区域网平差
待人工量测所需要的异源光学和SAR遥感影像之间的连接点和地面控制点对应像点的像素坐标后，可以基于类光学像方改正的异源光学和SAR遥感影像区域网平差，实现其精确联合定位。
对于光学遥感影像，直接采用像方仿射变换模型，即：
式中，(x,y)为物方点在光学遥感影像上对应像点的像素坐标，(Sample,Line)为物方点利用所在光学遥感影像的RPC参数投影到光学遥感影像上的像素坐标，(pxi,pyi)(i＝0,1,2)是待求光学遥感影像对应的仿射变换系数，px1、py2的初始值为1，其他系数初始值均为0。
对于SAR遥感影像，因将其侧视斜距投影的像素坐标进行了透视中心投影的转换，可采用光学遥感影像的像方仿射变换模型，即：
式中，(x′,y′)为对物方点在SAR遥感影像上对应像点的像素坐标进行透视中心投影转换后的像素坐标，(sample,line)为利用所在SAR遥感影像的RPC参数计算的地面点投影到等效光学影像面上的像素坐标，(px′i,py′i)(i＝0,1,2)是SAR遥感影像对应的仿射变换系数。
本实施例中，首先人工量测所需要的异源光学和SAR遥感影像之间的连接点和地面控制点对应像点的像素坐标，然后利用光学、SAR遥感影像的RPC参数，对每个连接点和控制点，以其对应的物方坐标，以及各图像对应的仿射变换系数为未知数列立误差方程。
对于光学遥感影像上的每个连接点和地面控制点的像素坐标，按公式(16)进行线性化，列立误差方程：
V1＝A1t+B1X-L1,P1 (18)
对于SAR遥感影像上的每个连接点和地面控制点的像素坐标，按公式(17)进行线性化，列立误差方程：
V2＝A2t+B2X-L2,P2 (19)
对于地面控制点，将其物方坐标视为带权观测值进行处理，还需列立如下误差方程：
V3＝EX-L3,P3 (20)
联立式(18)～式(20)，即为类光学像方改正的异源光学和SAR遥感影像区域网平差模型。式中，t为各影像对应的仿射变换系数改正数向量，X为连接点、地面控制点对应的物方坐标改正数向量，Ai和Bi(i＝1,2)分别为未知数t和X对应的系数矩阵，Li(i＝1,2,3)为对应误差方程的常数向量，Pi为(i＝1,2,3)为相应的权矩阵。利用最小二乘原理，可解算出未知参数的改正数。将与上次迭代的参数值相加，并利用更新后的参数值重新进行平差计算，直至达到收敛条件，即可完成异源光学和SAR遥感影像的联合定位。

Claims (4)

1.一种基于类光学像方改正的异源光学和SAR遥感影像联合定位方法，其特征在于：所述基于类光学像方改正的异源光学和SAR遥感影像联合定位方法包括以下步骤：
1)解算异源光学和SAR遥感影像的RPC参数：
1.1)构建SAR遥感影像在物方均匀分布的虚拟控制点格网：
1.1.1)根据SAR遥感影像的四个角点的像素坐标求解对应物方点的大地经纬度坐标，确定SAR遥感影像的覆盖范围：
其中求解对应物方点的大地经纬度坐标的具体实现方式是：
由某像点像素坐标(x,y)和人工设定的投影面高程h，利用SAR严格几何模型并引入地球椭球模型，计算该像点对应物方点的地心旋转坐标系坐标(X,Y,Z)，通过坐标转换可得到对应的大地经纬度坐标(lat,lon,h)；
所述SAR严格几何模型是根据距离条件和多普勒频率条件来表达雷达图像瞬时构像的几何关系：
(X-Xs)VSX+(Y-Ys)VSY+(Z-Zs)VSZ＝0
其中：
RSlant为卫星天线中心S到地面点P的斜距；
Rp以及Vp分别为地面点P的物方坐标(X,Y,Z)以及速度矢量(Vpx,Vpy,Vpz)；
Rs以及Vs分别为卫星天线中心S成像瞬间对应的物方空间坐标(XS,YS,ZS)以及速度矢量(VSX,VSY,VSZ)；
所述地球椭球模型的函数表达式是：
其中：
re,rp分别为地球椭球的长、短半轴；
h为人工设定的投影面高程；
在计算过程中，首先从SAR遥感影像头文件中读取所需的模型参数，然后以物方点的地心旋转坐标系坐标(X,Y,Z)为未知数，对上述公式线性化，建立误差方程的矩阵形式为：
V＝BX-L,P
式中：
为(X,Y,Z)的改正数对应的系数矩阵；
X为(X,Y,Z)的改正数向量；
L＝-[F1 F2 F3]T为常数向量；
P为权矩阵；
根据误差方程，按照最小二乘法迭代求解解算出(X,Y,Z)，经坐标转换即可得到大地经纬度坐标(lat,lon,h)；
1.1.2)读取SAR遥感影像的覆盖范围内的免费全球DEM数据，确定最大高程值以及最小高程值，并根据给定的高程分层数建立三维立体格网；
1.1.3)对步骤1.1.2)中所建立得到的三维立体格网中各格网点，由其大地经纬度坐标计算各格网点在SAR遥感影像上对应的像点坐标，并进一步对像点坐标进行从侧视斜距投影到光学透视投影的转换，建立物方均匀分布的虚拟控制点格网；
其中：由其大地经纬度坐标计算各格网点在SAR遥感影像上对应的像点坐标的具体实现方式是：
由某物方点的大地经纬度坐标(lat,lon,h)，转换为地心旋转坐标系坐标(X,Y,Z)，通过SAR严格几何模型的间接定位公式，计算该物方点在SAR遥感影像上对应像点的像素坐标(x,y)：
所述SAR严格几何模型是根据距离条件和多普勒频率条件来表达雷达图像瞬时构像的几何关系：
其中：
RSlant为卫星天线中心S到地面点P的斜距；
Rp以及Vp分别为地面点P的物方坐标(X,Y,Z)以及速度矢量(Vpx,Vpy,Vpz)；
Rs以及Vs分别为卫星天线中心S成像瞬间对应的物方空间坐标(XS,YS,ZS)以及速度矢量(VSX,VSY,VSZ)；
fD为多普勒频率；
λ为卫星天线所接收的雷达的波长；
其中，斜距RSlant为距离向坐标y的函数；
对于斜距影像，RSlant＝R0+yMy，R0为第一斜距，My为斜距影像的距离向分辨率，参数R0和My可以从SAR遥感影像头文件中直接读取，y为像点的像素坐标(x,y)在斜距向的值；
对于SAR遥感影像，方位向上每一行对应不同的卫星位置和速度矢量，表示为成像时刻t的函数：
t＝t0+(x-1)Δt
其中：
Δt为每行SAR遥感影像成像的时间间隔，Δt是根据SAR遥感影像的头文件提供的脉冲重复频率PRF得到，Δt＝1/PRF；
t为对应的成像时刻；
x为方位向坐标；
t0为参考时刻；
ai,bi,ci(i＝0,1,2)为多项式拟合系数，由SAR遥感影像的头文件中提供的若干个卫星状态向量拟合得到；
星载SAR遥感影像一般都进行了多普勒归零处理，且地物在地心旋转坐标系中速度为零，为此多普勒条件方程(2)可简化为：
(X-Xs)VSX+(Y-Ys)VSY+(Z-Zs)VSZ＝0 (4)
联立公式(1)和公式(4)即为距离多普勒数学模型；
对距离多普勒数学模型线性化，建立误差方程的矩阵形式为：
V＝AX-L,P (5)
其中：
为未知数t、y改正数的系数；
X为未知数t、y改正数的向量；
L＝-[F1 F2]T为常数向量，其中，T表示矩阵转置，F1以及F2分别代表多普勒对应两个方程计算得到的值；
P为权矩阵；
根据误差方程(5)，按照最小二乘法迭代求解可解算出未知数(t,y)，进而由公式t＝t0+(x-1)Δt可得到像点像素坐标(x,y)；
对像点坐标进行从侧视斜距投影到光学透视投影的转换的具体实现方式是：
假设O点为卫星天线中心所在的位置，卫星天线中心距离整个测区平均高程面的高度均为H，O在地面上的投影为E点；A点为SAR遥感影像的近地点所对应的地面点，而B点为远地点所对应的地面点；OO'为SAR遥感影像成像面，DD'为等效焦距为f时的光学影像成像面；
则近地点A对应侧视角为：
θ0＝arccos(H/R0)
其中：
R0为近地点斜距OA的距离；
H为卫星天线航高；
远地点B的斜距为：
OB＝R0+MxW
其中：
Mx为SAR遥感影像的斜距分辨率；
W为SAR遥感影像的宽度；
则有近地点A与远地点B之间的距离AB是：
令近地点A与远地点B在等效焦距为f时的光学影像成像面DD'上的投影长度ab为W，则等效焦距f为：
设物方点P高程为hp，物方点P在SAR遥感影像对应像点坐标为(x,y)，像点坐标(x,y)转换为透视中心投影坐标为(x',y')，则有：
1.2)根据步骤1.1)所建立的虚拟控制点格网求解SAR遥感影像有理函数多项式模型；
1.3)对步骤1.2)所建立得到的有理函数多项式模型进行线性化，同时采用谱修正迭代算法计算得到异源光学和SAR遥感影像的RPC参数；
2)利用影像匹配或者人工量测的方法，获取光学以及SAR遥感影像间连接点和地面控制点对应像点的像素坐标；
3)待完成步骤2)后，基于步骤1)所解算得到的异源光学和SAR遥感影像的RPC参数以及基于类光学像方改正的异源光学和SAR遥感影像区域网平差，实现精确联合定位。
2.根据权利要求1所述的基于类光学像方改正的异源光学和SAR遥感影像联合定位方法，其特征在于：所述步骤1.2)中有理函数多项式模型的具体定义是：
其中：
NumL(P,L,H)、DenL(P,L,H)、Nums(P,L,H)以及Dens(P,L,H)为三次多项式，三次多项式的系数为RPC参数；
(P,L,H)为正则化的地面坐标，(X,Y)为正则化的影像坐标，计算公式如下：
其中：
LAT_OFF、LAT_SCALE、LONG_OFF、LONG_SCALE、HEIGHT_OFF以及HEIGHT_SCALE为地面点坐标的正则化参数；
SAMP_OFF、SAMP_SCALE、LINE_OFF以及LINE_SCALE为影像坐标的正则化参数；
Lat表示纬度，Lon表示经度，Height表示高程，此三项即代表地面点的空间坐标；
Sample代表影像列坐标；Line代表影像行坐标。
3.根据权利要求2所述的基于类光学像方改正的异源光学和SAR遥感影像联合定位方法，其特征在于：步骤3)的具体实现方式是：
3.1)计算和给定各未知参数的初始值；
3.2)利用异源光学和SAR遥感影像的RPC参数，对步骤2)中所得到的异源光学和SAR遥感影像之间的连接点和地面控制点的每个连接点和地面控制点，以连接点和地面控制点对应的物方坐标，以及各图像对应的仿射变换系数为未知参数列立误差方程；
3.3)利用最小二乘原理，可解算出未知参数的改正数；3.4)将步骤3.3)解算得到的未知参数的改正数与上次迭代的未知参数值相加，得到本次迭代更新后的未知参数值；当第一次迭代时，未知参数值为初始值；
3.5)利用更新后的未知参数值按步骤3.2)～3.4)重新进行计算，直至达到收敛条件，即完成异源光学和SAR遥感影像的联合定位。
4.根据权利要求3所述的基于类光学像方改正的异源光学和SAR遥感影像联合定位方法，其特征在于：所述步骤3.2)的具体实现方式是：
对于光学遥感影像，直接采用像方仿射变换模型，即：
其中：
(x,y)为物方点在光学遥感影像上对应像点的像素坐标；
(Sample,Line)为物方点利用所在光学遥感影像的RPC参数投影到光学遥感影像上的像素坐标；
(pxi,pyi)(i＝0,1,2)是待求光学遥感影像对应的仿射变换系数；其中，px1、py2的初始值为1，其他系数初始值均为0；
对于SAR遥感影像，采用光学遥感影像的像方仿射变换模型，即：
其中：
(x′,y′)为对物方点在SAR遥感影像上对应像点的像素坐标进行透视中心投影转换后的像素坐标；
(sample,line)为利用所在SAR遥感影像的RPC参数计算的地面点投影到等效光学影像面上的像素坐标；
(px′i,py′i)(i＝0,1,2)是SAR遥感影像对应的仿射变换系数；
对于光学遥感影像上的每个连接点和地面控制点的像素坐标，按公式(6)进行线性化，列立误差方程：
V1＝A1t+B1X-L1,P1 (8)
对于SAR遥感影像上的每个连接点和地面控制点的像素坐标，按公式(7)进行线性化，列立误差方程：
V2＝A2t+B2X-L2,P2 (9)
对于地面控制点，将地面控制点的物方坐标视为带权观测值进行处理，还需列立如下误差方程：
V3＝EX-L3,P3 (10)
联立公式(8)～公式(10)，即为类光学像方改正的异源光学和SAR遥感影像区域网平差模型；
其中：
t为各影像对应的仿射变换系数改正数向量；
X为连接点以及地面控制点对应的物方坐标改正数向量；
Ai和Bi(i＝1,2)分别为未知数t和X对应的系数矩阵；
Li(i＝1,2,3)为对应误差方程的常数向量；
Pi为(i＝1,2,3)为相应的权矩阵。
CN201510767975.1A 2015-11-11 2015-11-11 基于类光学像方改正的异源光学和sar遥感影像联合定位方法 Expired - Fee Related CN105510913B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510767975.1A CN105510913B (zh) 2015-11-11 2015-11-11 基于类光学像方改正的异源光学和sar遥感影像联合定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510767975.1A CN105510913B (zh) 2015-11-11 2015-11-11 基于类光学像方改正的异源光学和sar遥感影像联合定位方法

Publications (2)

Publication Number Publication Date
CN105510913A CN105510913A (zh) 2016-04-20
CN105510913B true CN105510913B (zh) 2018-01-09

Family

ID=55719023

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510767975.1A Expired - Fee Related CN105510913B (zh) 2015-11-11 2015-11-11 基于类光学像方改正的异源光学和sar遥感影像联合定位方法

Country Status (1)

Country Link
CN (1) CN105510913B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110500995A (zh) * 2019-07-12 2019-11-26 武汉大学 利用rpc参数建立高分辨率卫星影像等效几何成像模型的方法

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108176854A (zh) * 2017-12-27 2018-06-19 深圳市圆梦精密技术研究院 电子束散焦的修正方法
CN108761453B (zh) * 2018-03-29 2021-02-05 航天东方红卫星有限公司 一种光学卫星与sar卫星图像融合的成像视角优化方法
CN108761444B (zh) * 2018-05-24 2021-12-21 中国科学院电子学研究所 联合星载sar和光学图像计算地面点高度的方法
CN108764161B (zh) * 2018-05-30 2022-05-03 北京大学 基于极坐标系的破解稀疏阵引发的病态奇异性的遥感影像处理方法和装置
CN109115218A (zh) * 2018-07-06 2019-01-01 航天星图科技（北京）有限公司 一种无地面控制点的遥感影像定位方法
CN109100719B (zh) * 2018-07-27 2023-01-10 国家测绘地理信息局第三航测遥感院 基于星载sar影像与光学影像的地形图联合测图方法
CN109919835B (zh) * 2019-03-20 2022-07-26 湖北省电力勘测设计院有限公司 基于多源卫星遥感影像联合平差的境外电力选线方法
CN110133653B (zh) * 2019-05-29 2020-12-08 中国空间技术研究院 一种基于dsm数据的星载sar图像快速间接定位方法
CN111197975A (zh) * 2019-12-09 2020-05-26 中国资源卫星应用中心 一种基于Rodrigues的影像姿态估计方法
CN111754458B (zh) * 2020-05-18 2023-09-15 北京吉威空间信息股份有限公司 面向几何精处理的卫星影像立体空间基准框架构建方法
CN111553986B (zh) * 2020-05-19 2021-03-05 北京数字绿土科技有限公司 三角网的构建方法、构建装置和数字表面模型的生成方法
CN113899387B (zh) * 2021-09-27 2023-09-22 武汉大学 一种基于验后补偿的光学卫星遥感影像区域网平差方法及系统
CN114565653B (zh) * 2022-03-02 2023-07-21 哈尔滨工业大学 一种存在旋转变化和尺度差异的异源遥感图像匹配方法
CN114562982B (zh) * 2022-03-09 2023-09-26 北京市遥感信息研究所 一种光学和sar异源卫星影像联合平差的定权方法和装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102213762A (zh) * 2011-04-12 2011-10-12 中交第二公路勘察设计研究院有限公司 基于rfm模型的多源星载sar影像自动匹配方法
CN102520406A (zh) * 2011-12-20 2012-06-27 中国测绘科学研究院 基于影像坐标映射的sar立体模型构建与量测方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9335409B2 (en) * 2013-03-20 2016-05-10 Raytheon Company Bistatic inverse synthetic aperture radar imaging

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102213762A (zh) * 2011-04-12 2011-10-12 中交第二公路勘察设计研究院有限公司 基于rfm模型的多源星载sar影像自动匹配方法
CN102520406A (zh) * 2011-12-20 2012-06-27 中国测绘科学研究院 基于影像坐标映射的sar立体模型构建与量测方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"基于岭估计的有理多项式参数求解方法";袁修孝等;《武汉大学学报信息科学版》;20081130;第33卷(第11期);第1130-1133页 *
"多源SAR影像区域网平差的比较研究";吴颖丹;《测绘科学》;20140930;第39卷(第9期);第24-28页 *
"星载SAR遥感影像的精确几何定位";吴颖丹;《中国博士学位论文全文数据库基础科学辑》;20150515;A008-37 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110500995A (zh) * 2019-07-12 2019-11-26 武汉大学 利用rpc参数建立高分辨率卫星影像等效几何成像模型的方法

Also Published As

Publication number Publication date
CN105510913A (zh) 2016-04-20

Similar Documents

Publication Publication Date Title
CN105510913B (zh) 基于类光学像方改正的异源光学和sar遥感影像联合定位方法
CN106526593B (zh) 基于sar严密成像模型的子像素级角反射器自动定位方法
Fraser et al. Insights into the affine model for high-resolution satellite sensor orientation
CN106960174A (zh) 高分影像激光雷达高程控制点提取及其辅助定位方法
CN106556822B (zh) 星载滑动聚束sar瞄准精度在轨测试方法
CN109613583A (zh) 基于单星与地面站测向及联合测时差的无源目标定位方法
Wang et al. Planar block adjustment and orthorectification of ZY-3 satellite images
Jiao et al. A generic framework for improving the geopositioning accuracy of multi-source optical and SAR imagery
CN104406770B (zh) 波像差测量模块的畸变测量装置和畸变校正方法
CN110030968A (zh) 一种基于星载立体光学影像的地面遮挡物仰角测量方法
CN102466799A (zh) 基于pos运动数据的干涉sar回波数据仿真方法
Jiao et al. Geolocation accuracy improvement of multiobserved GF-3 spaceborne SAR imagery
CN110160503A (zh) 一种顾及高程的无人机景观匹配定位方法
Jo et al. Dense thermal 3d point cloud generation of building envelope by drone-based photogrammetry
Zhang Photogrammetric processing of low altitude image sequences by unmanned airship
CN116045920A (zh) 一种基于gf7号立体像对提取dem的方法
Kartal et al. Comperative analysis of different geometric correction methods for very high resolution pleiades images
CN111505626B (zh) 一种利用底视差分干涉测量二维地形坡度的方法
CN108008367A (zh) 星载单航过InSAR系统电离层误差校正方法
Liu et al. Adaptive re-weighted block adjustment for multi-coverage satellite stereo images without ground control points
CN109828234B (zh) 一种异源光学平台对空间目标融合定位方法
CN114839632B (zh) 一种联合非摄影测量观测条件约束的sar无控几何定标方法及其系统
Murthy et al. Analysis of DEM generated using Cartosat-1 stereo data over Mausanne Les Alpiles–Cartosat scientific appraisal programme (CSAP TS–5)
Yao et al. Multi-observation block adjustment by rational function model without ground control points
CN108764161A (zh) 基于极坐标系的破解稀疏阵引发的病态奇异性的遥感影像处理方法和装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Wu Yingdan

Inventor after: Zhang Heng

Inventor after: Yan Jin

Inventor after: Ming Yang

Inventor after: Huang Chuyun

Inventor after: Lv Hui

Inventor after: Ye Zhiwei

Inventor after: Lei Yong

Inventor after: Zhu Ying

Inventor after: Yang Fei

Inventor after: Wang Xinying

Inventor before: Wu Yingdan

Inventor before: Ming Yang

Inventor before: Huang Chuyun

Inventor before: Lv Hui

Inventor before: Ye Zhiwei

Inventor before: Lei Yong

Inventor before: Zhu Ying

GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180109

Termination date: 20201111