[发明专利]一种基于SAR的地表监测图像处理方法

摘要

一种基于SAR的地表监测图像处理方法，包括分别获取同一时刻和同一地表监测区域的SAR图像，进行边缘轮廓提取，进行滤波预处理，像素处理后得到平均像素，利用最小二乘法选出预处理SAR图像，将预处理SAR图像中的所有像素对应求和后取平均值，作为处理后SAR图像中的对应像素，重构后得到处理后的SAR图像，再获取不同时刻和同一地表监测区域SAR图像，得到所有时刻的经过处理后的SAR图像，本方法处理方式简单，计算量小，匹配对准容易，并且经过处理后的地表监测图像更加精确，对于土地情况的地表监测更加准确。

主权项

一种基于SAR的地表监测图像处理方法，其特征在于，依次包括以下步骤：(1)分别获取同一时刻和同一地表监测区域的分辨率为N×M的5幅SAR图像，采用边缘检测算法分别对5幅SAR图像进行边缘轮廓提取，将边缘轮廓提取部分，分别进行滤波预处理；(2)将分别进行滤波预处理过的边缘轮廓提取部分中的每一个像素与其邻近的8个像素构成3×3的像素块，按照如下公式分别得到处理后的每一个像素：$P_{\mathrm{IJ}}=(\frac{\frac{P_{(i-1)(j-1)}+{2P}_{(i-1)j}+P_{(i-1)(j+1)}}{4}+\frac{P_{(i-1)(j-1)}+{2P}_{i(j-1)}+P_{(i+1)(j-1)}}{4}+\frac{P_{(i+1)(j-1)}+{2P}_{(i+1)j}+P_{(i+1)(j+1)}}{4}+\frac{P_{(i-1)(j+1)}+{2P}_{i(j+1)}+P_{(i+1)(j+1)}}{4}}{4}+P_{\mathrm{ij}})/2$其中，P_ij为对应处理前第i行、第j列的像素值，P_IJ为对应经过处理后的第i行、第j列的像素值，其中i＝1,2,…,N,j＝1,2,…,M；(3)将与5幅SAR图像对应的经过处理后的第i行、第j列的像素值P_IJ进行平均值计算，得到平均像素(4)利用最小二乘法，分别将与5幅SAR图像对应的经过处理后的第i行、第j列的像素值P_IJ与平均像素相减后作平方，求得每一幅SAR图像对应的误差平方和，选出误差平方和最小的2幅作为预处理SAR图像；(5)将所述预处理SAR图像中的所有像素对应求和后取平均值，作为处理后SAR图像中的对应像素，重构后得到处理后的SAR图像；(6)分别获取不同时刻和同一地表监测区域的分辨率为N×M的5幅SAR图像，重复步骤(1)‑(5)得到所有时刻的经过处理后的SAR图像。