[发明专利]平面地表环境中地下建筑的红外成像探测定位方法

权利要求书

1.一种平面地表环境中地下建筑的红外成像探测定位方法，其特征在于，所述方法根据地下建筑的能量扩散的高斯模型，对地下建筑经地层调制后所形成的原始红外图像进行反调制处理，得到地下建筑的目标图像，所述方法包括以下步骤：

（1）获取地下建筑经地层调制后所形成的原始红外图像g₀，并确定地下建筑在原始红外图像g₀中的大体位置的局部红外图像g；

（2）设置迭代终止条件，并设定高斯热扩散函数的初始值h₀；

（3）以所述局部红外图像g作为初始目标图像f₀，根据所述高斯热扩散函数的初始值h₀，利用最大似然估计算法迭代求解热扩展函数h_n和目标图像f_n；

（4）判断是否满足迭代终止条件，如果满足，则本次迭代求解得到的目标图像f_n即为最终的目标图像f；若不满足，则返回步骤（3），继续迭代计算。

2.如权利要求1所述的方法，其特征在于，所述步骤（1）具体为：

（1.1）将原始红外图像g₀分成N个m×m像素大小的红外图像区域；

（1.2）计算得到N个红外图像区域的平均灰度值，其中第i幅红外图像区域的平均灰度值为g_ij为第i幅红外图像区域中的第j个像素的灰度值；

（1.3）分别计算出相邻红外图像区域的平均灰度差值ΔG＝G_i+1-G_i；

（1.4）将平均灰度值大且相邻平均灰度值差值ΔG小的区域作为需要反调制处理的局部红外图像g。

3.如权利要求1或2所述的方法，其特征在于，所述步骤（3）中迭代求解热扩展函数h_n和目标图像f_n具体根据下面两式迭代计算：

hn+1(x,y)={[g(i,j)hn(x,y)*fn(i,j)]*fn(-i,-j)}*hn(x,y)]]>

fn+1(i,j)={[g(i,j)hn+1(x,y)*fn(i,j)}]*hn+1(-x,-y)}*fn(i,j)]]>

其中n为当前迭代次数,(x,y)属于热扩散函数支持域，(i,j)属于图像支持域，f(i,j)表示目标图像,g(i,j)表示局部红外图像,h(x,y)是热扩散函数，h(-x,-y)表示h(x,y)的共轭，f(-i,-j)表示f(i,j)的共轭，*为卷积运算符。

4.如权利要求3所述的方法，其特征在于，所述步骤（2）中的迭代终止条件为迭代终止次数n＞N₀或误差ε，所述步骤（4）中判断是否满足迭代条件具体为：

判断是否满足|g-hⁿ⁺¹*fⁿ⁺¹|＜ε或者n＞N₀，如果二者中任一个满足，则满足迭代终止条件，否则不满足。