模式识别-第八章-实例教学ppt课件.ppt

《模式识别-第八章-实例教学ppt课件.ppt》由会员分享，可在线阅读，更多相关《模式识别-第八章-实例教学ppt课件.ppt（58页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第八章第八章 实例教学实例教学8.1 应用统计决策的肤色提取应用统计决策的肤色提取 8.1.1图像基础知识 nHSI表色模型表色模型n面向彩色处理的最常用的模型，符合人类的视觉面向彩色处理的最常用的模型，符合人类的视觉感受。感受。H（Hue）色调（色相），表明物体是什么颜色。色调（色相），表明物体是什么颜色。S（Saturation）饱和度，表明该种颜色的浓淡。饱和度，表明该种颜色的浓淡。I（Intensity）亮度，与图像的彩色信息无关。亮度，与图像的彩色信息无关。nRGB模型与模型与HSI模型可以相互转换。模型可以相互转换。8.1.2样本获取 原始图像 H H分量分量 S S分量分量 I

2、I分量分量由于由于H H分量较好的区分了肤色与分量较好的区分了肤色与背景，因此选背景，因此选H H分量图像中的像分量图像中的像素点作为样本。素点作为样本。8.1.3样本学习样本学习 H分量 H分量的二值化图皮肤点数的先验概率：皮肤点数的先验概率：P P（ 1 1）= =是皮肤的像素点是皮肤的像素点数数 / / 所有像素点数所有像素点数 = 7.99%= 7.99%非皮肤点数的先验概率：非皮肤点数的先验概率：P P（ 2 2）= =不是皮肤的像不是皮肤的像素点数素点数 / / 所有像素点数所有像素点数= 92.01%= 92.01%（1 1）计算先验概率：）计算先验概率：图像二值化，为图像二值化

3、，为1 1的是的是皮肤，为皮肤，为0 0的是背景。的是背景。（2 2）计算类条件概率密度：）计算类条件概率密度：计算计算H H分量图像中皮肤区域的直方图，就是皮肤的类分量图像中皮肤区域的直方图，就是皮肤的类条件概率密度；计算条件概率密度；计算H H分量图像中非皮肤区域的直方分量图像中非皮肤区域的直方图，就是非皮肤的类条件概率密度。图，就是非皮肤的类条件概率密度。 非皮肤的类条件概率密度非皮肤的类条件概率密度 皮肤的类条件概率密度皮肤的类条件概率密度 横坐标是像素点可能的灰度值（将灰度量化成了横坐标是像素点可能的灰度值（将灰度量化成了8等分），纵坐标是相应区域中灰度值等分），纵坐标是相应区域中灰

4、度值为横坐标值的像素点数。为横坐标值的像素点数。 P(x| 1): 0.9855 0.0143 0 0 0 0 0 0.0002 P(x| 2): 0.0790 0.7820 0.0742 0.0192 0.0199 0.0074 0.0072 0.01088.1.4模式分类模式分类 输入下图，用最小误判概率准则和最小损失准输入下图，用最小误判概率准则和最小损失准则判决判断每个像素点是否是皮肤。则判决判断每个像素点是否是皮肤。 待分类的图像 H分量 a)最小误判概率准则最小误判概率准则判断每个像素点是否使皮肤。判断每个像素点是否使皮肤。两类问题的两类问题的Bayes最小误判概率准则为最小误判概

5、率准则为如果如果 ，则判，否则判，则判，否则判，x为为每个像素点的灰度值，若该点属于皮肤，置位每个像素点的灰度值，若该点属于皮肤，置位1；若不；若不是皮肤，置位是皮肤，置位0，得到下图。，得到下图。 最小误判概率准则的识别结果12121221( |)()( )( |)()p xPlxp xP1x2xb)b)最小损失准则最小损失准则 1111= 0,= 0, 1212=10=10, , 2121=1, =1, 2222=0=0 如果，则判，否则如果，则判，否则判。判。x x为每个像素点的灰度值，若该点为每个像素点的灰度值，若该点属于皮肤，置位属于皮肤，置位1 1；若不是皮肤，置位；若不是皮肤，置

