图像盲恢复的算法研究论文

张玉君

（地质矿产部航空物探遥感中心研究所北京）

摘要：本文介绍一种独特的航放数据图像复原方法。该方法的主要技术关键是：提出航放数据图像复原原理和理论基础；建立航放数据图像复原处理流程；制定重建数据网格文件的途径；进行航放数据图像复原效果及误差评价。

关键词：航放数据，大气本底，图像处理，图像复原技术。

一、引言

自20世纪70年代初期方柱形NaI晶体进入机载综合航空站以来，航放测量的灵敏度和有效性显著提高，地质及地球物理勘探界对于航空放射性测量的需求产生了根本变化。

在近20年的应用实际中，大气氡本底（简称大气本底）的改正很难准确，它始终是困扰该方法应用效果的主要难点；其后果是在图面上造成条带现象，从而严重地影响着图件的可用性及方法的效果。其原因[1]可概括为：空中所测放射性不仅仅来源于地下，而且受飞机硬件环境、宇宙射线、大气中氡及其子体的影响。后者称之为大气本底干扰，它又受气候、风力、风向、温度、季节及一天中何时测量等因素的影响。大气本底干扰的主要表现形式是架次与架次间本底水平不同。受干扰最大的是铀道，钾道次之，钍道和总道虽较小，但也不可忽视，（见彩版附图7中的图3、图4）。由于这种噪声的存在，来源于地质体的信息常常淹没于噪声之中。图3a（彩版附图7）为哈密土墩测区K（红）、Th（绿）、U（蓝）三元素复原图像，图3b为该测区航放原始数据合成图像，图4a（彩版附图7）表示各架次早、晚校准读数，图4b为总道原始数据图像。条带噪声的存在，可以形象地比喻为挂在有用信息图像前面的彩条窗帘，条带的严重性使得该工区原始航放数据无法绘制等值线图。

航放图面条带问题是一个“世界性”的问题[2]。解决得较好的是加拿大，靠星罗棋布的湖泊河流等水域上空测量结果改正本底，取得较好的效果，而且不使用向上探头[1]。美国Geometrics公司及其他航空物探公司则是靠向上探头测量，作为本底改正的依据[3]。1986年Grasty[4]提出当测区内没有湖泊时，可用测线上无异常区的平均值代替本底。

本文介绍的方法与国际上已采用的各种方法全然不同，该方法在数字图像处理学中可称之为航放图像复原技术。图像复原技术的主要目的，是要改善给定的图像。复原是一个过程，它试图利用蜕化现象的某种先验知识，把已经蜕化了的图像加以重建或恢复。因此，复原技术是把蜕化模型化，并运用相反的过程在某种程度上恢复原来的图像。

Cannon博士[5]研究了一种图像复原技术或称图案去除技术，它适用于：从规则图案（如纺织品）上提取指纹图形，改善散焦图像，消除卫片图像探测器与探测器间的噪声，使在曝光过程中相机或物体平移造成模糊的图像清晰化等。Srinivasan也报道了这类研究[6]。张玉君等研究了深海锰结核照片光照不均匀等蜕化现象的图像复原问题[7]。航放数据图像复原处理是数字图像复原技术在地学界成功应用的又一实例，但航放数据图像所存在的蜕化问题与上述各例均不相同。该方法研究成功后，曾在6个测区得到验证。

二、航放数据图像复原技术原理和理论基础

航放所测到的是一幅蜕化了的图像G（x,y）它可视为由真实图像F（x,y）与干扰图像η（x,y）叠加而成，简化了的蜕化过程，见图1。航放图像蜕化现象的先验知识来自对航放测量过程及原始图像的分析。在测量过程中，来自地质体的有用信息是不随时间为转移的。而干扰在本质上是随时间变化的，但在图像上干扰已变为（x,y）的函数，因为：

张玉君地质勘查新方法研究论文集

图1航放数据图像蜕化示意图

η的变化可分为架次之间的跳变及架次之内的渐变，见图4（彩版附图7），在每一测线上此干扰大致为一个常数，如果将x（即图像上的列）表示垂直测线的方向，则η（x,y）简化为η（x），则有

张玉君地质勘查新方法研究论文集

航放图像复原的目的，就在于设法近似地求出η（x），从而近似地得到F（x,y）。为此，沿测线方向对原始图像进行多次单列

多行窄长窗口褶积：

张玉君地质勘查新方法研究论文集

