1 前言
计算机技术不断发展，硬件性能不断提升，软件效能不断优化，在其带动下，数字图像处理技术也得到了迅速的发展。以往难以解决的技术性瓶颈，而今已经可以轻而易举的突破。但是在一些专业性较强的领域，对于数字图像处理有着一些特殊的要求，很多通用的数字图像处理技术、方法，不能充分发挥其优势，甚至无法使用。本文针对数字图像去噪这一典型问题，在简单回顾相关通用数字图像处理技术发展的基础上，针对医学PET图像的特殊性，提出一些自己的看法及思考。
虽然数字图像处理技术很早便被应用于医学相关领域，但是由于医学图像自身的分类多样、成像技术相对复杂，目前仍然有许多影响医学图像成像的诸多因素无法从理论层面解释，所以数字图像处理技术在医学领域的发展仍然相对滞后。
2 数字图像处理技术的回顾
数字图像处理作为一个学科，经过将近60年的发展，已经形成了完善的理论体系，并细化为多个专业方向。数字图像处理技术诞生之初，图像去噪就是其主要目标之一。在通用数字图像处理领域，常见的噪声主要包括乘性噪声和加性噪声。数字图像的噪声往往和图像的特征信息交织在一起，如何有效地区别并去除噪声，同时尽量保持图像的细节完整，是数字图像去除噪声要面对的难题[1].
2.1 空间域通过卷积处理图像噪声 对获得的图像在空间域直接进行处理，往往采用卷积的数学形式。例如常见的均值滤波、中值滤波、为纳滤波等。常见的中值滤波可定义为：
式子中{xij（i,j）∈I2} 表示数字图像个点的灰度值。根据图像处理的实际需要，目前已经有许多成熟的滤波器可供选择。该类方法的优势在于数学结构相对简单，运算量较小，但是图像处理结果相对较差，容易造成图像细节丢失或者噪声抑制不足等结果。
2.2 频率域通过相应转换处理图像噪声 将要处理的图像进行转换，根据实际情况选择适当的频带进行处理，然后经过反变换获得去噪后的图像[2].基于小波及其衍生分析方法的多尺度分析是该种方法中最具代表性的。因为该类方法理论较为成熟，数学机构灵活，所以至今仍然被学术界关注，并且每年都有一些新的方法被提出。虽然该类方法与前面提到的滤波器方法相比在运算复杂性上有着较大的劣势，但其在图像处理结果上的优势明显，并且随着计算机硬件效能的提升，在频率域的多尺度分析方法已经得到了较为普及的应用[3].
传统的小波去噪方法（wavelet-based denoising）：将含有噪声的图像进行小波变换，转化成小波系数Wf（j,k），选择合适的阈值 Tf,根据一定的处理规则，对小波系数进行处理，把小波系数中噪声的部分去除，最后经过小波反变换得到去噪后的图像。
根据实际的图像分析需求，在小波的基础上发展出了很多小波的衍生方法。这些方法又可以按照处理图像时基函数的变化与否分为自适应分析和非自适应分析。其中自适应分析比较常见的方法有：Brushlet、Wedgelet、Bandelet和Directionlet等。非自适应分析常见的方法有：Ridgelet、Curvelet、Contourlet、Shearlet、NSCT等[3].
3 数字图像技术在医学图像领域的应用
数字图像处理在医学领域有着广泛的应用，在这里只针对医疗中常见的大型设备来简单阐述。这里仅例举CT、MRI和PET（positron emissiontomography）三种具有代表性的影像设备来简单分析。
CT、MRI和PET都属于临床常见的影像设备，虽然成像原理各不相同，但却使用许多相同的数字成像技术。CT主要是运用X光穿透被照体，检测X光的衰减，从而反推出被照体的衰减系数，根据已知的衰减系数对照表重建出被照体的解剖图像。MRI的结构成像主要检测静态磁场中外加射频脉冲对人体内氢质子的影响，从而得到对应的MR信号，重建出对应的解剖图像。PET主要是检测注入人体内的放射性核医药物在人体内代谢过程中发生“湮灭”放射出的γ射线，通过图像重建得到对应的功能性信息[4].
以上三种设备普遍使用了空间域通过卷积去噪的方法，最为常见的就是Gaussion filter.几乎常见的品牌设备都将此种方法作为图像重建过程中的去噪标准步骤。一般作为补充，许多品牌的设备也内建了一些“更为复杂”和“有效”的去噪方法。但是实际使用中，去噪效果并不十分理想[5].
4 局限与对策
CT在图像重建过程中，如果只简单使用Gaussion filter,忽略X光的实际物理效应，如康普顿散射、X光的能谱分布等，得不到较为准确的图像信息，这也是目前临床使用的普通CT无法做到准确“定量”的原因之一。MRI成像相较CT更为复杂，而且截至目前尚无研究表明MRI对人体有任何直接危害，所以MRI设备生产商将更多资源投入到如何从设备的硬件改善来提高图像品质，而不是考虑如何从图像处理的角度来解决问题。因此Gaussion filter在MRI中使用也非常广泛。PET成像与CT原理不同，但是同样受限于放射剂量。与MRI相比PET的图像处理技术相对较为多样，因为PET图像自身的一些特点，也决定其需要配合更多的图像处理才能获得更多有价值的信息[6].而对于PET来说，虽然核医药物的γ射线产生符合泊松分布，但是在临床实际图像处理中往往并不针对泊松来处理，Guassion filter依然被使用。
由于早期计算机硬件的效能不高，许多优秀的数字图像处理方法无法应用于临床实践，当今计算机硬件效能已经有了大幅度提升，可以重新考虑在相关大型医疗影像设备中植入较为先进的算法，从而获得更为优良的图像质量。并且应当根据不同影像设备的成像特点，合理优化图像处理技术，做到通过软件提升，使设备硬件发挥最大效率，使病人在最小的辐射剂量下得到最优化的诊疗。
5 讨论与总结
数字图像处理经历了近60年的发展，涌现出了许多优秀的理论与研究成果，也在生产生活的诸多领域得参 考 文 献到了广泛应用。医疗领域对于图像处理有着强大的需求，但就目前来看，该领域的数字图像处理技术发展仍然滞后。对于科研工作者，医疗领域的图像处理有较高的研究价值和较大的研究空间。要根据临床实践，通过深入研究医学图像的成像原理，理清问题的研究重点，有针对性地设计实验方法，从而获得较好的研究结果。在社会发展，人们越来越重视自身健康和生活质量的大背景下，只有把研究工作与生产生活实践相结合，才能获得科研转化为生产力的良性循环。在科学技术飞速发展的大背景下，只有深入研究，细化问题，有针对性地解决具体问题，才是科研的最重要任务。
参 考 文 献
[1] 宁 媛 , 李 皖 . 图 像 去 噪 的 几 种 方 法 分析比较[J].贵州工业大学学报：自然科学版，2006（34）：63-66.
[2] 马国兵，肖培如。基于小波的图像去噪研究综述[J].工业控制计算机，2013（5）：91-92.
[3] 李彦，汪胜前，邓承志。多尺度几何分析的图像去噪方法综述[J].计算机工程与应用，2011（47）。
[4] 陈武凡，秦安，江少峰，等。医学图像分析的现状与展望[J].中国生物医学工程学报，2008（27）。
[5] Le Pogam A,Hanzouli H,Hatt M,et al.Denoising of PET images by combiningwavelets and curvelets for improvedpreservation of resolution and quantitation[J].Medical Image Analysis,2013（17）：877-891.
[6] Mejia JM,Ochoa Dominguez HDJ,VergaraVillegas OO,et al.Noise reduction in small-animal PET images using a multiresolutiontransform[J].IEEE Transactions on MedicalImaging,2014（33）：2010-2019.
