CN105631226B - 基于等位线滤波反投影算法的压力分布重建方法 - Google Patents

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

Info

Publication number
CN105631226B
CN105631226B CN201610016387.9A CN201610016387A CN105631226B CN 105631226 B CN105631226 B CN 105631226B CN 201610016387 A CN201610016387 A CN 201610016387A CN 105631226 B CN105631226 B CN 105631226B
Authority
CN
China
Prior art keywords
mrow
pressure
electrode
pressure pad
msub
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
Application number
CN201610016387.9A
Other languages
English (en)
Other versions
CN105631226A (zh
Inventor
李士强
王新立
刘国强
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Electrical Engineering of CAS
Original Assignee
Institute of Electrical Engineering of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Institute of Electrical Engineering of CAS filed Critical Institute of Electrical Engineering of CAS
Priority to CN201610016387.9A priority Critical patent/CN105631226B/zh
Publication of CN105631226A publication Critical patent/CN105631226A/zh
Application granted granted Critical
Publication of CN105631226B publication Critical patent/CN105631226B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/50—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • Biomedical Technology (AREA)
  • Databases & Information Systems (AREA)
  • Epidemiology (AREA)
  • Theoretical Computer Science (AREA)
  • Primary Health Care (AREA)
  • General Health & Medical Sciences (AREA)
  • Data Mining & Analysis (AREA)
  • Pathology (AREA)
  • Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Force Measurement Appropriate To Specific Purposes (AREA)

Abstract

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

Description

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

Claims (1)

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)施加压力的过程中，根据预设的激励电流注入方式，向压力垫(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)标准化处理检测电压的变化，得到矩阵Dij
<mrow> <msub> <mi>D</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>U</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>V</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> </mrow> <msub> <mi>V</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
式中，i，j分别为注入电极对和接收电极对编号，i，j＝1，2，…N，N为电极个数，Vij为所述第一步计算得到的的从i对电极注入电流时，在第j对电极间的响应电压，Uij为所述第二步检测得到的从i对电极注入电流时，在第j对电极间的响应电压；
2)根据式(3)求解第m个单元的电导率变化σm
<mrow> <msub> <mi>&amp;sigma;</mi> <mi>m</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>W</mi> <mi>i</mi> </msub> <mo>&amp;lsqb;</mo> <mi>m</mi> <mo>,</mo> <mi>j</mi> <mo>&amp;rsqb;</mo> <mo>&amp;CenterDot;</mo> <msub> <mi>B</mi> <mi>i</mi> </msub> <mo>&amp;lsqb;</mo> <mi>m</mi> <mo>,</mo> <mi>j</mi> <mo>&amp;rsqb;</mo> <mo>&amp;CenterDot;</mo> <msub> <mi>D</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
式中，i，j分别为注入电极对和接收电极对编号，i，j＝1，2，…N，N为电极个数，m为有限单元编号，m＝1，2，…M，Wi[m,j]为加权系数矩阵W中的对应加权系数；
3)根据式(4)获得新的电阻率分布变化矩阵ρ′
<mrow> <msubsup> <mi>&amp;rho;</mi> <mi>m</mi> <mo>&amp;prime;</mo> </msubsup> <mo>=</mo> <mfrac> <mn>1</mn> <msub> <mi>&amp;sigma;</mi> <mi>m</mi> </msub> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
式中，ρm′为第m个单元的电阻率变化，将每一个单元的电阻率变化组合在一起，其可得到整个压力垫(1)的电阻率分布变化矩阵ρ′；
第四步：利用压阻特性重建压力分布图像
根据第一步测量的压力垫(1)的压阻特性关系f(ρ，P)，获得与电阻率ρ′相对应的压力P的分布图像。
CN201610016387.9A 2016-01-12 2016-01-12 基于等位线滤波反投影算法的压力分布重建方法 Active CN105631226B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610016387.9A CN105631226B (zh) 2016-01-12 2016-01-12 基于等位线滤波反投影算法的压力分布重建方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610016387.9A CN105631226B (zh) 2016-01-12 2016-01-12 基于等位线滤波反投影算法的压力分布重建方法

