第3章MATLAB与图像处理幻灯片.pptVIP

* * 例3.1 figure, imhist(b,64); %作b的灰度直方图，64 可改为128, 256等 * * 例3.1 figure, imhist(b, 128); %作b的灰度直方图，64 可改为128, 256等 * * 例3.1 figure, imhist(b, 256); %作b的灰度直方图，64 可改为128, 256等 * * 例3图像增强 例3.2直方图均衡化 a=imread(lily.tif); b=a(:,:,1); imshow(b); * * 例3.2 c=histeq(b);　 %对b进行直方图均衡化 imshow(c) ; %(左边的图未均衡化,右边的图已均衡化) * * 例3.2 figure, imhist(b,64);　 %作b的灰度直方图，64 可改为128, 256等 * * 例3.2 figure, imhist(b, 128);　 %作b的灰度直方图，64 可改为128, 256等 * * 例3.2 figure, imhist(b, 256);　 %作b的灰度直方图，64 可改为128, 256等 * * 例3图像增强 例3.3邻域平均法线性平滑滤波 a=imread(lily.tif); b=a(:,:,1); imshow(b);%见左图 K1=filter2(fspecial(average,3),b)/255; %进行3*3均值滤波 imshow(K1);%见右图 * * 例3.3 K2=filter2(fspecial(average,5),b)/255; %进行5*5均值滤波 imshow(K2);%见右图,左图未滤波 * * 例3.3 K3=filter2(fspecial(average,7),b)/255; %进行7*7均值滤波 imshow(K3);%见右图,左图未滤波 * * 例3图像增强 例3.4中值滤波 a=imread(lily.tif); b=a(:,:,1); imshow(b);%见左图 d=medfilt2(b,[ 3 3]);% [3,3]为滤波模板的规格 imshow(d);%见右图 * * 例4 图像重建(参考) a=imread(1.jpg); b=a(:,:,1); r=radon(b,0:180); %从0度开始，对b进行180个角度扫描，步长为1度 ir=iradon(r,0:180); %对投影数据r进行重建 * * 例5 图像分析(参考) 例5.1 边缘提取 以下程序的功能为：仅根据Euler数这一个特征，判断手写数字0,1,8等。适当地多取几个特征，很容易将程序扩展为判断0，1，2，3，4，5，6，7，8，9等10个数字的场合。以下图像分析的实例非常简单，以该思想为基础，可进一步深化算法，设计出更复杂的高效率字符判读算法。 (使用wavread,wavplay两函数，其目的是为了将所判定结果，以声音形式输出。) * * 图像分析(参考) %y0= wavread(s0); %y1= wavread(s1); %y8 = wavread(s8); Fs = 11025; b=imread(b018.jpg); %用图画程序写一个0、1或8，保存为b018.jpg（图像获取） bb=b(:,:,1); %取b的一个分量 bw=bb100; %将b 二值化（相当于图像分割操作） c=bweuler(bw); %求bw的Euler数（相当于图像描述，特征提取操作） %以下部分为图像识别操作，（相当于分析与分类操作） switch c case -1 cc=8; % wavplay(y8,Fs); case 0 cc=0; %wavplay(y0,Fs); case 1 cc=1; %wavplay(y1,Fs); otherwise ccc=999; %ccc=999,不可识别 end cc ccc %显示判断结果 * * 作业 1.翻译第3章课后参考英文. * * 实验 见实验任务书 * * END 参考资料： 1.《数字图像处理与分析》，龚声蓉等编著，清华大学出版社 2.《MATLAB7.0图像处理》，贺兴华等编著，人民邮电出版社 3.数字图像处理第二版中文版（冈萨雷斯） 4. Digital Image Processing Using Matlab (Gonzalez) 5. MATLAB6.5辅助图像处理 * 逆矩阵的定义,对于方阵A,.若有方阵B,使得: AB=BA,则称A可逆,称