式中W为褶积模板，是由加权因子组成的矩阵。褶积过程是一种线性运算，其算子H不随空间变化。因为算子为线性的，则两个输入之和的响应等于两个响应之和。

张玉君地质勘查新方法研究论文集

由于假定了η仅与x有关，又由于褶积窗口为单列，则有：

张玉君地质勘查新方法研究论文集

现在分析HF（x,y）的性质，由于沿y方向的多次滑动平均，局部异常“淹没”于近区域特征之中，这种近区域特征表现为沿测线方向的低缓变化；如果用f（x,y）表示局部异常，用L（x,y）表示近区域场，于是：

张玉君地质勘查新方法研究论文集

再经如下处理

张玉君地质勘查新方法研究论文集

由（9）式可见，从原始图像中减去噪声图像后，所得复原图像f（x,y），它从局部异常角度是接近真实图像的，误差取决于所减去的“近区域背景值”在测线方向起伏变化的幅度。

三、航放数据图像复原处理过程

航放数据图像复原技术的研究是以多元统计为理论基础，以图像处理为工具来完成的，并体现了图像处理快速直观的特点，其流程如图2所示。

图2航放数据图像复原处理流程

此方法假定航放噪声本底沿测线方向不变或呈线性变化。通过沿测线方向的多次滑动平均，使局部异常逐渐淹没在噪声本底之中，得出一幅与噪声本底线性相关的噪声图像。噪声图像尚需进行边缘影响补偿；对去除噪声后的图像，经中值滤波和空间变量反差增强，达到最终复原的效果。这一复原过程集中表示于图2的左半部。

图2的右半部为数据网格文件的重建过程，它是实际应用所必不可少的。经过分类分区，求得各类别在复原前后的均值向量，经最小二乘拟合求出复原图像的元素含量或计数率值，重新建立为在主计算机上绘制等值线图用的网格文件。

本研究曾试验通过沿测线方向取平均值做为噪声水平，结果不及上述方法理想。

四、效果及误差评价

1.航放数据图像复原的效果

（1）图面直观效果的改善。

可以形象地说，航放图像复原好比揭去一层条带窗帘，使原来透过此窗帘隐约可见的图像显示出了真面目，见图3a（彩版附图7）。图面直观效果的改善还表现在由于定位问题所引起的岩体边界上的锯齿状噪声得以消除，见图5（彩版附图7）。图5为总道对比图像，5a为原始数据，5b表示噪声图像，5c为去噪声后的图像，5d为复原图像。

（2）用复原数据所做等值线图真实可信。

以哈密土墩测区为例，原始资料由于条带干扰，在主计算机上，钾、钍、铀道都无法绘等值线图，仅提供了平剖图；只有总道提供了等值线图，但仍可看到条带的影响。

经图像复原、重建网格文件，反馈回主计算机后绘制了TC、K、Th和U等值线图，现以K道经复原后数据等值线图为例示于图6（彩版附图7），与地质图对比，表明异常和地质体对应良好，各类岩性的放射性趋势也都吻合，证实这些等值线图的可靠性。利用复原图像所做分类图也证实了这一点，见图7（彩版附图7），图7中数字分别为：①超基性岩；②基性岩；③花岗岩；④闪长岩；⑤变质岩；⑥混合岩；⑦第四纪沉积；⑧第三、四纪沉积；⑨第三纪沉积。

（3）有用信息增加。

本研究利用多元统计的方法，对航放图像复原的效果给出定量评价。可用一幅图像有用信息构成的变异值的大小来对它做定量评价。为此应计算全图面总变异对于一个象素的平均值，即平均变异值。用C、C´和C＂分别表示原始图像中有用信息平均变异值、原始图像中干扰信息平均变异值和最终复原图像有用信息平均变异值。统计时以G´（x,y）近似代表η（x）；以[G（x,y）-G´（x,y）]近似代表F（x,y）；以P（x,y）表示最终复原图像，并假定它已无干扰存在。

张玉君地质勘查新方法研究论文集

式中，字母上加“—”表示平均值；M、N为图像的行、列数。

表1为哈密土墩测区航放数据图像按上述各式所做定量评价的统计结果。

表1

从表1可见，K、Th、U、TC经图像复原后，有用信息都有十分显著的增长；就此工区而论，TC和K原始图像相对质量较好，Th和U较差。

2.复原图像的准确度及误差评价

复原图像的主要误差来源是“近区域背景值”L（x,y），它是在多次滑动平均时形成的。通过对干扰图像剖面数据的统计，得到以下准确度评价：

