基于等位线滤波反投影算法的压力分布重建方法技术

一种基于等位线滤波反投影算法的压力分布重建方法，它通过采用压阻材料制作的压力垫(1)在受到压力时电阻率发生变化的特性，采用激励电源(2)向压力垫(1)注入激励电流，同时检测压力垫(1)表面边缘电压，利用电压信息，采用等位线滤波反投影算法反演计算得到压力垫电阻率的变化，再根据材料的压阻特性重建压力垫(1)的受力的分布情况。获得的压力分布图像可以用来分析患者坐、立、行、卧等姿势，调整康复辅具。

【技术实现步骤摘要】

本专利技术涉及一种压力分布重建方法，特别涉及一种基于等位线滤波反投影算法的压力分布重建方法。
技术介绍
中国已经进入了严重的老龄化社会，关注老年人健康是当前社会关注的热点问题。在现阶段，有大约5％的老年人身受偏瘫疾病的困扰，而针对偏瘫患者最有效的治疗手段是矫形康复治疗。在矫形治疗中，对人体压力分布的测量对于指导病人更好地康复至关重要。目前，针对人体压力分布的检测手段主要是采用压力传感器阵列检测。所用的压力传感器有电容式、电阻式、压电式、压阻式等等。它通过阵列分布的压力传感器检测到的压力变化，将检测到的传感器节点的压力值代替传感器所处区域的压力值，来得到整个界面的压力分布。这种重建压力分布的方法需要密集布置传感器阵列并精确知道每一个传感器的位置，重建方法简单，但成本很高，不适合经济基础薄弱的病患群体使用。2013年中国科学院电工研究所提出来了一种新的人体压力检测技术，专利申请号为201310006835.3，它通过向压力垫中注入200kHz以下频率的交流电流信号，利用压力垫在受力时电阻率发生变化的特性，通过检测压力垫表面边缘电压信息来重建压力垫所受压力的分布。这种技术由于不需要高精度、密集分布的压力传感器阵列，而仅仅通过几个检测电极就可以完成压力分布的检测，减少了电子线路和器件，极大的降低了成本。这种技术由于采用了全新的压力检测原理，因此在压力分布重建方面也需要一种全新的重建算法。
技术实现思路
本专利技术的目的是克服现有技术的的缺点，提出一种基于等位线滤波反投影算法的压力分布重建方法。本专利技术可以满足现有的压力分布重建算法的要求。本专利技术的原理为：采用压阻材料制作的压力垫在受到压力时，其受力部位会出现电阻率特性的变化，通过向压力垫注入激励电流，同时检测压力垫表面边缘电压，利用电压信息反演计算得到压力垫电阻率的变化，再根据材料的压阻特性重建压力垫的受力的分布情况。应用本专利技术的人体压力分布检测系统包括压力垫、激励电源、信号处理系统和计算机系统。压力垫表面边缘均匀布置电极，用于输入激励电流信号和检测电压信号。激励电源的输出端与压力垫相连，激励电源用于提供激励电流。信号处理系统的输入端与压力垫相连，用于对压力垫输出的检测信号进行放大滤波等处理。计算机系统的输入端与信号信号处理系统的输出端相连，用于处理信号，反演计算压力分布。本专利技术包括四个步骤：第一步构建无压力时电阻率分布矩阵、反投影矩阵以及加权系数矩阵；第二步获得电压检测信号；第三步采用等位线滤波反投影算法求解电阻率变化分布矩阵；第四步利用压阻特性重建压力分布图像。各步骤的具体方法如下：第一步：构建无压力时电阻率分布矩阵、反投影矩阵以及加权系数矩阵首先通过检测压力垫自身的压阻特性，构建压力垫的电阻率分布矩阵，然后根据压力垫的形状、尺寸、激励检测方式构建反投影矩阵以及加权系数矩阵，具体如下：1)使用电流表或电压表，采用四探针法或范德堡法(VanDerPauwMethod)测量压力垫的电阻率，得到压力垫的电阻率分布矩阵ρ，通过压力加载等方式获得压力垫的压阻特性关系f(ρ，P)，其中P为施加的压力；2)根据压力垫的形状和尺寸，采用计算仿真软件建立压力垫有限元模型，该有限元模型中的有限单元数为M，将该有限元模型有限单元的电阻率设置成电阻率分布矩阵ρ的元素。确定电极个数N、电极排列方式、激励电流注入方式、电压测量方式，以及加权系数矩阵W；3)采用所建立的压力垫有限元模型，计算压力垫表面边缘电极处的响应电压Vj1，Vj2，Vij，其中，Vj1，Vj2分别为第j对电极上的响应电压，Vij为从i对电极注入电流后，在第j对电极间的响应电压，Vij＝Vj1-Vj2，i，j＝1，2，…N，N为电极个数；4)计算有限元模型中有限单元的电位Vm，m为有限单元编号，m＝1，2，…M，Vm取第m个单元三个顶点电压的平均值；5)构建第i对电极激励时的反投影矩阵Bi，矩阵中的元素Bi[m，j]可由下式(1)决定：式中，m为有限单元编号，m＝1，2，…M；i为注入电极对编号，j为检测电极对编号，i，j＝1，2，…N，N为电极个数，Bi[m，j]表示从第i对电极注入电流时，第j对电极间的电压变化投影到第m个有限单元时的反投影矩阵元素；第二步：获得电压检测信号在给压力垫施加压力的过程中，根据预设的激励电流注入方式，向压力垫注入电流，信号为稳恒交流信号，频率在200kHz以下；通过布置在压力垫表面边缘的电极，检测电极处的电压信号Uj1，Uj2，Uij，其中，Uj1，Uj2分别为第j对电极上的响应电压，Uij为从i对电极注入电流后，在第j对电极间的响应电压，Uij＝Uj1-Uj2，i，j分别为注入电极对和接收电极对编号，i，j＝1，2，…N，N为电极个数；然后通过信号处理系统实现电压信号的前置放大，滤波，二级放大等处理后存储；第三步：采用等位线滤波反投影算法求解电阻率变化分布矩阵1)根据下式(2)标准化处理检测电压的变化，得到矩阵DijDij=Uij-VijVij---(2)]]>式中，i，j分别为注入电极对和接收电极对编号，i，j＝1，2，…N，N为电极个数，Vij为所述第一步计算得到的的从i对电极注入电流时，在第j对电极间的响应电压，Uij为所述第二步检测得到的从i对电极注入电流时，在第j对电极间的响应电压；2)根据式(3)求解第m个单元的电导率变化σmσm=1NΣi=1NΣj=1NWi[m,j]·Bi[m,j]·Dij---(3)]]>式中，i，j分别为注入电极对和接收电极对编号，i，j＝1，2，…N，N为电极个数，m为有限单元编号，m＝1，2，…M；3)根据式(4)获得新的电阻率分布变化矩阵ρ′ρm′=1σm---(4)]]>式中，ρm′为第m个单元的电阻率变化，将每一个单元的电阻率变化组合在一起，其可得到整个压力垫的电阻率分布变化矩阵ρ′；第四步：利用压阻特性重建压力分布图像根据第一步测量的压力垫的压阻特性关系f(ρ，P)，获得与电阻率ρ′相对应的压力P的分布。通过以上过程即可重建压力垫受力分布的情况，由获得的压力分布图像就可以分析患者坐、立、行、卧等姿势，更好地调整康复辅具，指导病人康复。本专利技术基于等位线滤波反投影算法的压力分布重建方法，可用于偏瘫、脑瘫、残疾等患者在矫形康复过程中对人体压力的检测，可以为病人的康复治疗提供指导或康复效果评价，有利于推动人体压力检测本文档来自技高网...

