該樓層疑似違規已被系統折疊?隱藏此樓查看此樓

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，請勿更改元素數目。

請問大神應該怎么改啊？

總結

以上是生活随笔為你收集整理的如何修改matlab中的语句,求大神帮忙看一下这个语句怎么改！！！！！的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。