用直接灰度变换改变图像（求反，增强对比度，动态范围压.doc

《用直接灰度变换改变图像（求反，增强对比度，动态范围压.doc》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、3.用直接灰度变换改变图像（求反，增强对比度，动态范围压缩，灰度切分）求反>>I=imread('moon.tif');>>imshow(I);>>level=graythresh(I);>>I1=im2bw(I,level);>>figure,imshow(I1);>>x1=128;>>y1=128;>>figure,plot([0,x1],[y1,0])>>axistight,xlabel('x'),ylabel('y')>>title('变换曲线')>>k=y1/x1;>>[m,n]=size(I);>>J=double(I);>>fori=1:mforj=1:nx=J(i,j

2、);y(i,j)=0;if(x>=0)&(x<=x1)y(i,j)=y1-k*x;elsey(i,j)=0;endendend>>figure,imshow(mat2gray(y))原图二值化后的图像求反后的图像增强对比度源代码：>>I=imread('pout.tif');>>imhist(I);>>I1=imadjust(I,[0.30.7],[]);>>figure,imshow(I1);>>imhist(I1);>>figure,imshow(I1);>>figure,imshow(I);原图原图的直方图原图增加对比度增加对比度后的直方图灰度切分:>>I=imread('pi

3、c1.jpg');>>J=roicolor(I,100,200);>>figure,imshow(I),figure,imshow(J);原图切分后的图动态范围压缩>>I=imread('tire.tif');>>imshow(I);>>c=64/log(256);>>x=0:1:255;>>y=c*log(1+x);>>figure,plot(x,y)>>axistight,xlabel('f'),ylabel('g')>>title('变换曲线')>>[m,n]=size(I);>>X2=double(I);>>fori=1:mforj=1:n;g(i,j)=c*log(X2(i

4、,j)+1);endend>>figure,imshow(mat2gray(g))原图动态范围压缩后的图