CN107505232B - 运动信息获取方法及装置 - Google Patents

运动信息获取方法及装置 Download PDF

Info

Publication number
CN107505232B
CN107505232B CN201710649554.8A CN201710649554A CN107505232B CN 107505232 B CN107505232 B CN 107505232B CN 201710649554 A CN201710649554 A CN 201710649554A CN 107505232 B CN107505232 B CN 107505232B
Authority
CN
China
Prior art keywords
signal
slope
time diagram
vibration
position time
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
CN201710649554.8A
Other languages
English (en)
Other versions
CN107505232A (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.)
Wuxi Hisky Medical Technologies Co Ltd
Original Assignee
Wuxi Hisky Medical Technologies Co Ltd
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 Wuxi Hisky Medical Technologies Co Ltd filed Critical Wuxi Hisky Medical Technologies Co Ltd
Priority to CN201710649554.8A priority Critical patent/CN107505232B/zh
Publication of CN107505232A publication Critical patent/CN107505232A/zh
Priority to AU2018301989A priority patent/AU2018301989A1/en
Priority to JP2020503053A priority patent/JP7403163B2/ja
Priority to PCT/CN2018/088406 priority patent/WO2019015398A1/zh
Priority to RU2020106886A priority patent/RU2020106886A/ru
Priority to EP18835507.7A priority patent/EP3656311A4/en
Priority to CA3065575A priority patent/CA3065575C/en
Priority to BR112019025264-9A priority patent/BR112019025264B1/pt
Priority to KR1020207004623A priority patent/KR102235737B1/ko
Priority to TW107121597A priority patent/TWI686180B/zh
Application granted granted Critical
Publication of CN107505232B publication Critical patent/CN107505232B/zh
Priority to US16/689,022 priority patent/US11517289B2/en
Priority to AU2022200463A priority patent/AU2022200463A1/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0833—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures
    • A61B8/085—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures for locating body or organic structures, e.g. tumours, calculi, blood vessels, nodules
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5223—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N11/00—Investigating flow properties of materials, e.g. viscosity, plasticity; Analysing materials by determining flow properties
    • G01N11/10—Investigating flow properties of materials, e.g. viscosity, plasticity; Analysing materials by determining flow properties by moving a body within the material
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/0048—Detecting, measuring or recording by applying mechanical forces or stimuli
    • A61B5/0051—Detecting, measuring or recording by applying mechanical forces or stimuli by applying vibrations
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/0093—Detecting, measuring or recording by applying one single type of energy and measuring its conversion into another type of energy
    • A61B5/0097—Detecting, measuring or recording by applying one single type of energy and measuring its conversion into another type of energy by applying acoustic waves and detecting light, i.e. acoustooptic measurements
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/42—Detecting, measuring or recording for evaluating the gastrointestinal, the endocrine or the exocrine systems
    • A61B5/4222—Evaluating particular parts, e.g. particular organs
    • A61B5/4244—Evaluating particular parts, e.g. particular organs liver
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48—Diagnostic techniques
    • A61B8/485—Diagnostic techniques involving measuring strain or elastic properties
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H5/00—Measuring propagation velocity of ultrasonic, sonic or infrasonic waves, e.g. of pressure waves
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H9/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H9/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means
    • G01H9/008—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means by using ultrasonic waves
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/20—Analysis of motion
    • G06T7/262—Analysis of motion using transform domain methods, e.g. Fourier domain methods
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2576/00—Medical imaging apparatus involving image processing or analysis
    • 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
    • G16H30/00—ICT specially adapted for the handling or processing of medical images
    • G16H30/40—ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Medical Informatics (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Molecular Biology (AREA)
  • Public Health (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Immunology (AREA)
  • Physiology (AREA)
  • Biochemistry (AREA)
  • Mathematical Analysis (AREA)
  • Data Mining & Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Gastroenterology & Hepatology (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • Operations Research (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • Endocrinology (AREA)

Abstract

本发明公开了运动信息获取方法及装置。对在介质中传播的振动的检测信号进行频域变换，得到频域信号，接着去除该频域信号中位于设定振动速度范围之外的信号，即只保留振动的信号，再沿设定的振动传播方向得到位置时间图。无需通过复杂的计算对振动的传播进行运动估计，只需通过在频域进行处理来判断振动的有无，即可得到位置时间图，是一种高效的运动信息获取方式。

Description

运动信息获取方法及装置
技术领域
本发明涉及测量技术领域，特别涉及运动信息获取方法及装置。
背景技术
对介质进行振动激励时，振动在介质中的传播特性与介质的粘弹性有关，通过测量振动的传播特性，可以对介质的粘弹性进行定量。获取振动的传播特性，又需要先利用针对该振动的检测信号，来获取该振动的运动信息。
上述原理目前已被应用至多个技术领域，以医学检测为例，检测肝脏、甲状腺和肌肉等器官或组织时，通过对介质的粘弹性进行定量，可以对病变进行定位。
因此，如何高效的获取在介质中传播的振动的运动信息，是一个需要解决的问题。
发明内容
本发明实施例提供了运动信息获取方法及装置。为了对披露的实施例的一些方面有一个基本的理解，下面给出了简单的概括。该概括部分不是泛泛评述，也不是要确定关键/重要组成元素或描绘这些实施例的保护范围。其唯一目的是用简单的形式呈现一些概念，以此作为后面的详细说明的序言。
根据本发明实施例的第一方面，提供了一种运动信息获取方法，所述方法包括：
对在介质中传播的振动的检测信号进行频域变换，得到频域信号；
去除所述频域信号中位于设定振动速度范围之外的信号，得到处理信号；
利用所述处理信号得到所述振动的位置时间图。
根据所述方法，作为可选的实施例一，所述去除所述频域信号中位于设定振动速度范围之外的信号，得到处理信号，包括：
对所述频域信号进行滤波或特征值选取，得到处理信号；
所述滤波的参数与所述设定振动速度范围相关；
所述特征值选取与所述设定振动速度范围相关。
根据所述方法，作为可选的实施例二，所述利用所述处理信号得到所述振动的位置时间图，包括：
根据设定的振动传播方向，利用所述处理信号得到所述振动的位置时间图。
根据所述方法，作为可选的实施例三，所述方法还包括：
对所述位置时间图进行图像分割；
提取图像特征；
利用所述图像特征进行线性拟合，得到所述位置时间图的斜率线的斜率；
根据所述斜率，计算所述介质的粘弹性参数。
根据所述方法，作为可选的实施例四，所述方法还包括：
沿所述位置时间图的预设角度范围内的各个角度进行角度投影，确定信号能量最大的角度所对应的所述位置时间图的斜率；
根据所述斜率，得到所述介质的粘弹性参数。
根据所述实施例四，作为可选的实施例五，所述沿所述位置时间图的预设角度范围内的各个角度进行角度投影，确定信号能量最大的角度所对应的所述位置时间图的斜率，包括：
对所述位置时间图进行沿预设角度范围内的各个角度的积分计算；
将积分值最大的角度，确定为所述位置时间图的斜率线的斜率角；
利用所述斜率角确定所述斜率线的斜率。
根据本发明实施例的第二方面，提供了一种运动信息获取装置，所述装置包括：
第一处理模块，用于对在介质中传播的振动的检测信号进行频域变换，得到频域信号；
第二处理模块，用于去除所述频域信号中位于设定振动速度范围之外的信号，得到处理信号；
获取模块，用于利用所述处理信号得到所述振动的位置时间图。
根据所述方法，作为可选的实施例一，所述第二处理模块对所述频域信号进行滤波或特征值选取，得到处理信号；
所述滤波的参数与所述设定振动速度范围相关；
所述特征值选取与所述设定振动速度范围相关。
根据所述装置，作为可选的实施例二，所述获取模块根据设定的振动传播方向，利用所述处理信号得到所述振动的位置时间图。
根据所述装置，作为可选的实施例三，所述装置还包括：
粘弹性定量模块，用于对所述位置时间图进行图像分割；提取图像特征；利用所述图像特征进行线性拟合，得到所述位置时间图的斜率线的斜率；
根据所述斜率，计算所述介质的粘弹性参数。
根据所述装置，作为可选的实施例四，所述装置还包括：
粘弹性定量模块，用于沿所述位置时间图的预设角度范围内的各个角度进行角度投影，确定信号能量最大的角度所对应的所述位置时间图的斜率；根据所述斜率，得到所述介质的粘弹性参数。
根据所述实施例四，作为可选的实施例五，所述粘弹性定量模块，包括：
计算子模块，用于对所述位置时间图进行沿预设角度范围内的各个角度的积分计算；
确定子模块，用于将所述计算子模块计算出的积分最大值的角度，确定为所述位置时间图的斜率线的斜率角；利用所述斜率角确定所述斜率线的斜率；
定量子模块，用于根据所述斜率，得到所述介质的粘弹性参数。
根据本发明实施例的第三方面，提供了一种运动信息获取装置，所述装置包括：
存储器，存储执行指令；
处理器，被配置于读取所述执行指令，完成如下操作：
用于对在介质中传播的振动的检测信号进行频域变换，得到频域信号；
去除所述频域信号中位于设定振动速度范围之外的信号，得到处理信号；
利用所述处理信号得到所述振动的位置时间图。
本发明实施例提供的技术方案可以包括以下有益效果：
对在介质中传播的振动的检测信号进行频域变换，得到频域信号，接着去除该频域信号中位于设定振动速度范围之外的信号，即只保留振动的信号，再得到振动的位置时间图。无需通过复杂的计算对振动的传播进行运动估计，只需通过在频域进行处理来判断振动的有无，即可得到位置时间图，是一种高效的运动信息获取方式。
应当理解的是，以上的一般描述和后文的细节描述仅是示例性和解释性的，并不能限制本发明。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分，示出了符合本发明的实施例，并与说明书一起用于解释本发明的原理。
图1是根据一示例性实施例示出的运动信息获取方法的流程图；
图2是根据一示例性实施例示出的介质粘弹性定量方法的流程图；
图3是根据一示例性实施例示出的介质粘弹性定量方法的流程图；
图4是根据一示例性实施例示出的介质粘弹性定量方法的流程图；
图5是根据一示例性实施例示出的运动信息获取方法的流程图；
图6是根据一示例性实施例示出的运动信息获取装置的框图；
图7是根据一示例性实施例示出的运动信息获取装置的框图；
图8是图7中所示粘弹性定量模块的框图；
图9是根据一示例性实施例示出的运动信息获取装置的框图；
图10是根据一示例性实施例示出的运动信息获取装置的框图。
具体实施方式
以下描述和附图充分地示出本发明的具体实施方案，以使本领域的技术人员能够实践它们。实施例仅代表可能的变化。除非明确要求，否则单独的部件和功能是可选的，并且操作的顺序可以变化。一些实施方案的部分和特征可以被包括在或替换其他实施方案的部分和特征。本发明的实施方案的范围包括权利要求书的整个范围，以及权利要求书的所有可获得的等同物。在本文中，各实施方案可以被单独地或总地用术语“发明”来表示，这仅仅是为了方便，并且如果事实上公开了超过一个的发明，不是要自动地限制该应用的范围为任何单个发明或发明构思。本文中，诸如第一和第二等之类的关系术语仅仅用于将一个实体或者操作与另一个实体或操作区分开来，而不要求或者暗示这些实体或操作之间存在任何实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素。本文中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的结构、产品等而言，由于其与实施例公开的部分相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
图1是根据一示例性实施例示出的运动信息获取方法的流程图。如图1所示，该方法包括如下步骤。
在步骤11中，对在介质中传播的振动的检测信号进行频域变换，得到频域信号。
在步骤12中，去除频域信号中位于设定振动速度范围之外的信号，得到处理信号。
在步骤13中，利用处理信号得到振动的位置时间图。
介质受到振动激励后，振动在介质中传播，在不同的时刻，波前会沿着传播方向到达不同的位置，这种位置和时刻的对应关系就是振动的运动信息。目前常用的运动信息获取方法，利用的是在介质振动时，针对该振动的检测信号会发生相位去相关，根据这种相位去相关的特性，可以通过互相关、自相关、光流等算法，得到该振动的位置时间图。上述算法可以选择任意基于块匹配的方法。这种常用的运动信息获取方法，相当于需要通过介质的位移和应变等信息先对振动的传播进行运动估计，然后才能得到该振动的位置时间图。
在本示例性实施例中，利用检测信号中包含振动信息产生多普勒效应的特性，对检测信号在成像时间维度进行频域变换，得到频域信号，接着去除该频域信号中位于设定振动速度范围之外的信号，即去除相对静止或振动速度低的信号，再得到振动的位置时间图。可见，本示例性实施例中的运动信息获取方法，不需要复杂的计算，通过在频域进行处理，就可以得到不以位移或应变为特征的位置时间图。这种方法相当于无需先对振动的传播进行运动估计，只需判断振动的有无即可得到位置时间图，是一种高效的运动信息的获取方法。
在一示例性实施例中，上述步骤11中，可以采用多种方式进行频域变换，例如采用傅里叶变换或奇异值分解来进行频域变换。
在一示例性实施例中，上述步骤12中，去除频域信号中位于设定振动速度范围之外的信号，得到处理信号，可以通过对频域信号滤波或进行特征值选取来实现。以采用滤波实现上述去除操作为例，为了只保留振动速度在设定范围之内的信号，考虑信号在空间、时间上的采样率，结合设定的振动速度范围，例如0.1米每秒至30米每秒，来设置滤波器的滤波参数，就可以基于该滤波去除频域信号中位于设定振动速度范围之外的信号。在采用特征值选取实现上述去除操作时，同样可以设置特征值的选取与设定的振动速度范围相关，从而实现去除频域信号中位于设定振动速度范围之外的信号。通过步骤12，可以只将振动的信号保留下来，提高后续形成位置时间图时的准确度。
在一示例性实施例中，通过机械振动、声辐射力或其他可以产生振动的方式，对介质进行振动激励后，介质产生振动，振动在介质中传播。由于上述振动在介质中的传播速度有限，因此可以利用检测信号对介质进行动态成像。上述检测波可以是光波、超声波等。上述动态成像可以是一维成像、二维成像或三维成像等。无论成像方式是哪种，都可以根据设定的振动传播方向，利用通过去除操作之后的处理信号得到振动的位置时间图。上述设定的振动传播方向，在振动只在一个传播方向上传播时，为振动的实际传播方向，在振动在多个传播方向上传播时，为选择出的某一个传播方向。例如，当介质为一均匀薄片时，对介质进行振动激励后，振动将沿该薄片的延展方向进行传播，此时设定的振动传播方向就是该振动的实际传播方向。又例如，当介质为立体不规则形状时，振动传播的波前为立体形状，例如振动传播的波前为椭球，则沿不同的振动传播方向得到的位置时间图是不同的，此时设定的振动传播方向就是选择出的某一个感兴趣的传播方向。上述感兴趣的传播方向根据实际要测量的方向来确定，例如可以是振动传播最快的方向、振动传播最慢的方向和振动传播速度为某一区间的方向中的至少一个。
在一示例性实施例中，图1所示的方法中还可以包括介质粘弹性定量的步骤。图2是根据一示例性实施例示出的介质粘弹性定量方法的流程图，该流程基于图1所示流程实现，包括如下步骤。
在步骤21中，沿位置时间图上的预设角度范围内的各个角度进行角度投影，确定信号能量最大的角度所对应的位置时间图的斜率。
预设角度范围，指根据实际情况，选取的进行角度投影的角度范围。作为一种可选的实施方式，预设角度范围可以是360度，此时需进行全角度的角度投影。作为另一种可选的实施方式，根据得到的位置时间图的特性，来选择进行角度投影的角度范围。步骤11中得到的位置时间图的横轴指示时间、纵轴指示位置，如果振动只由振动激励起点处向远处传播，则振动传播的速度无限大时，在位置时间图上接近于与纵轴平行的直线，而振动传播的速度无限小时，在位置时间图上接近于与横轴平行的直线，此时预设角度范围是90度即可满足需求，而无需做全角度投影，进而提高介质粘弹性定量的效率。如果振动除了由振动激励起点处向远处传播，还可沿相反方向继续传播，则预设角度范围可以是180度。至于预设角度范围的实际起始点和终止点，在直角坐标系保持不变的情况下，与0度的起始点及逆时针或顺时针的旋转方向相关，可以按需要选取，只要保证预设角度范围即可。
各个角度，指在预设角度范围内进行角度投影的各个角度。具体角度的选取根据时间精度需求和计算速度需求来确定，时间精度要求越高时，角度选取的精度要求越高，计算速度需求越高时，角度选取的精度要求越低。例如，可以从0.0.1度到1度的范围内选择。
角度投影，指针对设定角度进行图像特征识别或提取，以确定信号能量最大的角度。
在步骤22中，根据斜率，得到介质的粘弹性参数。
粘弹性参数包括粘性参数和弹性参数中的至少一项。
位置时间图的斜率由单位时间的振动传播的距离决定，即介质中振动传播的速度。在均匀的介质中，振动传播的速度和介质的粘弹性有关，在得到位置时间图的斜率之后，就可以定量计算出介质的粘弹性参数，因此如何高效准确的获得上述斜率成为介质粘弹性定量的关键。本示例性实施例利用角度投影确定信号能量最大的角度，由于该信号能量最大的角度对应位置时间图的斜率，即相当于获得了位置时间图的斜率。这种方法无需从位置时间图中选取波峰、波谷或振动的某一相位作为特征点来计算位置时间图的斜率，不会受到噪声干扰且计算量小，是一种高效准确的进行介质粘弹性定量的方法。
由于振动在介质中传播时，遇介质的边缘或异物时会产生反射波，为提高后续处理的精度，如图3所示，在进行角度投影之前，还可以包括步骤21’，即对位置时间图中的反射波进行滤除。滤除的方式可以有很多种，方向滤波是其中一种实现方式。
作为一种可选的实施方式，通过角度投影来确定信号能量最大的角度，进而得到位置时间图的斜率，可以通过积分计算来实现。例如，对位置时间图进行沿预设角度范围内的各个角度的积分计算，当积分角度和振动传播方向一致时能量聚集，此时得到的积分值最大，因此将积分值最大的角度为位置时间图的斜率线的斜率角。根据得到的斜率角，再结合位置和时间信息，即可以得到位置时间图的斜率线的斜率。上述积分计算又称为Radon变换。
作为另一种可选的实施方式，通过角度投影来确定信号能量最大的角度，进而得到位置时间图的斜率，还可以通过计算灰度共生矩阵来实现。由于通过计算灰度共生矩阵可以获得图像纹理特征，图像纹理特征又可以体现信号能量的大小，因此可以用灰度共生矩阵来获取信号能量最大的角度的信息。基于上述原理，通过角度投影来确定信号能量最大的角度，进而得到位置时间图的斜率，可以通过灰度共生矩阵的计算来实现。例如，针对位置时间图，先沿预设角度范围内的各个角度计算灰度共生矩阵。然后，利用灰度共生矩阵，得到各个角度的图像纹理特征。接着，利用图像纹理特征，将信号能量最大的角度，确定为位置时间图的斜率线的斜率角。最后，利用斜率角确定斜率线的斜率。
在一示例性实施例中，图1所示的方法中还可以包括介质粘弹性定量的步骤。图4是根据一示例性实施例示出的介质粘弹性定量方法的流程图，该流程基于图1所示流程实现，包括如下步骤。
在步骤41中，对位置时间图进行图像分割。
在步骤42中，提取图像特征。
图像特征可以是中轴线、峰值、谷值和过零点中的至少一个。上述中轴线指的是位置时间图上图纹的骨架，上述过零点指的是斜率最大值点或二阶导的最大值点。提取出的图像特征中均包含有振动的信息。
在步骤43中，利用图像特征进行线性拟合，得到位置时间图的斜率线的斜率。
在步骤44中，根据斜率，计算介质的粘弹性参数。
在前述两个介质粘弹性定量的示例性实施例中，根据力学原理，介质的粘弹性决定了振动在所述介质中的传播速度，因此通过得到位置时间图的斜率，可以获知振动在介质中传播的速度，进而根据力学原理，可以定量得出介质的粘弹性参数。这里的粘弹性参数可以包括剪切模量、杨氏模量、粘性模量、剪切粘弹性、剪切粘度、机械阻抗、机械松弛时间、各向异性等。
可选的，在使用线性拟合来实现介质粘弹性定量时，也可以先滤除位置时间图中的反射波，以达到更精确的定量效果。
下面以一种具体的应用场景给出本发明实施例中运动信息获取方法的应用。
在对人体肝脏等粘粘弹性介质进行无损粘弹性检测时，需要对介质粘弹性进行定量，定量之前就需要先获取运动信息。检测设备中包括激发装置和成像装置，其中激发装置对待检测介质进行振动激励，成像装置利用超声波对振动激励后的介质进行成像。振动在介质中传播时，在不同的时刻，波前会沿着传播方向到达不同的位置，形成位置时间图。上述波前可以是波峰、波谷、或振动的同一相位中的一种。
如图5所示，在这种具体应用场景的运动信息获取方法可以包括如下步骤。
在步骤51中，对介质进行振动激励。
在步骤52中，对在介质中传播的振动的超声检测信号进行频域变换，得到频域信号。
在步骤53中，对频域信号进行滤波或特征值选取，得到处理信号。
在步骤54中，根据设定的振动传播方向，利用处理信号得到振动的位置时间图。
在步骤55中，利用线性拟合或Radon变换得到位置时间图的斜率线的斜率。
在步骤56中，根据斜率，计算介质的粘弹性参数。
在以上运动信息获取方法的各个示例性实施例中，针对介质粘弹性定量的步骤，当设定的振动传播方向为至少两个时，每个设定的振动传播方向会对应得到一个位置时间图，进而会得到该位置时间图对应的介质的粘弹性参数。综合得到的至少两套粘弹性参数，可以更全面的评价介质的粘弹性。
上述给出的运动信息获取方法的各个示例性实施例，可以根据情况进行组合，这里并不限定各个示例性实施例之间的组合关系。
图6是一示例性实施例示出的运动信息获取装置的框图，该装置可以位于介质粘弹性检测设备的控制主机中，例如在医疗检测领域，可以位于肝脏无损检测设备的控制主机中。该装置还可以位于云端，介质粘弹性检测设备的检测数据需要在云端进行处理。
图6所示的装置包括：第一处理模块61、第二处理模块62和获取模块63。
第一处理模块61，用于对在介质中传播的振动的检测信号进行频域变换，得到频域信号。第一处理模块61在进行频域变换时，可以采用傅里叶变换、奇异值分解等多种方法。
第二处理模块62，用于去除频域信号中位于设定振动速度范围之外的信号，得到处理信号。
获取模块63，用于利用处理信号得到振动的位置时间图。
在一示例性实施例中，第二处理模块62对频域信号进行滤波或特征值选取，得到处理信号。其中，滤波的参数与设定振动速度范围相关，特征值选取与设定振动速度范围相关。
在一示例性实施例中，获取模块63根据设定的振动传播方向，利用处理信号得到振动的位置时间图。
在一示例性实施例中，如图7所示，图6所示的装置中还可以包括：粘弹性定量模块64，用于对位置时间图进行图像分割；提取图像特征；利用图像特征进行线性拟合，得到位置时间图的斜率线的斜率；根据斜率，计算介质的粘弹性参数。此处图像特征的含义与前文方法中所述的相同。
作为另一种可选的实施方式，粘弹性定量模块64还可以使用角度投影来实现相同功能。粘弹性定量模块64，用于沿位置时间图的预设角度范围内的各个角度进行角度投影，确定信号能量最大的角度所对应的位置时间图的斜率；根据斜率，得到介质的粘弹性参数。
进一步可选的，如图8所示，粘弹性定量模块64可以包括：计算子模块641和确定子模块642。
计算子模块641，用于对位置时间图进行沿预设角度范围内的各个角度的积分计算。
确定子模块642，用于将计算子模块641计算出的积分最大值的角度，确定为位置时间图的斜率线的斜率角。利用斜率角确定所述斜率线的斜率。
定量子模块643，用于根据斜率，得到介质的粘弹性参数。
作为另一种可选的实施方式，除了积分计算的方式之外，粘弹性定量模块64还可以通过计算灰度共生矩阵来确定斜率。此时，计算子模块641用于针对所述位置时间图，沿预设角度范围内的各个角度计算灰度共生矩阵。确定子模块642用于得到各个角度的图像纹理特征；利用图像纹理特征，将信号能量最大的角度，确定为位置时间图的斜率线的斜率角；利用斜率角确定斜率线的斜率。
进一步可选的，如图9所示，运动信息获取装置中还包括：滤波模块65，用于在粘弹性定量模块64进行角度投影之前，滤除位置时间图中的反射波。当然，在采用线性拟合来实现介质粘弹性定量时，也可以先由滤波模块65来滤除位置时间图中的反射波。
图10是根据一示例性实施例示出的运动信息获取装置的框图，该装置可以位于介质粘弹性检测设备的控制主机中，例如在医疗检测领域，可以位于肝脏无损检测设备的控制主机中。该装置还可以位于云端，介质粘弹性检测设备的检测数据需要在云端进行处理。
图10所示的装置包括：存储器101和处理器102。
存储器101中存储有执行指令。
处理器102，被配置于读取存储器101中的执行指令，执行前文所述运动信息获取方法的各示例性实施例中的部分或全部步骤。处理器102可以由芯片实现。
如果图10所示的运动信息获取装置位于介质粘弹性检测设备的控制主机中，可以通过总线、无线等方式与介质粘弹性定量设备中的激励装置、成像装置进行连接，此时该装置上具备与实现上述连接的接口及相应通信机制。
如果图10所示的运动信息获取装置位于云端，可以通过网络与介质粘弹性检测设备进行通信。
应当理解的是，本发明并不局限于上面已经描述并在附图中示出的流程及结构，并且可以在不脱离其范围进行各种修改和改变。本发明的范围仅由所附的权利要求来限制。

Claims (11)

1.一种运动信息获取方法，其特征在于，所述方法包括：
对在介质中传播的振动的检测信号进行频域变换，得到频域信号；
去除所述频域信号中位于设定振动速度范围之外的信号，得到处理信号；
利用所述处理信号得到所述振动的不以位移或应变为特征的位置时间图；
其中，所述去除所述频域信号中位于设定振动速度范围之外的信号，得到处理信号，包括：
对所述频域信号进行滤波或特征值选取，得到处理信号；
所述滤波的参数与所述设定振动速度范围相关；
所述特征值选取与所述设定振动速度范围相关。
2.如权利要求1所述的方法，其特征在于，所述利用所述处理信号得到所述振动的不以位移或应变为特征的位置时间图，包括：
根据设定的振动传播方向，利用所述处理信号得到所述振动的不以位移或应变为特征的位置时间图。
3.如权利要求1所述的方法，其特征在于，所述方法还包括：
对所述位置时间图进行图像分割；
提取图像特征；
利用所述图像特征进行线性拟合，得到所述位置时间图的斜率线的斜率；
根据所述斜率，计算所述介质的粘弹性参数。
4.如权利要求1所述的方法，其特征在于，所述方法还包括：
沿所述位置时间图的预设角度范围内的各个角度进行角度投影，确定信号能量最大的角度所对应的所述位置时间图的斜率；
根据所述斜率，得到所述介质的粘弹性参数。
5.如权利要求4所述的方法，其特征在于，所述沿所述位置时间图的预设角度范围内的各个角度进行角度投影，确定信号能量最大的角度所对应的所述位置时间图的斜率，包括：
对所述位置时间图进行沿预设角度范围内的各个角度的积分计算；
将积分值最大的角度，确定为所述位置时间图的斜率线的斜率角；
利用所述斜率角确定所述斜率线的斜率。
6.一种运动信息获取装置，其特征在于，所述装置包括：
第一处理模块，用于对在介质中传播的振动的检测信号进行频域变换，得到频域信号；
第二处理模块，用于去除所述频域信号中位于设定振动速度范围之外的信号，得到处理信号；
获取模块，用于利用所述处理信号得到所述振动的不以位移或应变为特征的位置时间图；
其中，所述第二处理模块对所述频域信号进行滤波或特征值选取，得到处理信号；
所述滤波的参数与所述设定振动速度范围相关；
所述特征值选取与所述设定振动速度范围相关。
7.如权利要求6所述的装置，其特征在于，所述获取模块根据设定的振动传播方向，利用所述处理信号得到所述振动的不以位移或应变为特征的位置时间图。
8.如权利要求6所述的装置，其特征在于，所述装置还包括：
粘弹性定量模块，用于对所述位置时间图进行图像分割；提取图像特征；利用所述图像特征进行线性拟合，得到所述位置时间图的斜率线的斜率；
根据所述斜率，计算所述介质的粘弹性参数。
9.如权利要求6所述的装置，其特征在于，所述装置还包括：
粘弹性定量模块，用于沿所述位置时间图的预设角度范围内的各个角度进行角度投影，确定信号能量最大的角度所对应的所述位置时间图的斜率；根据所述斜率，得到所述介质的粘弹性参数。
10.如权利要求9所述的装置，其特征在于，所述粘弹性定量模块，包括：
计算子模块，用于对所述位置时间图进行沿预设角度范围内的各个角度的积分计算；
确定子模块，用于将所述计算子模块计算出的积分最大值的角度，确定为所述位置时间图的斜率线的斜率角；利用所述斜率角确定所述斜率线的斜率；
定量子模块，用于根据所述斜率，得到所述介质的粘弹性参数。
11.一种运动信息获取装置，其特征在于，所述装置包括：
存储器，存储执行指令；
处理器，被配置于读取所述执行指令，完成如下操作：
用于对在介质中传播的振动的检测信号进行频域变换，得到频域信号；
去除所述频域信号中位于设定振动速度范围之外的信号，得到处理信号；
利用所述处理信号得到所述振动的不以位移或应变为特征的位置时间图；
其中，所述去除所述频域信号中位于设定振动速度范围之外的信号，得到处理信号，包括：
对所述频域信号进行滤波或特征值选取，得到处理信号；
所述滤波的参数与所述设定振动速度范围相关；
所述特征值选取与所述设定振动速度范围相关。
CN201710649554.8A 2017-07-21 2017-07-21 运动信息获取方法及装置 Active CN107505232B (zh)

Priority Applications (12)

Application Number Priority Date Filing Date Title
CN201710649554.8A CN107505232B (zh) 2017-07-21 2017-07-21 运动信息获取方法及装置
CA3065575A CA3065575C (en) 2017-07-21 2018-05-25 Method and apparatus for acquiring motion information
KR1020207004623A KR102235737B1 (ko) 2017-07-21 2018-05-25 운동 정보 획득 방법 및 장치
PCT/CN2018/088406 WO2019015398A1 (zh) 2017-07-21 2018-05-25 运动信息获取方法及装置
RU2020106886A RU2020106886A (ru) 2017-07-21 2018-05-25 Способ и устройство для получения информации о перемещении
EP18835507.7A EP3656311A4 (en) 2017-07-21 2018-05-25 METHOD AND DEVICE FOR DETECTING MOVEMENT INFORMATION
AU2018301989A AU2018301989A1 (en) 2017-07-21 2018-05-25 Method and device for acquiring motion information
BR112019025264-9A BR112019025264B1 (pt) 2017-07-21 2018-05-25 método e aparelho para adquirir informação de movimento
JP2020503053A JP7403163B2 (ja) 2017-07-21 2018-05-25 動き情報取得方法及び装置
TW107121597A TWI686180B (zh) 2017-07-21 2018-06-22 運動資訊獲取方法及裝置
US16/689,022 US11517289B2 (en) 2017-07-21 2019-11-19 Method and apparatus for acquiring motion information
AU2022200463A AU2022200463A1 (en) 2017-07-21 2022-01-25 Method and device for acquiring motion information

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710649554.8A CN107505232B (zh) 2017-07-21 2017-07-21 运动信息获取方法及装置

Publications (2)

Publication Number Publication Date
CN107505232A CN107505232A (zh) 2017-12-22
CN107505232B true CN107505232B (zh) 2019-09-03

Family

ID=60690286

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710649554.8A Active CN107505232B (zh) 2017-07-21 2017-07-21 运动信息获取方法及装置

Country Status (11)

Country Link
US (1) US11517289B2 (zh)
EP (1) EP3656311A4 (zh)
JP (1) JP7403163B2 (zh)
KR (1) KR102235737B1 (zh)
CN (1) CN107505232B (zh)
AU (2) AU2018301989A1 (zh)
BR (1) BR112019025264B1 (zh)
CA (1) CA3065575C (zh)
RU (1) RU2020106886A (zh)
TW (1) TWI686180B (zh)
WO (1) WO2019015398A1 (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107440740B (zh) 2017-07-21 2021-06-25 无锡海斯凯尔医学技术有限公司 介质粘弹性定量方法及装置
CN107505232B (zh) * 2017-07-21 2019-09-03 无锡海斯凯尔医学技术有限公司 运动信息获取方法及装置
CN113075120B (zh) * 2021-02-03 2022-08-30 湖南大学 一种实时土体类别识别方法、系统及一种盾构机

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101034004A (zh) * 2007-01-15 2007-09-12 国家海洋技术中心 视频波浪测量方法和测量系统
CN102641137A (zh) * 2011-02-17 2012-08-22 美国西门子医疗解决公司 使用幅度-相位调制超声波的粘弹性测量
CN103054552A (zh) * 2012-12-24 2013-04-24 深圳先进技术研究院 生物组织粘弹性测量方法和系统
CN103995053A (zh) * 2014-03-13 2014-08-20 天津大学 声表面波技术探测各向同性材料表层杨氏模量的检测方法
CN106037796A (zh) * 2015-04-16 2016-10-26 美国西门子医疗解决公司 量化的粘弹性超声成像
CN106419961A (zh) * 2015-08-11 2017-02-22 美国西门子医疗解决公司 声学辐射力成像中的自适应运动估计

Family Cites Families (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS55165587A (en) * 1979-06-11 1980-12-24 Hitachi Cable Method of forming rubber and plastic insulated cable connector
US5099848A (en) * 1990-11-02 1992-03-31 University Of Rochester Method and apparatus for breast imaging and tumor detection using modal vibration analysis
US5810731A (en) 1995-11-13 1998-09-22 Artann Laboratories Method and apparatus for elasticity imaging using remotely induced shear wave
EP1030592A4 (en) 1997-11-10 2004-04-07 Medacoustics Inc IMAGE GENERATION ARRANGEMENT FOR NON-INVASIVELY MEASURED TURBULENT BLOOD FLOWS
US20040034304A1 (en) * 2001-12-21 2004-02-19 Chikayoshi Sumi Displacement measurement method and apparatus, strain measurement method and apparatus elasticity and visco-elasticity constants measurement apparatus, and the elasticity and visco-elasticity constants measurement apparatus-based treatment apparatus
US8211019B2 (en) 2005-01-21 2012-07-03 Chikayoshi Sumi Clinical apparatuses
WO2006124603A2 (en) 2005-05-12 2006-11-23 The Trustees Of Columbia University In The City Of New York System and method for electromechanical wave imaging of body structures
US7464583B1 (en) 2005-06-10 2008-12-16 Carnegie Mellon University Methods and apparatuses using proximal probes
EP2310829B1 (en) * 2008-07-30 2016-04-06 Centre Hospitalier De L'Universite de Montreal A system and method for detection, characterization and imaging of heterogeneity using shear wave induced resonance
DE202008014792U1 (de) * 2008-11-07 2010-03-25 Qass Gmbh Vorrichtung zum Bewerten von Zerspanungsprozessen
CN101699280B (zh) 2009-10-15 2011-08-17 北京索瑞特医学技术有限公司 超声无损检测粘弹性介质弹性的方法及其装置
US9237878B2 (en) * 2011-04-22 2016-01-19 Mayo Foundation For Medical Education And Research Generation and assessment of shear waves in elasticity imaging
WO2014055410A1 (en) * 2012-10-01 2014-04-10 Mayo Foundation For Medical Education And Research Shear wave attenuation from k-space analysis system
KR101307526B1 (ko) * 2012-11-19 2013-09-12 (주)나다에스앤브이 실시간 진동 모니터링 장치
US9895137B2 (en) * 2013-03-05 2018-02-20 Hitachi, Ltd. Ultrasonic diagnosis device and transmission/reception method
EP4105640A1 (en) 2013-06-19 2022-12-21 The General Hospital Corporation Apparatus, devices and methods for obtaining omnidirectional viewing by a catheter
CN104640506B (zh) * 2013-09-18 2017-06-30 东芝医疗系统株式会社 超声波诊断装置、医用图像处理装置和医用图像处理方法
US10278671B2 (en) 2013-09-30 2019-05-07 Siemens Medical Solutions Usa, Inc. Shear wave detection in medical ultrasound imaging
JP6223169B2 (ja) * 2013-12-19 2017-11-01 キヤノン株式会社 情報処理装置、情報処理方法およびプログラム
US10722174B2 (en) 2014-07-11 2020-07-28 The Board Of Trustees Of The Leland Stanford Junior University Skin-conformal sensors
KR101620507B1 (ko) * 2015-10-07 2016-05-13 한국지질자원연구원 볼 밀 시스템의 진동 모니터링 장치 및 방법
CN105455851B (zh) * 2015-12-24 2018-03-13 无锡海斯凯尔医学技术有限公司 粘弹性介质的粘弹性参数检测方法和设备
CN107505232B (zh) * 2017-07-21 2019-09-03 无锡海斯凯尔医学技术有限公司 运动信息获取方法及装置
CN107505233A (zh) * 2017-07-21 2017-12-22 无锡海斯凯尔医学技术有限公司 介质粘弹性定量方法及装置
CN107440740B (zh) * 2017-07-21 2021-06-25 无锡海斯凯尔医学技术有限公司 介质粘弹性定量方法及装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101034004A (zh) * 2007-01-15 2007-09-12 国家海洋技术中心 视频波浪测量方法和测量系统
CN102641137A (zh) * 2011-02-17 2012-08-22 美国西门子医疗解决公司 使用幅度-相位调制超声波的粘弹性测量
CN103054552A (zh) * 2012-12-24 2013-04-24 深圳先进技术研究院 生物组织粘弹性测量方法和系统
CN103995053A (zh) * 2014-03-13 2014-08-20 天津大学 声表面波技术探测各向同性材料表层杨氏模量的检测方法
CN106037796A (zh) * 2015-04-16 2016-10-26 美国西门子医疗解决公司 量化的粘弹性超声成像
CN106419961A (zh) * 2015-08-11 2017-02-22 美国西门子医疗解决公司 声学辐射力成像中的自适应运动估计

Also Published As

Publication number Publication date
CA3065575C (en) 2023-08-01
BR112019025264B1 (pt) 2021-05-11
KR20200026302A (ko) 2020-03-10
TW201907869A (zh) 2019-03-01
US20200100769A1 (en) 2020-04-02
RU2020106886A (ru) 2021-08-23
BR112019025264A2 (pt) 2020-06-16
WO2019015398A1 (zh) 2019-01-24
EP3656311A4 (en) 2021-04-21
JP2020527429A (ja) 2020-09-10
TWI686180B (zh) 2020-03-01
KR102235737B1 (ko) 2021-04-01
CA3065575A1 (en) 2019-01-24
US11517289B2 (en) 2022-12-06
CN107505232A (zh) 2017-12-22
AU2018301989A1 (en) 2019-12-19
JP7403163B2 (ja) 2023-12-22
EP3656311A1 (en) 2020-05-27
AU2022200463A1 (en) 2022-02-17

Similar Documents

Publication Publication Date Title
CN102641137B (zh) 使用幅度-相位调制超声波的粘弹性测量
CN107505232B (zh) 运动信息获取方法及装置
US11719613B2 (en) Method and device for quantifying viscoelasticity of a medium
CN103857343A (zh) 使用薄体积的超声成像的弹性成像
Kijanka et al. Two-point frequency shift method for shear wave attenuation measurement
Kijanka et al. Robust phase velocity dispersion estimation of viscoelastic materials used for medical applications based on the multiple signal classification method
WO2018226688A1 (en) Estimating phase velocity dispersion in ultrasound elastography using a multiple signal classification
JP2018528013A5 (zh)
CN107505233A (zh) 介质粘弹性定量方法及装置
US12104999B2 (en) Method and device for quantifying viscoelasticity of a medium
Aristizabal et al. In-plane anisotropy method for the characterization of the elastic properties of anisotropic materials
Yuliuhina Synthesizing of brain MRE wave data
Lin et al. Determination of ultrasound phase velocity in trabecular bone using time dependent phase tracking technique
Zahiri-Azar et al. Soft Tissue Characterization Using One Dimensional Lumped Mass-Spring-Damper Model
Eriksson Dimension Estimation of an Impact Oscillator, an Experimental And Theoretical Approach

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
PE01 Entry into force of the registration of the contract for pledge of patent right
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: Method and device for obtaining motion information

Effective date of registration: 20230923

Granted publication date: 20190903

Pledgee: Bank of China Co.,Ltd. Wuxi High tech Industrial Development Zone Branch

Pledgor: WUXI HISKY MEDICAL TECHNOLOGIES Co.,Ltd.

Registration number: Y2023980058069