冈萨雷斯数字图像处理中文版_第二章讲述.ppt

* * * * * * 三、灰度的对数变换 c为尺度比例常数，s为源灰度值，t为变换后的目标灰度值。k为常数。灰度的对数变换可以增强一幅图像中较暗部分的细节，可用来扩展被压缩的高值图像中的较暗像素。广泛应用于频谱图像的显示中。 Warning：log函数会对输入图像矩阵s中的每个元素进行 操作，但仅能处理double类型的矩阵。而从图像文件中得到的 图像矩阵大多是uint8类型的，故需先进行im2double数据类型 转换。 * * I=imread(nir.bmp);%读入图像 F=fft2(im2double(I));%FFT F=fftshift(F);%FFT频谱平移 F=abs(F); T=log(F+1);%频谱对数变换 figure;imshow(F,[]);title(未经变换的频谱); figure;imshow(T,[]);title(对数变换后); * * 图像的几何变换 二、图像镜像 B=imtransform(A,TFORM,method); TFORM=makeform(transformtype,Matrix);%空间变换结构 Method合法值 含义 ‘bicubic’ 双三次插值 ‘bilinear’ 双线性插值 ‘nearest’ 最近邻插值 参数transformtype指定了变换的类型，常见的’affine’为二维或 多维仿射变换，包括平移、旋转、比例、拉伸和错切等。 Matrix为相应的仿射变换矩阵。 * * A=imread(nir.bmp); [height,width,dim]=size(A); tform=maketform(affine,[-1 0 0;0 1 0;width 0 1]); B=imtransform(A,tform,nearest); tform2=maketform(affine,[1 0 0;0 -1 0;0 height 1]); C=imtransform(A,tform2,nearest); figure;imshow(A); figure;imshow(B);imwrite(B,nir水平镜像.bmp); figure;imshow(C);imwrite(B,nir垂直镜像.bmp); 原图像 水平镜像图像 垂直镜像图像 * * A=imread(nir.bmp); tform=maketform(affine,[0 1 0;1 0 0;0 0 1]); B=imtransform(A,tform,nearest); figure;imshow(A); figure;imshow(B);imwrite(B,nir转置后图像.bmp); /3、图像的几何变换 三、图像转置 * * 四、图像中心旋转 B=imrotate(A,angle,method,’crop’); angle为旋转角度，正值为逆时针旋转。可选参数method为imrotate函数 指定插值方法。‘crop’选项会裁减旋转后增大的图像，保持和原图像 同样大小。 A=imread(nir.bmp); B=imrotate(A,30,nearest,crop); figure;imshow(B);imwrite(B,逆时针中心旋转30度.bmp); 逆时针30度 /3、图像的几何变换 * * 空间域图像增强 一、噪声添加 h=imnoise(I,type,parameters); type为噪声类型，合法值如下： ‘gaussian’ 高斯白噪声：幅度为高斯分布，功率谱均匀分布 ‘saltpepper’ 黑点如胡椒，白点如盐粒。由图像传感器、传输信道、解码处理、图像切割等产生的黑白相间的亮暗点噪声。 添加高斯白噪声 添加椒盐噪声 * * /4、空间域图像增强 二、空间域滤波 滤波过程就是在图像f(x,y)中逐点移动模板，使模板中心和点(x,y)重合，滤波器 在每一点(x,y)的响应是根据模板的具体内容并通过预先定义的关系来计算的。 W(-1,-1) W(-1,0) W(-1,1) W(0,-1) W(0,0) W(0,1) W(1,-1) W(1,0) W(1, 1) f(x-1,y-1) f(x-1,y) f(x-1,y+1) f(x,y-1) f(x,y) f(x,y+1) f(x+1,y-1) f(x+1,y) f(x+1,y+1) * * B=imfilter(f,w,option1,option2,…); f 为要进行滤波操作的图像。 w为滤波操作使用的模板，为一个二维数组，可自己定义。 option1……是可选项，包括： 1、边界选项（’symmetric’、’replicate’、’circular’） 2、尺寸选项（’sam