【技术保护点】
一种基于等位线滤波反投影算法的压力分布重建方法，其特征在于：所述的基于快速牛顿一步迭代算法的压力分布重建方法包括四个步骤：第一步：构建无压力时电阻率分布矩阵、反投影矩阵以及加权系数矩阵首先通过检测压力垫(1)自身的压阻特性，构建压力垫(1)的电阻率分布矩阵，然后根据压力垫(1)的形状、尺寸、激励检测方式构建反投影矩阵以及加权系数矩阵，具体如下：1)使用电流表或电压表，采用四探针法或范德堡法测量压力垫(1)的电阻率，得到压力垫(1)的电阻率分布矩阵ρ，通过压力加载获得压力垫的压阻特性关系f(ρ，P)，其中P为施加的压力；2)根据压力垫(1)的形状、尺寸，采用计算仿真软件建立压力垫有限元模型，该有限元模型中的有限单元数为M，将该有限元模型有限单元的电阻率设置成电阻率分布矩阵ρ的元素；确定电极个数N、电极排列方式、激励电流注入方式、电压测量方式，以及加权系数矩阵W；3)采用所建立的压力垫有限元模型，计算压力垫表面边缘电极处的响应电压Vj1，Vj2，Vij，其中，Vj1，Vj2分别为第j对电极上的响应电压，Vij为从i对电极注入电流后，在第j对电极间的响应电压，Vij＝Vj1‑Vj2，i，j＝1，2，…N，N为电极个数；4)计算有限元模型中有限单元的电位Vm，m为有限单元编号，m＝1，2，…M，Vm取第m个单元三个顶点电压的平均值；5)构建第i对电极激励时的反投影矩阵Bi，矩阵中的元素Bi[m，j]可由下式(1)决定：式中，m为有限单元编号，m＝1，2，…M；i为注入电极对编号，j为检测电极对编号，i，j＝1，2，…N，N为电极个数，Bi[m，j]表示从第i对电极注入电流时，第j对电极间的电压变化投影到第m个有限单元时的反投影矩阵元素；第二步：获得电压检测信号在给压力垫(1)施加压力的过程中，根据预设的激励电流注入方式，向压力垫(1)注入电流，信号为稳恒交流信号，频率在200kHz以下；通过布置在压力垫(1)表面边缘的电极，检测电极处的电压信号Uj1，Uj2，Uij，其中，Uj1，Uj2分别为第j对电极上的响应电压，Uij为从i对电极注入电流后，在第j对电极间的响应电压，Uij＝Uj1‑Uj2，i，j分别为注入电极对和接收电极对编号，i，j＝1，2，…N，N为电极个数；然后通过信号处理系统实现电压信号的前置放大，滤波，二级放大等处理后存储；第三步：采用等位线滤波反投影算法求解电阻率变化分布矩阵1)根据下式(2)标准化处理检测电压的变化，得到矩阵DijDij=Uij-VijVij---(2)]]>式中，i，j分别为注入电极对和接收电极对编号，i，j＝1，2，…N，N为电极个数，Vij为所述第一步计算得到的的从i对电极注入电流时，在第j对电极间的响应电压，Uij为所述第二步检测得到的从i对电极注入电流时，在第j对电极间的响应电压；2)根据式(3)求解第m个单元的电导率变化σmσm=1NΣi=1NΣj=1NWi[m,j]·Bi[m,j]·Dij---(3)]]>式中，i，j分别为注入电极对和接收电极对编号，i，j＝1，2，…N，N为电极个数，m为有限单元编号，m＝1，2，…M；3)根据式(4)获得新的电阻率分布变化矩阵ρ′ρm′=1σm---(4)]]>式中，ρm′为第m个单元的电阻率变化，将每一个单元的电阻率变化组合在一起，其可得到整个压力垫(1)的电阻率分布变化矩阵ρ′；第四步：利用压阻特性重建压力分布图像根据第一步测量的压力垫(1)的压阻特性关系f(ρ，P)，获得与电阻率ρ′相对应的压力P的分布图像。...

