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

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

Info

Publication number
WO2019015398A1
WO2019015398A1 PCT/CN2018/088406 CN2018088406W WO2019015398A1 WO 2019015398 A1 WO2019015398 A1 WO 2019015398A1 CN 2018088406 W CN2018088406 W CN 2018088406W WO 2019015398 A1 WO2019015398 A1 WO 2019015398A1
Authority
WO
WIPO (PCT)
Prior art keywords
slope
signal
time map
position time
vibration
Prior art date
Application number
PCT/CN2018/088406
Other languages
English (en)
French (fr)
Inventor
何琼
邵金华
孙锦
段后利
王强
Original Assignee
无锡海斯凯尔医学技术有限公司
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 无锡海斯凯尔医学技术有限公司 filed Critical 无锡海斯凯尔医学技术有限公司
Priority to AU2018301989A priority Critical patent/AU2018301989A1/en
Priority to KR1020207004623A priority patent/KR102235737B1/ko
Priority to CA3065575A priority patent/CA3065575C/en
Priority to EP18835507.7A priority patent/EP3656311A4/en
Priority to RU2020106886A priority patent/RU2020106886A/ru
Priority to JP2020503053A priority patent/JP7403163B2/ja
Priority to BR112019025264-9A priority patent/BR112019025264B1/pt
Publication of WO2019015398A1 publication Critical patent/WO2019015398A1/zh
Priority to US16/689,022 priority patent/US11517289B2/en
Priority to AU2022200463A priority patent/AU2022200463A1/en

Links

Images

Classifications

    • 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
    • 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
    • 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

