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

摘要

一种机载InSAR困难区域的DEM重建方法，包括：对各个方向照射的InSAR数据，分别进行初步干涉处理，生成主、辅图像、干涉相位图及相干系数图；将待测绘区域投影到高斯坐标系下，并划分像素网格；将粗DEM数据插值使其与上述像素网格一一对应；将粗DEM值为待测量高程的初值，设置该像素网格高程值的变化范围；设置每一次迭代处理的高程值；分别计算该像素网格投影到各个方向对应的主图像斜距平面中的坐标，读取对应坐标处的干涉相位和相干系数观测值，并计算迭代高程值对应的理想干涉相位；计算迭代高程值的联合似然函数；取使联合似然函数最大时对应的迭代高程值即为该像素网格的最终高程估计值，即重建得到了每个像素网格的DEM，完成机载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重建。