6、位0 0，得到下，得到下图。图。 可以看出最小损可以看出最小损失准则得到的皮肤更失准则得到的皮肤更完整，漏检率下降，完整，漏检率下降，但是将更多非皮肤点但是将更多非皮肤点误判为皮肤，虚警率误判为皮肤，虚警率增加。增加。最小损失准则的识别结果 1221221212211211( |)()()( )( |)( )()p xPlxp xP 1x2x8.2 车牌识别车牌识别 8.2.1车牌相关知识 英文字母中的英文字母中的I I和和O O一般避而不用，以免和数字中的一般避而不用，以免和数字中的1 1和和0 0混淆。混淆。 车牌大小为车牌大小为440440* *140mm140mm。汉字字体为黑体，数字

7、和字母字体为国家公安部门开汉字字体为黑体，数字和字母字体为国家公安部门开发的特殊字体。发的特殊字体。8.2.2 图像采集 8.2.3车牌识别预处理(1)字符分割 (2)字符细化 1.特征选择与提取（1 1）特征空间向量的选择）特征空间向量的选择 由于车牌中字母和数字都有特定的字体，且作为汽车的由于车牌中字母和数字都有特定的字体，且作为汽车的“身份证身份证”，车牌制作过程严谨，相同字符形态区别不大，他们的差异性和共同点都很车牌制作过程严谨，相同字符形态区别不大，他们的差异性和共同点都很明显。基于上节描述的设计准则，综合考虑了提取特征的编程复杂度和特明显。基于上节描述的设计准则，综合考虑了提取特征

8、的编程复杂度和特征对不同字符的区分度，本系统分别从字符的形态、结构、笔画特点出发征对不同字符的区分度，本系统分别从字符的形态、结构、笔画特点出发，采用了，采用了4 4个类别共个类别共1010个特征值作为判断的矢量依据（特征提取时均以细个特征值作为判断的矢量依据（特征提取时均以细化后图像为基础）。化后图像为基础）。字符的重心判断：字符的重心判断： 字符重心位于上方字符重心位于上方/ /下方，分别计算字符上下半区的像素数目，如果下方，分别计算字符上下半区的像素数目，如果其差值在一个固定的小范围内其差值在一个固定的小范围内( (试验得出此差值在试验得出此差值在0-200-20像素内为最佳像素内为最佳

9、) )，则该特征值记为则该特征值记为2 2，上方偏多记为，上方偏多记为1 1，下方偏多记为，下方偏多记为3 3； 字符重心位于左方字符重心位于左方/ /右方，分别计算字符左右半区的像素数目，如果右方，分别计算字符左右半区的像素数目，如果其差值在一个固定的小范围内其差值在一个固定的小范围内( (试验得出此差值在试验得出此差值在0-200-20像素内为最佳像素内为最佳) )，则该特征值记为则该特征值记为2 2，左方偏多记为，左方偏多记为1 1，右方偏多记为，右方偏多记为3 3。 特定位置的扫描线特征：特定位置的扫描线特征： 图像垂直图像垂直1/31/3处扫描线穿越的线条数目；处扫描线穿越的线条数目

10、； 图像垂直图像垂直1/21/2处扫描线穿越的线条数目；处扫描线穿越的线条数目； 图像垂直图像垂直2/32/3处扫描线穿越的线条数目；处扫描线穿越的线条数目； 图像水平图像水平1/31/3处扫描线穿越的线条数目；处扫描线穿越的线条数目； 图像水平图像水平2/32/3处扫描线穿越的线条数目；处扫描线穿越的线条数目； 为了消除可能的细化算法不完善的影响，以上在计算穿越次数时均为为了消除可能的细化算法不完善的影响，以上在计算穿越次数时均为0-10-1改变次数。改变次数。 笔画特征：笔画特征： 字符中字符中“竖竖”的数目，计算图像在水平方向上的投影，投影点的像素数目的数目，计算图像在水平方向上的投影，

11、投影点的像素数目累加值超过累加值超过1/31/3图像高度且无超过图像高度且无超过5 5像素的中断，记为像素的中断，记为“竖竖”的数目加的数目加1 1； 字符中字符中“横横”的数目，计算图像在垂直方向上的投影，投影点的像素数目的数目，计算图像在垂直方向上的投影，投影点的像素数目累加值超过累加值超过1/31/3图像宽度且无超过图像宽度且无超过5 5像素的中断，记为像素的中断，记为“横横”的数目加的数目加1 1； 在实际情况中，由于车牌图像并不总是标准的，当出现倾斜、几何失真在实际情况中，由于车牌图像并不总是标准的，当出现倾斜、几何失真等情况时，字符中的横或竖不可能全部投影在一点上，而是在一个小范围

