[发明专利]基于相位相关的高重叠度遥感影像满度连接点匹配方法

权利要求书

1.一种基于相位相关的高重叠度遥感影像满度连接点匹配方法，其特征在于所述方法包括：

步骤一、输入多张重叠影像，构建匹配模型，对基准影像进行分块和特征点提取；

步骤二、基于RPC成像模型和DEM获得特征点的经纬度坐标和高程；

步骤三、以特征点的经纬度坐标为中心进行物方空间重采样；

步骤四、相位相关计算匹配窗口的偏移量获得匹配点对；

步骤五、将匹配点的坐标通过RPC模型反算，计算匹配点对在像方的行列号坐标；

步骤六、剔除粗差点得到最终匹配结果；

其中，所述步骤一具体为选取其中一张作为基准影像，其余的影像作为目标影像；对基准影像进行分块，得到N×M个大小相同的影像区域；在每个影像块的中心选取一个S×S的窗口进行特征点提取，保留n个特征点作为候选点；

其中，所述步骤二具体为：设特征点在基准影像上的行列号坐标为(l,s)，由RPC模型和DEM数据，通过迭代计算可以得到特征点对应于物方的经纬度和高程(B,L,H)；

其中，RPC模型的表示如下：

其中，

Num_L(U,V,W)＝a₁+a₂V+a₃U+a₄W+a₅VU+a₆VW+a₇UW+a₈V²+a₉U²+a₁₀W²+a₁₁UVW+a₁₂V³+a₁₃VU²+a₁₄VW²+a₁₅V²U+a₁₆U³+a₁₇UW²+a₁₈V²W+a₁₉U²W+a₂₀W³

Den_L(U,V,W)＝b₁+b₂V+b₃U+b₄W+b₅VU+b₆VW+b₇UW+b₈V²+b₉U²+b₁₀W²+b₁₁UVWb₁₂V³+b₁₃VU²+b₁₄VW²+b₁₅V²U+b₁₆U³+b₁₇UW²+b₁₈V²W+b₁₉U²W+b₂₀W³

Num_S(U,V,W)＝c₁+c₂V+c₃U+c₄W+c₅VU+c₆VW+c₇UW+c₈V²+c₉U²+c₁₀W²+c₁₁UVW+c₁₂V³+c₁₃VU²+c₁₄VW²+c₁₅V²U+c₁₆U³+c₁₇UW²+c₁₈V²W+c₁₉U²W+c₂₀W³

Den_S(U,V,W)＝d₁+d₂V+d₃U+d₄W+d₅VU+d₆VW+d₇UW+d₈V²+d₉U²+d₁₀W²+d₁₁UVW+d₁₂V³+d₁₃VU²+d₁₄VW²+d₁₅V²U+d₁₆U³+d₁₇UW²+d₁₈V²W+d₁₉U²W+d₂₀W³

其中，a_i,b_i,c_i,d_i(i＝1,2…20)为有理多项式系数,LineOff,SampleOff,LonOff,LatOff,HeiOff,LineScale,SampleScale,LonScale,LatScale,HeiScale为模型已知的重心化系数；

其中，所述步骤三根据所述步骤二得到特征点在物方空间中的经纬度坐标(B，L)，对于每一幅输入的目标影像，由该目标影像与基准影像的最低分辨率及匹配窗口大小确定物方的经纬度范围BL，BR，LT，LB，由四个角点的经纬度坐标内插得到匹配窗口内每一个点的经纬度坐标，结合DEM和RPC模型计算匹配窗口每一个点的像方坐标，通过双线性内插获取匹配窗口像素的灰度。