第6次课 医学图像恢复技术

医学图像处理

MedicalImageProcessing陈家益(beyond38@163.com)广东医学院信息工程学院计算机科学教研室主要内容：5.2医学图像中的噪声及滤除5.3系统退化模型5.1图像恢复的基本概念第五章医学图像恢复技术（3）噪声图像的恢复3）顺序统计滤波器修正后阿尔法均值滤波器：中值滤波器：g(x,y)=median(f(x,y))最大值滤波器：g(x,y)=max(f(x,y))最小值滤波器：g(x,y)=min(f(x,y))中点滤波器：g(x,y)=(1/2)*(min(f(x,y))+max(f(x,y)))5.2医学图像中的噪声及滤除（3）噪声图像的恢复3）顺序统计滤波器修正后阿尔法均值滤波器：适用于混合噪声中值滤波器：适用于脉冲噪声最大值滤波器：适用于胡椒噪声，滤除暗点最小值滤波器：适用于盐噪声，滤除亮点中点滤波器：适用于高斯噪声和均匀随机分布的噪声5.2医学图像中的噪声及滤除（3）噪声图像的恢复①修正后阿尔法均值滤波设g(x,y)表示带噪图像，Sxy表示图像中点(x,y)的m*n邻域。对邻域Sxy的所有像素g(s,t)排序，去掉d/2个最低灰度值和d/2个最高灰度值。用gr(s,t)表示剩余的m*n-d个像素。由这些像素的平均值形成的滤波器称为修正后阿尔法均值滤波器。5.2医学图像中的噪声及滤除（3）噪声图像的恢复①修正后阿尔法均值滤波其中0<d<mn-1，当d=0,修正后阿尔法均值滤波退化为算术均值滤波；当d=mn-1，滤波退化为中值滤波。当d取中间的其他值时，修正后阿尔法均值滤波适用于处理包括多种噪声的混合情况。5.2医学图像中的噪声及滤除①修正后阿尔法均值滤波functionf=alphatrim(g,m,n,d)g=im2double(g);f=imfilter(g,ones(m,n),'symmetric');fork=1:d/2f=imsubtract(f,ordfilt2(g,round(k),ones(m,n),'symmetric'));endfork=(m*n-(d/2)+1):m*n

f=imsubtract(f,ordfilt2(g,round(k),ones(m,n),'symmetric'));endf=f/(m*n-d);

5.2医学图像中的噪声及滤除（3）噪声图像的恢复①修正后阿尔法均值滤波

若对修正后阿尔法均值滤波稍加改变，就可以得到K邻近均值（中值）滤波。它是在Sxy中选择K个与待处理像素的灰度差为最小的像素，再用这K个像素的均值（中值）替换原来的像素，因此有mn-d=K.5.2医学图像中的噪声及滤除（3）噪声图像的恢复①修正后阿尔法均值滤波为了减少平滑后图像的模糊，在对当前像素进行平滑时，应尽量保留边界点和细节点，而选择非边界点非细节点进行平滑。修正后阿尔法均值滤波和K邻近均值（中值）滤波都是对与当前像素的灰度差较小的邻域像素进行均值或中值处理，具有一定的边界点和细节点保持效果。5.2医学图像中的噪声及滤除（3）噪声图像的恢复②中值滤波器g=medfilt2(f,[mn],'symmetric');③最大值滤波g=ordfilt2(f,m*n,ones(m,n),'symmetric');④最小值滤波g=ordfilt2(f,1,ones(m,n),'symmetric');5.2医学图像中的噪声及滤除（3）噪声图像的恢复⑤中点滤波g1=ordfilt2(f,1,ones(m,n),'symmetric');g2=ordfilt2(f,m*n,ones(m,n),'symmetric');g=imlincomb(0.5,g1,0.5,g2);5.2医学图像中的噪声及滤除（3）噪声图像的恢复4）自适应滤波器

自适应中值滤波算法(AdaptiveMedianFilteringAlgorithm,AMF)是一种中值滤波，使用一个邻域区域作为滤波窗口，在滤波过程中会根据一定的设定条件改变滤波窗口的大小。5.2医学图像中的噪声及滤除（3）噪声图像的恢复4）自适应滤波器

另外，自适应中值滤波算法主要通过判断当前滤波窗口的中心像素点的性质来决定滤波的输出：当判断窗口中心的像素是脉冲噪声时，将中心像素用滤波窗口中的像素排序后的中值代替；否则不改变当前的像素值，直接输出。5.2医学图像中的噪声及滤除（3）噪声图像的恢复4）自适应滤波器

自适应中值滤波算法分为两步进行：噪声检测和噪声滤波，基本步骤如下：

Wij表示正被处理的、中心在(i,j)处的滤波窗口；

S0为默认的初始窗口大小，Smax为允许的最大窗口大小。

fij为(i,j)处的灰度值，fmin为Wij中的最小灰度值，fmax为Wij中的最大灰度值，fmed为Wij中的灰度中值。

5.2医学图像中的噪声及滤除（3）噪声图像的恢复4）自适应滤波器

