数字图像处理

1、空间域图像增强Matlab实现,一、MATLAB数字图像处理基础 （一）读取图像 1.语法：Imread (filename) imread(filename); 具体实现 imread( D: filename.jpg) 2.给出一幅图像的行数和列数 size(f); M N=size(f); 3.显示一个数组的附加信息 whos f,（二）显示图像 1.语法 imshow(f,G); imshow(f,low high); imshow(f, ); 2.光标上显示图像的像素值 pixval; 3.同时显示两个窗口 figure;,4.显示不同类型的图像 imshow(X,map)%显示索引

2、图像 imshow(I);%显示灰度图像 BW=imread(001.jpg); imshow(BW);%显示二值图像 RGB=imread(001.jpg); imshow(RGB);%显示RGB图像,（三）保存图像 imwrite(f,001,tif); imwrite(f,001.jpg,quality,q);,（四）图像类型之间的变换 1.将灰度值限定在0，1内 mat2gray(f); 2.将灰度值限定在0，255的范围内 im2unit8(f); 3.抖动函数 BW=dither(I);%将灰度图像I抖动成二值图像,4.将灰度图像转化索引图像 X,map=gray2ind(I,n)

3、; %例 I=imread(001.jpg); X,map=gray2ind(I,6); imshow(X,map); 5.通过设定域值将灰度图象转化成索引图像 X=grayslice(I,n); X=grayslice(I,v); 6.将灰度图像变成RGB图像 RGB=cat(3,I,I,I);,7.通过设置亮度域值将真彩图像、索引图像、灰度图像转化成二值图像 BW=im2bw(I,level); BW=im2bw(X,map,level); BW=im2bw(RGB,level); 8.将索引图像转化成灰度图像 I=ind2gray(X,map) 9.将索引图像转化成真彩色图像 RGB=i

4、nd2rgb(X,map); 10.将一个矩阵转换成灰度图像 I=mat2gray(A,amin amax); I=mat2gray(A);,（五）图像的运算 imadd(I,J); immultiply(I,J); imdivide(I,J); k=imdivide(imadd(I,J),2);,（六）图像的缩放 imresize(I,1.25); imresize(I,100 150); imresize(I,100 150,bilinear); （七）旋转图像 imrotate(I,35,bilinear)；,（八）裁切图像 imshow 001.jpg; I=imcrop; imsho

5、w(I),二、亮度变换函数,1.imadjust函数 g=imadjust(f,low_in high_in,low_out high_out,gamma); 2.图像的负片 g=imcomplement(f); 3.对数和对比度拉伸变换 g=c*log(1+double(f) 4.将灰度值限定在0，1内 mat2gray(f); 5.将灰度值限定在0，255的范围内 im2unit8(f);,三、直方图处理,1.直方图的绘制 （1）基本语法：imhist(f,b) b的默认值为256 （2）归一化：p=imhist(f,b)/numel(f) 例：具体操作。,2.直方图均衡化 （1）基本语法

6、 g=histeq(f,nelv) nelv为指定灰度等级（64） 例： 3.直方图规定化 基本语法：g=histeq(f,hspec),四、空间滤波 1.线性滤波 语法： g=imfilter(f,w,filtering_mode,boudary_options,size_options) 参数： corr 相关 conv 卷积,P 用p值来填充，默认值是0 replicate 通过复制外边界的值来扩展 symmetric通过镜像反射其边界值 Circular 将图像看做是一个二维周期函数的一个周期来扩展 Full 图像大小 与扩展后的图像大小相同 same 与输入图像大小相同 通用语法 g

7、=imfilter(f, w,replicate) 举例,2.非线性滤波器 语法： g=colfilt(f,m n,sliding,fun,parameters) 填充 fp=padarray(f,r c,method,direction) 方法； replicate，symmetric，circular 方向： pre 在每一维的第一个元素前填充 post 在每一维的第一个元素前填充 both 在每一维的前一个元素和最后一个元素填充 举例演示,五、标准处理工具箱的标准空间滤波器 1.线性空间滤波器 生成滤波器的函数fspecial 语法：w= fspecial(type,parameter)

8、 类型 (1)average 语法 fspecial(average, r c) 默认3*3，一个数代表方形 (2)disk 语法 fspecial(average,r) 半径默认为5 (3)gaussian 语法 fspecial(average,r s,sig) 标准偏差为正，默认值为3*3 0.5,(4)laplacian 语法 fspecial(laplacian,alpha) 0 1 默认值为0.5 (5)log 高斯-拉普拉斯滤波器 语法 fspecial(log,r c, sig) 5*5 0.5 (6)motion 语法 fspecial(motion,len,theta) 9 0 (7)prewitt 语法 fspecial(prewitt) 垂体梯度 (8)sobel 语法 fspecial(sobel) (9)unsharp 语法 fspecial(unsharp,alpha) 3*3 0.2,2.非线性滤波器 生成非线性滤波器的函数是 ordfilt2，生成统计排序滤波器。 (1)基本语法 g=ordfilt2(f,order,domain) (2)最小值滤波器 g=ordfi