CN107705329B - 基于几何约束的高分辨率光学卫星凝视影像配准方法 - Google Patents
基于几何约束的高分辨率光学卫星凝视影像配准方法 Download PDFInfo
- Publication number
- CN107705329B CN107705329B CN201711003290.5A CN201711003290A CN107705329B CN 107705329 B CN107705329 B CN 107705329B CN 201711003290 A CN201711003290 A CN 201711003290A CN 107705329 B CN107705329 B CN 107705329B
- Authority
- CN
- China
- Prior art keywords
- image
- dense
- point
- model
- coordinates
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/38—Registration of image sequences
-
- 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/10028—Range image; Depth image; 3D point clouds
-
- 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/30—Subject of image; Context of image processing
- G06T2207/30181—Earth observation
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Processing Or Creating Images (AREA)
Abstract
本发明公开了一种基于几何约束的高分辨率光学卫星凝视影像配准方法,包括:S100利用附加像方误差补偿的有理函数模型构建影像的自由网平差定向模型;S200基于自由网平差定向模型,对影像中最优交会条件的立体影像对进行密集匹配,获得密集同名点坐标;S300生成密集点云数据,对密集点云数据进行粗差点剔除,并生成数字表面模型数据;S400基于自由网平差定向模型和DSM数据,对影像依次进行几何校正。本发明利用几何定位一致性,实现了高分辨率光学卫星多角度凝视序列影像的高精度配准,解决了卫星平台控制不准、姿态测量不准引起的多角度凝视序列影像配准问题,可为后续应用提供高精度数据基础。
Description
技术领域
本发明属于遥感卫星数据处理技术领域,更确切的说,本发明是在光学遥感卫星地面数据处理程中,采用的一种基于几何约束的高分辨率光学卫星凝视影像配准方法。
背景技术
随着分辨率的逐步提高,为了弥补成像幅宽减小带来的成像效率降低的不足,高分辨率遥感卫星对平台敏捷性能的要求不断提高,使其能快速机动成像满足对目标区域大范围或多角度等成像需求。凝视成像模式是一种光学卫星敏捷成像模式,通过调整卫星平台姿态,卫星在飞行轨道上的不同位置可从不同角度对同一地区成像,获取同一地区的多角度凝视序列影像,为动目标检测、三维建模等应用提供数据来源。然而,在不同角度成像时,由于卫星对地成像存在指向误差,不同视角成像的姿态观测误差也不尽相同,导致获取的多角度序列影像间存在配准误差,无法直接使用,因此需对多角度凝视序列影像进行高精度配准处理,以满足后续应用的对影像间配准精度的需求。
现有的序列影像配准方法多采用像方配准方法。通过匹配影像间的同名点,然后建立像方仿射变换模型,补偿影像间相对误差,最后选取基准影像,对其他影像进行校正来实现序列影像间的配准。但由于高分辨率影像尺寸大,行列数可达数万以上,单一二维像方仿射变换无法补偿影像间的变形。同时考虑中心投影成像特性,对于地形起伏较大、建筑物较高的地区,由于不同角度成像导致影像间的投影差也是凝视成像的一个特点,像方配准方法可能会改变影像投影特性。
在高精度在轨几何定标、高精度定轨、定姿等技术的支持下,可建立严密几何成像模型,继而生成其替代模型—有理函数模型(rational function model,RFM)。高分辨率卫星影像通过RFM即可直接建立像点与地面点的坐标关系。因此可以基于凝视影像间几何定位一致性的约束来实现多角度凝视序列影像的配准。
发明内容
针对凝视成像模式获取的多角度序列影像配准问题,本发明提供了一种基于几何约束的高分辨率光学卫星凝视影像配准方法。
本发明提供的基于几何约束的高分辨率光学卫星凝视影像配准方法,包括步骤:
S100利用附加像方误差补偿的有理函数模型对影像进行自由网平差定向,获得各影像对的自由网平差定向模型;所述影像为多角度凝视序列影像;
S200基于自由网平差定向模型,对影像中最优交会条件的立体影像对进行密集匹配,获得密集同名点坐标;所述最优交会条件的立体影像对同时满足原则(1)和原则(2):
原则(1)为:根据地形情况立体影像对的交会角为10度~45度,对平坦地区,立体影像对的交会角为30度~45度,对于地形起伏大的地区,立体影像对的交会角为10度~30度;
原则(2)为:立体影像对为影像序列中成像俯仰角大小绝对值最接近的一组影像;
S300基于自由网平差定向模型和密集同名点坐标生成密集点云数据,对密集点云数据进行基于几何约束的粗差点剔除,根据粗差点剔除后的密集点云数据生成数字表面模型数据;
所述对密集点云数据进行基于几何约束的粗差点剔除具体为:
对每一对密集同名点对,将密集同名点的物方坐标带入所在影像对的自由网平差定向模型,求得像点像素坐标,记为(lnew,snew),计算(lnew,snew)与密集同名点匹配的像素坐标(l,s)的坐标差,若坐标差大于2倍中误差,即剔除该密集同名点;所述中误差为所有密集同名点对的坐标差的中误差;所述密集同名点对的坐标差指密集同名点对的(lnew,snew)与像素坐标(l,s)的坐标差;
S400基于自由网平差定向模型和DSM数据,对影像依次进行几何校正;
所述对影像依次进行几何校正,具体为:
根据自由网平差定向模型和DSM数据计算影像四角点的经纬度及高程坐标,求四角点的正北外接矩形,即几何校正后影像的地面覆盖范围;
将地面覆盖范围划分格网块;
对各格网块分别进行几何校正,具体为:
采用直接线性变换模型建立格网块中点的像点坐标和物方点坐标的关联;
根据数字表面模型数据生成地面覆盖范围内的高程最大值Hmax和高程最小值Hmin,并获得高程中值Hmid=0.5*(Hmax+Hmin);
分别以Hmax、Hmin和Hmid为高程建立三层高程面,各高程面上均均匀选取数量相等、分布相同的若干虚拟控制点;
利用自由网平差定向模型反投影计算虚拟控制点对应的像点坐标;
利用虚拟控制点,采用个最小二乘法估计直接线性变换模型中的直接线性变换系数;
利用直接线性变化模型将格网块内各点逐一反投影到原始影像,得到校正后影像每一个像素点在原始影像中的坐标;
根据校正后影像每一个像素点在原始影像中的坐标和灰度进行灰度重采样获得各像素点对应的灰度值,将灰度值赋给校正后影像中对应的像素点。
进一步的,步骤S100进一步包括子步骤:
S101影像两两进行组合获得影像对,对影像对进行匹配获得连接点;
S102利用附加像方误差补偿的有理函数模型,根据连接点构建自由网平差定向模型;
S103采用最小二乘法解算自由网平差定向模型的误差补偿系数。
进一步的,步骤S200中,所述对影像中最优交会条件的立体影像对进行密集匹配,进一步包括:
根据立体影像对的自由网平差定向模型建立影像金字塔,由粗到精逐层对立体影像对进行像素级匹配,得到像素级匹配结果;
以像素级匹配结果为初值,通过最小二乘法对立体影像对进行匹配,得到子像素匹配结果。
进一步的,步骤S300中,所述基于自由网平差定向模型和密集同名点坐标生成密集点云数据具体为:
利用自由网平差定向模型和密集同名点坐标,通过前方交会原理计算密集同名点的物方坐标(B,L,H),从而获得具有地理坐标的点云数据,即密集点云数据。
进一步的,步骤S300中,所述根据粗差点剔除后的密集点云数据生成数字表面模型数据,具体为:
根据粗差点剔除后的密集点云数据获得不规则三角网形式的数字表面模型,将三角网栅格化为带地理坐标的图像格式,即获得数字表面模型数据。
本发明提供的基于几何约束的高分辨率光学卫星凝视影像配准系统,包括:
第一模块,用来利用附加像方误差补偿的有理函数模型对影像进行自由网平差定向,获得各影像对的自由网平差定向模型;所述影像为多角度凝视序列影像;
第二模块,用来基于自由网平差定向模型,对影像中最优交会条件的立体影像对进行密集匹配,获得密集同名点坐标;所述最优交会条件的立体影像对同时满足原则(1)和原则(2):
原则(1)为:根据地形情况立体影像对的交会角为10度~45度,对平坦地区,立体影像对的交会角为30度~45度,对于地形起伏大的地区,立体影像对的交会角为10度~30度;
原则(2)为:立体影像对为影像序列中成像俯仰角大小绝对值最接近的一组影像;
第三模块,用来基于自由网平差定向模型和密集同名点坐标生成密集点云数据,对密集点云数据进行基于几何约束的粗差点剔除,根据粗差点剔除后的密集点云数据生成数字表面模型数据;
所述对密集点云数据进行基于几何约束的粗差点剔除具体为:
对每一对密集同名点对,将密集同名点的物方坐标带入所在影像对的自由网平差定向模型,求得像点像素坐标,记为(lnew,snew),计算(lnew,snew)与密集同名点匹配的像素坐标(l,s)的坐标差,若坐标差大于2倍中误差,即剔除该密集同名点;所述中误差为所有密集同名点对的坐标差的中误差;所述密集同名点对的坐标差指密集同名点对的(lnew,snew)与像素坐标(l,s)的坐标差;
第四模块,用来基于自由网平差定向模型和DSM数据,对影像依次进行几何校正;
所述对影像依次进行几何校正,具体为:
根据自由网平差定向模型和DSM数据计算影像四角点的经纬度及高程坐标,求四角点的正北外接矩形,即几何校正后影像的地面覆盖范围;
将地面覆盖范围划分格网块;
对各格网块分别进行几何校正,具体为:
采用直接线性变换模型建立格网块中点的像点坐标和物方点坐标的关联;
根据数字表面模型数据生成地面覆盖范围内的高程最大值Hmax和高程最小值Hmin,并获得高程中值Hmid=0.5*(Hmax+Hmin);
分别以Hmax、Hmin和Hmid为高程建立三层高程面,各高程面上均均匀选取数量相等、分布相同的若干虚拟控制点;
利用自由网平差定向模型反投影计算虚拟控制点对应的像点坐标;
利用虚拟控制点,采用个最小二乘法估计直接线性变换模型中的直接线性变换系数;
利用直接线性变化模型将格网块内各点逐一反投影到原始影像,得到校正后影像每一个像素点在原始影像中的坐标;
根据校正后影像每一个像素点在原始影像中的坐标和灰度进行灰度重采样获得各像素点对应的灰度值,将灰度值赋给校正后影像中对应的像素点。
和现有技术相比,本发明具有如下优点和有益效果:
本发明利用几何定位一致性,实现了高分辨率光学卫星多角度凝视序列影像的高精度配准,解决了卫星平台控制不准、姿态测量不准引起的多角度凝视序列影像配准问题。
本发明依靠多角度凝视序列影像的几何成像模型,利用自由网平差方法来精化几何定位一致性,同时利用交会角良好的立体影像对生成凝视区域的地表模型,以解决在几何校正过程中由于地形起伏高程误差引起的定位偏差,保证影像配准效为后续应用提供高精度数据基础。
附图说明
图1为本发明基于几何约束的高分辨率光学卫星凝视影像配准方法的具体流程图;
图2为本发明多角度凝视成像的原理示意图;
图3为基于自由网平差定向模型生成虚拟控制点的示意图。
具体实施方式
为了更清楚地说明本发明实施例和/或现有技术中的技术方案,下面将对照附图说明本发明的具体实施方式。显而易见地,下面描述中的附图仅仅是本发明的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图,并获得其他的实施方式。
图1表示的是本发明基于几何约束的高分辨率光学卫星凝视影像配准方法的具体流程图,以下将针对该具体流程中各步骤,对本发明方法的具体实施方式做详细描述。
本发明方法依靠多角度凝视序列影像的几何成像模型,利用自由网平差方法来精化几何定位一致性,同时利用交会角良好的立体影像对生成凝视区域的地表模型,以解决在几何校正过程中由于地形起伏高程误差引起的定位偏差,保证影像配准效果。所述方法具体步骤如下:
步骤100,利用附加像方误差补偿的有理函数模型对多角度凝视序列影像进行自由网平差定向,获得自由网平差定向模型。本步骤用来使凝视影像间具备几何定位一致性,从而作为影像配准的几何约束。
本步骤进一步包括:
步骤101,对凝视序列影像两两进行组合,获得影像对,利用SIFT算子对各影像对匹配一定数量的连接点。
步骤102,对连接点进行自由网平差定向模型构建。
影像对的相对定向基于误差补偿的有理多项式模型进行。公式(1)为有理函数模型(rational function model,RFM),公式(2)为附加像方误差补偿的RFM,即精化RFM。
RFM如下:
式(1)中:
(ln,sn)是像点的像素坐标(l,s)对应的像方归一化坐标;
(U,V,W)是物方坐标(B,L,H)的归一化坐标,物方坐标(B,L,H)即像点对应物方点的经度、纬度和高程;
NumL(U,V,W)=a1+a2V+a3U+a4W+a5VU+a6VW+a7UW+a8V2+a9U2+
a10W2+a11UVW+a12V3+a13VU2+a14VW2+a15V2U+a16U3+
a17UW2+a18V2W+a19U2W+a20W3;
DenL(U,V,W)=b1+b2V+b3U+b4W+b5VU+b6VW+b7UW+b8V2+b9U2+
b10W2+b11UVWb12V3+b13VU2+b14VW2+b15V2U+b16U3+
b17UW2+b18V2W+b19U2W+b20W3;
NumS(U,V,W)=c1+c2V+c3U+c4W+c5VU+c6VW+c7UW+c8V2+c9U2+
c10W2+c11UVW+c12V3+c13VU2+c14VW2+c15V2U+c16U3+
c17UW2+c18V2W+c19U2W+c20W3;
DenS(U,V,W)=d1+d2V+d3U+d4W+d5VU+d6VW+d7UW+d8V2+d9U2+
d10W2+d11UVW+d12V3+d13VU2+d14VW2+d15V2U+d16U3+;
d17UW2+d18V2W+d19U2W+d20W3
ai,bi,ci,di(i=1,2,...20)为有理多项式系数。
令LineOff、SampleOff为像方坐标的平移值,LineScale、SampleScale为像方坐标的缩放值,LonOff、LatOff、HeiOff为物方坐标的平移值,LonScale、LatScale、HeiScale为物方坐标的缩放值,则有
精化RFM如下:
其中,A0、A1、A2、B0、B1、B2为误差补偿参数。
利用匹配的连接点分别对各影像对建立基于精化RFM的自由网平差定向模型。对于八视凝视成像获取的序列影像,见图2,假设每一组连接点均为八度重叠同名点,每一组连接点即可建立16个方程,由于未知参数包括8组像方误差补偿参数和连接点对应的物方坐标,至少需要4组连接点建立64个方程可解算8组像方误差补偿参数和所有连接点对应的物方坐标。
步骤103,根据所建立的自由网平差定向模型,根据最小二乘法解算自由网平差定向模型的误差补偿系数,消除影像对间的相对误差,得到精化RFM,保证序列影像的几何定位具有一致性。
步骤2,基于自由网平差定向模型,对多角度凝视序列影像中最优交会条件的立体影像对进行密集匹配,获得密集同名点坐标。
本步骤进一步包括:
步骤201,从多角度凝视序列影像中选取最优交会条件的影像对作为立体影像对,所述立体影像从多角度凝视序列影像中的影像对中选取,同时满足原则(1)和原则(2):
(1)根据地形情况,立体影像对的交会角为10度~45度,对平坦地区,立体影像对的交会角应大于30度,这样可提高交会精度;对于地形起伏大地区,立体影像对的交会角则不应超过30度,避免地形起伏引起的变形导致后续密集匹配结果不理想问题;
(2)立体影像对尽量成对称成像,即俯仰方向相反,更具体的,即立体影像对为多角度凝视影像序列中成像俯仰角大小绝对值最接近的一组影像。
步骤202,根据立体影像对的自由王平差定向模型建立影像金字塔,由粗到精逐层对立体影像对进行像素级匹配,得到像素级匹配结果。
立体影像对前后片分别建立3层金字塔,每一级影像的像素总数对于前一级影像以3×3的倍数缩小。从金字塔顶端开始,以相关系数为匹配测度,逐层进行匹配,直到底层原始影像结束,得到同名点的像素级坐标,即像素级匹配结果。
步骤203,以像素级匹配结果为初值,通过最小二乘法对立体影像对进行匹配,得到子像素匹配结果。所述子像素匹配结果即同名点的子像素级坐标,也即密集同名点坐标。
步骤3,基于步骤1得到的自由网平差定向模型和步骤2得到的密集同名点坐标生成密集点云数据,对密集点云数据进行基于几何约束的粗差点剔除,根据粗差点剔除后的密集点云数据生成数字表面模型(digital surface model,DSM)数据。
本步骤进一步包括:
步骤301,利用步骤S100获得的自由网平差定向模型,通过前方交会原理计算步骤S200中密集同名点的物方坐标(B,L,H),得到具有地理坐标的点云数据,即密集点云数据。
步骤302,对密集点云数据进行基于几何约束的粗差点剔除,具体为:对每一对密集同名点对,其在影像上的像素坐标记为(l,s),将密集同名点的物方坐标(B,L,H)带入影像对的自由网平差定向模型,求得像点像素坐标,记为(lnew,snew),计算(lnew,snew)与密集同名点匹配的像素坐标(l,s)的坐标差,若该坐标差大于2倍中误差,即剔除该密集同名点。
所述中误差为所有密集同名点对的坐标差的中误差,所述密集同名点对的坐标差指密集同名点对的(lnew,snew)与像素坐标(l,s)的坐标差。
步骤303,根据粗差点剔除后的密集点云数据获得不规则三角网(TIN)形式的DSM,设置DSM的采样间隔,将TIN栅格化为带地理坐标的图像格式,图像灰度值即为高程值,即栅格化的DSM。
步骤4,基于步骤1得到的自由网平差定向模型和步骤3得到的DSM数据,依次对多角度凝视序列影像进行几何校正。
本步骤进一步包括子步骤:
步骤401,根据步骤1得到的自由网平差定向模型和步骤3得到的DSM数据,计算影像四角点的经纬度及高程坐标,求四角点的正北外接矩形,即几何校正后影像的地面覆盖范围,根据地面采样间隔得到校正后影像的列数和行数。
步骤402,将地面覆盖范围划分格网块,所划分的格网块有M行N列。
步骤403,为了简化计算,对于每一个地面格网块,将直接线性变换模型替代有理函数模型建立像点坐标和物方点坐标的关联。
直接线性变换模型如下式所示:
式(3)中,lj(j=1,2,...,11)为直接线性变换系数。
如图3所示,根据步骤3生成的DSM数据得到地面覆盖范围内的高程最大值Hmax和高程最小值Hmin,分别以高程最大值Hmax、高程最小值Hmin和高程中值Hmid为高程建立三层高程面,其中,Hmid=0.5*(Hmax+Hmin)。在每一层高程面上均均匀选取m×n个虚拟控制点,利用步骤1得到的自由网平差定向模型反投影计算虚拟控制点对应的像点坐标。利用虚拟控制点,采用个最小二乘法估计直接线性变换系数,由此得到直接线性变换模型,从而建立格网块内像点坐标与物方点坐标映射关系。最后通过直接线性变化模型将格网内每一个点逐一反投影到原始影像上,得到校正后影像每一个像素点在原始影像中的坐标。由于反投影计算得到的像点坐标通常不是整数,则需通过灰度重采样,例如采用双线性内插模型,得到影像灰度值,并赋给对应校正后影像的像素点。
步骤404,重复步骤403对各格网块分别进行几何校正。
步骤405,重复步骤401-404,对所有多角度凝视序列影像分别进行几何校正,得到配准后影像。
上述实施例用来解释说明本发明,而不是对本发明进行限制,在本发明的精神和权利要求的保护范围内,对本发明做出任何的修改和改变,都落入本发明的保护范围。
Claims (6)
1.基于几何约束的高分辨率光学卫星凝视影像配准方法,其特征是,包括:
S100利用附加像方误差补偿的有理函数模型对影像进行自由网平差定向,获得各影像对的自由网平差定向模型;所述影像为多角度凝视序列影像;
S200基于自由网平差定向模型,对影像中最优交会条件的立体影像对进行密集匹配,获得密集同名点坐标;所述最优交会条件的立体影像对同时满足原则(1)和原则(2):
原则(1)为:根据地形情况立体影像对的交会角为10度~45度,对平坦地区,立体影像对的交会角为30度~45度,对于地形起伏大的地区,立体影像对的交会角为10度~30度;
原则(2)为:立体影像对为影像序列中成像俯仰角大小绝对值最接近的一组影像;
S300基于自由网平差定向模型和密集同名点坐标生成密集点云数据,对密集点云数据进行基于几何约束的粗差点剔除,根据粗差点剔除后的密集点云数据生成数字表面模型数据;
所述对密集点云数据进行基于几何约束的粗差点剔除具体为:
对每一对密集同名点对,将密集同名点的物方坐标带入所在影像对的自由网平差定向模型,求得像点像素坐标,记为(lnew,snew),计算(lnew,snew)与密集同名点匹配的像素坐标(l,s)的坐标差,若坐标差大于2倍中误差,即剔除该密集同名点;所述中误差为所有密集同名点对的坐标差的中误差;所述密集同名点对的坐标差指密集同名点对的(lnew,snew)与像素坐标(l,s)的坐标差;
S400基于自由网平差定向模型和DSM数据,对影像依次进行几何校正;
所述对影像依次进行几何校正,具体为:
根据自由网平差定向模型和DSM数据计算影像四角点的经纬度及高程坐标,求四角点的正北外接矩形,即几何校正后影像的地面覆盖范围;
将地面覆盖范围划分格网块;
对各格网块分别进行几何校正,具体为:
采用直接线性变换模型建立格网块中点的像点坐标和物方点坐标的关联;
根据数字表面模型数据生成地面覆盖范围内的高程最大值Hmax和高程最小值Hmin,并获得高程中值Hmid=0.5*(Hmax+Hmin);
分别以Hmax、Hmin和Hmid为高程建立三层高程面,各高程面上均均匀选取数量相等、分布相同的若干虚拟控制点;
利用自由网平差定向模型反投影计算虚拟控制点对应的像点坐标;
利用虚拟控制点,采用最小二乘法估计直接线性变换模型中的直接线性变换系数;
利用直接线性变化模型将格网块内各点逐一反投影到原始影像,得到校正后影像每一个像素点在原始影像中的坐标;
根据校正后影像每一个像素点在原始影像中的坐标和灰度进行灰度重采样获得各像素点对应的灰度值,将灰度值赋给校正后影像中对应的像素点。
2.如权利要求1所述的基于几何约束的高分辨率光学卫星凝视影像配准方法,其特征是:
步骤S100进一步包括子步骤:
S101影像两两进行组合获得影像对,对影像对进行匹配获得连接点;
S102利用附加像方误差补偿的有理函数模型,根据连接点构建自由网平差定向模型;
S103采用最小二乘法解算自由网平差定向模型的误差补偿系数。
3.如权利要求1所述的基于几何约束的高分辨率光学卫星凝视影像配准方法,其特征是:
步骤S200中,所述对影像中最优交会条件的立体影像对进行密集匹配,进一步包括:
根据立体影像对的自由网平差定向模型建立影像金字塔,由粗到精逐层对立体影像对进行像素级匹配,得到像素级匹配结果;
以像素级匹配结果为初值,通过最小二乘法对立体影像对进行匹配,得到子像素匹配结果。
4.如权利要求1所述的基于几何约束的高分辨率光学卫星凝视影像配准方法,其特征是:
步骤S300中,所述基于自由网平差定向模型和密集同名点坐标生成密集点云数据具体为:
利用自由网平差定向模型和密集同名点坐标,通过前方交会原理计算密集同名点的物方坐标(B,L,H),从而获得具有地理坐标的点云数据,即密集点云数据。
5.如权利要求1所述的基于几何约束的高分辨率光学卫星凝视影像配准方法,其特征是:
步骤S300中,所述根据粗差点剔除后的密集点云数据生成数字表面模型数据,具体为:
根据粗差点剔除后的密集点云数据获得不规则三角网形式的数字表面模型,将三角网栅格化为带地理坐标的图像格式,即获得数字表面模型数据。
6.基于几何约束的高分辨率光学卫星凝视影像配准系统,其特征是,包括:
第一模块,用来利用附加像方误差补偿的有理函数模型对影像进行自由网平差定向,获得各影像对的自由网平差定向模型;所述影像为多角度凝视序列影像;
第二模块,用来基于自由网平差定向模型,对影像中最优交会条件的立体影像对进行密集匹配,获得密集同名点坐标;所述最优交会条件的立体影像对同时满足原则(1)和原则(2):
原则(1)为:根据地形情况立体影像对的交会角为10度~45度,对平坦地区,立体影像对的交会角为30度~45度,对于地形起伏大的地区,立体影像对的交会角为10度~30度;
原则(2)为:立体影像对为影像序列中成像俯仰角大小绝对值最接近的一组影像;
第三模块,用来基于自由网平差定向模型和密集同名点坐标生成密集点云数据,对密集点云数据进行基于几何约束的粗差点剔除,根据粗差点剔除后的密集点云数据生成数字表面模型数据;
所述对密集点云数据进行基于几何约束的粗差点剔除具体为:
对每一对密集同名点对,将密集同名点的物方坐标带入所在影像对的自由网平差定向模型,求得像点像素坐标,记为(lnew,snew),计算(lnew,snew)与密集同名点匹配的像素坐标(l,s)的坐标差,若坐标差大于2倍中误差,即剔除该密集同名点;所述中误差为所有密集同名点对的坐标差的中误差;所述密集同名点对的坐标差指密集同名点对的(lnew,snew)与像素坐标(l,s)的坐标差;
第四模块,用来基于自由网平差定向模型和DSM数据,对影像依次进行几何校正;
所述对影像依次进行几何校正,具体为:
根据自由网平差定向模型和DSM数据计算影像四角点的经纬度及高程坐标,求四角点的正北外接矩形,即几何校正后影像的地面覆盖范围;
将地面覆盖范围划分格网块;
对各格网块分别进行几何校正,具体为:
采用直接线性变换模型建立格网块中点的像点坐标和物方点坐标的关联;
根据数字表面模型数据生成地面覆盖范围内的高程最大值Hmax和高程最小值Hmin,并获得高程中值Hmid=0.5*(Hmax+Hmin);
分别以Hmax、Hmin和Hmid为高程建立三层高程面,各高程面上均均匀选取数量相等、分布相同的若干虚拟控制点;
利用自由网平差定向模型反投影计算虚拟控制点对应的像点坐标;
利用虚拟控制点,采用最小二乘法估计直接线性变换模型中的直接线性变换系数;
利用直接线性变化模型将格网块内各点逐一反投影到原始影像,得到校正后影像每一个像素点在原始影像中的坐标;
根据校正后影像每一个像素点在原始影像中的坐标和灰度进行灰度重采样获得各像素点对应的灰度值,将灰度值赋给校正后影像中对应的像素点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711003290.5A CN107705329B (zh) | 2017-10-24 | 2017-10-24 | 基于几何约束的高分辨率光学卫星凝视影像配准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711003290.5A CN107705329B (zh) | 2017-10-24 | 2017-10-24 | 基于几何约束的高分辨率光学卫星凝视影像配准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107705329A CN107705329A (zh) | 2018-02-16 |
CN107705329B true CN107705329B (zh) | 2021-03-16 |
Family
ID=61182254
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711003290.5A Active CN107705329B (zh) | 2017-10-24 | 2017-10-24 | 基于几何约束的高分辨率光学卫星凝视影像配准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107705329B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10733470B2 (en) * | 2018-01-25 | 2020-08-04 | Geomni, Inc. | Systems and methods for rapid alignment of digital imagery datasets to models of structures |
CN108761271A (zh) * | 2018-03-30 | 2018-11-06 | 广州中科云图智能科技有限公司 | 一种电网树障检测方法及系统 |
CN110134914B (zh) * | 2019-05-15 | 2020-06-02 | 哈尔滨工业大学 | 一种规划敏捷卫星二维姿态机动成像任务的方法、装置及计算机存储介质 |
CN110595442A (zh) * | 2019-08-13 | 2019-12-20 | 中国南方电网有限责任公司超高压输电公司昆明局 | 输电线路通道树障检测方法、储存介质和计算机设备 |
CN111010494B (zh) * | 2019-10-28 | 2020-11-03 | 武汉大学 | 一种带有地理编码的光学卫星视频稳像方法及系统 |
CN111126148B (zh) * | 2019-11-25 | 2022-11-29 | 长光卫星技术股份有限公司 | 一种基于视频卫星影像的dsm生成方法 |
CN111044037B (zh) * | 2019-12-26 | 2021-06-18 | 中国人民解放军战略支援部队信息工程大学 | 一种光学卫星影像的几何定位方法及装置 |
CN112258422B (zh) * | 2020-08-17 | 2023-04-28 | 中国人民解放军61540部队 | 立体影像有理多项式参数(rpc)自动精化方法 |
CN112419380B (zh) * | 2020-11-25 | 2023-08-15 | 湖北工业大学 | 一种基于云掩膜的静止轨道卫星序列影像高精度配准方法 |
CN113255740B (zh) * | 2021-05-07 | 2024-04-19 | 北京市遥感信息研究所 | 一种多源遥感影像平差定位精度分析方法 |
CN112991525B (zh) * | 2021-05-07 | 2021-09-24 | 北京道达天际科技有限公司 | 像方和物方混合匹配基元的数字表面模型生成方法 |
CN114782525B (zh) * | 2022-06-22 | 2022-09-20 | 中国地质大学(武汉) | 基于全局几何约束的无人机影像定位定向方法及设备 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1987896A (zh) * | 2005-12-23 | 2007-06-27 | 中国科学院中国遥感卫星地面站 | 高分辨率sar影像配准处理方法及系统 |
CN103823981A (zh) * | 2014-02-28 | 2014-05-28 | 武汉大学 | 一种数字高程模型辅助的卫星影像区域网平差方法 |
CN104361590A (zh) * | 2014-11-12 | 2015-02-18 | 河海大学 | 一种控制点自适应分布的高分辨率遥感影像配准方法 |
CN105761248A (zh) * | 2016-01-30 | 2016-07-13 | 武汉大学 | 超大规模无控制区域网稳健平差方法及系统 |
CN106887016A (zh) * | 2017-03-09 | 2017-06-23 | 中国科学院遥感与数字地球研究所 | 一种gf‑4卫星序列图像自动相对配准方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8694258B2 (en) * | 2010-02-14 | 2014-04-08 | Vermeer Manufacturing Company | Derivative imaging for subsurface object detection |
US20150356341A1 (en) * | 2013-01-07 | 2015-12-10 | Bae Systems Plc | Fusion of multi-spectral and range image data |
-
2017
- 2017-10-24 CN CN201711003290.5A patent/CN107705329B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1987896A (zh) * | 2005-12-23 | 2007-06-27 | 中国科学院中国遥感卫星地面站 | 高分辨率sar影像配准处理方法及系统 |
CN103823981A (zh) * | 2014-02-28 | 2014-05-28 | 武汉大学 | 一种数字高程模型辅助的卫星影像区域网平差方法 |
CN104361590A (zh) * | 2014-11-12 | 2015-02-18 | 河海大学 | 一种控制点自适应分布的高分辨率遥感影像配准方法 |
CN105761248A (zh) * | 2016-01-30 | 2016-07-13 | 武汉大学 | 超大规模无控制区域网稳健平差方法及系统 |
CN106887016A (zh) * | 2017-03-09 | 2017-06-23 | 中国科学院遥感与数字地球研究所 | 一种gf‑4卫星序列图像自动相对配准方法 |
Non-Patent Citations (3)
Title |
---|
On-Orbit Geometric Calibration Using a Cross-Image Pair for the Linear Sensor Aboard the Agile Optical Satellite;Ying-Dong Pi等;《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》;20170731;第14卷(第7期);第1176-1180页 * |
TDI CCD 相机动态成像对几何质量的影响研究;岳庆兴等;《测绘科学》;20120531;第37 卷(第3 期);第14-17页 * |
一种光学遥感卫星多相机成像系统的高精度影像拼接方法;程宇峰等;《光学学报》;20170831;第37卷(第8期);第1-10页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107705329A (zh) | 2018-02-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107705329B (zh) | 基于几何约束的高分辨率光学卫星凝视影像配准方法 | |
CN109903352B (zh) | 一种卫星遥感影像大区域无缝正射影像制作方法 | |
CN110388898B (zh) | 构建虚拟控制点约束的多源多重覆盖遥感影像平差方法 | |
CN107767440B (zh) | 基于三角网内插及约束的文物序列影像精细三维重建方法 | |
Oh et al. | A piecewise approach to epipolar resampling of pushbroom satellite images based on RPC | |
CN108592884B (zh) | 一种通用线阵卫星核线影像生成方法 | |
CN106780321B (zh) | 一种cbers-02卫星hr传感器影像整体严密定向与纠正拼接方法 | |
CN113358091B (zh) | 一种利用三线阵立体卫星影像生产数字高程模型dem的方法 | |
Li et al. | A new analytical method for estimating Antarctic ice flow in the 1960s from historical optical satellite imagery | |
CN112070891B (zh) | 数字地面模型作为三维控制的影像区域网平差方法及系统 | |
CN113899387A (zh) | 一种基于验后补偿的光学卫星遥感影像区域网平差方法及系统 | |
CN108562900B (zh) | 一种基于高程校正的sar图像几何配准方法 | |
CN108594255B (zh) | 一种激光测距辅助光学影像联合平差方法及系统 | |
Jacobsen | Problems and limitations of satellite image orientation for determination of height models | |
CN110853140A (zh) | 一种dem辅助的光学视频卫星影像稳像方法 | |
CN108447100B (zh) | 一种机载三线阵ccd相机的偏心矢量和视轴偏心角标定方法 | |
CN111508028A (zh) | 光学立体测绘卫星相机的自主在轨几何定标方法及系统 | |
CN108830921A (zh) | 基于入射角度的激光点云反射强度纠正方法 | |
CN110631555A (zh) | 基于二阶多项式无控制点区域网平差的历史影像正射纠方法 | |
KR100870894B1 (ko) | 선형 스캐닝 영상의 자동 영상기하 보정 방법 | |
CN109029379B (zh) | 一种高精度小基高比立体测绘方法 | |
CN109579796B (zh) | 一种投影后影像的区域网平差方法 | |
CN115326025B (zh) | 一种用于海浪的双目影像测量与预测方法 | |
Takaku et al. | RPC generations on ALOS prism and AVNIR-2 | |
CN107146281B (zh) | 一种月球表面高分辨率dem提取方法 |
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 |