[发明专利]一种适用于中子成像系统图像的基于高斯‑泊松混合噪声模型的图像复原方法

说明书

技术领域

本发明属于信息处理技术领域，涉及一种图像降噪方法。

背景技术

中子成像技术是用于非破坏性测试及其他许多基础研究领域的一种有效工具。在可移动中子成像系统中，由于高达10⁹n/s的中子发生器，小于15m的校直系统，慢中子流和受限的长度/直径比值，所得图像通常伴有严重的高斯泊松混合噪声。图像噪声会造成：①图像主观视觉效果差；②图像特征被掩盖，降低图像质量和精度，影响对图像进行特征提取和分析等后期工作。故对于中子成像系统，图像降噪是一个至关重要的后期处理过程。

图像的信息主要分为高频和低频两个部分，高频部分是图像中像素亮度值变换较快的部分，主要包含图像中的轮廓信息，表现图像的主体特征。低频部分是图像中较为平滑的部分，表现了图像的细节特征。对于噪声而言，也属于图像中高频的部分。由于图像边缘特征反映图像主题内容，直接冲击人的视觉感受，所以图像降噪处理过程中应尽量保持图像的边缘特征。而在传统的BM3D降噪方法中，由于使用全局相似块匹配，相似块匹配效果较差，在较高方差的高斯噪声干扰下，滤波过程常造成图像中物体边缘和噪声一起作为高频信息被削减或去除，造成物体边缘模糊严重，导致最终降噪效果不良。

发明内容

本发明的目的是提出一种适用于中子成像系统图像的基于高斯-泊松混合噪声模型的图像复原方法。

本发明是通过以下技术方案实现的。

本发明所述的一种适用于中子成像系统图像的基于高斯-泊松混合噪声模型的图像复原方法，其特征是包括以下步骤：

对中子成像图像的降噪问题进行数学建模，设中子成像系统直接得到的图像为z_i＝p_i(y_i)+n_i，其中n_i为加性的方差为σ²的高斯白噪声，p_i(y_i)为乘性的参数为y_i的泊松噪声项。其中泊松噪声项的参数y_i即为无噪图像i位置的像素亮度值。中子成像图像的降噪问题就是从z_i中尽可能恢复出y_i。降噪质量可使用降噪后图像与无噪图像的PSNR值进行评估，此值越高则降噪效果越优秀。

S1：对要处理的原始图像进行GAT变换；

S2：对变换后的图像使用BM3D算法Stage1处理，得到基本估计图像；

S3：对S2所得基本估计图像使用k-均值聚类算法划分区域；

S4：在各区域内使用BM3D算法Stage2处理(即将Stage2的相似块匹配限制在同一图像区域内)，得到最终估计图像；

S5：对最终估计图像做无偏GAT逆变换得到复原后的图像。

本发明步骤S3所述的区域划分，有多种方法，本发明优选以下方法：

一种用于BM3D Stage1与Stage2阶段之间的，基于物体边缘特征的图像区域划分方法，也即一种k-均值聚类算法。k-均值聚类算法使用图像中的随机像素点作为初始值进行迭代。迭代过程使用Lloyd算法进行：

①在图像中选取随机分布的k个像素作为聚类中心点初始值；

②计算每个聚类中心点周围N/k×N/k范围内的像素亮度中位值，以这些中位值作为新的聚类中心点；

③检查迭代过程中心点移动距离是否小于给定阈值，若不小于则回到步骤②，否则认为迭代过程收敛，进入下一步；

④以收敛后得到的中心点对图像进行聚类，对某个图像中的像素，计算其与每个中心点的亮度差异值，将此像素归入与之亮度最接近的中心点。

如此将图像划分为以物体边缘划分的k个区域。通常k值与画面中物体数量成正相关。

本发明主要目的是解决传统BM3D算法在强噪声下边缘信息丢失严重的技术问题，提出一种改进的中子成像图像降噪算法，有效降低图像噪声，提高中子成像系统的图像质量。

更具体地说，本发明通过以下技术方案来实现。

本发明所述的一种适用于中子成像系统图像的基于高斯-泊松混合噪声模型的图像复原方法，利用GAT变换将高斯泊松混合噪声稳定化为高斯噪声，并利用区域限制BM3D方法减少变换后图像的噪声和模糊，最后使用GAT的无偏逆变换得到最终的图像。具体实施方法包括如下步骤：

S1：对要处理的原始图像z进行GAT变换

其中σ²为高斯噪声方差，以将高斯泊松混合噪声变换化为高斯白噪声，便于后续BM3D算法处理。