deconvlucy
使用Lucy-Richardson方法解模糊图像
语法
描述
例子
去除高斯模糊使用deconvlucy
读取和显示原始图像没有模糊或噪音。本例中选择作物图像的大小为256 -,- 256左上的(x, y)坐标(50)。
我= imread (“board.tif”);我= imcrop(我[2 50 255 255]);imshow (I)标题(原始图像的)
创建一个代表一个高斯模糊的PSF 5-by-5标准差5和过滤器的大小。
PSF = fspecial (“高斯”5、5);
模拟模糊图像中。
模糊= imfilter (PSF,我“对称”,“conv”);
添加模拟零均值高斯噪声。
V = 0.002;blurred_noisy = imnoise(模糊,“高斯”0 V);imshow (blurred_noisy)标题(模糊和噪声图像)
使用deconvlucy恢复和嘈杂的图像模糊。指定PSF用来创建模糊和迭代的数量减少到5。
luc1 = deconvlucy (blurred_noisy PSF 5);imshow (luc1)标题(恢复图像的)
删除模糊使用几个deconvlucy可选的语法
创建一个示例图像和模糊。
我=棋盘(8);PSF = fspecial (“高斯”7、10);V =。;BlurredNoisy = imnoise (imfilter (PSF),“高斯”0 V);
创建一个权重数组和电话deconvlucy使用几个可选参数。
WT = 0(大小(I));WT (5: end-4, 5: end-4) = 1;j - 1 = deconvlucy (BlurredNoisy PSF);J2 = deconvlucy BlurredNoisy PSF, 20日√(V));J3 = deconvlucy (BlurredNoisy PSF 20, sqrt (V)、WT);
显示结果。
次要情节(221);imshow (BlurredNoisy);标题(“=模糊和噪声”);次要情节(222);imshow (j - 1);标题(“deconvlucy (PSF)”);次要情节(223);imshow (J2);标题(“deconvlucy (PSF,镍、DP) ');:次要情节(224);imshow (J3);标题(“deconvlucy (PSF,倪,DP, WT)”);
输入参数
输出参数
提示
您可以使用
deconvlucy进行反褶积,开始之前的反褶积停了下来。要使用该功能,通过输入图像我作为一个单元阵列,{我}。当你做什么,deconvlucy函数返回的输出图像J作为一个单元阵列,然后你可以通过数组作为输入到下一个deconvlucy调用。输出单元阵列J包含四个要素:J {1}包含我原始图像。J {2}包含最后一次迭代的结果。J {3}包含倒数第二迭代的结果。J {4}是一个数组生成的迭代算法。输出图像
J可能引入的展览响离散傅里叶变换中使用的算法。为了减少振铃,使用I = edgetaper (psf)在调用之前deconvlucy。deconvlucy转换的PSF双没有标准化。deconvlucy可以在输出图像,返回值范围之外的输入图像。
引用
[1]D.S.C.比格斯和m·安德鲁斯加速迭代图像重建算法应用光学,36卷,8号,1997年。
[2]R.J. Hanisch R.L.白色,R.L. Gilliland,反褶积的哈勃太空望远镜的图像和光谱反褶积的图像和光谱,艾德。P.A.简颂,第二版,学术出版社,1997。
