CN114609630B 一种基于导航卫星双基地InSAR的DEM成像方法及装置 （北京理工大学）.docxVIP

(19)国家知识产权局

(12)发明专利

(10)授权公告号CN114609630B(45)授权公告日2025.07.08

(21)申请号202210120140.7

(56)对比文件

(22)申请日2022.01.30

CN108681621CN109035312

A,2018.10.19

A,2018.12.18

(65)同一申请的已公布的文献号

申请公布号CN114609630A

审查员曾思萌

(43)申请公布日2022.06.10

(73)专利权人北京理工大学

地址100081北京市海淀区中关村南大街5

号

(72)发明人刘飞峰高检王战泽徐智祥周经天

(74)专利代理机构北京理工大学专利中心

11120

专利代理师温子云李爱英

(51)Int.CI.

G01S13/90(2006.01)

权利要求书3页

说明书8页附图7页

(54)发明名称

一种基于导航卫星双基地InSAR的DEM成像方法及装置

(57)摘要

CN114609630B本发明公开了一种基于导航卫星双基地InSAR的DEM成像的方法及装置，能够提高PS点的时间相关性，提升了系统的配准精度。所述方法包括：使用滑动最小二乘拟合方法拟合本地坐标系下的DEM;对拟合后的DEM的各项进行二次差分，找到所述拟合后的DEM存在的所有不连续单元；得到拟合后的不连续单元相对应的待平滑核；对每个待平滑核，生成与所述待平滑核对应的、以所述待平滑核为基准进行中心对称变换后的高斯平滑核；对每个待平滑核，将其与对应的高斯平滑核进行加权求和运算，进而得到每个待平滑核对应的不连续单元的拟合平滑结果；将所有不连续单元的拟合平滑结果作为最终平滑后的DEM,将所述最终平滑后的

CN114609630B

使用滑动最小二乘拟合方法拟合本地坐标系下的

使用滑动最小二乘拟合方法拟合本地坐标系下的DE

对拟合后的DEM的各项进行二次差分，找到所述拟合后的DEM存在的所有不连续单元

对每个拟合后的不连续单元，以所述拟合后的不连续单元为中心进行边界扩展，扩展后的结果作为与所述拟合后的不连续单元相对应的待平滑核

对每个待平滑核，生成与所述待平滑核对应的、以所述待平滑核为基准进行中心对称变换后的高斯平滑核

对每个待平滑核，将其与对应的高斯平滑核进行加权求和运算，进而得到每个待平滑核对应的不连续单元的拟合平滑结果

将所有不连续单元的拟合平滑结果作为最终平滑后、的DEM,将所述最终平滑后的DEM作为后向投影成像平面进行成像，得到最终的成像结果

结束

CN114609630B权利要求书1/3页

2

1.一种基于导航卫星双基地InSAR的DEM成像方法，其特征在于，所述方法包括：

步骤S1:使用滑动最小二乘拟合方法拟合本地坐标系下的DEM;

步骤S2:对拟合后的DEM的各项进行二次差分，找到所述拟合后的DEM存在的所有不连续单元；

步骤S3:对每个拟合后的不连续单元，以所述拟合后的不连续单元为中心进行边界扩展，扩展后的结果作为与所述拟合后的不连续单元相对应的待平滑核；

步骤S4:对每个待平滑核，生成与所述待平滑核对应的、以所述待平滑核为基准进行中心对称变换后的高斯平滑核；

步骤S5:对每个待平滑核，将其与对应的高斯平滑核进行加权求和运算，进而得到每个待平滑核对应的不连续单元的拟合平滑结果；

步骤S6:将所有不连续单元的拟合平滑结果作为最终平滑后的DEM,将所述最终平滑后的DEM作为后向投影成像平面进行成像，得到最终的成像结果；

所述步骤S2,包括：对DEM全局拟合结果进行二次差分处理，并以所述DEM全局拟合结果的均值作为阈值；对DEM的N个单元进行遍历，找到二次差分结果大于所述阈值的单元，确定为拟合不连续单元；

所述二次差分处理为：求所述DEM全局拟合结果的二阶偏导数，对x方向求偏导或对y方向求偏导。

2.如权利要求1所述的方法，其特征在于，所述步骤S1:使用滑动最小二乘拟合方法拟合本地坐标系下的DEM,包括：

步骤S11:DEM的N个单元，记为F=[f(x?,y?)f(x?,y)…f(x,y)],

其中，f(x:,y)是DEM的第i单元，i=1,2,…,N;x:为第i单元的横坐标，y;为第i单元的纵坐标；

步骤S12:对于DEM的N个单元中的每一个f(x,,y),采用二阶加权最小二乘拟合，包括：

对于第i单元f(x