【技术特征摘要】
1.一种基于等位线滤波反投影算法的压力分布重建方法，其特征在于：所述的基于快
速牛顿一步迭代算法的压力分布重建方法包括四个步骤：
第一步：构建无压力时电阻率分布矩阵、反投影矩阵以及加权系数矩阵
首先通过检测压力垫(1)自身的压阻特性，构建压力垫(1)的电阻率分布矩阵，然后
根据压力垫(1)的形状、尺寸、激励检测方式构建反投影矩阵以及加权系数矩阵，具体如
下：
1)使用电流表或电压表，采用四探针法或范德堡法测量压力垫(1)的电阻率，得到压
力垫(1)的电阻率分布矩阵ρ，通过压力加载获得压力垫的压阻特性关系f(ρ，P)，其中P
为施加的压力；
2)根据压力垫(1)的形状、尺寸，采用计算仿真软件建立压力垫有限元模型，该有限
元模型中的有限单元数为M，将该有限元模型有限单元的电阻率设置成电阻率分布矩阵ρ的
元素；确定电极个数N、电极排列方式、激励电流注入方式、电压测量方式，以及加权系数
矩阵W；
3)采用所建立的压力垫有限元模型，计算压力垫表面边缘电极处的响应电压Vj1，Vj2，
Vij，其中，Vj1，Vj2分别为第j对电极上的响应电压，Vij为从i对电极注入电流后，在第j
对电极间的响应电压，Vij＝Vj1-Vj2，i，j＝1，2，…N，N为电极个数；
4)计算有限元模型中有限单元的电位Vm，m为有限单元编号，m＝1，2，…M，Vm取第m个单元三个顶点电压的平均值；
5)构建第i对电极激励时的反投影矩阵Bi，矩阵中的元素Bi[m，j]可由下式(1)决定：
式中，m为有限单元编号，m＝1，2，…M；i为注入电极对编号，j为检测电极对编号，i，
j＝1，2，…N，N为电极个数，Bi[m，j]表示从第i对电极注入电流时，第j对电极间的电压
变化投影到第m个有限单元时的反投影矩阵元素；
第二步：获得电压检测信号
在给压力垫(1)施加压力的过程中，根据预设的激励电流注入方...

【专利技术属性】
技术研发人员：李士强，王新立，刘国强，
申请(专利权)人：中国科学院电工研究所，
类型：发明
国别省市：北京;11