[发明专利]一种平面地形中地下管状设施红外成像探测定位方法

权利要求书

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

（1）获取地下管状设施经地层调制后所形成的原始红外图像g；

（2）根据所述原始红外图像g的热辐射能量带的宽度，设定高斯热扩散函数的初始值h₀，并设置迭代终止条件；

（3）以所述原始红外图像g作为初始目标图像f₀，根据所述高斯热扩散函数的初始值h₀，利用基于贝叶斯理论的单帧图像盲目去卷积方法迭代求解热扩展函数h_n和目标图像f_n；

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

2.如权利要求1所述的方法，其特征在于，所述步骤（3）中利用基于贝叶斯理论的单帧图像盲目去卷积方法迭代求解热扩展函数h_n和目标图像f_n具体包括：

在第n次盲目迭代时首先按照下面两式进行Richardson－Lucy迭代计算，

hm+1n(x)={[g(x)hmn(x)*fn-1(x)]*fn-1(-x)}hmn(x),]]>

fm+1n(x)={[g(x)hn(x)*fmn(x)]*hn(-x)}fmn(x),]]>

其中，m表示当前的Richardson－Lucy迭代次数，m＝1□M，M表示总的Richardson－Lucy迭代次数；

经过M次Richardson－Lucy迭代后，得到和再将和作为下式中的hⁿ(x)和fⁿ(x)，求解hⁿ⁺¹(x)和fⁿ⁺¹(x)：

hn+1(x)={[g(x)hn(x)*fn(x)]*fn(-x)}nhn(x)]]>

fn+1(x)={[g(x)hn+1(x)*fn(x)]*hn+1(-x)}nfn(x)]]>

其中n为当前盲目迭代次数,f(x)表示目标图像,g(x)表示原始红外图像,h(x)是热扩散函数，h(-x)表示h(x)的共轭，f(-x)表示f(x)的共轭，*为卷积运算符。

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

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