数据图像处理.docVIP

例4.1 代码及结果如下： I=imread(l2.jpg); imshow(I); I=double(I); [M,N]=size(I); for i=1:M for j=1:N if I(i,j)=30 I(i,j)=I(i,j); elseif I(i,j)=150 I(i,j)=(200-30)/(150-30)*(I(i,j)-30)+30; else I(i,j)=(255-200)/(220-150)*(I(i,j)-150)+220; end end end figure(2);imshow(uint8(I)); 截图如下： 例4.2 代码及结果如下： I=imread(l3.jpg) figure;imshow(I); I=double(I); I2=41*log(1+I); I2=uint8(I2); figure;imshow(I2); 截图如下： 例4.3 代码及结果如下： I=imread(l2.jpg); figure;imshow(I); I=double(I) [M,N]=size(I); for i=1:M for j=1:N if I(i,j)=50 I(i,j)=40; elseif I(i,j)=180 I(i,j)=220; else I(i,j)=40; end end end I=uint8(I); figure;imshow(I); 截图如下： 例4.4 代码及结果如下： I=imread(l2.jpg); imshow(I); I=double(I); [M,N]=size(I); for k=1:8 J=zeros(M,N); for i=1:M for j=1:N temp=I(i,j); s1=0;s2=0; range=[k:-1:1]; for d=range s1=2^(8-d)+s1;s2=2^(8-d+1); if temp=s1temps2; J(i,j)=255;break; end end end end J=uint8(J); figure;imshow(J); end 截图如下： 例4.6 代码及结果如下： I=imread(cameraman.tif); figure subplot(221);imshow(I); subplot(222);imhist(I); I1=histeq(I); figure; subplot(221);imshow(I1) subplot(222);imhist(I1) 截图如下： 例4.8 代码及结果如下： I=imread(cameraman.tif); [M,N]=size(I); for i=1:8:257 couts(i)=i; end Q=imread(cameraman.tif); N=histeq(Q,couts); figure subplot(221);imshow(N); subplot(222);imhist(N); axis([0 260 0 5000]) 截图如下： 例4.9 代码及结果如下： I=imread(circuit.tif) [M,N]=size(I); II1=zeros(M,N); for i=1:16 II(:,:,i)=imnoise(I,gaussian,0,0.01); II1=II1+double(II(:,:,i)); if or(or(i==1,i==4),or(i==8,i==16)); figure; imshow(uint8(II1/i)); end end 例4.10 代码及结果如下： I=imread(concordorthophoto.png ); J=imnoise(I,salt pepper,0.02); subplot(231),imshow(I);