自适应中值滤波算法工作在两个层面，表示为LevelA和LevelB。LevelA:

如果fmin<fmed<fmax，转向LevelB，否则，增大窗口尺寸；如果窗口Wij≤Wmax，重复LevelA，否则输出fmed。LevelB:

如果fmin<fij<fmax，输出fij，否则输出fmed。5.2医学图像中的噪声及滤除（3）噪声图像的恢复functionf=adpmedian(g,Smax)if(Smax<=1)||(Smax/2==round(Smax/2))||(Smax~=round(Smax))error('SMAXmustbeanoddinteger>1.')endf=g;f(:)=0;alreadyProcessed=false(size(g));fork=3:2:Smaxzmin=ordfilt2(g,1,ones(k,k),'symmetric');zmax=ordfilt2(g,k*k,ones(k,k),'symmetric');zmed=medfilt2(g,[kk],'symmetric');processUsingLevelB=(zmed>zmin)&(zmax>zmed)&~alreadyProcessed;zB=(g>zmin)&(zmax>g);5.2医学图像中的噪声及滤除（3）噪声图像的恢复

outputZxy=processUsingLevelB&zB;outputZmed=processUsingLevelB&~zB;f(outputZxy)=g(outputZxy);

f(outputZmed)=zmed(outputZmed);alreadyProcessed=alreadyProcessed|processUsingLevelB;ifall(alreadyProcessed(:))break;endendf(~alreadyProcessed)=zmed(~alreadyProcessed);5.2医学图像中的噪声及滤除（1）常见的系统退化函数实际获得的图像不仅包含噪声，还包含各种因素引起的图像退化，如相机聚焦不佳、系统散射等。对于线性移不变成像系统，退化因素可由退化函数h(x,y)反映。5.3系统退化复原（1）常见的系统退化函数成像系统的系统函数或退化函数h(x,y)又称点扩散函数(PSF)。输出的退化图像可以表示为原图像与点扩散函数的卷积，再加上系统的噪声。5.3系统退化复原（1）常见的系统退化函数图像恢复的效果取决于对图像退化过程的先验知识所掌握的准备程度，知道关于退化函数h和噪声n的信息越多，对原图像的估计就越准确。

对于退化图像的恢复，关键是确定或得到系统的点扩散函数h。5.3系统退化复原（1）常见的系统退化函数成像系统的退化函数一般可以通过对成像系统退化过程的先验知识得到。

常用的估计方法有：模型估计法图像观察估计法实验估计法5.3系统退化复原（1）常见的系统退化函数①图像观察估计法利用退化图像自身携带的信息对退化模型进行估计，它通过选择图像中具有简单结构或较强特征的局部区域（如点源的模糊图像），并构造一个具有相同大小和特征、但没有退化的估计图像（点源）通过解卷积来得到退化函数的近似估计。5.3系统退化复原（1）常见的系统退化函数①图像观察估计法如果能确定原图像的某部位有一个清晰的点，假定噪声可忽略，则使得该点退化的模糊图像就是h(x,y)。从原图像中找到一块信号内容较强的区域，构建出该区域的估计图像，进而大致推出退化函数h(x,y):5.3系统退化复原（1）常见的系统退化函数①图像观察估计法图像观察估计法特点和应用场合：没有任何关于退化函数h的知识，只能从图像本身来收集信息。选择强信号区域（如高对比度区域）进行处理，是为了降低噪声的影响。处理过程复杂，仅在特殊环境下使用，如恢复具有历史价值的老照片。5.3系统退化复原（1）常见的系统退化函数②实验估计法实验估计法通过使用或设计一个与原图像退化过程相似的成像系统装置或过程，并利用这个系统或过程对点源成像，从而得到原系统点扩散函数的较准确估计。5.3系统退化复原（1）常见的系统退化函数②实验估计法实验估计法特点和应用场合：可以使用与获取退化图像设备一致或相似的装置时使用。用上述系统对一个冲激成像，冲激可由一个亮点来模拟，该点应尽可能亮，以便将噪声的影响降低到可以忽略的程度。对冲激的成像结果即为PSF函数。5.3系统退化复原（1）常见的系统退化函数③模型估计法模型估计法基于对不同成像过程基本原理和物理特性的了解和建模。模型估计法特点和应用场合：一般需要从基本物理原理推到数学模型。模型建立后应用范围广。5.3系统退化复原（1）常见的系统退化函数可通过调用函数fspecial创建常见的退化滤波器。h=fspecial('type'，parameter);常用的滤波器：typeparametergaussianhsize，sigmasobelprewitt5.3系统退化复原（1）常见的系统退化函数常用的滤波器：typeparameterlaplacianalphaloghsize，sigmaaveragehsizediskradiusMotionlen,theta5.3系统退化复原（1）常见的系统退化函数用退化滤波器退化图像：h=fspecial('motion',10,60);g=imfilter(f,h);h=fspecial('disk',10);g=imfilter(f,h);5.3系统退化复原（2）逆滤波

对于线性移不变系统，有退化模型为了方便讨论和求解，用顺序堆叠的方式将二维图形表示为一维，上述退化模型可转化为：原