K±（绝对含量）；Th+ pp m；

U± ppm；TC±计数。

五、结论

（1）本文介绍的方法是在国内外首次提出的方法独特的航放数据图像的复原技术，并在多个工区验证了其可靠性和实用性。

（2）本技术可以基本上消除由于大气本底及阈值变动所造成的图面条带现象，基本复原航放图像的真面目，为进一步图像处理（诸如：求导、增强、分类、逻辑运算等）做了准备，因此本技术也是一种快速预处理方法。

（3）本方法改善了由于飞行往返定位位移所引起的某些地质体边缘呈锯齿状的图像噪声问题。

（4）本研究建立了“有用信息平均变异值”做为定量评价航放数据图像复原效果的尺度。还讨论了图像复原做为一种预处理过程，对于元素含量值可能导入的绝对误差或称为方法的准确度。

参考文献

[1]Grasty, ., Gamma ray spectrometric methods in uranium exploration—Thcory and operational procedures, Geophysics and Geochemistry in the Search for Metallic Ores,GSC,Ottawa,147-162,1977.

[2]Creen, airborne gamma-radiation data using between-channel correlation information,Geophysics,52,1557-1562,1987.

[3]Foote, ., Improvement in airborne gamma-radiation data analysis by removal of environmental and pedologic radiation changes, in the Use of Nuclear Techniques in Prospecting and Developmcnt of Mineral Resources: Energy Agency Mtg.,Buenos Aires,187-196,1968.

[4]Grasty, system for computing on-line atmospheric backgrounds,GSC paper,1-52,1987.

[5]Cannon,M.,Lehar, Preston,F.,Background pattern removal by power spectral filtering,Applied Optics,22,777-779,1983.

[6]Srinivasan,R.,Software image restoration techniques,Digital Design, 16,4,27-34, 1986.

[7]张玉君，史鉴文.深海多金属结核照片的图像复原和图像处理技术研究.物探与化探，1989,（13）:435～441.

致谢林振民同志对本文提出了宝贵的意见，史鉴文同志参加了重复工区试验，张志民和谢欣同志分别编制了网格文件转换和最小二乘拟合程序，杨星虹同志拍摄了屏幕图片，水恩海同志搜集了试验工区校准资料，在此一并致谢。

A STUDY ON IMAGE RESTORATION TECHNIQUES FOR AERORADIOMETRIC DATA

Zhang Yu jun

（Research Institute, Center of Aero-Gcophysics and Remote Sensing,Ministry of Geology and Mineral Resources, Beijing）

Abstract

This paper represents a specific methodfor restoration of images of airborne radiometric main technical keys involved in this study are；the advancementof the principles and theory；the establishment of the flow-diagram for processing；the formulation of the means for reestablishment of the gridded data file；the evaluation of the restoration results and the errors, involved by the restoration processing.

Key words Aeroradiometric data, Atmospheric background, Image processing, Image restoration techniques.

原载《地球物理学报》，1990,。

​

论文简介： 利用图像传输理论测量海水的点扩散函数和调制传递函数并且使用维纳滤波器复原模糊的图像。退化方程H(u,v)在水槽中测量得到。在实验中利用狭缝图像和光源，第一步:一维光照射到水中从而得到不同距离下的狭缝图像数据，这样一维的海水点扩散函数就可以通过去卷积得到。又因为点扩散函数的对称性二维的函数模型也可以通过数学方法得到。利用相似的方法调制传递函数也可以得到。这样传输方程便可以得到:

图像可以由下式获得:

​

论文简介： 论文中提出自然光照下的水下图像退化效果与光偏振相关,而场景有效箱射则与光偏振无关。在相机镜头端安装可调偏振器,使用不同偏振角度对同一场景成两幅图像,所得到的图像中的背景光会有明显不同。通过对成像物理模型的分析,利用这两幅图像和估计出的偏振度,就能恢复出有效场景辐射。他还提出了一个计算机视觉方法水下视频中的退化效应。分析清晰度退化的物理原因发现主要与光的部分偏振有关。然后提出一个逆成像方法来复原能见度。该方法基于几张通过不同偏振方向的偏振片采集图像。

