Xmean=mean(mean(I))*ones(size(I))
segma2=sum(sum((I-Xmean).^2))/prod(size(I))
D=sum(sum((I-I2).^2))/prod(size(I2))
SNR=10*log10(segma2/D)%%%%% 信噪比 (SNR:Signal to Noise Ratio)
PSNR=10*log10(255^2/D)%%%%% 峰值信噪比(PSNR:Peak Signal to Noise Ratio)
你这个就是简单的加权滤波算法(模板是3*3),不过加权系数不对,正常应该是中心位置(也就是a(5))的权重最高。改成春知中腊丛值滤波算法很简单,将a(1)~a(9)递增或递减排序(用sort()函数,去help里查轮森樱),然后取出中间元素,赋值给X2(i,j)就哦了~
/埋尺/ 中值滤波对椒盐噪声肢皮RGB=imread('历液差peppers','PNG')
I=rgb2gray(RGB)
J=imnoise(I,'salt &pepper',0.02)
subplot(121),imshow(J)
L=medfilt2(J,[3 3])
subplot(122),imshow(L)
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)