基于数字图像处理技术的答题卡识别方法(图像处理课程练习)

XX大学2013-2014学年第二学期课程考核

《图像处理》综合设计报告

基于数字图像处理技术的答题卡识别方法

学号 ___________________

姓名 ___________________

班级 ________________

日期 ______________________

本人郑重声明：本人认真、独立完成了查找资料、编写程序、撰写报告等考核任务。

签字:

日期:

随着科技的发展，电子与计算机技术的进步，答题卡的出现大大减轻教学工作者们批改试卷的工作量。答题卡是光标阅读机输入信息的载体，是配套光标阅读机的各种信息录入表格的总称。答题卡将用户

需要的信息转化为可选择的选项，供用户涂写。OMR是用光学扫描

的方法来识别按一定格式印刷或书写的标记，并将其转换为计算机能接受的电信号的设备，并根据信息点的涂与未涂和格式文件设置将信息还原。因此，如何将答题卡填涂的黑色区域识别出来并使用计算机进行处理是极为关键的。本论文探索了有效识别答题卡的方法，以matlab为工具，基于数字图像处理技术对答题卡填涂区域进行了识别，并对识别的结果进行了处理，得到了结果。本论文利用Hough

变换的直线检测技术检测图像的倾斜度，判断图像是否倾斜，对存在倾斜的图像进行旋转校正。最终实现答题卡答案的定位和检测。论文使用像素统计方法进行识别，利用黑白颜色灰度值的巨大差异对二值图像的灰度值进行累加并进行阈值判定，识别错误效率极低，能够准确的识别答题卡的涂卡标记。

关键词：Hough变换，答题卡识别，matlab,

一、设计任务、目的和要求

本设计以matlab为工具利用数字图像处理技术对答题卡进行了识别，并对识别结果进行了处理。注意到答题卡在采集图像的过程中由于各种原因可能会产生图像倾斜、水平或垂直错位，要进行正确识别首先必须对其进行校正，再依据像素

检索技术进行识别。

二、总体方案设计

说明系统运行环境，编程软件平台，编码算法原理，算法流程图设计

本系统运行在当今主流的Win dows7系统，使用较新的MATLAB2012b进行设计。在进行旋转校正时，先使用Hough变换检测出答题卡边缘直线，调用lines函数的参数

得到倾斜角度并进行旋转变换从而消除答题卡倾斜状态。将扫描部分分为

7个扫描区域，分别使用像素检索技术，对各个检索区域内的灰度值进行累加，并对累加和进行阈值判定，从而得出检索区域颜色，完成判定。最后再与标准答案进行比对，得出成绩。

流程图:

阈值判定，得出所检测区域颜色并赋值

三、程序设计

clc;clear;close;

A=imread('原图.jpg');%读取图像

imshow(A),title('原图');

gray=rgb2gray(A);

bw=edge(gray,'ca nn y');%ca nny算子边缘检测得到二值边缘图像

[h,t,r]=hough(bw,'RhoResolutio n',0.5,'ThetaResolutio n',0.5); %Ho ugh变换

figure,imshow(imadjust(mat2gray(h)),'XData',t,'YData',r,'l ni tialMag ni ficatio n','fit'),titl e('Hough变换矩阵');%显示Hough变换矩阵

xlabel('\theta'),ylabel('\rho');

axis on, axis no rmal,hold on;

P=houghpeaks(h,2);

x=t(P(:,2));y=r(P(:,1));

plot(x,y,'s','color','r'),title('');%获取并标出参数平面的峰值点

lin es=houghli nes(bw,t,r,P,'FillGap',5,'Mi nlen gth',7);%检测图像中的直线段

figure,imshow(gray);

hold on;

max_le n=0;

for i=1:le ngth(li nes)

xy=[li nes(i).po in t1;li nes(i).po in t2];

plot(xy(:,1),xy(:,2),'LineWidth',2,'Color','g');% 用绿色线段标注直线段

plot(xy(:,1),xy(1,2),'x','L in eWidth',2,'Color','y');

plot(xy(:,1),xy(2,2),'x','LineWidth',2,'Color','r');% 标注直线段端点

end

x=li nes.theta;%由与图像边缘平行的直线段的斜率得到整个图像旋转的角度

B=imrotate(gray,x);% 图像修正figure,imshow(B),title('旋转后图像');

T=graythresh(B);%使用OTSU方法获得阀值T result=im2bw(B,T);% 二值化图像figure,imshow(result),title('二值化后结果');

%扫描区域：在本程序中，有7个扫描区域，分别对应题号为

1-5,6-10,11-15,16-20,21-25,26-30,31-35这几个区域

n=52;m=244;%首先测得答题区域第一题A选项左上角的坐标(52,244) s=60;t=90;%s 为上下两个扫描区域的距离，t为左右两个扫描区域的距离。P=4;Q=5;%P为字母项的个数(这里有ABCD共有4项)，Q为每个小区域选项数。

in terval_le ngth=14;%左右相邻两题填涂区域左边线的距离

in terval_width=9;%上下相邻两题填涂区域上边线的距离

len gth=11;%填涂区域的长度

width=6;%填涂区域的宽度

a2=zeros(4,5);%初始化灰度值统计矩阵

%第1-5题区域的检测

for P=1:4

for Q=1:5

%m=244;n=52;%第1-5题区域的检测的起始坐标(即第一题A选项左上角的坐标)

a仁result(m+(P-1)*(interval_width+1) :m+(P-1)*(interval_width+1)+width,n+(Q-1)* (in terval_le ngth+1): n+(Q-1)*(i nterval_le ngth+1)+le ngth);

-%得到矩形区域内各像素的灰度值

sum仁sum(sum(a1));%计算统计的矩形区域内像素灰度值之和

a2(P,Q)=sum1;%多次循环后得到20个统计区域的灰度值，并依次放在a2矩阵中

end

end

a2(a2<80)=1;

a2(a2>=81)=0;%对灰度值统计矩阵里的数值进行处理，大于某一阀值的值赋1, 其余的赋0。

%在此程序中，涂黑则相应数值为1。

result1_5=a2;%存储1-5题的结果

%第6-10题区域的检测

for P=1:4

for Q=1:5

a仁result(m+s*1+(P-1)*(i nterval_width+1) :m+s*1+(P-1)*(i nterval_width+1)+width, n+(Q-1)*(i nterval」e ngth+1): n+(Q-1)*(i nterval」e ngth+1)+le ngth);

sum仁sum(sum(a1)); a2(P,Q)=sum1;

end

end

a2(a2<80)=1;

a2(a2>=81)=0; result6_10=a2; %存储6-10 题的结果

%第11-15题区域的检测

for P=1:4

for Q=1:5

a仁result(m+(P-1)*(interval_width+1):m+(P-1)*(interval_width+1)+width,n+t*1+(Q -