12、内等情况时，字符中的横或竖不可能全部投影在一点上，而是在一个小范围内密集分布。经过多次尝试，结合车牌字符本身特点，本文将判定存在竖密集分布。经过多次尝试，结合车牌字符本身特点，本文将判定存在竖/ /横笔横笔画的门限值设为图像高度画的门限值设为图像高度/ /宽度的宽度的1/31/3左右。左右。 为了消除可能的图像几何失真带来的影响，计算笔画数目时均不重复计为了消除可能的图像几何失真带来的影响，计算笔画数目时均不重复计算算2020像素点以内的笔画，既若该点存在横像素点以内的笔画，既若该点存在横/ /竖，则在该投影方向上向下竖，则在该投影方向上向下/ /后移后移动动2020像素继续判断。像素继续判断

13、。 结构特征：结构特征： 字符中存在闭合回路的数目。利用递归实现区域生长算法，对所有独立连字符中存在闭合回路的数目。利用递归实现区域生长算法，对所有独立连通区域进行标号处理，最后将连通区域个数减一通区域进行标号处理，最后将连通区域个数减一( (外围背景区域也作为一个连外围背景区域也作为一个连通区域被标号，需要在特征值中减去通区域被标号，需要在特征值中减去) )记为该特征值。记为该特征值。（2 2）模板数据的生成）模板数据的生成考虑到车牌字符的特殊性，无法使用Windows系统自带的普通字体来作为标准模板。我们对拍摄到的车牌图像进行如下处理，形成了标准的模版数据。为了杜绝可能的系统图像处理缺陷，

14、使用PhotoShop等商用软件将其调整为二值化图像；为了减少系统缩放函数可能存在的不完善对画面失真的影响，将不同字符做成140*440的标准图像，供系统提取特征值；使用程序中的细化模块细化图像，并以特征提取模块计算标准图像的特征值，将得出的数据结果作为模板数据，以矩阵的形式保存在template.ini文件中。 2.2.字符识别字符识别采用最小距离准则进行模板匹配的识别。判决的标准采用欧式距离最小准则，计算当前图形的特征向量与34个特征向量(除字母I和O以外的所有数字和字母，共34个)的欧式距离，取其中最小的为判决结果。于是待识别样品与模板k之间的欧式距离dk为： 10122|iiikyxd

15、33010122|minkiiiyxd其中距离最短者为：对每一个样品的判决过程为：遍历所有模板，计算其模板与样品间的欧式距离，取最小值作为最终判决结果。如果最短距离小于某个规定的阈值，则结果为这个最小值相对应的类别；否则判定为无法识别。识别结果如图所示 在本系统中，使用一个识别模块来识别当前样品，其算法流程图如下图：8.3 纸币识别器纸币识别器 、数据采集数据采集、特征提取特征提取、训练学习训练学习、分类识别分类识别8.3.1数据采集 纸币面额：5元，10元，20元，50元，100元8.3.2特征提取与选择 特征提取：长度、宽度、磁性、磁性特征提取：长度、宽度、磁性、磁性的位置，光反射亮度、光

16、透射亮度等的位置，光反射亮度、光透射亮度等8.3.38.3.3训练学习训练学习训练集：是一个已知样本集，在监督学习方法中，训练集：是一个已知样本集，在监督学习方法中，用它来开发出模式分类器。图中给出了对一个样本用它来开发出模式分类器。图中给出了对一个样本进行特征提取、分类的学习过程进行特征提取、分类的学习过程8.3.4分类识别分类识别对测试集来进行分类识别，确定纸币的面额及真伪对测试集来进行分类识别，确定纸币的面额及真伪。为了更好地对模式识别系统性能进行评价，必须。为了更好地对模式识别系统性能进行评价，必须使用一组独立于训练集的测试集对系统进行测试。使用一组独立于训练集的测试集对系统进行测试。

17、图中是对图中是对1010元纸币分类识别的界面。元纸币分类识别的界面。8.4最近邻方法用于图像拼接最近邻方法用于图像拼接8.4.18.4.1图像获取图像获取获取有重叠区域的两幅图像，如下图所示： （a a） (b)(b)8.4.2 特征提取特征提取n分别在两副图中用分别在两副图中用SIFT算法提取特征点，算法提取特征点，(a)图中的特图中的特征点数为征点数为5793，(b)图中的特征点数为图中的特征点数为3199。nScale Invariant Feature Transform ，尺度不变量特，尺度不变量特征提取算法，简称征提取算法，简称 SIFT 算法，是一种提取控制点特征算法，是一种提取

18、控制点特征的经典算法（以下所说的的经典算法（以下所说的“关键点关键点”（Keypoint）是）是文献文献1中提到的说法，中提到的说法， 可以理解为控制点）。可以理解为控制点）。nSIFT算法提取的算法提取的SIFT特征向量具体具有如下特性：特征向量具体具有如下特性：SIFT特征是图像的局部特征，其对旋转、尺度缩放、特征是图像的局部特征，其对旋转、尺度缩放、亮度变化保持不变性，对视角变化、仿射变换、噪声也亮度变化保持不变性，对视角变化、仿射变换、噪声也保持一定程度的稳定；独特性保持一定程度的稳定；独特性(Distinctiveness)好，信好，信息量丰富，适用于在海量特征数据库中进行快速、准确

19、息量丰富，适用于在海量特征数据库中进行快速、准确的匹配；多量性，即使少数的几个物体也可以产生大量的匹配；多量性，即使少数的几个物体也可以产生大量SIFT特征向量；高速性，经优化的特征向量；高速性，经优化的SIFT匹配算法甚至匹配算法甚至可以达到实时的要求；可扩展性，可以很方便的与其他可以达到实时的要求；可扩展性，可以很方便的与其他形式的特征向量进行联合。形式的特征向量进行联合。nSIFT SIFT 算法首先在尺度空间进行特征检测，并确定关算法首先在尺度空间进行特征检测，并确定关键点的位置和关键点所处的尺度，然后使用关键点键点的位置和关键点所处的尺度，然后使用关键点邻域梯度的主方向作为该点的方向

20、特征，以实现算邻域梯度的主方向作为该点的方向特征，以实现算子对尺度和方向的无关性。主要计算步骤如下：子对尺度和方向的无关性。主要计算步骤如下： n尺度空间极值检测：搜索整个尺度和图像位置。通尺度空间极值检测：搜索整个尺度和图像位置。通过使用高斯差分函数（过使用高斯差分函数（Difference of GaussianDifference of Gaussian）确定对尺度和方向具有不变性的兴趣点。以初步确确定对尺度和方向具有不变性的兴趣点。以初步确定关键点位置和所在尺度。定关键点位置和所在尺度。n关键点定位：在每一个候选位置，详细地模型拟合关键点定位：在每一个候选位置，详细地模型拟合以确定位置

21、和尺度。根据稳定性选择关键点。通过以确定位置和尺度。根据稳定性选择关键点。通过拟合三维二次函数以精确确定关键点的位置和尺度拟合三维二次函数以精确确定关键点的位置和尺度，同时消除低对比度的关键点和不稳定的边缘响应，同时消除低对比度的关键点和不稳定的边缘响应点点( (因为高斯差分算子会产生较强的边缘响应因为高斯差分算子会产生较强的边缘响应) )，以，以增强匹配稳定性、提高抗噪声能力。增强匹配稳定性、提高抗噪声能力。n方向分配：每一个关键点根据局部图像梯度方向分配一个或多个方向分配：每一个关键点根据局部图像梯度方向分配一个或多个方向。其后所有操作都根据分配的方向、尺度、位置。利用关键方向。其后所有操

22、作都根据分配的方向、尺度、位置。利用关键点邻域像素的梯度方向分布特性为每个关键点指定方向参数，使点邻域像素的梯度方向分布特性为每个关键点指定方向参数，使算子具备旋转不变性。算子具备旋转不变性。n关键点描述符：在每个关键点周围在选定的尺度下测量局部图像关键点描述符：在每个关键点周围在选定的尺度下测量局部图像梯度。为了增强匹配的稳健性，对每个关键点使用梯度。为了增强匹配的稳健性，对每个关键点使用44共共16个个种子点来描述，这样对于一个关键点就可以产生种子点来描述，这样对于一个关键点就可以产生128个数据，最个数据，最终形成终形成128维的维的SIFT特征向量。特征向量。n其中第一步骤里高斯差分函

23、数是这样定义的：其中第一步骤里高斯差分函数是这样定义的：n高斯卷积核是实现尺度变换的唯一线性核，于是一副二维图像的高斯卷积核是实现尺度变换的唯一线性核，于是一副二维图像的尺度空间定义为：尺度空间定义为：n n其中其中 是尺度可变高斯函数，是尺度可变高斯函数，n n（x，y）是空间坐标，是尺度坐标。）是空间坐标，是尺度坐标。),(),(),(yxIyxGyxL2)(22/21),(22yxeyxGn为了有效的在尺度空间检测到稳定的关键点，提出了为了有效的在尺度空间检测到稳定的关键点，提出了高斯差分尺度空间（高斯差分尺度空间（DOG scale-space）。利用不同）。利用不同尺度的高斯差分核与

24、图像卷积生成。尺度的高斯差分核与图像卷积生成。n以上就是高斯差分函数计算式。高斯差分函数（以上就是高斯差分函数计算式。高斯差分函数（Difference of Gaussian）计算比较简单，是尺度归）计算比较简单，是尺度归一化的拉普拉斯一化的拉普拉斯-高斯算子的近似。图像经过高斯差分高斯算子的近似。图像经过高斯差分算子计算，可以像用算子计算，可以像用Sobel 算子计算时一样，提取局算子计算时一样，提取局部极值点作为控制点。部极值点作为控制点。n下面用下图形象地描述下面用下图形象地描述SIFT 算法中第三步和第四步的算法中第三步和第四步的过程。过程。),(),(),(),(),(),(yxL

25、kyxLyxIyxGkyxGyxDnSIFT算法形象表示图算法形象表示图 如图所示，在控制点（左图中心点）周围如图所示，在控制点（左图中心点）周围8 88 8邻域内邻域内，每个点的梯度大小和方向均为已知，中央点为当前控，每个点的梯度大小和方向均为已知，中央点为当前控制点的位置，每个小方格代表控制点邻域所在尺度空间制点的位置，每个小方格代表控制点邻域所在尺度空间的一个像素，箭头方向代表该像素的梯度方向，箭头长的一个像素，箭头方向代表该像素的梯度方向，箭头长度代表梯度模值（注：实际应用中，中心点也是像素点度代表梯度模值（注：实际应用中，中心点也是像素点，因此可将左图中靠近中心的四个像素点选出一个作

26、为，因此可将左图中靠近中心的四个像素点选出一个作为中心点）。然后根据像素点和中心点的远近距离来进行中心点）。然后根据像素点和中心点的远近距离来进行高斯函数加权，左图圆圈表示高斯函数加权的范围。高斯函数加权，左图圆圈表示高斯函数加权的范围。2222212xye 3230,1N 加权方法如下：对圆圈内的每个像素点对应的梯度，令其方向不变，使其梯度模值乘以加权系数。n其中用来确定高斯函数加权的范围，是因为 所以在上图中，可以近似，即。其中和是所考察像素点与中心点X方向和Y方向的坐标差值。n加权完毕后，将88的邻域分块进行梯度直方图统计。梯度直方图统计方法如下：如左图所示粗线将邻域分成4块，每一块里都

27、含有16个像素点。接下来给定梯度8个方向，如右图所示，将左图每一块里的各个像素点所对应的梯度，归类到8个给定的方向。再将归为一类的梯度模值求和，就得到如右图所示的8个方向长短不一的梯度直方图了。n如此进行得到控制点周围4个梯度直方图，将这4个直方图化为向量，再将其首尾相连，就成为控制点的32位SIFT特征向量。n以上就是SIFT 算法的大致思路。这种邻域方向性信息联合的思想增强了算法抗噪声的能力，同时对于含有定位误差的特征匹配也提供了较好的容错性。8.4.3 8.4.3 特征匹配特征匹配 对图对图(a)(a)中的中的57935793个点，依次用最近邻方法个点，依次用最近邻方法K-NNK-NN算

28、法，在图算法，在图(b)(b)中的中的31993199个点中寻找与个点中寻找与(a)(a)每每个点的最佳配对点，即属于同一类的点，没有匹个点的最佳配对点，即属于同一类的点，没有匹配上的点删除。实验结果为找到了配上的点删除。实验结果为找到了5656对匹配点，对匹配点，如下图所示。如下图所示。8.4.4 图像拼接图像拼接通过56对匹配点，将两幅图像拼接成下图。8.5从红外图像序列中识别目标从红外图像序列中识别目标n、跑步人识别跑步人识别n、船只识别船只识别n、车辆识别车辆识别8.5.8.5.1 1 跑步人识别跑步人识别红外图像序列是由红外图像序列是由ECCO ENGNEERING CO.ECCO

29、ENGNEERING CO.生产的生产的FLIRFLIR摄像机采集的。监视状态下视场角为摄像机采集的。监视状态下视场角为2013.32013.3度，分度，分辨率辨率800800线线 450 450线，帧频线，帧频30 Hz30 Hz。检测算法中的运。检测算法中的运动强度系数取动强度系数取4 4，位移系数取，位移系数取1010。在视场中能否检测。在视场中能否检测出运动目标，与目标的距离、运动速度和目标面积出运动目标，与目标的距离、运动速度和目标面积密切相关。在一定的距离上能检测的目标的运动速密切相关。在一定的距离上能检测的目标的运动速度和大小也有一个确定的范围。测试结果如下：度和大小也有一个确定

30、的范围。测试结果如下：8.5.2 8.5.2 船只识别船只识别 船只在海面上行驶时，海水、船只和天空在红外图像中特性各不一样，海水温度较低，灰度值较低，天空次之，灰度较高，船只温度最高，灰度最亮，因此在理想情况下，只要从红外图像中检测出最亮的部分即可识别船只。然而，在实际图像中，海水、天空和船只在红外图像中的灰度并不是均一的，而是受强烈的噪声干扰，噪声的方差已达到几十个灰度级，因此简单的门限分割是不可能解决问题，整个问题的核心是如何从强噪声环境中检测出分割门限。321,321,www vp 3122221imviiiiewvp2550 v下面我们来分析问题的数学模型：我们把噪声视为高斯分布，设

31、海水、天空和船只的平均灰度分别为，对应方差为它们在整个图像所占面积的比例为则图像灰度分布的概率密度函数 若根据输入图像的灰度直方图得到灰度分布概率密度函数，可以依据贝叶斯准则得到海水、船只和天空之间的分割门限。红外图像船只识别系统界面如下图所示：321,mmm8.5.3 车辆识别车辆识别n由红外传感器所获取的运动车辆的图像序列，具由红外传感器所获取的运动车辆的图像序列，具有强噪声、目标边缘模糊以及同一目标由于温度有强噪声、目标边缘模糊以及同一目标由于温度分布不均匀造成目标内部灰度变化剧烈等特点，分布不均匀造成目标内部灰度变化剧烈等特点，这些特点使得许多常规图像分析和处理方法失效这些特点使得许多

32、常规图像分析和处理方法失效或效果不佳。另一方面，由于传感器与目标之间或效果不佳。另一方面，由于传感器与目标之间的距离较近，运动目标在序列的各幅图像中的姿的距离较近，运动目标在序列的各幅图像中的姿态变化比较剧烈，这也对目标的检测识别和跟踪态变化比较剧烈，这也对目标的检测识别和跟踪工作提出了挑战。我们从运动目标的检测和跟踪工作提出了挑战。我们从运动目标的检测和跟踪两方面着手，解决了在红外图像序列中车辆的捕两方面着手，解决了在红外图像序列中车辆的捕获及其姿态判定问题。获及其姿态判定问题。n上面的上面的4幅图像是我们从序列中抽取出来的反映车辆在幅图像是我们从序列中抽取出来的反映车辆在行驶过程中不同的姿

33、态：行驶过程中不同的姿态：n 本系统能够从序列的任意时刻开始进行运动目标的检本系统能够从序列的任意时刻开始进行运动目标的检测和跟踪，准确地报告目标的位置和姿态，并且在传感测和跟踪，准确地报告目标的位置和姿态，并且在传感器指向发生切换时，系统能够及时停止跟踪，自动完成器指向发生切换时，系统能够及时停止跟踪，自动完成目标的重新检测识别，然后继续进行跟踪，相应的测试目标的重新检测识别，然后继续进行跟踪，相应的测试结果如下：结果如下： 8.6心电图自动诊断心电图自动诊断、心电自动诊断介绍、心电自动诊断介绍、预处理、预处理、特征提取与选择、特征提取与选择、心脏病诊断、心脏病诊断8.6.1 心电自动诊断介

34、绍 心脏病是威胁人类生命的主要疾病之一，长期以来，对心脏病的研究一直是医学界的重要课题。人体心电图(Electrocardiograph，ECG，如图1所示) 作为心脏电活动在体表的综合表现，蕴涵着丰富的反映心脏节律及其电传导的生理和病理信息，在一定程度上可以客观反映心脏各部位的生理状况，是诊断心脏疾病、评价心脏功能的重要依据之一，尤其对于各种心律失常和传导障碍的诊断分析具有重大诊断价值，是目前分析与鉴别各种心律失常最精确的方法。图1 正常的典型心电波形图n由于病态心电图种类繁多、变异极大，同种病理不同患者的心电图、甚至同一患者不同时刻的心电图都存在着较大的差异；若医师长期从事大量图形的识别工

35、作，极易疲劳，容易漏检、出错。而且，人工进行心电信号分析都是事后分析，分析速度达不到准实时性要求。为了把医师从繁琐的图形识别工作中解脱出来，将精力集中在异常心电波形的识别和判断上，从而提高医疗诊断效率、缩短诊断时间，研究心电图的计算机辅助分析算法，开发出一种功能强大的心电信号自动分析系统，特别是可以准确、快速地对心律失常病类进行判别的系统，近年来越来越受到医务人员和科技人员的重视。n心电信号自动分析是指在采集到ECG信号的基础上，通过对其处理提取表征心脏状态的波形信息和特征参数，获取心脏工作状态的相关信息，然后利用这些特征信息分析、判别心电信号类型及所对应的疾病类型或健康水平，进而对心脏状态和

36、健康状况进行预测。整个心电信号自动分析的流程如图2所示。ECG预处理特征提取与选择波形分类疾病诊断结果图2 心电信号自动分析流程图n由图2可知，心电信号自动分析研究的内容包括以下四个主要方面：一是心电 信号预处理技术的研究；二是波形检测与特征点定位算法的研究；三是心电 特征提取与选择方法的研究；四是心电波形自动分类与疾病诊断方法的研究。 8.6.2 预处理 n临床采集的心电信号，尤其是动态心电信号，由于受试者所在环境、运动量和活动的剧烈程度不一样，系统所采集的ECG质量也不一样。心电信号是mV级的微弱信号，在采集过程中，受人体运动、采集仪器、电磁环境等影响，不可避免地会引入许多噪声干扰，这些噪

37、声主要有以下几种表现形式：工频干扰，基线漂移，肌电干扰，电极接触噪声，运动伪迹等。来源于心脏以外的这些噪声信号以加性或乘性方式与心电信号相混杂，会引起ECG信号的各种畸变，给准确检测和判断ECG信号带来极大的困难。在以上各种噪声干扰中，工频干扰和基线漂移是最为普遍的，也是心电信号预处理中必须首要考虑对其消除或有效抑制的主要噪声成份。小波滤波结果：小波滤波结果：图3 ECG信号的小波阈值滤波结果(a)含高频噪声的心电信号及其小波滤波结果(b)含基线漂移的心电信号及其小波滤波结果形态学滤波结果：形态学滤波结果： 图4 数学形态学滤波器对ECG信号的滤波结果(a)含高频噪声的ECG信号及其滤波结果

38、(b)含基线漂移干扰的ECG信号及其滤波结果8.6.3 特征提取与选择 图5 特征检测界面1RR2RRwQRShRppQRS均值方差均值方差均值方差均值方差均值方差N-类型0.75590.05350.75910.06020.04050.00991.24850.15941.40260.1690R-类型0.70520.03790.70450.04280.03920.00771.26380.11032.27150.1719A-类型0.58720.22070.53570.17820.03500.00740.92700.72401.66370.4228E-类型1.00250.14570.99180.15

39、080.02970.00940.80580.15611.14890.1655表1 时域特征参数的均值和方差统计结果图6 不同类型心电信号在小波系数分布上的差异 妈妈新开了个淘宝店，欢迎前来捧场妈妈新开了个淘宝店，欢迎前来捧场 妈妈的淘宝点开了快半年了，主要卖的是毛绒玩具、坐垫、抱枕之类的，妈妈的淘宝点开了快半年了，主要卖的是毛绒玩具、坐垫、抱枕之类的，但生意一直不是很好，感觉妈妈还是很用心的，花了不少功夫，但是就是没但生意一直不是很好，感觉妈妈还是很用心的，花了不少功夫，但是就是没有人气，所以我也来出自己的一份力，帮忙宣传一下。有人气，所以我也来出自己的一份力，帮忙宣传一下。 并且妈妈总是去五

40、亭龙挑最好的玩具整理、发货，质量绝对有保证。并且妈妈总是去五亭龙挑最好的玩具整理、发货，质量绝对有保证。 另外我家就在扬州五亭龙玩具城旁边，货源丰富，质量可靠，价格便宜。另外我家就在扬州五亭龙玩具城旁边，货源丰富，质量可靠，价格便宜。 欢迎大家来逛逛欢迎大家来逛逛【扬州五亭龙玩具总动员扬州五亭龙玩具总动员】 个人小广告：个人小广告：原始QRS波形二阶累积量三阶累积量四阶累积量最大幅度方差最大幅度方差最大幅度方差最大幅度方差N-类型14270.0126225000.00425.32641060.00521.17431090.0063V-类型15240.0194783810.00381.13911

41、070.02965.13881090.0095A-类型13110.0151233010.01172.94131060.02611.02061090.0157表2 三类ECG信号QRS波形及其高阶累积量特征信息统计表图7 特征提取结果显示界面8.6.4 心脏病诊断 图8 心脏病智能诊断专家系统界面图9 诊断结果的解释说明类型训练样本集1测试样本集1全体样本集样本数错检数错检率(%)样本数错检数错检率(%)样本数错检数错检率(%)误判漏判误判漏判误判漏判N6242000.0062422230.407489834570.12L807000.00807220.5080645350.50R725000.

42、00725381.5272503731.05A1267000.00126724304.26253483304.46a75000.007481936.49149101919.46J41000.00405320.00817312.35V2374000.00237458233.417122138512.65F399000.003995216.5279810213.88I236000.00236453.81472552.12j114000.0011321212.392272126.17E53000.0052227.69105425.71P3508000.003508940.3770162140.36f

43、491000.004912143.269824141.83小计16332000.00163281461461.791096983263260.59表3 SVM分类器的心电信号分类结果8.7遥感图像中的目标识别遥感图像中的目标识别、飞机识别、飞机识别、机场识别、机场识别、桥梁识别、桥梁识别 8.7.18.7.1 飞机识别飞机识别n算法是基于部分的（算法是基于部分的（Part-basedPart-based）和层次的（）和层次的（hierarchicalhierarchical），），包括两个部分：学习和识别，算法框图如下图所示。包括两个部分：学习和识别，算法框图如下图所示。n（1 1）学习阶段。

44、通过对属于同一类别的学习样本的学习，得到）学习阶段。通过对属于同一类别的学习样本的学习，得到基于部分的结构模型，并针对问题特点用一个属性关系图来表示基于部分的结构模型，并针对问题特点用一个属性关系图来表示。然后由这些图匹配得到一个模糊属性关系图，以此作为形状模。然后由这些图匹配得到一个模糊属性关系图，以此作为形状模型。型。n（2 2）识别阶段。用上述形状模型来检测识别形状。首先用柔性）识别阶段。用上述形状模型来检测识别形状。首先用柔性形态学凸集检测器检测出最显著的子部分，得到候选区域，然后形态学凸集检测器检测出最显著的子部分，得到候选区域，然后在其邻域检测识别其它子部分，分析这些子部分的几何结构关系在其邻域检测识别其它子部分，分析这些子部分的几何结构关系，以此作进一步证实。，以此作进一步证实。n下面就对各个部分作详细论述。下面就对各个部分作详细论述。训练训练结构模型输入图像凸子部分预处理预处理形状分解形状分解图表示图表示建建 模模识别子部分识别子部分关系证实关系证实描述图候选区 训练样本识别识别一组由40幅图像（其中包含97架飞机）组成的图像库被用来实验，其中之一示于下图。 (a)训练样本 (b) 识别结果 8.7.28.7.2、机场识别、机场识别机场识别系统界面：机场识别过程： 8.7.38.7.3 桥梁识别桥梁识别桥梁识别系统界面：桥梁识别过程：