演示文稿在图像处理中的应用.pptVIP

subplot(221);imshow(BlurredNoisy); title(A = Blurred and Noisy); subplot(222);imshow(J1); title(deconvlucy(A,PSF)); subplot(223);imshow(J2); title(deconvlucy(A,PSF,NI,DP)); subplot(224);imshow(J3); title(deconvlucy(A,PSF,NI,DP,WT)); 现在是63页\一共有107页\编辑于星期六 图45 用Lucy-Richardson 算法进行恢复的图像比较 现在是64页\一共有107页\编辑于星期六 5 用盲去卷积算法进行恢复 使用 deconvblind 函数，采用盲去卷积算法恢复图像。该算法用在不知道失真信息的情况下比较有效。deconvblind 函数使用与加速衰减 Lucy-Richardson 算法近似的迭代处理来同时恢复图像和 PSF。 与 deconvlucy 函数类似，deconvblind 函数也实现了几个原 Lucy-Richardson 最大似然算法的改进版本来完成复杂的图像恢复任务。 现在是65页\一共有107页\编辑于星期六 【例13】Deconvblind函数应用举例： 参考程序： I = checkerboard(8); PSF = fspecial(gaussian,7,10); V = .0001; BlurredNoisy = imnoise(imfilter(I,PSF),gaussian,0,V); WT = zeros(size(I)); WT(5:end-4,5:end-4) = 1; INITPSF = ones(size(PSF)); [J P] = deconvblind(BlurredNoisy,INITPSF,20,10*sqrt(V),WT); 现在是66页\一共有107页\编辑于星期六 subplot(221);imshow(BlurredNoisy); title(A = Blurred and Noisy); subplot(222);imshow(PSF,[]); title(True PSF); subplot(223);imshow(J); title(Deblurred Image); subplot(224);imshow(P,[]); title(Recovered PSF); 现在是67页\一共有107页\编辑于星期六 图46 用盲卷积算法进行恢复的图像比较 现在是68页\一共有107页\编辑于星期六 6 图像的重建操作 图像重建操作指通过对离散图像进行线性空间内插或线性空间滤波来重新获得连续图像的方法。本章节主要以线性空间滤波对二维静态图像进行重建。 MATLAB图像处理工具箱中的Phantom函数，可以产生一幅Shepp-Logan人脑模型图，该图作为一个二维平面测试图，对其进行重建操作可以模拟出许多人脑横切面的性质和特征。 现在是69页\一共有107页\编辑于星期六 【例14】调用Phantom函数，产生Shepp-Logan人脑模型图，利用Radon Transform和Inverse Radon Transform函数实现图像的重建。 参考程序： angularvector=0:179;%设置角度从0到179. P = phantom(Modified Shepp-Logan,200);%调用Phantom函数产生重建原始模型. imshow(P);%显示原始图像. angularvector1=0:10:170;%设置角度的起止点和步长. 现在是70页\一共有107页\编辑于星期六 [R1,xp]=radon(P,angularvector1);%在设置角度情况下，对原始图像调用Radon变换. angularvector2=0:5:175; %设置角度新的起止点和步长，赋予新的变量. [R2,xp]=radon(P,angularvector2);%对应新的角度，对原始图像调用Radon变换. angularvector3=0:2:178;%原理同上. [R3,xp]=radon(P,angularvector3);%综上为R1,R2,R3三个不同角度下的Radon变换. 现在是71页\一共有107页\编辑于星期六 figure,imagesc(angularvector,xp,R1);%独立显示 imagesc绘制Radon变换后的离散点阵. title(radon transform of the Sh