matlab实现对机读卡识别

%function testsegment clc;clear; I = imread('E:/8.jpg'); if size(I,3)>1 % 转灰度图 I = rgb2gray(I); end I = im2double(I); H=fspecial('average',floor(size(I)/8)); %% 用滤波函数求图片背景 avgI=imfilter(I,H,'replicate'); I_s = 1 - avgI + I; %%用原图减去背景去阴影 I_s = I_s.^16; bws = im2bw(I_s,graythresh(I_s)); %二值化 bws = ~bws; %%图像取反 figure; imshow(bws); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Ie=edge(I,'canny'); S1 = regionprops(Ie,'BoundingBox','PixelIdxList'); max_area = 0; for i = 1:length(S1) %%%找图像最大连通域ROI area = S1(i).BoundingBox(3)*S1(i).BoundingBox(4); if area>max_area max_area = area; pos = i; end end bbox = S1(pos).BoundingBox; figure; imshow(I(bbox(2):bbox(2)+bbox(4),bbox(1):bbox(1)+bbox(3))) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% ROIbw=zeros(size(I)); ROIbw(S1(pos).PixelIdxList)=Ie(S1(pos).PixelIdxList); corners = detectHarrisFeatures(ROIbw); %%%角点检测 [c_number,ww]=size(corners); imshow(I); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% hold on; plot(corners.selectStrongest(c_number));%%%%%角点标记 hold off hc = corners.selectStrongest(c_number); hcxy = hc.Location; loccir = mean(hcxy); %%%%图像中点 quadrant=[]; %%%%象限 center_x=loccir(:,1); center_y=loccir(:,2); for i=1:c_number quadrant(i)=1; end for i=1:c_number %%%%%四个象限初始化 if hcxy(i,1)-center_x>0 && hcxy(i,2)-center_y<=0 quadrant(i)=1; end if hcxy(i,1)-center_x<=0 && hcxy(i,2)-center_y<0 quadrant(i)=2; end if hcxy(i,1)-center_x<0 && hcxy(i,2)-center_y>=0 quadrant(i)=3; end if hcxy(i,1)-center_x>=0 && hcxy(i,2)-center_y>0 quadrant(i)=4; end end point_len=0; %四个点 max_len=0; quadrant1=0; quadrant2=0; quadrant3=0; quadrant4=0; for i=1:length(quadrant) %%计算四个象限距离中心点最远的四个点 if quadrant(i)==1 %%第一象限 point_len=sqrt((hcxy(i,1)-center_x)^2+(hcxy(i,2)-center_y)^2); if max_len<point_len max_len=point_len; quadrant1=i; end end end point_len=0; max_len=0; for i=1:length(quadrant) if quadrant(i)==2 %%第二象限 point_len=sqrt((hcxy(i,1)-center_x)^2+(hcxy(i,2)-center_y)^2); if max_len<point_len max_len=point_len; quadrant2=i; end end end point_len=0; max_len=0; for i=1:length(quadrant) if quadrant(i)==3 %%第三象限 point_len=sqrt((hcxy(i,1)-center_x)^2+(hcxy(i,2)-center_y)^2); if max_len<point_len max_len=point_len; quadrant3=i; end end end point_len=0; max_len=0; for i=1:length(quadrant) if quadrant(i)==4 %%第四象限 point_len=sqrt((hcxy(i,1)-center_x)^2+(hcxy(i,2)-center_y)^2); if max_len<point_len max_len=point_len; quadrant4=i; end end end figure; imshow(I); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% hold on plot( center_x,center_y ,'r+'); %%标记四个角点 plot( hcxy(quadrant1,1),hcxy(quadrant1,2) ,'r+'); plot( hcxy(quadrant2,1),hcxy(quadrant2,2) ,'g+'); plot( hcxy(quadrant3,1),hcxy(quadrant3,2) ,'b+'); plot( hcxy(quadrant4,1),hcxy(quadrant4,2) ,'g+'); ROIpoints(1,1)= hcxy(quadrant1,1); ROIpoints(1,2)= hcxy(quadrant1,2); ROIpoints(2,1)= hcxy(quadrant2,1); ROIpoints(2,2)= hcxy(quadrant2,2); ROIpoints(3,1)= hcxy(quadrant3,1); ROIpoints(3,2)= hcxy(quadrant3,2); ROIpoints(4,1)= hcxy(quadrant4,1); ROIpoints(4,2)= hcxy(quadrant4,2); X=[510 10 10 510]; Y=[10 10 510 510]; %%%变换后的图像模板大小 tform = fitgeotrans(ROIpoints,[X' Y'],'Projective'); %%图像变换矫正 outputView = imref2d(size(I)); Is = imwarp(I,tform, 'OutputView', outputView); Is=Is(X(2):X(1),Y(2):Y(3)); figure; imshow(Is); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% H=fspecial('average',floor(size(Is)/8)); %%矫正后二值化 avgI=imfilter(Is,H,'replicate'); I_s = 1 - avgI + Is;% I_s = I_s.^16; bws = im2bw(I_s,graythresh(I_s)); bws=~bws; figure; imshow(bws); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% I_tem=bws(172:361,6:497); %%选择题区域提取 % I_show=I_tem;%%%%%%%%%%%%%%%%%% I_tem= bwareaopen(I_tem,50); %%过滤噪声点 figure; imshow(I_tem); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%% figure; imshow(I_tem); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%% choose %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% answer=cell(1,50); choose={'A','B','C','D'}; %%%%%%选择题识别，遍历像素点 for i=1:3 for j=1:4 for k=1:5 for p=(j-1)*124+32:26:(j-1)*124+32+(26*3) is_done=I_tem((i-1)*66+4+(k-1)*13,p); if is_done==1 a_num=(i-1)*20+(j-1)*5+k; answer(a_num)= choose(floor((p-((j-1)*124+32))/26)+1); end end end end end bws=bwareaopen(bws,10); %%%试卷号和身份证号区域提取 S2 = regionprops(bws,'BoundingBox','Area'); bwArea=cat(1,S2.Area); [bwAreasort,idx]=sort(bwArea); bbx1=round(S2(idx(end-2)).BoundingBox); bbx2=round(S2(idx(end-1)).BoundingBox); bwss1=bws(bbx1(2):bbx1(2)+bbx1(4),bbx1(1):bbx1(1)+bbx1(3)); bwss2=bws(bbx2(2):bbx2(2)+bbx2(4),bbx2(1):bbx2(1)+bbx2(3)); figure(7); subplot(121);imshow(bwss1);title('paper number'); subplot(122);imshow(bwss2);title('ID'); L1 = bwlabel(bwss1); L2 = bwlabel(bwss2); S1 = regionprops(L1,'BoundingBox'); S2 = regionprops(L2,'BoundingBox'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%% paper num %%%%%%%%%%%%%%%%%%%%%%%% %文件操作 cd('E:\'); % 设置当前目录 if ~exist('paper') rmdir('paper','s'); % 删除文件夹 end folderName = 'paper'; mkdir(folderName); % 新建一个文件夹 for i=2:length(S1) %%数字分割 s1 = sprintf('E:/paper/num(%d).jpg',i-1); tem_img = imcrop(bwss1,S1(i).BoundingBox); imwrite(~tem_img,s1,'jpg'); end paper_value=[]; for id_index=1:6 %%数字分割 s1 = sprintf('E:/paper/num(%d).jpg',id_index); tem_img=imread(s1); [m,n]=size(tem_img); if n<7 paper_value(id_index)=1; continue; end tem_img=imresize(tem_img,[25,15]);%图像大小调整 thresh = graythresh(tem_img); %二值化 bw = im2bw(tem_img,thresh); flag1=0;%%%%分别检测数字的