Definitions

  • the present invention relates to the field of measurement technologies, and in particular, to a method and apparatus for acquiring motion information.
  • the propagation characteristics of the vibration in the medium are related to the viscoelasticity of the medium.
  • the viscoelasticity of the medium can be quantified.
  • the lesions can be positioned by quantifying the viscoelasticity of the medium.
  • Embodiments of the present invention provide a method and an apparatus for acquiring motion information.
  • a brief summary is given below. This generalization is not a general comment, nor is it intended to identify key/critical constituent elements or to describe the scope of protection of these embodiments. Its sole purpose is to present some concepts in a simplified form as a prelude to the following detailed description.
  • a method for acquiring motion information comprising:
  • Frequency domain transforming the detection signal of the vibration propagating in the medium to obtain a frequency domain signal
  • a position time map of the vibration is obtained using the processed signal.
  • the removing the signal in the frequency domain signal that is outside the set vibration speed range, and obtaining the processed signal includes:
  • the filtered parameter is related to the set vibration speed range
  • the feature value selection is related to the set vibration speed range.
  • the obtaining a position time map of the vibration by using the processing signal includes:
  • a position time map of the vibration is obtained by the processing signal according to the set vibration propagation direction.
  • the method further includes:
  • the viscoelastic parameters of the medium are calculated.
  • the method further includes:
  • the angle projection is performed at each angle within a preset angle range of the position time map, and the position time map corresponding to the angle at which the signal energy is maximum is determined.
  • the slope including:
  • the slope of the slope line is determined using the slope angle.
  • an apparatus for acquiring motion information comprising:
  • a first processing module configured to perform frequency domain transformation on the detection signal of the vibration propagating in the medium to obtain a frequency domain signal
  • a second processing module configured to remove a signal in the frequency domain signal that is outside a set vibration speed range, to obtain a processed signal
  • an acquisition module configured to obtain a position time map of the vibration by using the processing signal.
  • the second processing module performs filtering or feature value selection on the frequency domain signal to obtain a processed signal.
  • the filtered parameter is related to the set vibration speed range
  • the feature value selection is related to the set vibration speed range.
  • the acquiring module obtains a position time map of the vibration by using the processing signal according to the set vibration propagation direction.
  • the device further includes:
  • a viscoelasticity quantification module configured to perform image segmentation on the position time map; extract image features; perform linear fitting using the image features to obtain a slope of a slope line of the position time map;
  • the viscoelastic parameters of the medium are calculated.
  • the apparatus further includes:
  • a viscoelasticity quantification module configured to perform an angular projection along each angle within a preset angle range of the position time map, and determine a slope of the position time map corresponding to an angle at which the signal energy is maximum; according to the slope, obtain a The viscoelastic parameters of the medium.
  • the viscoelasticity quantification module comprises:
  • a calculation submodule configured to perform integral calculation on each position angle within a preset angle range on the position time map
  • Determining a sub-module configured to determine an angle of an integral maximum value calculated by the calculation sub-module as a slope angle of a slope line of the position time map; determining a slope of the slope line by using the slope angle;
  • an apparatus for acquiring motion information comprising:
  • the processor is configured to read the execution instruction and complete the following operations:
  • Frequency domain transform for detecting a vibration of a vibration propagating in a medium to obtain a frequency domain signal
  • a position time map of the vibration is obtained using the processed signal.
  • Frequency domain transform is performed on the detection signal of the vibration propagating in the medium to obtain a frequency domain signal, and then the signal in the frequency domain signal outside the set vibration speed range is removed, that is, only the vibration signal is retained, and the vibration position is obtained.
  • Time map It is not necessary to perform motion estimation on the propagation of vibration through complicated calculations. It is only necessary to determine the presence or absence of vibration by processing in the frequency domain to obtain the position time map, which is an efficient way to acquire motion information.
  • FIG. 1 is a flowchart of a motion information acquisition method according to an exemplary embodiment
  • FIG. 2 is a flow chart showing a method for quantifying a medium viscoelasticity according to an exemplary embodiment
  • FIG. 3 is a flow chart showing a method of quantifying a medium viscoelasticity according to an exemplary embodiment
  • FIG. 4 is a flow chart showing a method for quantifying a medium viscoelasticity according to an exemplary embodiment
  • FIG. 5 is a flowchart of a motion information acquisition method according to an exemplary embodiment
  • FIG. 6 is a block diagram of an exercise information acquiring apparatus according to an exemplary embodiment
  • FIG. 7 is a block diagram of an exercise information acquiring apparatus according to an exemplary embodiment
  • Figure 8 is a block diagram of the viscoelasticity quantification module shown in Figure 7;
  • FIG. 9 is a block diagram of an exercise information acquiring apparatus according to an exemplary embodiment.
  • FIG. 10 is a block diagram of an exercise information acquiring apparatus according to an exemplary embodiment.
  • FIG. 1 is a flowchart of a motion information acquisition method according to an exemplary embodiment. As shown in FIG. 1, the method includes the following steps.
  • step 11 the detection signal of the vibration propagating in the medium is subjected to frequency domain transformation to obtain a frequency domain signal.
  • step 12 the signal outside the set vibration speed range of the frequency domain signal is removed to obtain a processed signal.
  • step 13 a positional time map of the vibration is obtained using the processed signal.
  • the vibration After the medium is excited by vibration, the vibration propagates in the medium. At different moments, the wavefront will reach different positions along the propagation direction. The correspondence between the position and the moment is the motion information of the vibration.
  • the commonly used motion information acquisition method utilizes a phase de-correlation of a detection signal for the vibration when the medium vibrates, and according to the phase de-correlation property, an algorithm such as cross-correlation, autocorrelation, and optical flow can be used. A position time map of the vibration is obtained.
  • the above algorithm can choose any method based on block matching. This commonly used method of acquiring motion information is equivalent to first estimating the motion of the vibration through the information such as the displacement and strain of the medium, and then obtaining the position time map of the vibration.
  • the detection signal is subjected to frequency domain transformation in the imaging time dimension to obtain a frequency domain signal, and then the frequency domain signal is removed and set in the frequency domain signal.
  • a signal outside the vibration speed range that is, a signal with a relatively low static or vibration speed is removed, and a positional time map of the vibration is obtained.
  • the motion information acquisition method in the exemplary embodiment does not require complicated calculation, and by processing in the frequency domain, a position time map not characterized by displacement or strain can be obtained. This method is equivalent to no need to first estimate the motion of the vibration, and it is only necessary to determine the presence or absence of the vibration to obtain the position time map, which is an efficient method for acquiring motion information.
  • the frequency domain transform may be performed in various manners, for example, using a Fourier transform or a singular value decomposition to perform frequency domain transform.
  • the signal outside the set vibration speed range in the frequency domain signal is removed, and the processed signal is obtained, which can be implemented by filtering the frequency domain signal or selecting the feature value.
  • the filter parameters of the filter in order to retain only the signal whose vibration speed is within the set range, consider the sampling rate of the signal in space and time, combined with the set vibration speed range, for example, 0.1 m to 30 In meters per second, to set the filter parameters of the filter, the signal in the frequency domain signal outside the set vibration speed range can be removed based on the filtering.
  • the selection of the feature value can also be set in relation to the set vibration velocity range, thereby removing the signal in the frequency domain signal that is outside the set vibration velocity range.
  • step 12 only the vibration signal can be retained, and the accuracy of the subsequent formation of the position time map can be improved.
  • the medium vibrates and the vibration propagates in the medium. Since the above-mentioned vibration has a limited propagation velocity in the medium, the medium can be dynamically imaged using the detection signal.
  • the detection wave may be a light wave, an ultrasonic wave or the like.
  • the above dynamic imaging may be one-dimensional imaging, two-dimensional imaging, three-dimensional imaging, or the like. Regardless of the imaging mode, the position time map of the vibration can be obtained by the processing signal after the removal operation according to the set vibration propagation direction.
  • the vibration propagation direction set as described above is the actual propagation direction of the vibration when the vibration propagates only in one propagation direction, and is a selected one of the propagation directions when the vibration propagates in a plurality of propagation directions.
  • the medium is a uniform sheet
  • the vibration will propagate along the extending direction of the sheet
  • the set vibration propagation direction is the actual propagation direction of the vibration.
  • the wavefront of the vibration propagation is a three-dimensional shape.
  • the wavefront of the vibration propagation is an ellipsoid, and the position time map obtained along different vibration propagation directions is different.
  • the direction of the vibration propagation is the selected one of the propagation directions of interest.
  • the above-described propagation direction of interest is determined according to the direction actually to be measured, and may be, for example, at least one of the direction in which the vibration propagation is the fastest, the direction in which the vibration propagation is the slowest, and the direction in which the vibration propagation velocity is a certain interval.
  • the method of Figure 1 may also include the step of quantifying the viscoelasticity of the medium.
  • 2 is a flow chart showing a method for quantifying a medium viscoelasticity, which is implemented based on the flow shown in FIG. 1 and includes the following steps, according to an exemplary embodiment.
  • step 21 the angle projection is performed along each angle within the preset angle range on the position time map to determine the slope of the position time map corresponding to the angle at which the signal energy is maximum.
  • the preset angle range refers to the range of angles at which the angle projection is selected according to the actual situation.
  • the preset angle range may be 360 degrees, and an angle projection of the full angle is required.
  • the angular range of the angular projection is selected according to the characteristics of the obtained position time map.
  • the horizontal axis of the position time map obtained in step 11 indicates the time and the vertical axis indicates the position. If the vibration propagates only from the vibration excitation starting point to the far side, the vibration propagation speed is infinite, and the position and time map are close to and vertical. A straight line parallel to the axis, and the speed of vibration propagation is infinitely small. It is close to a straight line parallel to the horizontal axis on the position time map.
  • the preset angle range is 90 degrees to meet the demand, without having to make a full-angle projection, thereby improving The efficiency of the quasi-elasticity of the medium.
  • the preset angle range may be 180 degrees.
  • the actual starting point and ending point of the preset angle range in the case where the Cartesian coordinate system remains unchanged, it is related to the starting point of 0 degree and the direction of rotation counterclockwise or clockwise, and can be selected as needed, as long as the preset is guaranteed.
  • the angle range is OK.
  • Each angle refers to each angle of angle projection within a preset angle range.
  • the selection of the specific angle is determined according to the time precision requirement and the calculation speed requirement.
  • the time precision requirement is higher, the accuracy requirement of the angle selection is higher, and the higher the calculation speed requirement, the lower the accuracy requirement of the angle selection. For example, it can be selected from the range of 0.0.1 degrees to 1 degree.
  • the pointer performs image feature recognition or extraction on the set angle to determine the angle at which the signal energy is maximum.
  • step 22 the viscoelastic parameters of the medium are obtained based on the slope.
  • the viscoelastic parameter includes at least one of a viscous parameter and an elastic parameter.
  • the slope of the position time map is determined by the distance traveled by the vibration per unit time, ie the speed at which the vibration propagates in the medium. In a uniform medium, the speed of vibration propagation is related to the viscoelasticity of the medium. After obtaining the slope of the position time map, the viscoelastic parameters of the medium can be quantitatively calculated. Therefore, how to obtain the above-mentioned slope efficiently and accurately becomes the viscoelasticity of the medium. key.
  • the present exemplary embodiment utilizes angular projection to determine the angle at which the signal energy is greatest, since the angle at which the signal energy is greatest corresponds to the slope of the position time map, that is, the slope of the position time map is obtained.
  • This method does not need to select a certain peak of the peak, trough or vibration as the feature point from the position time map to calculate the slope of the position time map. It is not subject to noise interference and the calculation amount is small, which is an efficient and accurate medium viscoelasticity. Quantitative method.
  • the step 21' may be included, that is, the position time
  • the reflected waves in the figure are filtered out.
  • filter There are many ways to filter, and directional filtering is one of the implementations.
  • the angle of the signal energy is determined by the angle projection to obtain the slope of the position time map, which can be realized by integral calculation.
  • the position time map is calculated by integrating the angles along the preset angle range. When the integration angle and the vibration propagation direction are the same, the energy is concentrated, and the integrated value obtained at this time is the largest, so the angle at which the integral value is the largest is the position time.
  • the slope angle of the slope line of the graph Based on the obtained slope angle, combined with the position and time information, the slope of the slope line of the position time map can be obtained.
  • the above integral calculation is also called Radon transform.
  • the angle of the signal energy is determined by the angle projection to obtain the slope of the position time map, and can also be realized by calculating the gray level co-occurrence matrix. Since the image texture feature can be obtained by calculating the gray level co-occurrence matrix, the image texture feature can reflect the signal energy, so the gray level co-occurrence matrix can be used to obtain the information of the angle of the signal energy. Based on the above principle, the angle at which the signal energy is maximum is determined by angular projection, and the slope of the position time map is obtained, which can be realized by calculation of the gray level co-occurrence matrix. For example, for a position time map, the gray level co-occurrence matrix is first calculated along various angles within a preset angle range.
  • image texture features of various angles are obtained.
  • the angle at which the signal energy is maximized is determined as the slope angle of the slope line of the position time map.
  • the slope angle is used to determine the slope of the slope line.
  • the method of Figure 1 may also include the step of quantifying the viscoelasticity of the medium.
  • FIG. 4 is a flow chart showing a method for quantifying a medium viscoelasticity, which is implemented based on the flow shown in FIG. 1 and includes the following steps, according to an exemplary embodiment.
  • step 41 image segmentation is performed on the position time map.
  • step 42 image features are extracted.
  • the image feature may be at least one of a center axis, a peak, a valley, and a zero crossing.
  • the above-mentioned central axis refers to the skeleton of the pattern on the position time map, and the above-mentioned zero-crossing point refers to the maximum point of the slope or the maximum point of the second-order guide.
  • the extracted image features contain vibration information.
  • step 43 a linear fit is performed using the image features to obtain the slope of the slope line of the position time map.
  • step 44 the viscoelastic parameters of the medium are calculated based on the slope.
  • the viscoelasticity of the medium determines the propagation speed of the vibration in the medium, so by obtaining the slope of the position time map, it can be known that the vibration is in the medium.
  • the speed of propagation and then according to the mechanics principle, can be used to quantitatively determine the viscoelastic parameters of the medium.
  • the viscoelastic parameters herein may include shear modulus, Young's modulus, viscous modulus, shear viscoelasticity, shear viscosity, mechanical resistance, mechanical relaxation time, anisotropy, and the like.
  • the reflected waves in the position time map can also be filtered out to achieve a more accurate quantitative effect.
  • the application of the motion information acquisition method in the embodiment of the present invention is given in a specific application scenario.
  • the detecting device includes an excitation device and an imaging device, wherein the excitation device vibrates the detection medium, and the imaging device images the vibration-excited medium by using the ultrasonic wave.
  • the wavefront will reach different positions along the propagation direction, forming a position time map.
  • the aforementioned wavefront may be one of the same phase of a crest, a trough, or a vibration.
  • the motion information acquisition method in this specific application scenario may include the following steps.
  • step 51 the medium is subjected to vibration excitation.
  • step 52 the ultrasonic detection signal of the vibration propagating in the medium is subjected to frequency domain transformation to obtain a frequency domain signal.
  • step 53 the frequency domain signal is filtered or the feature value is selected to obtain a processed signal.
  • step 54 a positional time map of the vibration is obtained from the processed signal based on the set vibration propagation direction.
  • step 55 the slope of the slope line of the position time plot is obtained using a linear fit or a Radon transform.
  • step 56 the viscoelastic parameters of the medium are calculated based on the slope.
  • each set vibration propagation direction corresponds to a position time map. Then, the viscoelastic parameters of the medium corresponding to the position time map are obtained.
  • the comprehensive set of at least two sets of viscoelastic parameters can more fully evaluate the viscoelasticity of the medium.
  • FIG. 6 is a block diagram of an exercise information acquiring apparatus shown in an exemplary embodiment, which may be located in a control host of a medium viscoelasticity detecting apparatus, for example, in the field of medical detection, may be located in a control host of a liver non-destructive detecting apparatus.
  • the device can also be located in the cloud, and the detection data of the medium viscoelasticity detecting device needs to be processed in the cloud.
  • the apparatus shown in FIG. 6 includes a first processing module 61, a second processing module 62, and an acquisition module 63.
  • the first processing module 61 is configured to perform frequency domain transformation on the detection signal of the vibration propagating in the medium to obtain a frequency domain signal.
  • the first processing module 61 may perform various methods such as Fourier transform and singular value decomposition when performing frequency domain transform.
  • the second processing module 62 is configured to remove signals in the frequency domain signal that are outside the set vibration speed range to obtain a processed signal.
  • the obtaining module 63 is configured to obtain a position time map of the vibration by using the processing signal.
  • the second processing module 62 filters or extracts the frequency domain signal to obtain a processed signal.
  • the filtered parameters are related to the set vibration speed range, and the feature value selection is related to the set vibration speed range.
  • the acquisition module 63 obtains a positional time map of the vibration using the processed signal according to the set vibration propagation direction.
  • the apparatus shown in FIG. 6 may further include: a viscoelasticity quantification module 64 for performing image segmentation on the position time map; extracting image features; and linearizing the image features Fit, to obtain the slope of the slope line of the position time map; calculate the viscoelastic parameters of the medium according to the slope.
  • a viscoelasticity quantification module 64 for performing image segmentation on the position time map; extracting image features; and linearizing the image features Fit, to obtain the slope of the slope line of the position time map; calculate the viscoelastic parameters of the medium according to the slope.
  • the viscoelastic quantitative module 64 can also use angular projection to achieve the same function.
  • the viscoelasticity quantification module 64 is configured to perform angular projection along various angles within a preset angle range of the position time map to determine a slope of a position time map corresponding to an angle at which the signal energy is maximum; and obtain a viscoelastic parameter of the medium according to the slope.
  • the viscoelasticity quantification module 64 may include a calculation sub-module 641 and a determination sub-module 642 .
  • the calculation sub-module 641 is configured to perform integral calculation on the position time map at various angles within a preset angle range.
  • the determining sub-module 642 is configured to determine the angle of the integral maximum value calculated by the calculation sub-module 641 as the slope angle of the slope line of the position time map.
  • the slope of the slope line is determined using the slope angle.
  • the quantitative sub-module 643 is configured to obtain a viscoelastic parameter of the medium according to the slope.
  • the viscoelasticity quantification module 64 may also determine the slope by calculating a gray level co-occurrence matrix.
  • the calculation sub-module 641 is configured to calculate a gray level co-occurrence matrix along each angle within a preset angle range for the position time map.
  • the determining sub-module 642 is configured to obtain image texture features at various angles; using the image texture feature, the angle at which the signal energy is maximized is determined as the slope angle of the slope line of the position time map; and the slope of the slope line is determined using the slope angle.
  • the motion information acquiring apparatus further includes: a filtering module 65, configured to filter the reflected wave in the position time map before the viscoelasticity quantitative module 64 performs angle projection.
  • a filtering module 65 configured to filter the reflected wave in the position time map before the viscoelasticity quantitative module 64 performs angle projection.
  • the filtering module 65 may first filter out the reflected waves in the position time map.
  • Figure 10 is a block diagram of an athletic information acquisition device, which may be located in a control host of a media viscoelasticity detection device, such as in the field of medical detection, may be located in a control host of a liver non-destructive testing device, in accordance with an exemplary embodiment.
  • the device can also be located in the cloud, and the detection data of the medium viscoelasticity detecting device needs to be processed in the cloud.
  • the apparatus shown in FIG. 10 includes a memory 101 and a processor 102.
  • An execution instruction is stored in the memory 101.
  • the processor 102 configured to read an execution instruction in the memory 101, performs some or all of the steps in the exemplary embodiments of the motion information acquisition method described above.
  • Processor 102 can be implemented by a chip.
  • the device can be connected to the excitation device and the imaging device in the medium viscoelasticity quantifying device by means of bus, wireless, etc., at this time, the device has An interface to implement the above connection and a corresponding communication mechanism.
  • the motion information acquiring device shown in FIG. 10 can communicate with the medium viscoelasticity detecting device through the network.

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Animal Behavior & Ethology (AREA)
  • Molecular Biology (AREA)
  • Public Health (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • Veterinary Medicine (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Analytical Chemistry (AREA)
  • Physiology (AREA)
  • Chemical & Material Sciences (AREA)
  • Biochemistry (AREA)
  • Immunology (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)
  • Operations Research (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (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 (13)

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

Priority Applications (9)

Application Number Priority Date Filing Date Title
AU2018301989A AU2018301989A1 (en) 2017-07-21 2018-05-25 Method and device for acquiring motion information
KR1020207004623A KR102235737B1 (ko) 2017-07-21 2018-05-25 운동 정보 획득 방법 및 장치
CA3065575A CA3065575C (en) 2017-07-21 2018-05-25 Method and apparatus for acquiring motion information
EP18835507.7A EP3656311A4 (en) 2017-07-21 2018-05-25 METHOD AND DEVICE FOR DETECTING MOVEMENT INFORMATION
RU2020106886A RU2020106886A (ru) 2017-07-21 2018-05-25 Способ и устройство для получения информации о перемещении
JP2020503053A JP7403163B2 (ja) 2017-07-21 2018-05-25 動き情報取得方法及び装置
BR112019025264-9A BR112019025264B1 (pt) 2017-07-21 2018-05-25 método e aparelho para adquirir informação de movimento
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 (2)

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

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US16/689,022 Continuation US11517289B2 (en) 2017-07-21 2019-11-19 Method and apparatus for acquiring motion information

Publications (1)

Publication Number Publication Date
WO2019015398A1 true WO2019015398A1 (zh) 2019-01-24

Family

ID=60690286

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2018/088406 WO2019015398A1 (zh) 2017-07-21 2018-05-25 运动信息获取方法及装置

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
CN107505232B (zh) 2017-07-21 2019-09-03 无锡海斯凯尔医学技术有限公司 运动信息获取方法及装置
CN107440740B (zh) * 2017-07-21 2021-06-25 无锡海斯凯尔医学技术有限公司 介质粘弹性定量方法及装置
CN113075120B (zh) * 2021-02-03 2022-08-30 湖南大学 一种实时土体类别识别方法、***及一种盾构机

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5099848A (en) * 1990-11-02 1992-03-31 University Of Rochester Method and apparatus for breast imaging and tumor detection using modal vibration analysis
CN102641137A (zh) * 2011-02-17 2012-08-22 美国西门子医疗解决公司 使用幅度-相位调制超声波的粘弹性测量
WO2014055410A1 (en) * 2012-10-01 2014-04-10 Mayo Foundation For Medical Education And Research Shear wave attenuation from k-space analysis system
CN105455851A (zh) * 2015-12-24 2016-04-06 无锡海斯凯尔医学技术有限公司 粘弹性介质的粘弹性参数检测方法和设备
CN106037796A (zh) * 2015-04-16 2016-10-26 美国西门子医疗解决公司 量化的粘弹性超声成像
CN107440740A (zh) * 2017-07-21 2017-12-08 无锡海斯凯尔医学技术有限公司 介质粘弹性定量方法及装置
CN107505233A (zh) * 2017-07-21 2017-12-22 无锡海斯凯尔医学技术有限公司 介质粘弹性定量方法及装置
CN107505232A (zh) * 2017-07-21 2017-12-22 无锡海斯凯尔医学技术有限公司 运动信息获取方法及装置

Family Cites Families (23)

* 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
US5810731A (en) 1995-11-13 1998-09-22 Artann Laboratories Method and apparatus for elasticity imaging using remotely induced shear wave
WO1999023940A1 (en) 1997-11-10 1999-05-20 Medacoustics, Inc. Non-invasive turbulent blood flow imaging system
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
US8858441B2 (en) 2005-05-12 2014-10-14 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
CN100470215C (zh) * 2007-01-15 2009-03-18 国家海洋技术中心 视频波浪测量方法和测量***
WO2010012092A1 (en) * 2008-07-30 2010-02-04 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
KR101307526B1 (ko) * 2012-11-19 2013-09-12 (주)나다에스앤브이 실시간 진동 모니터링 장치
CN103054552B (zh) * 2012-12-24 2014-12-10 深圳先进技术研究院 生物组织粘弹性测量方法和***
WO2014136502A1 (ja) * 2013-03-05 2014-09-12 日立アロカメディカル株式会社 超音波診断装置、及び送受信方法
EP4105640A1 (en) 2013-06-19 2022-12-21 The General Hospital Corporation Apparatus, devices and methods for obtaining omnidirectional viewing by a catheter
WO2015040710A1 (ja) * 2013-09-18 2015-03-26 株式会社 東芝 超音波診断装置、医用画像処理装置および医用画像処理方法
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 キヤノン株式会社 情報処理装置、情報処理方法およびプログラム
CN103995053A (zh) * 2014-03-13 2014-08-20 天津大学 声表面波技术探测各向同性材料表层杨氏模量的检测方法
US10722174B2 (en) 2014-07-11 2020-07-28 The Board Of Trustees Of The Leland Stanford Junior University Skin-conformal sensors
US10582911B2 (en) * 2015-08-11 2020-03-10 Siemens Medical Solutions Usa, Inc. Adaptive motion estimation in acoustic radiation force imaging
KR101620507B1 (ko) * 2015-10-07 2016-05-13 한국지질자원연구원 볼 밀 시스템의 진동 모니터링 장치 및 방법

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5099848A (en) * 1990-11-02 1992-03-31 University Of Rochester Method and apparatus for breast imaging and tumor detection using modal vibration analysis
CN102641137A (zh) * 2011-02-17 2012-08-22 美国西门子医疗解决公司 使用幅度-相位调制超声波的粘弹性测量
WO2014055410A1 (en) * 2012-10-01 2014-04-10 Mayo Foundation For Medical Education And Research Shear wave attenuation from k-space analysis system
CN106037796A (zh) * 2015-04-16 2016-10-26 美国西门子医疗解决公司 量化的粘弹性超声成像
CN105455851A (zh) * 2015-12-24 2016-04-06 无锡海斯凯尔医学技术有限公司 粘弹性介质的粘弹性参数检测方法和设备
CN107440740A (zh) * 2017-07-21 2017-12-08 无锡海斯凯尔医学技术有限公司 介质粘弹性定量方法及装置
CN107505233A (zh) * 2017-07-21 2017-12-22 无锡海斯凯尔医学技术有限公司 介质粘弹性定量方法及装置
CN107505232A (zh) * 2017-07-21 2017-12-22 无锡海斯凯尔医学技术有限公司 运动信息获取方法及装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See also references of EP3656311A4 *

Also Published As

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

Similar Documents

Publication Publication Date Title
Palmeri et al. Guidelines for finite-element modeling of acoustic radiation force-induced shear wave propagation in tissue-mimicking media
US11719613B2 (en) Method and device for quantifying viscoelasticity of a medium
CN107049360A (zh) 剪切波弹性成像方法和装置
US11517289B2 (en) Method and apparatus for acquiring motion information
CN106618635A (zh) 剪切波弹性成像方法和装置
Rose et al. Mindlin plate theory for damage detection: imaging of flexural inhomogeneities
Kijanka et al. Robust phase velocity dispersion estimation of viscoelastic materials used for medical applications based on the multiple signal classification method
Jarvis et al. Application of the distributed point source method to rough surface scattering and ultrasonic wall thickness measurement
Kijanka et al. Two-point frequency shift method for shear wave attenuation measurement
WO2018226688A1 (en) Estimating phase velocity dispersion in ultrasound elastography using a multiple signal classification
Ghosh et al. Modified error in constitutive equations (MECE) approach for ultrasound elastography
Sabbadini et al. Tapering of the interventricular septum can affect ultrasound shear wave elastography: An in vitro and in silico study
Yengul et al. Application of a forward model of axisymmetric shear wave propagation in viscoelastic media to shear wave elastography
Quang-Huy et al. Two-dimensional viscoelastic imaging using an enhanced FDTD-AHI approach
CN116784876A (zh) 可穿戴实时监测的新型便携式剪切波超声检测***及方法

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 18835507

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 3065575

Country of ref document: CA

REG Reference to national code

Ref country code: BR

Ref legal event code: B01A

Ref document number: 112019025264

Country of ref document: BR

ENP Entry into the national phase

Ref document number: 2018301989

Country of ref document: AU

Date of ref document: 20180525

Kind code of ref document: A

ENP Entry into the national phase

Ref document number: 2020503053

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

ENP Entry into the national phase

Ref document number: 20207004623

Country of ref document: KR

Kind code of ref document: A

ENP Entry into the national phase

Ref document number: 2018835507

Country of ref document: EP

Effective date: 20200221

ENP Entry into the national phase

Ref document number: 112019025264

Country of ref document: BR

Kind code of ref document: A2

Effective date: 20191129