论文简介： 论文提出了一种自适应滤波的水下图像复原方法。通过最优化图像局部对比度质量判决函数,可以估计出滤波器中所使用的参数值。 论文提出一种基于简化的Jaffe-McGlamery水下成像模型的自调谐图像复原滤波器。滤波器的最优参数值是针对每幅图像通过优化一个基于全局对比度的质量准则自动估算的。(对一幅图像滤波器能根据全局对比度自动估计最优参数值),简化的模型理想地适合后向散射较少的漫射光成像.1.首先简化Jaffe-McGlamery水下成像模型：假设光照均匀(浅水区阳光直射)，并且忽略后向散射部分.然后基于简化后的成像模型设计一个简单的反滤波器2.将滤波器设计成自适应滤波器。

论文简介： 论文对于调制传递函数给出了详细准确的系统函数信息,水下图像可以用它或点扩散函数进行复原.作者进行实验测量了水质参数得出了这些函数，并用得出的函数进行了图像复原。同时他还建立了一个框架来最大限度复原水下图像，在这个框架下传统的图像复原方法得到了拓展，水下光学参数被包含了进去，尤其时域的点扩散函数和频域的调制传递函数。设计了一个根据环境光学特性进行调整的客观图像质量度量标准来测量复原的有效性。

论文简介： 调制传递函数给出了详细准确的系统函数信息,水下图像可以用它或点扩散函数进行复原.作者进行实验测量了水质参数得出了这些函数，并用得出的函数进行了图像复原。(这一部分在王子韬的论文中有比较详细介绍)

论文简介： 在散射媒介中的正则化图像复原。论文在基于物理原因的复原方法难以去除噪声以及透射率低的基础上，提出一种自适应的过滤方法，即能明显的改善可见性，又能抑制噪声放大。本质上，恢复方法的正规化，是适合变化媒介的透射率，因此这个正则化不会模糊近距离的目标。

论文简介： 论文提出一种基于对边缘进行GSA（灰度规范角度）加权的测量图像清晰度的方法。图像首先被小波变换分解，去除部分随机噪声，增加真实边缘检测的可能性。每个边缘锐度由回归分析方法基于灰度的一个角的正切来确定边缘像素的灰度值之间的斜率和位置。整个图像的清晰度是平均每个测量的GSA的比例加权的第一级分解细节的量，作为图像的总功率，最后通过图像噪声方差自适应的边缘宽度。

论文简介： 论文提出了基于主动偏振的人工光照下水下图像处理技术。在宽场人工光照下的水下成像中,在光源端或相机端安装可调偏振器。通过调整光源或相机端的偏振器,同时拍摄两幅或多幅同一场景的图像,从两幅图像中可估计出背景光的偏振度。结合水下成像物理模型,就可以进行图像复原和场景3D信息估计。该方法操作简单,设备筒易,适用于水下画定目标的成像。 大范围人工照明条件下研究成像过程,基于该成像模型,提出一种恢复object signal的方法,同时能获得粗糙的3D scene structure.相机配备检偏振器,瞬间获取同一场景的两帧图片with different states of the analyzer or light-source polarizer,然后用算法处理获取的图片.它统一并推广了以前提出的基于偏振的方法.后向散射可以用偏振技术降低,作者在此基础上又用图像后处理去除剩余的后向散射,同时粗糙估测出3D场景结构.创新：之前的方法有的认为目标物反射光的偏振度可以忽略(即认为只有后向散射是偏振的)；另外还有的认为后向散射的偏振度可以忽略(即认为只有目标物反射光是偏振的)。本文作者认为两者都是部分偏振光。

论文简介： 论文在没有应用任何标准模式、图像先验、多视点或主动照明的条件下同时估算了水面形状和恢复水下二维场景。重点是应用水面波动方程建立紧凑的空间扭曲模型，基于这个模型，提出一个新的跟踪技术，该技术主要是解决对象模型的缺失以及水的波动存在的复杂的外观变化。在模拟的和真实的场景中，文本和纹理信息得到了有效的复原。

论文简介： 论文提出暗通道先验算法复原有雾图像。暗通道先验是一系列户外无雾图像的数理统计，基于观察户外无雾图像的大部分补丁补丁中包含至少一个颜色通道中低强度的像素点。在有雾图像中应用这些先验，我们可以直接的估算雾的厚度，复原成高质量的无雾图像，同时还能获得高质量的深度图。

论文简介： 论文比较研究了盲反卷积算法中的：R-L算法（Richardson-Lucy）、最小二乘法以及乘法迭代法。并且应用了水下图像去噪和威尔斯小角度近似理论推导出点分布函数。通过执行威尔斯的小角度散射理论和模糊度量方法对三种盲反卷积算法进行比较，确定总迭代次数和最佳图像复原结果。通过比较得出：最小二乘算法的复原率最高，但是乘法迭代的速度最好。

