该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
clc
clear
W=imread('e:\matlab\aaa\1.jpg');
imshow (W);
I = rgb2gray(W); % 灰度处理
imshow (I);
title('灰度图像');
J=imnoise(I,'salt & pepper',0.2); %椒盐噪声
figure
imshow (J) ;
title('椒盐噪声');
Q=imnoise(I,'gaussian',0,0.01); %高斯噪声
figure
imshow (Q);
title('高斯噪声');
E=imnoise(I,'speckle',0.2);
figure
imshow(E);
title('斑点噪声'); %斑点噪声
K1= medfilt2(J);
figure %中值滤波椒盐
imshow (K1)
title('中值滤波椒盐');
K2= medfilt2(Q); %中值滤波高斯
figure
imshow (K2)
title('中值滤波高斯');
K3= medfilt2(E); %中值滤波斑点
figure
imshow (K3)
title('中值滤波斑点');
F=double(J); %平均值滤波椒盐
[m,n]=size(F);
for i=2:m-1;
for j=2:n-1;
g1(i,j)=F(i-1,j-1)+F(i,j-1)+F(i+1,j-1)+F(i-1,j)+F(i,j)+F(i+1,j)+F(i-1,j+1)+F(i,j+1)+F(i+1,j+1);
h2(i,j)=g1(i,j)./9;
end
end
F=uint8(h2) ;
figure,
imshow (F)
title('平均值滤波椒盐');
D=double(Q); %平均值滤波高斯
[m,n]=size(D);
for i=2:m-1;
for j=2:n-1;
g(i,j)=D(i-1,j-1)+D(i,j-1)+D(i+1,j-1)+D(i-1,j)+D(i,j)+D(i+1,j)+D(i-1,j+1)+D(i,j+1)+D(i+1,j+1);
h1(i,j)=g(i,j)./9;
end
end
D=uint8(h1);
figure
imshow(D);
title('平均值滤波高斯');
Y=double(E); %平均值滤波斑点
[m,n]=size(Y);
for i=2:m-1;
for j=2:n-1;
g1(i,j)=Y(i-1,j-1)+Y(i,j-1)+Y(i+1,j-1)+Y(i-1,j)+Y(i,j)+Y(i+1,j)+Y(i-1,j+1)+Y(i,j+1)+Y(i+1,j+1);
h(i,j)=g1(i,j)./9;
end
end
Y=uint8(h) ;
figure,
imshow (Y)
title('平均值滤波斑点');
Z=wiener2(J,[5,5]); %维纳滤波椒盐
figure,
imshow (Z)
title('维纳滤波椒盐');
G=wiener2(Q,[5,5]); %维纳滤波高斯
figure
imshow (G)
title('维纳滤波高斯');
T=wiener2(Q,[5,5]); %维纳滤波斑点
figure
imshow (T)
title('维纳滤波斑点');
Imsz = size(I);
ngrid = Imsz(1)*Imsz(2); %求出原图像的大小
I=double(reshape(I,1,ngrid));
F=double(reshape(F,1,ngrid));
g_mean = mean(I); %求出图像的平均值
g_max = max(I); %求出图像的最大值
sqr_err = (I-F)*(I-F)';
MSE = sqr_err/ngrid %均方误差
SNR = 10.0*log10((I-g_mean)*(I-g_mean)'/sqr_err)%信噪比
PSNR = 10.0*log10(g_max*g_max*ngrid/sqr_err)%峰值信噪比
matlab提示错误说，
错误使用 reshape
要执行 RESHAPE，请勿更改元素数目。
请问大神应该怎么改啊？
)
