基于MATLAB的数字图像运算.doc

点运算： 线性点运算： a=imread(liftingbody.png); subplot(3,2,1);imshow(a); b1=a+50; subplot(3,2,2);imshow(b1); b2=1.2*a; subplot(3,2,3);imshow(b2); b3=0.6*a; subplot(3,2,4);imshow(b3); b4=-double(a)+255; subplot(3,2,5);imshow(b4); 非线性点运算; a=imread(cameraman.tif); subplot(1,3,1);imshow(a); x=1:255; y=x+x.*(255-x)/255; subplot(1,3,2);plot(x,y); b1=double(a)+0.006*double(a).*(255-double(a)); subplot(1,3,3);imshow(uint8(b1)); 图像代数运算： 加法运算： i=imread(rice.png); subplot(2,2,1);imshow(i); j=imread(cameraman.tif); subplot(2,2,2);imshow(j); k=imadd(i,j,uint16); subplot(2,2,3);imshow(k,[]); f=imadd(i,45); subplot(2,2,4);imshow(f); 利用imnoise（）对噪声进行相加运算： a=imread(eight.tif); a1=imnoise(a,gaussian,0,0.006); a2=imnoise(a,gaussian,0,0.006); a3=imnoise(a,gaussian,0,0.006); a4=imnoise(a,gaussian,0,0.006); k=imlincomb(0.25,a1,0.25,a2,0.25,a3,0.25,a4); subplot(1,3,1);imshow(a); subplot(1,3,2);imshow(a1); subplot(1,3,3);imshow(k); 图像减法运算： i=imread(cameraman.tif); subplot(2,2,1);imshow(i); back=imopen(i,strel(disk,15)); subplot(2,2,2);imshow(back); i1=imsubtract(i,back); subplot(2,2,3);imshow(i1); i2=imsubtract(i,45); subplot(2,2,4);imshow(i2); 图像乘法运算： i=imread(autumn.tif); subplot(2,2,1); imshow(i); i16=uint16(i); j=immultiply(i16,i16); subplot(2,2,2);imshow(j); j2=immultiply(i,0.5); subplot(2,2,3);imshow(j2); 图像除法运算： i=imread(autumn.tif); subplot(2,2,1);imshow(i); back=imopen(i,strel(disk,15)); i1=imdivide(i,back); subplot(2,2,2);imshow(i1,[]); i2=imdivide(i,2); subplot(2,2,3);imshow(i2); 绝对值运算： i=imread(cameraman.tif); subplot(1,3,1);imshow(i); i1=filter2(fspecial(gaussian),i); j=uint8(i1); subplot(1,3,2);imshow(j); k=imabsdiff(i,j); subplot(1,3,3);imshow(k,[]); 线性组合运算： i=imread(cameraman.tif); j=imlincomb(2.5,i); subplot(2,2,1);imshow(i); subplot(2,2,2);imshow(j); j1=uint8(filter2(fspecial(gaussian),i)); k=imlincomb(1,i,-1,j1,128); subplot(2,2,3);imshow(k); j2=imread(cameraman.tif); k2=imlincomb(1,i,1,j2,ui