论文简介： 论文提出点扩算函数（PSF）和调制解调函数（MFT）的方法用于水下图像复原，应用基于威尔斯小角度近似理论来进行图像增强。在本文中作者分析了水下图像退化的原因，在强化超快激光成像系统中采用了距离选通脉冲的方法，降低了反向散射中的加性噪声。本文对图像的基本噪声模式进行了分析，并使用算术平均滤波首先对图像进行去噪，然后，使用执行迭代盲反褶积方法的去噪图像的初始点扩散函数的理想值，来获得更好的恢复结果。本文通过比较得出，盲反褶积算法中，正确使用点扩散函数和调制解调函数对于水下图像复原的重要性。

论文简介： 本文提出一种图像复原的新方法，该方法不需要专门的硬件、水下条件或现在知识结构只是一个与小波变换的融合框架支持相邻帧之间的时间相干性进行一个有效的边缘保留噪声的方法。该图像增强的特点是降低噪声水平、更好的暴露黑暗区域、改善全局对比、增强细节和边缘显著性。此算法不使用补充信息，只处理未去噪的输入退化图像，三个输入主要来源于计算输入图像的白平衡和min-max增强版本。结论证明，融合和小波变换方法的复原结果优于直接对水下退化图像进行去雾得到的结果。

论文简介： 本文是一篇综述性质的论文。介绍了：1、水下光学成像系统 2、图像复原的方法（对各种图像复原方法的总结） 3、图像增强和颜色校正的方法总结 4、光学问题总结。

论文简介： 论文针对普通水下图像处理的方法不适用于水下非均匀光场中的问题，提出一种基于专业区域的水下非均匀光场图像复原方法，在该算法中，考虑去除噪声和颜色补偿，相对于普通的水下图像复原和增强算法，该方法获得的复原复原的清晰度和色彩保真度通过视觉评估，质量评估的分数也很高。

论文简介： 论文基于水下图像的衰减与光的波长的关系，提出一种R通道复原方法，复原与短波长的颜色，作为水下图像的预期，可以对低对比度进行复原。这个R通道复原的方法可以看做大气中有雾图像的暗通道先验方法的变体。实验表明，该方法在人工照明领域应用良好，颜色校正和可见性得到提高。

论文简介： 作者对各种水下图像增强和复原的算法做了调查和综述，然后对自己的提高水下质量的方法做了介绍。作者依次用到了过滤技术中的同态滤波、小波去噪、双边过滤和对比度均衡。相比于其他方法，该方法有效的提高了水下目标物的可见性。

论文简介： 论文应用湍流退化模型以质量标准为导向复原因水下湍流退化的图像。参考大气湍流图像复原的算法，省略了盐分的影响，只考虑水中波动引起的湍流对水下成像的影响，应用一种自适应的平均各向异性的度量标准进行水下图像复原。经过验证，使用STOIQ的方法优于双频谱的复原方法。

论文简介： 本文提出了一种新的方法来提高对比度和降低图像噪声，该方法将修改后的图像直方图合并入RGB和HSV颜色模型。在RGB通道中，占主导地位的直方图中的蓝色通道以95%的最大限度延伸向低水平通道，RGB通道中的低水平通道即红色通道以5%的最低限度向上层延伸且RGB颜色模型中的所有处理都满足瑞利分布。将RGB颜色模型转化为HSV颜色模型，S和V的参数以最大限度和最小限度的1%进行修改。这种方法降低了输出图像的欠拟合和过拟合，提高了水下图像的对比度。

论文简介： 论文根据简化的J-M模型提出一种水下图像复原的有效算法。在论文中定义了R通道，推导估算得到背景光和变换。场景可见度被深度补偿，背景与目标物之间的颜色得到恢复。通过分析PSF的物理特性，提出一种简单、有效的低通滤波器来去模糊。论文框架如下：1.重新定义暗通道先验，来估算背景光和变化，在RGB的每个通道中通过标准化变换来复原扭曲颜色。2.根据PSF的性能，选择没有被散射的光，用低通滤波器进行处理来提高图片的对比度和可见度。

论文简介： 论文中对当代水下图像处理的复原与增强做了综述，作者阐明了两种方法的模型的假设和分类，同时分析了优缺点以及适用的场景。

参考：