用matlab程序实现同时对比度实验、空间分辨率变化效果、实验幅度分辨率变化效果

1、作业11.用matlab程序实现同时对比度实验。程序代码如下：clearm=256;n=256;groun d_gray1=200;X1=on es( m,n )*gro un d_gray1; groun d_gray2=100;X2=on es(m, n)*gro un d_gray2; before_gray=150;m0=100;n 0=100;step=50;X1(m0:m0+step ,n0:n0+step)=before_gray; figure;imshow(ui nt8(X1);X2(m0:m0+step ,n0:n0+step)=before_gray; figure;im

2、show(ui nt8(X2);得到的图像如下:2.用matlab程序实验空间分辨率变化效果。clearX=imread( picture1.jpg);figure;imshow(ui nt8(X);m, n=size(X);X1=X(1:2:m,1:2: n);figure;imshow(ui nt8(X1);imwrite(X1, picture1.jpg);得到的图形如下3.用matlab程序实验幅度分辨率变化效果程序代码如下:clearX=imread( picture1.jpg ); mm,nn=size(X);for i=1:mm for j=1:nn X128(i,j)=2*(c

3、eil(X(i,j)/2);end end for i=1:mmfor j=1:nnX64(i,j)=4*(ceil(X(i,j)/4); end end for i=1:mmfor j=1:nnX32(i,j)=8*(ceil(X(i,j)/8); end endfiguresubplot(2,2,1); imshow(uint8(X);title( 256 级 );subplot(2,2,2); imshow(uint8(X128);title( 128 级 );subplot(2,2,3); imshow(uint8(X64);title( 64 级 );subplot(2,2,4);

4、imshow(uint8(X32);title( 32 级 );得到的图像如下);none);作业21.用matlab作出图像的直方图程序代码如下:clear;A = imread( picture2.jpg I = rgb2gray(A);add =;tab1 = zeros(1,256);for n = 0:255;X = I = n;add = add; sum(sum(X); end;a, b = size(I); final = add / (a*b); figure;stem(fi nal,Markerfigure(2);imshow(I);得到的图像如下070 6、0.5040.

5、302QL0501001502002503002.用matlab实现图像的直方图均衡化程序代码如下：clear allm=imread( picture2.jpg);A = rgb2gray(m);m, n=size(A); B=zeros(size(A); l=m*n; r=zeros(1,256); y=zeros(1,256); A=double(A); for i=1:m for j=1:n r(A(i,j)+1)=r(A(i,j)+1)+1;endendr=r./l;for k=1:256for j=1:k y(k)=y(k)+r(j);endendfor k=1:256y(k)=f

6、loor(255*y(k)+0.5); endfor i=1:mfor j=1:nB(i,j)=y(A(i,j)+1);endendA=uint8(A);B=uint8(B); subplot(2,2,1) imshow(A)title(原图像)subplot(2,2,2)imhist(A) title(原图像的直方图) subplot(2,2,3) imshow(B)title( 经过均衡处理的图像 )subplot(2,2,4)imhist(B)title( 经过均衡处理后该图像的直方图 )得到的图像如下凹& !?i-j ?坯过均flsih宦怡曲 ms作业31.用matlab实现多幅图像平

7、均去高斯白噪声。程序代码如下：clear;);A=imread( picture2.jpg l=rgb2gray(A);imshow(l);a,b=size(I);no ise=ra ndn( a,b)*30; X=double(l)+no ise; Y=zeros(a,b);figure(2); imshow(ui nt8(X);for n=1:30no ise=ra ndn( a,b)*30; X=double(I)+no ise; Y=Y+double(X);end;Y=Y/n; figure(3);imshow(ui nt8(Y)得到的图像如下:原图加了高斯白噪声的图像:去噪后的图像：

8、2.用matlab实现均值滤波去除高斯白噪声.程序代码如下：clear;A=imread( picture2.jpg);l=rgb2gray(A);figure;imshow(l);a,b=size(I);no ise=ra ndn( a,b)*30;I=double(I)+no ise;I2=zeros(a+2,b+2);I3=zeros(a,b);for n=1:afor m=1:bI2( n+1,m+1)=l( n,m);end ;end;figure(2);imshow(ui nt8(I2);for n=2:afor m=2:bI3( n-1,m-1)=I2( n-1,m-1)+I2(

9、 n-1,m)+I2( n-1,m+1)+I2( n,m- 1)+I2( n,m)+I2( n,m+1)+I2( n+1,m-1)+I2( n+1,m)+I2( n+1,m+1) /9;end ;end;figure (3); imshow(ui nt8(l3);得到的图像如下:原图加高斯白噪声图像均值滤波去噪后图像3.用matlab实现中值滤波去除脉冲噪声。 程序代码如下：clear;A=imread( picture2.jpg);l=rgb2gray(A);figure(1);imshow(l);X=zeros(1:9);a,b=size(I);Y=zeros(a+2,b+2);Z=zeros(a,b);for n=1:afor m=1:bY( n+1,m+1)=I( n,m);end ;end;for n=0:100Y(ro un d(ra nd*a),ro un d(ra nd*b)=0;end;for n=0:200Y(ro un d(ra nd*a),ro un d(ra nd*b)=255;end;figure(2)imshow(ui nt8(Y);Y=double(Y);for n=2:afor m=2:bX=Y( n-1,m-1),Y( n-1,m),Y( n-1,m+1)