[发明专利]机载InSAR困难区域的DEM重建方法

权利要求书

1.一种机载InSAR困难区域的DEM重建方法，包括：

步骤A：对各个方向照射的InSAR数据，分别进行初步的干涉处理，生成主图像和辅图像、干涉相位图以及相干系数图；

步骤B：将待测绘DEM的场景区域四角点的地理坐标，利用高斯投影的方法投影到高斯坐标系下，并将高斯投影后高斯坐标系下的矩形区域按照生成DEM产品要求的分辨率划分为等间隔的像素网格；

步骤C：选择待测绘DEM的场景区域对应的外部粗DEM数据，根据所选粗DEM数据的格式将粗DEM数据转换到高斯坐标系下，并将所述粗DEM数据插值到生成产品要求的分辨率，使粗DEM数据与步骤B所确定的等间隔的像素网格一一对应；

步骤D：对步骤B所划分的每一个像素网格，以该像素网格对应的粗DEM值为待测量高程的初值，设置该像素网格高程值的变化范围；

步骤E：对步骤D所处理后的每一个像素网格，设置每一次迭代处理的高程值；

步骤F：对于步骤E所述的每次迭代的高程值，根据各个方向多次照射的InSAR数据主图像的成像几何，分别计算该像素网格投影到各个方向对应的主图像斜距平面中的距离向和方位向坐标，并读取步骤A所生成的干涉相位图中该像素网格对应坐标处的干涉相位观测值，以及步骤A所生成的相干系数图中该像素网格对应坐标处的相干系数观测值；

步骤G：对于步骤E所述的每次迭代的高程值，根据各个方向多次照射的InSAR数据的成像几何，分别计算所述高程值对应的理想的干涉相位；

步骤H：对于步骤E所述的每次迭代的高程值，针对各个方向照射的InSAR数据，分别计算在干涉相位观测值一定的前提下关于所述高程值的似然函数；

步骤I：对于步骤E所述的每次迭代的高程值，计算各个方向多次照射的InSAR数据在联合的干涉相位观测值一定的前提下关于迭代高程值的联合似然函数；以及

步骤J：迭代次数每次增加1，重复步骤E至步骤I，取使联合似然函数最大时对应的迭代高程值即为该像素网格的最终高程估计值，即重建得到了步骤E所述的每个像素网格的DEM，完成机载InSAR困难区域的DEM重建。

2.根据权利要求1所述的机载InSAR困难区域的DEM重建方法，其中，所述步骤A，包括：

步骤A1：对各个方向照射的InSAR数据生成的主、辅图像进行复图像配准；

步骤A2：对步骤A1配准后的主、辅图像进行共轭相乘，对共轭相乘的结果取相位，生成干涉相位图；以及

步骤A3：对步骤A1配准后的主、辅图像计算每一个像素的相干系数，生成相干系数图。

3.如权利要求1所述的机载InSAR困难区域的DEM重建方法，其中，步骤C所述粗DEM数据包括：SRTM DEM数据或ASTER DEM数据。

4.根据权利要求1所述的机载InSAR困难区域的DEM重建方法，其中，所述步骤E在[h_init-Δh_range，h_init+Δh_range]范围内每次取高程值为h_j＝h_init-Δh_range+(j-1)·δh，其中，步骤B所划分的每一个像素网格均有一个高斯坐标(X_gauss，Y_gauss)，每一个像素网格的高斯坐标(X_gauss，Y_gauss)处对应的粗DEM值为Z_coarse，Z_coarse为待测量高程的初值h_init，Δh_range的取值为所选择的粗DEM数据的精度，δh为高程值的变化步长，j为迭代次数，表示向上取整。