Publications (2)

Publication Number Publication Date
CN105631226A CN105631226A (zh) 2016-06-01
CN105631226B true CN105631226B (zh) 2018-04-06

Family

ID=56046155

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610016387.9A Active CN105631226B (zh) 2016-01-12 2016-01-12 基于等位线滤波反投影算法的压力分布重建方法

Country Status (1)

Country Link
CN (1) CN105631226B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102908141A (zh) * 2012-10-29 2013-02-06 中国人民解放军第四军医大学 一种反正切压缩的等位线反投影电阻抗成像图像重建方法
CN103267597A (zh) * 2013-01-09 2013-08-28 中国科学院电工研究所 基于压阻材料的电阻率成像柔性压力检测系统及检测方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102908141A (zh) * 2012-10-29 2013-02-06 中国人民解放军第四军医大学 一种反正切压缩的等位线反投影电阻抗成像图像重建方法
CN103267597A (zh) * 2013-01-09 2013-08-28 中国科学院电工研究所 基于压阻材料的电阻率成像柔性压力检测系统及检测方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
MIT中反投影矩阵的计算与数据处理方法;柯丽 等;《仪器仪表学报》;20141031;第35卷(第10期);第2256-2262页 *
基于足底压力分布时空HOG特征的步态识别方法;夏懿 等;《模式识别与人工智能》;20130630;第26卷(第6期);第529-536页 *
足底压力分布测量系统及临床应用;袁刚 等;《中国康复》;20030228;第18卷(第1期);第22-25页 *

Also Published As

Publication number Publication date
CN105631226A (zh) 2016-06-01

Similar Documents

Publication Publication Date Title
CN103505203B (zh) 基于脑电源定位方法的人体精神状态检测方法
US20170010130A1 (en) Pliable capacitive structure apparatus and methods
JP5075966B2 (ja) 排泄検知装置
CN103610446B (zh) 一种足部运动状况检测装置
CN104605850A (zh) 一种基于电阻抗断层成像技术的膀胱尿量实时监测装置和方法
CN105631225B (zh) 基于快速牛顿一步迭代算法的压力分布重建方法
JP6297385B2 (ja) 就寝者の心拍測定方法および心拍測定装置
CN105030237A (zh) 一种脊髓损伤定位装置
CN105631226B (zh) 基于等位线滤波反投影算法的压力分布重建方法
GB2563644A (en) Matron
Corbellini et al. Low-cost wearable measurement system for continuous real-time pedobarography
CN103860172B (zh) 体表电导分布测试仪
Shi et al. Design and implementation of a high-precision electrical impedance tomography data acquisition system for brain imaging
CN111938642A (zh) 基于电阻抗成像技术的足底压力面参数测量系统及方法
KR20200051350A (ko) 다층 구조 압력 센서를 이용한 압력 분포 영상화 장치 및 방법
Kukucka et al. The influence of the amplitude of driving signal in measurement of the active points on human skin
CN107543600A (zh) 带有脂肪测量功能的电子秤
He et al. A combined regularization algorithm for electrical impedance tomography system using rectangular electrodes array
CN203988052U (zh) 用于测量生物电阻抗的装置
Kukučka et al. Various effects in voltage chart measurement of active points on skin
CN201260667Y (zh) 电子舌诊仪
James et al. Detection of posture and motion by accelerometer sensors
CN203749407U (zh) 无创伤亚健康检测仪
Kukučka et al. THE INFLUENCE OF THE SKIN FATIGUE, ITS PERSPIRATION AND THE TIME OF STIMULATION IN MEASUREMENT OF THE ACTIVE POINTS ON HUMAN SKIN
CN211911605U (zh) 一种肌电检测电路和肌电检测装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant