CN106887019B - 一种足迹压力分布特征表示方法 - Google Patents
一种足迹压力分布特征表示方法 Download PDFInfo
- Publication number
- CN106887019B CN106887019B CN201710099973.9A CN201710099973A CN106887019B CN 106887019 B CN106887019 B CN 106887019B CN 201710099973 A CN201710099973 A CN 201710099973A CN 106887019 B CN106887019 B CN 106887019B
- Authority
- CN
- China
- Prior art keywords
- footprint
- image
- pressure distribution
- pressure
- hindfoot
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30196—Human being; Person
Landscapes
- Image Analysis (AREA)
Abstract
本发明公开了一种足迹压力分布特征表示方法,包括以下步骤:S1:采用足迹采集仪采集足迹压力图像信息,对足迹压力图像进行预处理;S2:提取足迹前足、后足压力分布特征,构建足迹压力分布特征图;S3:将前足、后足区域压力分布特征图沿着各方向投影,寻找投影值(压力值)最大方向作为主方向,进行压力分布特征图主方向规范化;S4:将得到的足迹压力分布特征进行特征选择。本方法提取的足迹压力分布特征很好地体现了人体行走时脚底压力的分布特点,无需进行赤足和穿袜足迹的区分,适用性更加广泛。
Description
技术领域
本发明涉及图像处理与生物特征识别技术领域,尤其涉及一种足迹压力分布特征表示方法。
背景技术
目前足迹的特征表示主要分为三类:几何形态特征、压力形态特征和纹理特征。几何形态特征是建立统一的足迹图像坐标系,进行几何度量,对足长、掌宽、跟宽以及从足跟到每个脚趾的边缘,足跟中心到每个脚趾的中心进行测量。足迹压力形态特征主要包括足底与地面之间的相互作用力信息和足底压力分布信息。主要方法有:1)利用足底压力分布信息构建足底平均、最大、最小压力时间曲线,直接利用曲线的数据点来构建特征向量;2)提取足迹压力中心与区域几何中心的位置关系,足部重压和足部重压分布(COP)特征进行实时身份认定;3)基于足底压力引入时空HOG(方向梯度直方图)算子,用于刻画足底动态压力分布的时空特征实现分类识别。足迹纹理特征通过基于变换的方法被提取,如傅里叶变换,离散余弦变换,小波变换,以及主成分分析和独立主成分分析的方法获取纹理特征进行足迹识别。
目前,足迹生物特征分析和身份识别应用技术的研究,尚缺乏系统和深入的应用研究基础。足迹特征通常用足型特征来表示,例如足长、足宽和跟宽,很难精确提取。由于采集设备的长期使用,足迹压力图像容易受灰尘等固定噪声和随机噪声的影响,很难在降低压力图像噪声水平的基础上,同时最大限度的保留压力图像所携带的信息。足迹图像携带大量信息,存在高维小样本的现象,使整个过程很难达到准确、高效的识别目的。
发明内容
根据现有技术存在的问题,本发明公开了一种足迹压力分布特征表示方法,其特征还在于:包括以下步骤:
S1:采用足迹采集仪采集足迹压力图像信息,对足迹压力图像进行预处理;
S2:提取足迹前足、后足压力分布特征,构建足迹压力分布特征图;
S3:将前足、后足区域压力分布特征图沿着各方向投影,寻找投影值(压力值)最大方向作为主方向,进行压力分布特征图主方向规范化;
S4:将得到的足迹压力分布特征进行特征选择。
所述S1具体还包括以下步骤:
S11:采用足迹采集仪采集足迹压力图像信息,对足迹压力图像进行预处理:采用最大类间方差法设定阈值,使灰度值小于等于该阈值的像素灰度值设置为0,其余部分保持不变,消除图像中背景噪声对足迹压力分布特征的干扰;
S12:通过寻找整个图像的连通区域,求取包含足迹区域的最大连通区域;
S13:逐行扫描计算足迹区域坐标的最小和最大值,获取包含足迹区域的最小外接矩形;
S14:判断包含足迹最小外接矩形的高度和宽度,如果高度大于宽度,按照3:2的高度比例对足迹图像进行分区,从而将足迹图像分为前足区域A和后足区域B,否则按照3:2的宽度比例对足迹图像进行前足、后足的划分;
S15:根据S14得到了足迹图像后足区域B,将后足区域B在任意角度射线方向上做投影变换得到后足区域在某一角度方向上的特性,寻找投影值最大的角度作为后足区域的主方向,依据主方向的与水平方向的倾斜角对整幅足迹图像进行倾斜校正;
S16:将校正后的图像再一次进行步骤S13获取校正后足迹图像最小外接矩形并且进行前足、后足划分的操作。
S2采用如下方式:
S21:将前足区域A和后足区域B的图像质心近似为着力点cA,cB,分别以cA,cB为中心点,将压力图像从笛卡尔坐标(x,y)空间映射到极坐标(ρ,θ)空间,变换公式为:
式中:ρ表示点(x,y)与中心点的距离,最大值ρmax为前足区域最小外接矩形对角线长度的一半;θ表示该向量与x轴正向的夹角,最大值θmax为360度,映射后的图像列方向表示角度θ,行方向表示到中心点的距离ρ;
S22:在变换过程中,将ρmax分成W份,从而在前足或后足区域形成W个同心圆,从中心点出发做M条射线,将θmax均分成M等份,这M条射线与W个同心圆相交形成MW个互不重叠的区域,计算每个区域中像素点灰度值的均值,将各区域压力的均值作为特征值来描述足迹图像的压力分布。
S3步骤如下:
S31:分别寻找前足、后足压力分布特征图的主方向,即将特征图沿着各方向进行投影,投影值(压力值)最大的方向即为主方向,将主方向定为起始角度,即将主方向所在的列移到特征矩阵的首列,将特征图进行循环左移,从而使特征图具有旋转不变性;
S32:通过上述步骤得到了足迹压力图像的压力分布特征F。
S4采用如下方式:
S41:足迹压力图像特征集为FT={F1,F2,F3,...,FN}T,N为训练压力图像特征集的个数,Fi是足迹压力分布特征的列向量表示,其维数为:MW×1,求取特征均值和协方差:
S42:对协方差矩阵G进行特征值分解,求前d个最大特征值对应的特征向量:x1,x2,...,xd作为投影向量,对于特征矩阵FT,得到一组投影特征向量:Y1,Y2,...,Yd,令BMW×d=[Y1,Y2,...,Yd]MW×d,BMW×d为F的特征矩阵;
由于采用了上述技术方案,本发明提供的一种足迹压力分布特征表示方法,具有如下有益效果:1)、提出的足迹压力分布特征与人走路方式存在相关性,特征更加稳定,区分度更高,从而提高了识别的准确度。2)、提取的足迹压力分布特征具有旋转不变性。3)、提取的足迹压力分布特征很好地体现了人体行走时脚底压力的分布特点,无需进行赤足和穿袜足迹的区分,适用性更加广泛。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明的足迹压力分布特征图;
图2为本发明的足迹压力图像的最大连通区域示意图;
图3为本发明的校正规范足迹压力图像示意图
具体实施方式
为使本发明的技术方案和优点更加清楚,下面结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚完整的描述:
如图1-图3所示的一种足迹压力分布特征表示方法,它表示足迹在承受人体重力作用下,足底各部分相对于着力点的压力分布情况。该特征表示方法可表示为:F=(Ap,Bp)
Ap表示前足区域的压力分布特征图,Bp表示后足区域的压力分布特征图。压力分布特征图中的每一个元素值表示的是压力图对应区域内各点灰度值(压力)的均值,其定义如下:
式中,fA,fB分别表示足迹前足区域压力图像、后足区域压力图像,θi表示压力图像中对应区域上边界相对于水平线的夹角,同理θi-1表示下边界的夹角;ρi,j表示区域后边界相对于着力点的距离,ρi,j-1表示前边界相对于着力点的距离,N(θi,ρi,j)表示区域内的像素数。将足迹图像分为前足区和后足区,在前后足区,分别以着力点为中心做W个同心圆,这就产生了W个圆环。再以着力点为起始点,做M条射线,且各条射线之间的夹角相同,这些射线与圆环就将前足区或后足区各分成了MW份。由此,我们以各区域内压力的均值作为特征值来表征足迹,这些特征值按角度及其到着力点的距离顺序排列便形成足迹压力分布图。包括以下步骤:
S1:采用足迹采集仪采集足迹压力图像信息,对足迹压力图像进行预处理;
S2:提取足迹前足、后足压力分布特征,构建足迹压力分布特征图;
S3:将前足、后足区域压力分布特征图沿着各方向投影,寻找投影值(压力值)最大方向作为主方向,进行压力分布特征图主方向规范化;
S4:将得到的足迹压力分布特征进行特征选择。
所述S1具体还包括以下步骤:
S11:采用足迹采集仪采集足迹压力图像信息,对足迹压力图像进行预处理:采用最大类间方差法设定阈值,使灰度值小于等于该阈值的像素灰度值设置为0,其余部分保持不变,消除图像中背景噪声对足迹压力分布特征的干扰;
S12:通过寻找整个图像的连通区域,求取包含足迹区域的最大连通区域;
S13:逐行扫描计算足迹区域坐标的最小和最大值,获取包含足迹区域的最小外接矩形;
S14:判断包含足迹最小外接矩形的高度和宽度,如果高度大于宽度,按照3:2的高度比例对足迹图像进行分区,从而将足迹图像分为前足区域A和后足区域B,否则按照3:2的宽度比例对足迹图像进行前足、后足的划分;
S15:根据S14得到了足迹图像后足区域B,将后足区域B在任意角度射线方向上做投影变换得到后足区域在某一角度方向上的特性,寻找投影值最大的角度作为后足区域的主方向,依据主方向的与水平方向的倾斜角对整幅足迹图像进行倾斜校正;
S16:将校正后的图像再一次进行步骤S13获取校正后足迹图像最小外接矩形并且进行前足、后足划分的操作。
S2采用如下方式:
S21:将前足区域A和后足区域B的图像质心近似为着力点cA,cB,分别以cA,cB为中心点,将压力图像从笛卡尔坐标(x,y)空间映射到极坐标(ρ,θ)空间,变换公式为:
式中:ρ表示点(x,y)与中心点的距离,最大值ρmax为前足区域最小外接矩形对角线长度的一半;θ表示该向量与x轴正向的夹角,最大值θmax为360度,映射后的图像列方向表示角度θ,行方向表示到中心点的距离ρ;
S22:在变换过程中,将ρmax分成W份,从而在前足或后足区域形成W个同心圆,从中心点出发做M条射线,将θmax均分成M等份,这M条射线与W个同心圆相交形成MW个互不重叠的区域,计算每个区域中像素点灰度值的均值,将各区域压力的均值作为特征值来描述足迹图像的压力分布。
S3步骤如下:
S31:分别寻找前足、后足压力分布特征图的主方向,即将特征图沿着各方向进行投影,投影值(压力值)最大的方向即为主方向,将主方向定为起始角度,即将主方向所在的列移到特征矩阵的首列,将特征图进行循环左移,从而使特征图具有旋转不变性;
S32:通过上述步骤得到了足迹压力图像的压力分布特征F。
S4采用如下方式:
S41:足迹压力图像特征集为FT={F1,F2,F3,...,FN}T,N为训练压力图像特征集的个数,Fi是足迹压力分布特征的列向量表示,其维数为:MW×1,求取特征均值和协方差:
S42:对协方差矩阵G进行特征值分解,求前d个最大特征值对应的特征向量:x1,x2,...,xd作为投影向量,对于特征矩阵FT,得到一组投影特征向量:Y1,Y2,...,Yd,令BMW×d=[Y1,Y2,...,Yd]MW×d,BMW×d为F的特征矩阵;
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都应涵盖在本发明的保护范围之内。
Claims (1)
1.一种足迹压力分布特征表示方法,其特征还在于:包括以下步骤:
S1:采用足迹采集仪采集足迹压力图像信息,对足迹压力图像进行预处理;
S2:提取足迹前足、后足压力分布特征,构建足迹压力分布特征图;
S3:将前足、后足区域压力分布特征图沿着各方向投影,寻找投影值,即 压力值, 最大方向作为主方向,进行压力分布特征图主方向规范化;
S4:将得到的足迹压力分布特征进行特征选择;
所述S1具体还包括以下步骤:
S11:采用足迹采集仪采集足迹压力图像信息,对足迹压力图像进行预处理:采用最大类间方差法设定阈值,使灰度值小于等于该阈值的像素灰度值设置为0,其余部分保持不变,消除图像中背景噪声对足迹压力分布特征的干扰;
S12:通过寻找整个图像的连通区域,求取包含足迹区域的最大连通区域;
S13:逐行扫描计算足迹区域坐标的最小和最大值,获取包含足迹区域的最小外接矩形;
S14:判断包含足迹最小外接矩形的高度和宽度,如果高度大于宽度,按照3:2的高度比例对足迹图像进行分区,从而将足迹图像分为前足区域A和后足区域B,否则按照3:2的宽度比例对足迹图像进行前足、后足的划分;
S15:根据S14得到了足迹图像后足区域B,将后足区域B在任意角度射线方向上做投影变换得到后足区域在某一角度方向上的特性,寻找投影值最大的角度作为后足区域的主方向,依据主方向的与水平方向的倾斜角对整幅足迹图像进行倾斜校正;
S16:将校正后的图像再一次进行步骤S13获取校正后足迹图像最小外接矩形并且进行前足、后足划分的操作;
S2采用如下方式:
S21:将前足区域A和后足区域B的图像质心近似为着力点cA,cB,分别以cA,cB为中心点,将压力图像从笛卡尔坐标(x,y)空间映射到极坐标(ρ,θ)空间,变换公式为:
式中:ρ表示点(x,y)与中心点的距离,最大值ρmax为前足区域最小外接矩形对角线长度的一半;θ表示该向量与x轴正向的夹角,最大值θmax为360度,映射后的图像列方向表示角度θ,行方向表示到中心点的距离ρ;
S22:在变换过程中,将ρmax分成W份,从而在前足或后足区域形成W个同心圆,从中心点出发做M条射线,将θmax均分成M等份,这M条射线与W个同心圆相交形成MW个互不重叠的区域,计算每个区域中像素点灰度值的均值,将各区域压力的均值作为特征值来描述足迹图像的压力分布;
S3步骤如下:
S31:分别寻找前足、后足压力分布特征图的主方向,即将特征图沿着各方向进行投影,投影值,即 压力值, 最大的方向即为主方向,将主方向定为起始角度,即将主方向所在的列移到特征矩阵的首列,将特征图进行循环左移,从而使特征图具有旋转不变性;
S32:通过上述步骤得到了足迹压力图像的压力分布特征F;
所述压力分布特征的选择S4采用如下方式:
S41:足迹压力图像特征集为FT={F1,F2,F3,...,FN}T,N为训练压力图像特征集的个数,Fi是足迹压力分布特征的列向量表示,其维数为:MW×1,求取特征均值和协方差:
S42:对协方差矩阵G进行特征值分解,求前d个最大特征值对应的特征向量:x1,x2,...,xd作为投影向量,对于特征矩阵FT,得到一组投影特征向量:Y1,Y2,...,Yd,令BMW×d=[Y1,Y2,...,Yd]MW×d,BMW×d为F的特征矩阵;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710099973.9A CN106887019B (zh) | 2017-02-23 | 2017-02-23 | 一种足迹压力分布特征表示方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710099973.9A CN106887019B (zh) | 2017-02-23 | 2017-02-23 | 一种足迹压力分布特征表示方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106887019A CN106887019A (zh) | 2017-06-23 |
CN106887019B true CN106887019B (zh) | 2020-04-21 |
Family
ID=59179431
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710099973.9A Active CN106887019B (zh) | 2017-02-23 | 2017-02-23 | 一种足迹压力分布特征表示方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106887019B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107423715B (zh) * | 2017-07-31 | 2020-05-19 | 大连海事大学 | 一种基于多特征联合决策的脚印自动识别方法 |
CN109871827B (zh) * | 2019-03-14 | 2022-10-25 | 大连海事大学 | 一种结合区域置信度和压力分布方向强度的足迹表达方法 |
CN110188694B (zh) * | 2019-05-30 | 2022-10-21 | 大连海事大学 | 一种基于压力特征的穿鞋足迹序列识别方法 |
CN110292384A (zh) * | 2019-06-26 | 2019-10-01 | 浙江大学 | 一种基于足底压力数据分布的智能足弓指数测量方法 |
CN111358469B (zh) * | 2020-04-15 | 2024-01-02 | 青岛一小步科技有限公司 | 一种足底压力检测装置及足部健康检测装置和检测方法 |
CN112535472A (zh) * | 2020-12-18 | 2021-03-23 | 深圳行正科技有限公司 | 基于体重自校准的足底动态步态分析系统及分析方法 |
CN112766142B (zh) * | 2021-01-15 | 2022-12-20 | 天津大学 | 足底压力图像处理方法、识别方法及步态分析系统 |
CN113361342B (zh) * | 2021-05-20 | 2022-09-20 | 杭州好学童科技有限公司 | 一种基于多模态的人体坐姿检测方法及装置 |
CN113569872B (zh) * | 2021-08-10 | 2024-03-19 | 大连海事大学 | 一种基于压力显著性的多分辨穿鞋足迹序列识别方法 |
CN116793548A (zh) * | 2023-05-11 | 2023-09-22 | 西南交通大学 | 一种隧道锚杆压力数据表示预警方法及相关设备 |
CN117179741B (zh) * | 2023-07-24 | 2024-05-28 | 北京华益精点生物技术有限公司 | 基于压力标记足底位置的方法及相关设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101327126A (zh) * | 2008-07-23 | 2008-12-24 | 天津大学 | 人体赤足迹形态学特征提取方法 |
CN104156914A (zh) * | 2014-07-18 | 2014-11-19 | 北京海鑫科金高科技股份有限公司 | 一种足迹图像的自动处理方法 |
CN104573678A (zh) * | 2015-02-03 | 2015-04-29 | 大连恒锐科技股份有限公司 | 基于足迹图像的性别特征分析的方法及装置 |
CN104574426A (zh) * | 2015-02-03 | 2015-04-29 | 大连恒锐科技股份有限公司 | 基于赤足或穿袜足迹图像的人身特征分析的方法及装置 |
-
2017
- 2017-02-23 CN CN201710099973.9A patent/CN106887019B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101327126A (zh) * | 2008-07-23 | 2008-12-24 | 天津大学 | 人体赤足迹形态学特征提取方法 |
CN104156914A (zh) * | 2014-07-18 | 2014-11-19 | 北京海鑫科金高科技股份有限公司 | 一种足迹图像的自动处理方法 |
CN104573678A (zh) * | 2015-02-03 | 2015-04-29 | 大连恒锐科技股份有限公司 | 基于足迹图像的性别特征分析的方法及装置 |
CN104574426A (zh) * | 2015-02-03 | 2015-04-29 | 大连恒锐科技股份有限公司 | 基于赤足或穿袜足迹图像的人身特征分析的方法及装置 |
Non-Patent Citations (3)
Title |
---|
商标图像检索研究与系统实现;张玲;《中国优秀博硕士学位论文全文数据库(硕士) 信息科技辑》;20061215(第12期);论文第3章 * |
基于步态触觉信息的足迹提取与动态识别方法研究;汪飞跃;《中国优秀硕士学位论文全文数据库 信息科技辑》;20151015(第10期);论文摘要、第3章 * |
基于类内类间离散度的分类器设计方法;郭亚琴 等;《信息技术》;20150504(第5期);论文第2节 * |
Also Published As
Publication number | Publication date |
---|---|
CN106887019A (zh) | 2017-06-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106887019B (zh) | 一种足迹压力分布特征表示方法 | |
US11631192B2 (en) | Robot climbing control method and device and storage medium and robot | |
Jiang et al. | Automatic detection of crop rows based on multi-ROIs | |
Vidović et al. | Crop row detection by global energy minimization | |
CN105046271B (zh) | 基于模板匹配的melf元件定位与检测方法 | |
US20180114313A1 (en) | Medical Image Segmentation Method and Apparatus | |
WO2017059591A1 (zh) | 手指静脉识别方法及装置 | |
CN103218605B (zh) | 一种基于积分投影与边缘检测的快速人眼定位方法 | |
WO2019102608A1 (en) | Image processing device, image processing method, and image processing program | |
CN106203375A (zh) | 一种基于人脸图像中人脸和人眼检测的瞳孔定位方法 | |
Esteves et al. | Gradient convergence filters and a phase congruency approach for in vivo cell nuclei detection | |
Khotimah et al. | Iris recognition using feature extraction of box counting fractal dimension | |
CN113536958A (zh) | 导航路径提取方法、装置、农业机器人和存储介质 | |
CN108629297A (zh) | 一种基于空域自然场景统计的遥感图像云检测方法 | |
CN111879244B (zh) | 综采工作面液压支架支撑高度和顶梁倾角测量方法 | |
CN106127735A (zh) | 一种设施蔬菜边缘清晰类叶面病斑分割方法及装置 | |
CN115690086A (zh) | 一种基于对象的高分辨率遥感影像变化检测方法及系统 | |
CN104408721B (zh) | 基于背景密度估计的印章图像提取方法 | |
CN109871827B (zh) | 一种结合区域置信度和压力分布方向强度的足迹表达方法 | |
CN104036245A (zh) | 一种基于在线特征点匹配的生物特征识别方法 | |
CN116778537B (zh) | 一种指背静脉图像识别过程中指甲的检测与处理方法 | |
Purnama et al. | Follicle detection on the usg images to support determination of polycystic ovary syndrome | |
CN111402202B (zh) | 足底压力图像相似度评价方法 | |
Devasia et al. | Automatic Optic Disc Localization and Segmentation using Swarm Intelligence. | |
Sui et al. | Automatic image registration based on shape features and multi-scale image segmentation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |