实验六 图像处理实际应用
一、实验目的
1、熟悉并掌握MATLAB图像处理工具箱的使用； 2、理解并掌握常用的图像处理技术。 二、实验环境
MATLAB 6.5以上版本、WIN XP或WIN2000计算机 三、实验内容
调试运行下列程序,分析程序,对每条语句给出注释,并显示最终执行结果。总结算法思想及优缺点.
I=imread('Car.jpg'); [y,x,z]=size(I); myI=double(I);
%%%%%%%%%%% RGB to HSI %%%%%%%% tic
%%%%%%%% Y 方向 %%%%%%%%%% Blue_y=zeros(y,1); for i=1:y for j=1:x
if((myI(i,j,1)<=30)&&((myI(i,j,2)<=62)&&(myI(i,j,2)>=51))&&((myI(i,j,3)<=142)&&(myI(i,j,3)>=119))) % 蓝色RGB的灰度范围
Blue_y(i,1)= Blue_y(i,1)+1; end end end
[temp MaxY]=max(Blue_y); % Y方向车牌区域确定 PY1=MaxY;
while ((Blue_y(PY1,1)>=5)&&(PY1>1)) PY1=PY1-1; end PY2=MaxY;
while ((Blue_y(PY2,1)>=5)&&(PY2 IY=I(PY1:PY2,:,:); %%%%%%%% X 方向 %%%%%%%%%% Blue_x=zeros(1,x); for j=1:x for i=PY1:PY2 if((myI(i,j,1)<=30)&&((myI(i,j,2)<=62)&&(myI(i,j,2)>=51))&&((myI(i,j,3)<=142)&&(myI(i,j,3)>=119))) Blue_x(1,j)= Blue_x(1,j)+1; end end end PX1=1; while ((Blue_x(1,PX1)<3)&&(PX1 while ((Blue_x(1,PX2)<3)&&(PX2>PX1)) PX2=PX2-1; end PX1=PX1-2; % 对车牌区域的修正 PX2=PX2+2; Plate=I(PY1:PY2,PX1-2:PX2,:); t=toc % 读取计时 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% figure,imshow(I); figure,plot(Blue_y);grid figure,plot(Blue_x);grid figure,imshow(IY); 添加注释和改正后的程序： I=imread('f:\\Car.jpg');%读取图片 [y,x,z]=size(I);%给定图片大小 myI=double(I);%返回双精度值 %%%%%%%%%%% RGB to HSI %%%%%%%% tic %计时开始 %%%%%%%% Y 方向 %%%%%%%%%% Blue_y=zeros(y,1);%一列全零矩阵 for i=1:y%给定i的范围 for j=1:x%给定j的范围 if((myI(i,j,1)<=30)&&((myI(i,j,2)<=62)&&(myI(i,j,2)>=51))&&((myI(i,j,3 )<=142)&&(myI(i,j,3)>=119))) % 蓝色RGB的灰度范围 Blue_y(i,1)= Blue_y(i,1)+1; % y矩阵加一 end end end [temp MaxY]=max(Blue_y); % Y方向车牌区域确定 PY1=MaxY;%y矩阵的最大元素 while ((Blue_y(PY1,1)>=5)&&(PY1>1))% 确定蓝色RGB Blue_y所在位置 PY1=PY1-1; %对车牌区域的修正，向上调整 end PY2=MaxY; %y矩阵的最大元素 while ((Blue_y(PY2,1)>=5)&&(PY2 PY2=PY2+1; %对车牌区域的修正向下调整 end IY=I(PY1:PY2,:,:);%行的范围 %%%%%%%% X 方向 %%%%%%%%%% Blue_x=zeros(1,x); %一行全一矩阵 for j=1:x%给定j的取值范围 for i=PY1:PY2%给定i的取值范围 if((myI(i,j,1)<=30)&&((myI(i,j,2)<=62)&&(myI(i,j,2)>=51))&&((myI(i,j,3)<=142)&&(myI(i,j,3)>=119)))%调整I中的范围 Blue_x(1,j)= Blue_x(1,j)+1; % Blue_x 的矩阵加一 end end end PX1=1;%当PX1等于1时 while ((Blue_x(1,PX1)<3)&&(PX1 PX1=PX1+1; %对车牌区域的修正向下调整 end PX2=x; %当PX2等于x时 while ((Blue_x(1,PX2)<3)&&(PX2>PX1)) %确定蓝色RGB Blue_x的位置 PX2=PX2-1; %对车牌区域的修正向下调整 end PX1=PX1-2; % 对车牌区域的修正 PX2=PX2+2; %对车牌区域的修正 Plate=I(PY1:PY2,PX1-2:PX2,:);%矩阵行列的范围 t=toc % 读取计时 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% figure,imshow(I);%创建图像窗口，显示图像I figure,plot(Blue_y);grid%创建图像窗口，绘制Blue_y图像，画出网格线 figure,plot(Blue_x);grid%创建图像窗口，绘制Blue_x图像，画出网格线 figure,imshow(Plate); %创建图像窗口，显示图像Plate t = 0.5630