杨宇明《数字图像处理》《数字图像处理》书中代码.docVIP

例1-1 m=uint8(zeros(128,128,1,27)); for i=1:27 [m(:,:,:,i),map]=imread(mri.tif,i); end montage(m,map) 例2-1 a=imread(cameraman.tif); %读入图像 figure(1);imshow(); b1=a+50; %b1=a+45图像灰度值增加figure(2);imshow(b1); b2=1.2*a; %b=1.2*a图像对比度增大 figure(3);imshow(2) b3=0.65*a; %b=0.65*a图像对比度减少 figure(4);imshow(3); b4=-double(a)+255; %b4=-1*a+255，图像求补，注意把的类型转换为double figure(5);imshow(uint8(4)); %再把double类型转换为unit8 %读入lena图像 =imread(cameraman.tif); %读取原始图像 figure(1);imshow(a); %显示原始图像 xlabel ((a)原始图像); %显示函数的曲线图 x=1:255; y=x+x.*(255-x)/255; figure(2);plot(x,y);%绘制的曲线图的曲线图1=double(a)+0.006*double(a) .*(255-double(a)); figure(3);imshow(uint8(b1)); %显示非线性处理图像 xlabel((c)非线性处理效果); 例2-3 clear; %%读入图像 a=imread(cameraman.tif); imhist(a); title(原始cameraman图像的直方图); %%b=f(a) b1=1.25*double(a)+45; figure(2); imhist(uint8(b1)); title (变换后的直方图); 例2-4 clear; histgram=zeros(1,256); %生成直方图数组 cdf=zeros(1,256); [cm,map]=imread(cameraman.tif); [a,b]=size(cm); for i=1:a for j=1:b k=cm(i,j); histgram(k)=histgram(k)+1; end end %得到直方图 cdf(1)=histgram(1); for i=2:256 cdf(i)=cdf(i-1)+histgram(i); end for i=1:a %点运算 for i=1:b k=cm(i,j); cm_equ(i,j)=cdf(k)*256/(a*b); end end imshow(uint8(cm_equ)); figure(2);imhist(uint8(cm_equ)); 例2-5 I=imread(tire.tif); J=histeq(I); H=adapthisteq(I); figure(1),imshow(I); xlabel(原始图像); figure(2),imshow(J); xlabel(histeq均衡化); figure(3),imshow(H); xlabel(adapthisteq均衡化); 例2-6 a=imread(hill.jpg); s=size(a); b=double(a); c(:,:,1)=b(:,:,1)+b(:,:,2); c(:,:,2)=b(:,:,2); c(:,:,3)=b(:,:,3)-b(:,:,2); for i=1:s(1) for j=1:s(2) for k=1:s(3) if c(i,j,k)0 c(i,j,k)=0; end if c(i,j,k)255 c(i,j,k)=255; end end end end c=uint8(c); subplot(121);imshow(a);xlabel(原始图像); subplot