✅博主简介：热爱科研的Matlab仿真开发者，修心和技术同步精进，Matlab项目合作可私信。
🍎个人主页：海神之光
🏆代码获取方式：
海神之光Matlab王者学习之路—代码获取方式
⛳️座右铭：行百里者，半于九十。
更多Matlab仿真内容点击👇
Matlab图像处理（进阶版）
路径规划（Matlab）
神经网络预测与分类（Matlab）
优化求解（Matlab）
语音处理（Matlab）
信号处理（Matlab）
车间调度（Matlab）
⛄一、手写大写字母识别技术简介
0 引言
在高校教学过程中，考试是最为普遍的一种教学评估、综合练习的教学手段，随着科技进步，考试阅卷的方式也发生了巨大的变革。传统的阅卷方式主要以人工阅卷为主， 存在效率低下等不足； 现代的阅卷方式采用了光学标记识别(Optica Mark Recognition， OMR) 技术， 考生只需在答题卡上填涂， 计算机会通过对答题卡进行处理从而实现自动阅卷， 但这种方式需使用特别设计的答题卡与铅笔，并且遵循一定的填涂规范。以上两种方法均给老师、考生带来了一定的限制，耗费了较多的时间；因此为了让被试者灵活、高效的进行答题，对手写体的字符识别进行研究十分有必要。
手写体的字符识别按照采取的技术路线大体上分为两类：一类是特征提取法；另一类是神经网络法。特征提取法按照提取字符特征的不同又可分为两类：一类是基于轮廓、结构特征的字符识别算法，这类算法通过识别字符图像的轮廓特征、端点特征、凹凸结构特征等，对字符进行细化处理，采用模板匹配的方式实现手写字符的自动识别。这类方式虽然能够对字符的结构进行直观地描述，但是对于噪声及产生形变的字符不能很好的识别，缺乏鲁棒性。
另一类是基于统计特征的手写体字符识别算法，这类算法通过对大量的样本训练相应的分类器，并利用这些分类器对待识别字符进行分类。优点是当采样样本足够充分时，这类方法能够具有较好的识别效果；缺点就是往往不能构造足够丰富的样本。神经网络法得到较高的识别率，但需要对大量样本进行训练，耗费较多时间，降低了识别的实时性；即使对神经网络进行改进、优化，也没有从根本上弥补其不足。
在考试中由于被试者的书写风格不一造成了手写字符的形态各异的特点，这就决定了在手写字符的识别中单一方案不会得到很好的识别效果。试卷客观题的评阅中，大多只包含A、B、C、D四个字符，字符个数少，仅对A~D四个字符进行识别能够得到较好的阅卷效率及较高的正确识别率。针对手写英文字母的特点及应用场景，本文提出一种基于组合特征的手写英文字母识别方法。该方法在轮廓特征提取的基础上加入形状特征提取，提取特征信息简单，同时不需要样本训练，因此提高了手写英文字母的识别成功率与识别速度。
1 图像预处理
由于手写字符形状存在风格不一的特点，且图像拍摄时的光照变化、纸张整洁程度等都会使得字符图像存在一定程度的噪声，因此必须对图像进行预处理，这样便于对后面的字母切分和特征提取，提高最终的字母识别准确率。图像预处理包括使用高斯滤波去除拍照设备产生的高斯噪声；对图像进行二值化，使得图像去掉背景；最后使用轮廓提取去除较大的噪点。
1.1 高斯滤波
数字图像在获取、传输的过程中都可能受到噪声的污染，而摄像机传感器内部元器件最容易产生的噪声是高斯噪声。高斯噪声是指它的概率密度函数服从高斯分布的一类噪声，对于此类噪声常用的方法为高斯滤波。高斯滤波是一种线性平滑滤波，适用于消除高斯噪声，广泛应用于图像处理的减噪过程。根据对测试样本分别使用1×1、3×3、5×5的高斯核进行实验发现，使用3×3高斯核的样本实验结果表现最好，因此本文选用3*3的高斯核对原始图像进行高斯滤波。
1.2 图像二值化
对图像进行二值化处理是利用手写字符区域与其背景存在灰度差异的特点，将手写字符与背景以灰度区分开，便于后期的处理。常用的二值化算法有最大类间方差法(OTSU法) 、迭代法、Bers en、Niblack、简单统计法、Niblack与模糊C均值法(Niblack and Fuzzy C-Means， NFC M) 等。OTSU法是一种自适应的阈值确定的方法， 其基本思想是求取最佳门限阈值，此阈值将图像灰度直方图分割成黑白两部分，使两部分类间方差取得最大值，并使类内方差值最小，因此OTSU法也称为最大类间方差法。本文图像中的手写英文字母与背景具有明显的灰度差异， 因此采用OTSU法较为合适。图1为图像二值化过程。
图1 图像二值化
1.3 轮廓去噪
经过高斯滤波后，可以去除或者减少部分由摄像头传感器引起的噪声，但是对于纸张本身存在的噪点处理效果达不到要求。在实际应用过程中，除字母以外的背景区域有可能会出现分布不均、不规则噪点，因此本文采用寻找连通域轮廓的方式，找到所有连通域的轮廓，并计算它们的轮廓面积，设定一个面积阈值，轮廓面积小于阈值的连通域可以判断为噪点，将其剔除。经过大量实验，本文设定将像素面积小于12的像素点去除。具体的轮廓提取方法在后文详细叙述。图2是轮廓去噪后的图像。
2 字符分割
图像的投影特征在图像处理中是重要技术之一，一般是通过计算图像在X轴或者Y轴的每个位置上的像素点(黑色或白色)总数量，并描绘出相应的投影曲线图来分析图像特征，该技术常用于图像分割、字符检测及提取等。公式(1)(2)分别是水平投影与垂直投影的计算公式。其中，h，w分别是二值化图像的高度和宽度；f(i，j)为图像第i行第j列元素的值1或0。
2.1 行分割
行分割就是对图像X轴方向进行水平投影，得到X轴方向上每个位置的白点总数，若某个位置上的白点总数为0，则说明该位置没有手写痕迹。因为人们的常用书写习惯是一行一行写，且行与行之间存在一定距离，所以通过行分割就能知道哪些位置是没有书写过的，这样就能大致将手写区域按行进行分割。
图3 行分割过程
2.2 列分割
列分割就是对图像Y轴方向进行垂直投影，得到Y轴方向上每个位置的白点总数，若某个位置上的白点总数为0，则说明该位置为字母间的间隙。通过列分割很容易地将每个字母提取出来进行后面的识别。
图4 列分割过程
2.3 最小图像分割
在进行过行分割与列分割后，基本将字母分割成单独的图像，但是此时的图像往往尺寸较大，字母在整个图像中所占比例较小，为了减少图像处理时间，方便后面的字母识别，需要对分割后的字母再次进行切割，使得字母尽量占满整个图片。图5中的黑色边框为图像边界。
图5 最小图像示意图
3 分类识别
对于字母A~D很容易发现它们的轮廓特征，字母A、D包含一个闭合曲线；字母B包含2个闭合曲线；字母C包含0个闭合曲线，因此可以通过提取字母的闭合曲线数量对字母进行粗分，然后再对A、D进行区分。提取轮廓特征后接着提取形状特征，将两种特征的识别结果进行融合即可得到最后的识别结果
3.1 轮廓特征提取
轮廓提取在图像处理中是常用方法之一，轮廓可以简单认为成将连续的点连在一起的曲线，具有相同的颜色或者灰度。通过对二值图像进行拓扑分析，对图像进行行扫描，为不同的轮廓赋予不同的整数值，从而确定外轮廓、内轮廓以及它们的层级关系。如图6有2个黑色轮廓，1a为第一个黑色轮廓的外边界，1b为第一个黑色轮廓的孔边界，2a为第二个黑色轮廓的外边界，2b为第二个黑色轮廓的孔边界。
图6 轮廓示意图
具体的操作是先对二值图像进行行扫描，判断像素点的像素值，用式(3)表示图像的像素值，i是图像的行数，j是图像的列数。
当行扫描到轮廓的外边界与孔边界时，扫描终止，式(4)为外边界终止条件，式(5)为孔边界终止条件。以扫描终止时的f(i，j)为起点，标记边界上的像素，给边界赋予一个数值(如1)，以后每次发现一个新边界即在数值上加1，在图6中有两条轮廓，第一条轮廓的标记值为1，第二条轮廓标记值为2，这样即完成了轮廓的检测。
本文中，通过对分割后的手写英文字母进行轮廓提取可以很容易地将A~D分为3类：轮廓数为0的有C；轮廓数为1的有
A、D；轮廓数为2的有B。
由于字母A、D具有相同的轮廓数，因此需要对两个字母进行区分。通过对字母A、D的轮廓特征进行分析可以发现，对于同样尺寸的字母A和字母D，A字母轮廓面积比D字母轮廓面积小；对于已是最小图象的字母A和字母D，经过对182份手写英文字母数据进行统计分析得到字母A轮廓面积占图像面积的比例为0.075194884，字母D轮廓面积占图像面积比例为0.321412412。由此可见，A字母轮廓面积占尺寸面积的比例小于字母D。因此在王识别中只需对待识别的字母轮廓面积占比分别与字母A的占比系数和字母D的占比系数求差值，相比较即可判别A、D。图7为字母A与字母D提取轮廓面积示意图。
图7 字母A与字母D轮廓面积提取示意图
3.2形状特征提取
通过对英文字母A~D的形状进行分析可以找到各个字母的形状特征，具体操作为：将预处理后得到的二值图像(如图8(a))进行填充，得到图8(b)，可以看出字母A的形状近似三角形；字母B的形状包含两个凸起部分和一个凹进部分；字母C同样包含两个凸起部分和一个凹进部分；字母D只有一个凸起部分。然后对图8(b)进行水平投影将字母每行的像素和显示在投影直方图上(图8©)，对各个字母的像素和数据进行分析，可以得到字母A的数据曲线呈现单调递增；字母B的数据曲线包含2个波峰一个波谷；字母C的数据曲线同样包含2个波峰一个波谷；字母D的数据曲线包含1个波峰。
图8 形状特征提取过程
得到各个字母的数据后, 对数据曲线进行差分运算, 曲线的峰值点满足一阶导数为0, 二阶导数为负;波谷点满足一阶导数为0, 二阶导数为正;单调性曲线满足一阶导数不为0, 由此可以将字母区分为A、B和C、D三类。算法流程如下:数据曲线可以看作是一维向量如下:
其中, vi (i∈[1, 2, …, N]) 代表图像在第i行上的白色像素和。计算V的一阶差分向量Diffv:
得到差分向量后进行符号函数运算, 如式 (8) :
此时可以得到各点的变化情况，但是有些点值相同，却不是波峰波谷，那么再进行一次一阶差分即可得到数据曲线的波峰或波谷点。
3.3特殊处理
由于在考试中难免出现答案填写错误的情况，因此对于这种情况需要作特殊处理。本文在对字母识别的基础上设计一种错误答案标注方式，对于填写错误的答案，可在答案上从左至右划2笔不重复的斜杠，从右至左划2笔不重复的斜杠，如图9所示，经过标注后的错误答案轮廓数会超出正常字母的轮廓数，这样即可区分出填写错误的答案。
图9 特殊情况处理
⛄二、部分源代码
function varargout = untitled(varargin)
% UNTITLED MATLAB code for untitled.fig
% UNTITLED, by itself, creates a new UNTITLED or raises the existing
% singleton*.
%
% H = UNTITLED returns the handle to a new UNTITLED or the handle to
% the existing singleton*.
%
% UNTITLED(‘CALLBACK’,hObject,eventData,handles,…) calls the local
% function named CALLBACK in UNTITLED.M with the given input arguments.
%
% UNTITLED(‘Property’,‘Value’,…) creates a new UNTITLED or raises the
% existing singleton*. Starting from the left, property value pairs are
% applied to the GUI before untitled_OpeningFcn gets called. An
% unrecognized property name or invalid value makes property application
% stop. All inputs are passed to untitled_OpeningFcn via varargin.
%
% *See GUI Options on GUIDE’s Tools menu. Choose “GUI allows only one
% instance to run (singleton)”.
%
% See also: GUIDE, GUIDATA, GUIHANDLES
% Edit the above text to modify the response to help untitled
% Last Modified by GUIDE v2.5 05-Jan-2022 19:52:04
% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct(‘gui_Name’, mfilename, …
‘gui_Singleton’, gui_Singleton, …
‘gui_OpeningFcn’, @untitled_OpeningFcn, …
‘gui_OutputFcn’, @untitled_OutputFcn, …
‘gui_LayoutFcn’, [] , …
‘gui_Callback’, []);
if nargin && ischar(varargin{1})
gui_State.gui_Callback = str2func(varargin{1});
end
if nargout
[varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
else
gui_mainfcn(gui_State, varargin{:});
end
% End initialization code - DO NOT EDIT
% — Executes just before untitled is made visible.
function untitled_OpeningFcn(hObject, eventdata, handles, varargin)
% This function has no output args, see OutputFcn.
% hObject handle to figure
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% varargin command line arguments to untitled (see VARARGIN)
% Choose default command line output for untitled
handles.output = hObject;
% Update handles structure
guidata(hObject, handles);
% UIWAIT makes untitled wait for user response (see UIRESUME)
% uiwait(handles.figure1);
% — Outputs from this function are returned to the command line.
function varargout = untitled_OutputFcn(hObject, eventdata, handles)
% varargout cell array for returning output args (see VARARGOUT);
% hObject handle to figure
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Get default command line output from handles structure
varargout{1} = handles.output;
% — Executes on button press in pushbutton1.
function pushbutton1_Callback(hObject, eventdata, handles)
% hObject handle to pushbutton1 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
[filename, pathname] = uigetfile({‘.jpg;.tif;.png;.gif;.bmp’,‘All Image Files’;…
'.*’,‘All Files’ },‘载入图像’,…
fullfile(pwd, ‘images’));
if isequal(filename, 0) || isequal(pathname, 0)
return;
end
% [FileName,Pathname]=uigetfile(‘.jpg;.jpeg;.bmp;.png;*.gif’);
F_P=fullfile(pathname,filename);
img=imread(F_P);
handles.path=F_P;
handles.img=img;
axes(handles.axes1);imshow(img);
guidata(hObject,handles);
% — Executes on button press in pushbutton2.
function pushbutton2_Callback(hObject, eventdata, handles)
% hObject handle to pushbutton2 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
[out1,out2]=main_(handles.img);
axes(handles.axes3);imshow(out2);
set(handles.edit1,‘string’,out1);
function edit1_Callback(hObject, eventdata, handles)
% hObject handle to edit1 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hints: get(hObject,‘String’) returns contents of edit1 as text
% str2double(get(hObject,‘String’)) returns contents of edit1 as a double
% — Executes during object creation, after setting all properties.
function edit1_CreateFcn(hObject, eventdata, handles)
% hObject handle to edit1 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called
% Hint: edit controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispc && isequal(get(hObject,‘BackgroundColor’), get(0,‘defaultUicontrolBackgroundColor’))
set(hObject,‘BackgroundColor’,‘white’);
end
⛄三、运行结果
⛄四、matlab版本及参考文献
1 matlab版本
2014a
2 参考文献
[1]和小娟,罗勇.边界链码在字母与数字混合识别中的应用[J].郑州大学学报(理学版). 2011,43(03)
3 备注
简介此部分摘自互联网，仅供参考，若侵权，联系删除
🍅 仿真咨询
1 各类智能优化算法改进及应用
生产调度、经济调度、装配线调度、充电优化、车间调度、发车优化、水库调度、三维装箱、物流选址、货位优化、公交排班优化、充电桩布局优化、车间布局优化、集装箱船配载优化、水泵组合优化、解医疗资源分配优化、设施布局优化、可视域基站和无人机选址优化
2 机器学习和深度学习方面
卷积神经网络（CNN）、LSTM、支持向量机（SVM）、最小二乘支持向量机（LSSVM）、极限学习机（ELM）、核极限学习机（KELM）、BP、RBF、宽度学习、DBN、RF、RBF、DELM、XGBOOST、TCN实现风电预测、光伏预测、电池寿命预测、辐射源识别、交通流预测、负荷预测、股价预测、PM2.5浓度预测、电池健康状态预测、水体光学参数反演、NLOS信号识别、地铁停车精准预测、变压器故障诊断
3 图像处理方面
图像识别、图像分割、图像检测、图像隐藏、图像配准、图像拼接、图像融合、图像增强、图像压缩感知
4 路径规划方面
旅行商问题（TSP）、车辆路径问题（VRP、MVRP、CVRP、VRPTW等）、无人机三维路径规划、无人机协同、无人机编队、机器人路径规划、栅格地图路径规划、多式联运运输问题、车辆协同无人机路径规划、天线线性阵列分布优化、车间布局优化
5 无人机应用方面
无人机路径规划、无人机控制、无人机编队、无人机协同、无人机任务分配
6 无线传感器定位及布局方面
传感器部署优化、通信协议优化、路由优化、目标定位优化、Dv-Hop定位优化、Leach协议优化、WSN覆盖优化、组播优化、RSSI定位优化
7 信号处理方面
信号识别、信号加密、信号去噪、信号增强、雷达信号处理、信号水印嵌入提取、肌电信号、脑电信号、信号配时优化
8 电力系统方面
微电网优化、无功优化、配电网重构、储能配置
9 元胞自动机方面
交通流 人群疏散 病毒扩散 晶体生长
10 雷达方面
卡尔曼滤波跟踪、航迹关联、航迹融合
2419
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
