如何修改matlab中的语句,求大神帮忙看一下这个语句怎么改!!!!!

这篇博客探讨了图像处理中的噪声类型,如椒盐噪声、高斯噪声和斑点噪声,并展示了如何使用MATLAB进行中值滤波和平均值滤波来去除这些噪声。此外,还应用了维纳滤波器进行进一步的噪声抑制。最后,博客计算了滤波后的均方误差、信噪比和峰值信噪比等质量指标。
摘要由CSDN通过智能技术生成

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

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,请勿更改元素数目。

请问大神应该怎么改啊?

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值