当前位置:太阳集团娱乐网址 > 太阳集团娱乐网址 >

一种高精度小基高比立体测绘方法与流程

一种高精度小基高比立体测绘方法与流程

分类:

  小基高比立体测绘方法与传统大基高比立体测绘方法完全不同。大基高比 是通过图像匹配在两幅图像中找到同名点,进而将同名控制点坐标代入共线方 程求得共线方程参数,将待求点像方坐标代入已求解好的共线方程中,求得待 求点的物方三维坐标。而小基高比立体测绘方法的核心计算公式为像方相对高 程等于两幅图像的平面视差除以基高比数值。早在2002年,法国人就利用 SPOT5卫星上的全色影像和多光谱影像组成立体像对,对小基高比情况下获取 数字高程模型进行了试验,结果表明在立体交会角仅为0.02时,仍能获取一定 精度的高程信息,从而验证了小基高比立体测绘技术的可行性。

  传统的大基高比立体测绘方法是逐点代入共线方程,逐点获取相应点处的 高程值;单台相机的大基高比立体测绘对平台的姿控提出了更高的要求,要卫 星平台灵活、敏捷度高、稳定性强、姿态机动能力强,实现难度更大;大基高 比条件下,拍摄城市区域(高大建筑物密集)时,丢失高程信息的遮挡区域将 增大,不利于城市立体测绘。

  2007年第28卷《Journal of Mathematical Imaging and Vision》上,法国 空间研究中心的Julie Delon等人发表的《Small Baseline Stereovision》论文 通过理论推导和仿真分析,讨论了小基高比立体测绘方法的可行性及相关仿真 实验验证,其不足之处在于:该方法仅通过两幅图像的平面视差除以基高比数 值求得像方相对高程,没有对相机的内外方位元素进行校正,没有对原始采集 到的图像进行合理的MTF补偿,不能计算出物方绝对高程。

  本发明解决的技术问题是:克服现有技术的不足,提供了一种高精度小基 高比立体测绘方法,该方法提出了必须对采集到的图像进行MTF补偿预处理, 并校正相机的内外方位元素,最后通过控制点计算参考面的绝对高程来计算整 个重叠区域的物方绝对高程值,实现最终的地物高精度立体测绘。

  1)相机在轨道高度为H,基线长度为B的两个位置对地进行拍摄,采集 获取两幅被噪声污染后的退化图像,并建立所采集图像的退化模型 其中,g为采集到的被噪声污染后的退化图像,为傅里叶 逆变换算子,f为理想图像,n为噪声,H′为频域晶胞上的调制传递函数,*为 卷积运算;

  2)对两幅被噪声污染后的退化图像分别进行基于频域晶胞的总变分正则化 MTFC预处理,得到放大后的两幅复原图像f1、f2;

  5)采用基于解最优化的相位相关法对校正了外方位元素畸变的两幅图像 fext1、fext2进行亚像素级匹配,得到匹配后的视差图V;

  其中,s(a1)、s(c1)分别为a1、c1点视差值,a1、c1为的地面任意两点A和C在 第一幅图像上的成像点,P为探测器像元大小,f′为光学系统焦距,为物方地面采样距离;

  7)通过约束K个控制点计算的绝对高程与真实高程的误差平方和最小,求 解参考面高度h,即地面点C处的绝对高程为

  所述步骤2)进行基于频域晶胞的总变分正则化MTFC预处理的具体方法 为:

  其中,α为正则化参数,β为可调参数,Df是f的支持域, 即成像系统的频域晶胞;表示散度,对于二维向量f=[fx,fy]T,

  其中,为点(ν,ω)处的系统调制传递函数数值,为系统调制传递函数矩 阵,其大小为K×L,K、L均为正整数;comba(v,w)表示采样边带;θalias为阈值, 其取值范围为[1,0);

  所述步骤3)计算得到校正内方位元素畸变后的两幅图像fcal1、fcal2的具体 过程为:

  32)对两幅复原后的图像f1、f2进行校正,计算每幅图像校正后的图像 坐标(xcal,ycal)对应的原坐标(x,y),即

  33)经过双线性内插算法求得校正后坐标处(xcal,ycal)的像素值;

  34)由上面的过程最终实现对复原后图像f1、f2的内方位元素校正,得 到校正畸变后的两幅图像fcal1、fcal2。

  所述步骤4)计算得到校正外方位元素畸变后的两幅图像fext1、fext2的具体 方法为:

  其中,x、y分别为内方位元素校正后图像fcal1上像点的横、纵坐标; a1,a2,a3,b1,b2,b3,c1,c2,c3为第一幅图像的方向余弦;取xk、yk为整数,求得一系列 的像点坐标(x,y);

  对于图像fcal2,将y′k=yk代入下面的第一幅图像共线方程,求得对应外方位 元素校正后坐标(x′k,y′k)的原坐标值(x′,y′):

  其中,x′k为整数,x′、y′为内方位元素校正后图像fcal2上像点的横、纵坐标;a′1, b′1,…,c′3为第一幅图像的方向余弦,分别是第二幅图像相对于摄影基线的角 方位元素的函数;

  42)经过双线性内插算法求得校正后坐标处(xk,yk)、(x′k,y′k)的像素值;

  43)将校正过内方位元素的图像fcal1、fcal2上的所有点进行步骤41)、42) 的核线重采样操作,得到校正了外方位元素畸变的两幅图像fext1、fext2。

  51)通过求解下面的最优化问题来得到目标区域中心点的相对平移量 (Δx,Δy),即

  其中,Q为含有噪声的局部互相关功率谱,表示理论局部互相关功率 谱,W是权重矩阵;wx、wy分别为x方向和y方向的权重系数,φ(Δx,Δy)局部互 相关功率谱误差函数;

  判断是否小于条件阈值,如果小于某一阈值则迭代结束,得到不同迭 代次数对应的平移量(Δxi,Δyi),若不收敛,调整阈值再进行步骤54)迭代计算;

  所述步骤31)计算得到相机镜头x、y两个方向上的几何畸变量Δx,Δy的 具体方法为:

  其中,x、y为退化图像上像点的横、纵坐标;X、Y、Z为地面坐标系, a1,a2,a3,b1,b2,b3,c1,c2,c3为通过相机的姿态构建的旋转矩阵,x0、y0为x、y两个方 向上的主点偏移,XS,YS,ZS为相机的轨道参数;Δx、Δy为x、y两个方向上的几 何畸变量。

  所述步骤33)经过双线性内插算法求得校正后坐标处(xcal,ycal)的像素值 的具体方法为:

  其中,Δd为原始变形图像的采样间距,g′为原坐标(x,y)处待求像素灰度值, g1,g2,g3,g4分别为其周围4个像元的灰度值。

  所述步骤42)经过双线性内插算法求得校正后坐标处(xk,yk)、(x′k,y′k)的像 素值的具体方法为:

  其中,Δd′为原始变形图像的采样间距,g″为原坐标(x′,y′)处待求像素灰度 值,g′1,g′2,g′3,g′4分别为原坐标(x′,y′)周围4个像元的灰度值;同理求出第二幅图 像坐标处(x′k,y′k)的像素值。

  所述步骤1)中相机在轨道高度为H,基线长度为B的两个位置对地进行 拍摄,两个位置关于光轴方向对称。

  现有方法只是核线重采样后,通过两幅图像的平面视差除以基高比数值求 得像方相对高程。而本发明通过限制复原后的混叠效应,构建了频域晶胞模型, 利用频域晶胞上降晰函数进行反卷积,实现了采集图像无虚假添加的MTF补 偿,弥补相机自身的不足;通过内外方位元素重采样进行系统几何误差校正, 提高小基高比立体测绘精度;通过基于解最优化的相位相关匹配法,实现亚像 素级的平面视差计算,提高最终的测绘精度;通过计算控制点计算参考面的绝 对高程来计算整个重叠区域的物方绝对高程值,实现最终的地物高精度立体测 绘。

  小基高比立体测绘图像预处理主要进行了MTF补偿,本发明提出了基于频 域晶胞模型的总变分正则化MTFC算法。首先,整个采样系统的模型可以表示 为

  其中,g为采集到的被噪声污染后的退化图像矩阵,为傅里叶逆变换算 子,f为理想图像矩阵,n为噪声矩阵,H′为降晰函数傅里叶变换的模,*为 卷积运算。基于频域晶胞模型的总变分正则化MTFC算法的技术途径如下:

  其中,β0是可调参数。利用 该总变分最小化模型进行去噪,该模型迭代计算公式为

  其中,为系统MTF,comba(v,w)表示采样边带,θalias为阈值,其取值 范围为[1,0)。

  c)小基高比立体测绘中需要利用频域晶胞进行H′的计算,由步骤b)中频 域晶胞的计算公式可以得到与相应成像系统MTF最为匹配的频域晶胞,利用频 域晶胞上的MTF进行图像复原,则有

  基于频域晶胞的图像复原方法与传统方法的最大不同在于式(1)中的H′由 式(3)计算求得。

  由步骤a)的迭代计算公式可得到的值,将上式计算的H′代人,太阳集团娱乐网址 进行反卷积计算,求解复原图像f1、f1。

  物理验证样机的试验数据在MTFC前最优高程精度为5.66个GSD,利用 基于频域晶胞的总变分正则化MTFC方法的最优高程精度为4.76个GSD,基 于频域晶胞的MTFC方法处理后的高程精度,相比于MTFC前提高了16.1%。

  在小基高比立体测绘中,内方位元素校正主要是对由相机自身特性引起的 几何变形进行校正。考虑改正项后的共线方程如式:

  其中,x,y为像点坐标,X,Y,Z为地面点坐标,a1,a2,a3,b1,b2,b3,c1,c2,c3为 通过相机的姿态构建的旋转矩阵,x0,y0为主点偏移,f′为光学系统焦距, XS,YS,ZS为相机的轨道参数。Δx,Δy为附加参数,主要由相机镜径向畸变系数和 切向畸变系数、CCD安装误差系数等组成,主要依据实际定标过程中考虑哪些 误差项。

  通过精确测量墙面214个控制点的物方精确坐标,与图像上像方坐标匹配 成控制点对,M个控制点对就有关于共线方程系数的M对方程组,通过对方程 组求解得到修正后的共线方程参数,进而得到内方位元素。如图所示,为通过 计算墙面靶标得到的相机几何畸变图。

  由通过共线方程计算出的几何畸变量Δx,Δy,对两幅复原后的图像f1、f2进 行校正,计算每幅图像校正后的图像坐标(xcal,ycal)对应的原坐标(x,y),即

  变形图形坐标不一定为整数,必须进行灰度重采样得到变形坐标点(x,y)的 灰度值作为对应校正后坐标的灰度值。如图3所示,Δd为原始变形图像的采样 间距,待求像素灰度值为g′,可由坐标点(x,y)周围4个像元的灰度值g1,g2,g3,g4, 经双线)核线重采样

  核线在摄影测量中是一个重要概念。两个相机光心的连线称为基线,通过 基线的平面(即核面)与像片的交线为核线,同一核面与左右像片的交线称为 同名核线。由此可知任意物点一定位于通过该点的核面与像片的交线上,影像 校正可以使匹配点搜索位于同名核线上,若如图所示基线平行于行方向,则平 行核线重采样后,行就是核线,通过摄影基线OO和物点A的平面为通过像点 a的核平面,图4中I代表位于左方的航摄像片,k代表相应的平行于摄影基线 的水平像片,ak为A点在左水平像片上的相应像点。设a、ak在各自的像片坐 标系中的坐标分别为(x,y)和(xk,yk),则有

  其中,a1、a2、…、c3为第一幅图像的方向余弦,它们是这张像片相对于 摄影基线的角方位元素的函数。

  在相对水平像片上,当yk=α(α为正整数),则为核线。将yk=α代入式 (7),且xk值取正整数,即求得一系列的像点坐标(x1,y1)、(x2,y2)、(x3,y3)、…, 这些像点就位于倾斜像片I的核线上,将这些像点的灰度g′(x1y1),g′(x2y2),…… 直接赋给相对水平像片上相应的像点,就能获得相对水平像片上的核线影像。

  由于在相对水平像片上,同名核线yt的坐标值相等,因此将同样的y′t=α代 入第二幅图像共线方程:

  其中,a′1,b′1,…,c′3为第二幅图像的方向余弦,分别是第二幅图像相对 于摄影基线的角方位元素的函数。由此即得第二幅图像上的同名核线。

  图像是规则排列的灰度格网序列,为了获取核线的灰度序列,必须对原始 图像灰度进行重采样。

  用小基高比(小于0.1)立体测绘可以克服大基高比遮挡和辐射差异大的 缺点。在小基高比条件下,两幅图像的差异性小,图像匹配的相关性变大,可 以显著提高同名点配准精度。但为了获得同样的高程精度,要求同名点配准精 度有一定数量级的提高(一般要求优于1/20像元以上)。太阳集团娱乐网址因此,研究小基高比 成像条件下亚像素级的图像处理匹配方法是一项关键技术。本项目提出并采用 了基于解最优化的相位相关法进行图像匹配,基本原理是构建优化目标函数

  其中,Q为含有噪声的局部互相关功率谱,表示理论局部互相关 功率谱,W是权重矩阵。通过求解下面的最优化问题来得到目标区域中心点的 相对平移量(Δx,Δy),即

  式中,∑1=diag(σ1,σ2,…,σr),r=rank(Q)。用最大奇异值和相应 的奇异量来估计相位相关矩阵Q,进而通过奇异值分解得到qx,qy,进而求解得 到相应的相位px,py。用最小二乘法拟合px,py,得px=kxx+bx,py=kyy+by。则

  ),那么可以从下面的定义计算Q1(wx,wy),即c)由下面的迭代计算公式计算Qi+1(wx,wy)和

  判断是否小于某一阈值,如果小于某一阈值则迭代结束,否则返回步 骤a);d)最终的平移量由下式计算

  全局匹配的仿真试验结果表明,全局像差估计比较精确,可以达到1/100 个像素精度。但用于局部像素点平移量预估时会受到局部噪声的干扰比较大, 因此在应用该方法进行亚像素级匹配时,一方面要对遥感器的噪声进行控制, 另一方面还需对匹配数据进行抑噪预处理。

  由上述算法进行图像匹配的结果如图5(b)所示。视差图中颜色越浅代表 高度越高。图5(a)为原参考图。

  上述基于解最优化的相位相关匹配算法,可以实现小基高比成像条件下的 亚像素级图像匹配。然而还需要针对遥感图像进行亚像素级匹配验证,加快算 法的处理速度,使其适用于航天立体测绘。

  本节主要工作是基于上节匹配的视差图求取相对高程值。视差图上的像素 点a1表征参考图上点a1与待匹配图像上对应匹配点a2的像方像素差。假设a1点 视差值为s(a1),探测器像元大小为P,A点像方坐标差为

  如图6所示为将实验所得视差图进行相对像方高程计算,得到的三维重建 效果图。

  本节主要阐述基于控制点的相对高程值计算绝对高程的方法。假设已知K 个控制点的真实高程值相同参考点的高度为h,则计算绝对 高程的模型为

  通过约束K个控制点计算的绝对高程与真实高程的误差平方和最小求解参 考点高度,即

文章标签: 太阳集团娱乐网址 ,基一高比

上一篇:大连化物所MXene基高比能超级电容器研究取得新进展

下一篇:绵竹沃尔沃挖掘机维修油温高比憋机