WO2016033752A1 - Elasticity measurement and detection method and system - Google Patents

Elasticity measurement and detection method and system Download PDF

Info

Publication number
WO2016033752A1
WO2016033752A1 PCT/CN2014/085821 CN2014085821W WO2016033752A1 WO 2016033752 A1 WO2016033752 A1 WO 2016033752A1 CN 2014085821 W CN2014085821 W CN 2014085821W WO 2016033752 A1 WO2016033752 A1 WO 2016033752A1
Authority
WO
WIPO (PCT)
Prior art keywords
data
calculation
time
echo
measurement
Prior art date
Application number
PCT/CN2014/085821
Other languages
French (fr)
Chinese (zh)
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 PCT/CN2014/085821 priority Critical patent/WO2016033752A1/en
Priority to CN201910741683.9A priority patent/CN110432926B/en
Priority to CN201480081082.5A priority patent/CN106572838B/en
Priority to CN201910741685.8A priority patent/CN110368031B/en
Publication of WO2016033752A1 publication Critical patent/WO2016033752A1/en

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52023Details of receivers
    • G01S7/52036Details of receivers using analysis of echo signal for target characterisation
    • G01S7/52042Details of receivers using analysis of echo signal for target characterisation determining elastic properties of the propagation medium or of the reflective target
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52019Details of transmitters
    • G01S7/5202Details of transmitters for pulse systems
    • G01S7/52022Details of transmitters for pulse systems using a sequence of pulses, at least one pulse manipulating the transmissivity or reflexivity of the medium

Definitions

  • the invention relates to the field of medical detection, and in particular to an elastic measurement and detection method and system. Background technique
  • Ultrasound elastography is used to reflect the elasticity or softness of tissue. It is used more and more in the auxiliary detection, benign and malignant discrimination, and prognosis evaluation of tissue cancer.
  • the conventional superacoustic imaging method mainly produces a certain deformation by pressing the tissue by the probe, and then calculates and images the parameters related to the tissue elasticity such as strain and strain rate, and indirectly reflects the difference in elasticity between different tissues in the region of interest.
  • the strain parameters and the elastic parameters is affected by the magnitude of the pressure, even the same tissue, or different tissues with the same elastic parameters, if the applied pressure is different, the strain generated is different, which may cause Misjudgment of the difference in elasticity.
  • the operator needs to control the pressure of each operation consistently and maintain the stability of the hook, which makes the image repeatability and stability challenge. Summary of the invention
  • An elastic measurement detection method comprising the following steps,
  • the measurement calculation result is displayed.
  • the elastic measurement detecting method further comprises the following steps of generating a shear wavefront in the target region,
  • the data calculation acquires the displacement data, and when the displacement data is acquired, the following steps are included,
  • the reference echo data uses the echo data of the tracking pulse at any time or the reference at any time. Pulse echo data.
  • the particle velocity at different times is obtained by data calculation, and the following steps are included when acquiring the particle velocity
  • the displacement data of the tissue of the target area After obtaining the displacement data of the tissue of the target area, the displacement data of the same position is taken at different times; the gradient is obtained along the time direction to obtain the particle velocity curve at different times at the position.
  • the propagation distance of the shear wave in a period of time is obtained by data calculation, and the following steps are included when the propagation distance is obtained for calculation.
  • Obtaining a displacement-transverse position curve or a particle velocity-transverse position curve in a period of time before and after the time corresponding to the peak in the target area, a displacement-lateral position curve or a particle velocity-horizontal position for each time in the time period The curve performs cross-correlation judgment to obtain the propagation distance of the shear wave in the time period.
  • ⁇ 7 represent the propagation distance of the shear wave for a different target region within a period of time 7 ⁇
  • c represents the propagation velocity of the shear wave.
  • An elastic measurement and detection system includes an ultrasonic probe, a control module, a signal processing module, a calculation module and a display device, wherein the ultrasonic probe is provided with a transceiver module, a transceiver module, a signal processing module, a calculation module and a display device of the ultrasonic probe Connected in sequence, the control module is connected to the transceiver module,
  • the transceiver module is configured to transmit a push pulse, a tracking pulse and a reference pulse, and receive echo data of the tracking pulse and the reference pulse;
  • the control module is configured to control the transceiver module to transmit a push pulse, a tracking pulse, and a reference pulse;
  • the signal processing module is configured to perform signal preprocessing on the echo data
  • the calculation module is configured to process and calculate a signal output by the beam synthesis
  • the display device is configured to display an elastic image generated by the calculation module.
  • the calculating module includes:
  • a particle velocity calculation unit for calculating a particle velocity of the target tissue
  • a propagation distance calculation unit for calculating a propagation distance of the shear wave
  • the Young's modulus ratio calculation unit is used to calculate the Young's modulus ratio Young's modulus ratio of the target tissue.
  • the method of the invention can calculate and acquire related parameters reflecting the elasticity of the tissue without the operator controlling the strength of pressing the probe, so that the repeatability and stability of the measurement are ensured.
  • FIG. 1 is a schematic flow chart of an elastic measurement detecting method provided by a first preferred embodiment
  • FIG. 2 to FIG. 4 are shear wave transmissions generated by the elastic measurement detecting method provided by the first preferred embodiment.
  • FIG. 5 is a schematic diagram of displacement calculation in the elastic measurement detecting method provided by the first preferred embodiment
  • FIG. 6 is a schematic diagram of a propagation large off-time curve in different target regions in the elastic measurement detecting method provided by the first preferred embodiment
  • FIG. 7 is a schematic flow chart of an elastic measurement detecting method provided by a second preferred embodiment
  • Fig. 8 is a schematic view showing the configuration of an elastic measurement detecting system provided by the present invention. detailed description
  • a first preferred embodiment of the present invention provides an elastic measurement and detection method, including the following steps:
  • Step S101 generating a shear wave in the target area.
  • the target area can be determined according to the demand of the elastic measurement, and the target area is rearward to transmit a push pulse (push pulse) to the target area to generate a shear wave.
  • the target area can be selected and determined by any suitable method, such as a conventional two-dimensional B mode imaging, a conventional elastography E mode, and the like, and various types of applicable imaging detection methods are used for preliminary detection, and may also be selected according to detection requirements.
  • the number of target areas may be one or more. When the number of target areas is plural, the respective longitudinal depths or lateral positions of the plurality of target areas may be different. It can be understood that when measuring a plurality of target regions at the same time, an average distance ratio between the plurality of target regions can be obtained by a subsequent step, thereby reflecting a difference in elasticity between the plurality of target regions.
  • a specific pulse can be transmitted through the probe near the target area, ie, the artery is pushed.
  • the push pulse emission duration is longer than a conventional ultrasonic emission pulse, on the order of tens to hundreds of microseconds.
  • the push pulse may generate a shear wave source at the emission position, the shear wave is emitted from the shear wave source and propagates in a direction different from the direction in which the push pulse is emitted.
  • the driving pulse may be transmitted by focusing emission or non-focusing, and the driving pulse may also be used for transmitting once or continuously for multiple times, and may also be adjusted.
  • the depth of focus or/and the lateral position of the multiple successive shots are formed to form a specific shear wave propagation direction and broaden the propagation range of the shear wave.
  • each of the target areas may emit a plurality of sets of push pulses having the same or different emission focus modes, consecutive number of shots, and focus positions, each set of push pulses being transmitted for the respective target areas and in the target area Shear waves are generated nearby to analyze and calculate the echo data of the corresponding shear wave. Since the push pulse emission requires a large sound field energy to increase the shear wave intensity, it may take a certain time between multiple sets of push pulse emission to ensure energy safety, the intensity of each group of push pulses and the emission of each group of push pulses. The time interval is predetermined by the system.
  • Multiple target areas can also share the same set of push pulse transmissions, that is, after a set of push pulses are transmitted, shear wave propagation occurs, and shear wave propagation passes through multiple target areas, and then echoes of the plurality of target regions passing through are acquired. The data is analyzed and calculated. The system can determine whether a common push pulse emission is required depending on the depth of the target area and the distance of the lateral position.
  • each set of shots may include multiple consecutive shots, each of which may have a different focus position, longitudinal depth, and lateral position.
  • the direction of propagation and the propagation width of the shear wave generated by each of the shots can be controlled by adjusting the position of each shot and/or the interval between shots.
  • the calculation results of multiple sets of shear wave propagation in different directions or different positions can be combined and weighted to increase the accuracy of the results.
  • Step S102 transmitting a tracking pulse (tracking wave) to the target area and receiving its echo information, and acquiring echo data.
  • the echo data refers to echo data of a tracking pulse of a tracking pulse.
  • the tracking pulse is transmitted and its echo information is received, thereby obtaining echo data of the tracking pulse in a range of propagation within a target area for a period of time.
  • the transmission interval of the tracking pulse can be predetermined.
  • the echo data of the tracking pulse records the tissue information at each position within the above-mentioned propagation range during the shear wave propagation.
  • Step S103 Perform data calculation according to the echo data, and obtain a measurement calculation result.
  • the data can be calculated by tracking the echo data of the pulse, and the calculation results of the displacement, the particle velocity, the propagation distance of the shear wave at each moment in a period of time, and the propagation time of the shear wave within a certain distance can be obtained, and Show it.
  • the displacement calculation refers to calculating the tissue displacement of the target area and acquiring displacement data.
  • the displacement calculation can further include the following steps:
  • step S1031 echo data at any time in any position in the target area is taken as the reference echo data.
  • Step S1032 Segment the reference echo data, each segment is set as a reference echo core data (kernel), and set a reference echo core data center of the reference echo core data. It will be appreciated that the interval between the length of each of the reference echo core data and the reference echo core data center may be predetermined by the system.
  • Step S1033 Obtain a position where the cross-correlation between the echo data of each time and the reference echo core data center is the largest, and obtain the difference between the position with the largest cross-correlation and the center position of the reference echo core data. Displacement data for each moment.
  • the difference between the position of the echo data at each time and the position of the reference echo core data center and the center position of the reference echo core data is the reference echo core data position.
  • Corresponding displacement the echo data of the tracking pulse at different times in a certain horizontal position in the target area may be compared with the reference echo data, and the displacement data of the tissue at different depth positions may be obtained. This displacement reflects the relative longitudinal displacement of the tissue endoplasmic point relative to the above propagation moment caused by shear wave propagation.
  • the reference echo data and the echo data of the push pulse at different times may be cross-correlated by block-matching to obtain the displacement value of the position change.
  • the matching criteria of block matching algorithm and block matching displacement can also be set by themselves.
  • the segment data of the same position on the echo data can be phase-shifted, and the displacement value is directly calculated by the phase shift, and the position value is the displacement data. .
  • the above displacement calculation process can be performed for one or more locations within the target area.
  • a certain filtering process can be performed in the spatial or temporal dimension to reduce noise.
  • the particle velocity meter includes the following steps:
  • Step S1034 After acquiring the displacement data of the tissue of the target area, taking the displacement data of the same position at different times and obtaining the gradient along the time direction, the particle velocity curve at different times at the position can be obtained.
  • the particle velocity curve reflects the instantaneous velocity of the particle at various points in the shear wave propagation process at different times.
  • the above particle velocity calculation process can be performed for one or more locations within the target area.
  • a certain filtering process can be performed in the spatial or temporal dimension to reduce noise.
  • the propagation distance calculation includes the following steps:
  • Step S1035 Take displacement data of each time corresponding to any position in the target area, form a displacement-time curve, and find a time corresponding to the peak on the curve.
  • the time corresponding to the peak reflects the moment when the shear wave passes through the position. It can be understood that the time when the shear wave passes through the position is also the time corresponding to the peak of the particle velocity-time curve, and the particle velocity-time curve of the position can also be used to find the time when the corresponding shear wave passes the position.
  • the displacement or velocity data within a certain period of time ⁇ ⁇ near a specific time can be directly taken out, and the above propagation distance calculation is performed. The only difference is that the calculated distance is difficult to accurately determine which position is within the target area, but it can be judged after how long the displacement shear wave is transmitted.
  • Step S1036 Obtain a displacement-transverse position curve or a particle velocity-transverse position curve in a period of time ⁇ ⁇ before and after the peak corresponding to the peak in the target region, and a displacement-transverse position at each moment in the period ⁇ ⁇
  • the curve or the particle velocity-transverse position curve performs a cross-correlation judgment to obtain a difference in the lateral position of the shear wave in the time period, the difference reflecting the propagation distance of the shear wave within the above-mentioned period ⁇ ⁇ .
  • the above calculation processing is performed on all the positions in the target area, and the propagation distance of the shear wave within a period of time ⁇ ⁇ when the shear wave passes near each position in the target area can be obtained.
  • the length ⁇ ⁇ of the time period may be determined by itself, and the time period may be before or after the above time, before the above time or after the above time.
  • step S1035 the displacement or velocity data within a period of time ⁇ T near a certain specific time can be directly taken out, and the above-described propagation distance calculation is performed.
  • the average propagation distance within a period of time ⁇ ⁇ in the target region or the average propagation distance within a period of time ⁇ ⁇ during shear wave propagation can be obtained.
  • the propagation distances in ⁇ ⁇ corresponding to a plurality of consecutive moments are taken out. As shown in FIG. 6, the propagation distance-time curve can be obtained, or the total propagation distance within 7 ⁇ for a long period of time can be accumulated. It is also possible to take out the displacement data of each time corresponding to two different positions, and perform cross-correlation comparison to find out The time difference between the two can obtain the propagation time of the shear wave between these two positions.
  • the elastic ratio calculation may be included.
  • the elastic ratio mainly refers to a ratio of Young's modulus of each target region when the system has multiple target regions, to reflect each The degree of hardness difference between the target areas.
  • the propagation velocity of the shear wave has a nearly fixed relationship with the tissue hardness:
  • E ⁇ l E 2 c l I c 2 can be used for a period of time 7. Calculate the squared ratio of the shear wave propagation distance within ⁇ , ie: d x 2 ld 2 2
  • ⁇ 7 , ⁇ 7 represent the shear wave propagation distances of different target regions over a period of time T long .
  • T long the propagation velocity in each target region
  • the propagation speed calculation may be included.
  • the particle velocity-time curve or displacement-time curve of each lateral position is taken out, and the time shift between the curves is found by cross-correlation comparison, which represents the shear wave passing through the above-mentioned mutual participation.
  • the time difference between the two lateral positions of the correlation comparison represents the propagation speed at the corresponding position.
  • the approximate calculation formula can be directly derived by using the wave propagation equation as follows:
  • c is the propagation velocity and can be regarded as longitudinal displacement data or longitudinal velocity data
  • X represents the lateral coordinate
  • z represents the ⁇ slave coordinate
  • step S104 the measurement calculation result is displayed. After the above elastic measurement calculation, there are many ways The result is displayed.
  • the displacement data or the particle velocity data at a certain time is taken out to form a distribution map of the current time at each position of the target region, and the distribution map can reflect the propagation position of the shear wave at the current time.
  • the displacement data or the particle velocity data distribution map of multiple moments is continuously played, the shear wave wavefront propagating film map is formed, and the propagation process of the shear wave is visually represented.
  • the propagation distance distribution map is used to image the propagation distance of ⁇ ⁇ at various positions in the target area.
  • the propagation time distribution map For example, if the propagation time distribution map is used, the propagation time between a fixed distance near each position in the target area is imaged and displayed. Of course, the propagation time can also be displayed in a curved form.
  • the propagation distance curve For example, if the propagation distance curve is used, the propagation distance-time curve in the target area from a certain moment to a later period is drawn, and multiple curves can be displayed in multiple target areas.
  • the schematic diagram of the average propagation distance at each position in the target area is directly displayed, or the average propagation distance of the shear wave at a certain time in a certain period of time is directly displayed.
  • a schematic diagram of the elastic ratio in multiple target regions or a schematic diagram of the square of the propagation velocity, or a square of the propagation distance is displayed.
  • grayscale or color coding may be used, or the display effect may be enhanced by superimposing and blending with other modes of images.
  • a second preferred embodiment of the present invention provides an elastic measurement and detection method.
  • the elastic measurement and detection method of the present embodiment is substantially the same as the first preferred embodiment, and includes the following steps:
  • Step S201 transmitting a reference pulse (reference wave) to the target area, and receiving echo information of the reference pulse in the target area, acquiring echo data of the reference pulse, and taking echo data of the reference pulse as a reference echo. data.
  • the reference pulse is transmitted and the echo information in the target region is received before each group of push pulses is transmitted, and the echo data of the reference pulse is acquired and set as the reference echo data.
  • the echo information of the reference pulse records the tissue information of the target area before the shear wave is generated and starts to propagate.
  • Step S202 generating a shear wave in the target area.
  • Step S203 transmitting a tracking pulse to the target area and receiving its echo information, and acquiring echo data.
  • the echo data refers to echo data of a tracking pulse of a tracking pulse.
  • Step S204 performing data calculation according to the echo data, and obtaining a measurement calculation result.
  • the echo data of the reference pulse is taken as the reference echo data
  • the data is calculated by the echo data of the tracking pulse and the echo data of the reference pulse as the reference echo data, and the displacement, the particle velocity, and the time are obtained.
  • the echo data of the tracking pulses at different times in a certain horizontal position in the target region are compared with the reference echo data at the lateral position, and the displacement data of the tissue at different depth positions can be obtained.
  • This displacement mainly reflects the longitudinal displacement of the intragranular mass in the tissue caused by shear wave propagation relative to the intragranular mass before the shear wave propagates.
  • the calculation process such as the particle velocity, the propagation distance of the shear wave at each time in a period of time, and the propagation time of the shear wave within a certain distance are substantially the same in the first preferred embodiment, and will not be described herein.
  • Step S205 displaying the measurement calculation result.
  • the reference pulse before the shearing wave is generated in the target region, the reference pulse is first transmitted to the target region, and the echo data of the reference pulse is acquired as the reference echo data, thereby acquiring the intragranular particle relative to the shear caused by the shear wave propagation.
  • the calculation results such as the longitudinal displacement of the mass point in the tissue before the wave propagation, the corresponding particle velocity, the propagation distance of the shear wave at each moment in a period of time, and the propagation time of the shear wave in a distance.
  • the other steps in this embodiment are substantially the same as those in the first preferred embodiment, and are not described herein again.
  • the present invention further provides an ultrasonic elastic measuring system, which includes an ultrasonic probe 11 , a control module 12 , a signal processing module 13 , a calculation module 15 , and a display device 17 .
  • the ultrasonic probe 11 is provided with a transceiver module 110 .
  • the transceiver module 110, the signal processing module 13, the calculation module 15, and the display device 17 of the ultrasound probe 11 are sequentially connected, and the control module 12 is connected to the transceiver module 110. among them:
  • the transceiver module 110 is configured to transmit a push pulse, a tracking pulse, and a reference pulse, and receive echo data of the tracking pulse and the reference pulse.
  • the control module 12 is configured to control the transceiver module 110 to transmit a push pulse, a tracking pulse, and a reference pulse.
  • the control module 12 transmits a specific ultrasound sequence consisting of a push pulse, a reference pulse, and a tracking pulse according to a preset time interval, thereby facilitating generation of a shear wave in the target region and providing the transceiver module of the ultrasound probe 11 .
  • 110 receives the corresponding echo data.
  • the signal processing module 13 is configured to perform signal preprocessing on the echo data, thereby facilitating computational mode Block 15 performs subsequent calculations, which may include beamforming processing, and may also include, for example, signal amplification, analog to digital conversion, orthogonal decomposition, and the like.
  • the calculation module 15 is configured to perform processing calculation on the signal of the beam synthesis output.
  • the calculation module 15 can be used to calculate various types of propagation parameters of the shear wave.
  • the calculation module 15 includes:
  • a particle velocity calculation unit 151 configured to calculate a particle velocity of the target tissue
  • a propagation distance calculation unit 152 configured to calculate a propagation distance of the shear wave
  • the Young's modulus ratio calculation unit 153 is configured to calculate the Young's modulus ratio Young's modulus ratio of the target tissue.
  • the calculation module 15 calculates the shear wave propagation distance at each moment from a certain position, or the average distance within a certain period of time, or the average distance ratio of a certain period of time, or the average propagation speed in several target regions.
  • the ratio, etc. reflects the physical characteristics of the target area and generates a corresponding elastic image.
  • the display device 17 is for displaying the elasticity image generated by the calculation module 15.
  • the elastic measurement detecting method and system of the present invention generates a shear wave propagation by transmitting a series of specific pulses in a target region inside a tissue to be detected, and then continuously transmits and receives an echo signal within a range of the shear wave propagation path, from The echo distance of the shear wave is extracted from the echo signal for a period of time, and the average distance is calculated and finally displayed. If the shear wave propagates at the same time under the same emission conditions, the higher the hardness of the tissue, otherwise the softer the tissue.
  • the method of the invention can calculate and acquire related parameters reflecting the elasticity of the tissue without the operator controlling the strength of pressing the probe, so that the repeatability and stability of the measurement are ensured.
  • the present invention can obtain quantitative parameter values, and is particularly suitable for elastic measurement of uniform hook structures.
  • the present invention can also provide elastic difference parameters between different tissues to reflect the degree of elasticity difference between different tissues.

Abstract

An elasticity measurement and detection method and an elasticity measurement and detection system. The method comprises the following steps: generating a shear wave in a target region (S101); emitting a tracking pulse to the target region and acquiring echo data of the tracking pulse (S102); performing data calculation according to the echo data, so as to acquire a measurement calculation result (S103); and displaying the measurement calculation result (S104). In the method, an operator can acquire through calculation related parameters for reflecting the elasticity of a tissue, without the need of controlling the strength of pressing a probe, thereby ensuring the repeatability and stability of the measurement.

Description

弹性测量检测方法及系统 技术领域  Elastic measurement detection method and system
本发明涉及医疗检测领域, 尤其涉及一种弹性测量检测方法及系统。 背景技术  The invention relates to the field of medical detection, and in particular to an elastic measurement and detection method and system. Background technique
超声弹性成像用于反映组织 (Tissue ) 的弹性或软硬程度, 在组织癌症病 变的辅助检测、 良恶性判别、 预后恢复评价等方面得到越来越多应用。 常规超 声弹性成像方法主要通过探头按压组织产生一定的形变, 再将应变量、应变率 等与组织弹性相关的参数计算出来并成像,间接反映感兴趣区域中不同组织之 间的弹性差异。但是, 由于应变参数与弹性参数之间的对应关系受到压力大小 的影响, 即使是相同的组织, 或者具有相同的弹性参数的不同组织, 如果施加 的压力不同, 产生的应变也不同, 可能造成对弹性差异的误判。 而操作者使用 时需要人为控制每次操作的压力都保持一致,且维持均勾稳定, 这使得图像的 重复性和稳定性得到挑战。 发明内容  Ultrasound elastography is used to reflect the elasticity or softness of tissue. It is used more and more in the auxiliary detection, benign and malignant discrimination, and prognosis evaluation of tissue cancer. The conventional superacoustic imaging method mainly produces a certain deformation by pressing the tissue by the probe, and then calculates and images the parameters related to the tissue elasticity such as strain and strain rate, and indirectly reflects the difference in elasticity between different tissues in the region of interest. However, since the correspondence between the strain parameters and the elastic parameters is affected by the magnitude of the pressure, even the same tissue, or different tissues with the same elastic parameters, if the applied pressure is different, the strain generated is different, which may cause Misjudgment of the difference in elasticity. The operator needs to control the pressure of each operation consistently and maintain the stability of the hook, which makes the image repeatability and stability challenge. Summary of the invention
提供一种不受操作要求限制,具有稳定测量结果的弹性测量检测方法及系 统。  An elastic measurement detection method and system with stable measurement results that are not limited by operational requirements are provided.
一种弹性测量检测方法, 包括以下步骤,  An elastic measurement detection method, comprising the following steps,
于目标区域产生剪切波;  Producing shear waves in the target area;
向所述目标区域发射追踪脉冲并接收所述追踪脉冲的回波数据; 根据所述回波数据进行数据计算以获取测量计算结果;  Transmitting a tracking pulse to the target area and receiving echo data of the tracking pulse; performing data calculation according to the echo data to obtain a measurement calculation result;
显示所述测量计算结果。  The measurement calculation result is displayed.
进一步的,所述弹性测量检测方法于目标区域产生剪切波前还包括以下步 骤,  Further, the elastic measurement detecting method further comprises the following steps of generating a shear wavefront in the target region,
向目标区域发射参考脉冲,接收所述参考脉冲于所述目标区域内的回波信 息, 获取所述参考脉冲的回波数据。  And transmitting a reference pulse to the target area, receiving echo information of the reference pulse in the target area, and acquiring echo data of the reference pulse.
进一步的, 当根据所述回波数据进行数据计算以获取测量计算结果时, 通 过数据计算获取位移数据, 当获取位移数据时包括以下步骤, Further, when performing data calculation according to the echo data to obtain a measurement calculation result, The data calculation acquires the displacement data, and when the displacement data is acquired, the following steps are included,
取目标区域内任一位置的任一时刻的回波数据作为基准回波数据; 将基准回波数据分为多段基准回波核数据,设置所述基准回波核数据的基 准回波核数据中心;  Taking echo data at any time in any position in the target area as reference echo data; dividing the reference echo data into multi-segment reference echo core data, and setting a reference echo core data center of the reference echo core data ;
获取各个时刻的回波数据上与所述基准回波核数据中心之间互相关性最 大的位置,取该互相关性最大的位置与所述基准回波核数据的中心位置之差为 各个时刻的位移数据。  Obtaining a position where the cross-correlation between the echo data of each time and the reference echo core data center is the largest, and taking the difference between the position with the largest cross-correlation and the center position of the reference echo core data for each moment. Displacement data.
进一步的,当取目标区域内任一位置的任一时刻的回波数据作为基准回波 数据时,所述基准回波数据釆用任一时刻的追踪脉冲的回波数据或任一时刻的 参考脉冲的回波数据。  Further, when the echo data at any time in any position in the target area is taken as the reference echo data, the reference echo data uses the echo data of the tracking pulse at any time or the reference at any time. Pulse echo data.
进一步的, 当根据所述回波数据进行数据计算以获取测量计算结果时, 通 过数据计算获取不同时刻的质点速度, 当获取质点速度时包括以下步骤,  Further, when data calculation is performed according to the echo data to obtain a measurement calculation result, the particle velocity at different times is obtained by data calculation, and the following steps are included when acquiring the particle velocity,
获取目标区域的组织的位移数据后, 取同一个位置不同时刻的位移数据; 沿时间方向求梯度以获取该位置处不同时刻的质点速度曲线。  After obtaining the displacement data of the tissue of the target area, the displacement data of the same position is taken at different times; the gradient is obtained along the time direction to obtain the particle velocity curve at different times at the position.
进一步的, 当根据所述回波数据进行数据计算以获取测量计算结果时, 通 过数据计算获取一段时间段内的剪切波的传播距离,当获取所述传播距离进行 计算时包括以下步骤,  Further, when data calculation is performed according to the echo data to obtain a measurement calculation result, the propagation distance of the shear wave in a period of time is obtained by data calculation, and the following steps are included when the propagation distance is obtained for calculation.
取目标区域中任一位置对应的各个时刻的位移数据,形成位移 -时间曲线, 找出该曲线上峰值对应的时刻;  Taking displacement data of each moment corresponding to any position in the target area, forming a displacement-time curve, and finding a time corresponding to the peak on the curve;
获取目标区域内的所述峰值对应的时刻前后一段时间段内的位移-横向位 置曲线或者质点速度-横向位置曲线, 对所述时间段内各个时刻的位移 -横向位 置曲线或者质点速度 -横向位置曲线进行互相关判断, 获取所述时间段内的剪 切波的传播距离。  Obtaining a displacement-transverse position curve or a particle velocity-transverse position curve in a period of time before and after the time corresponding to the peak in the target area, a displacement-lateral position curve or a particle velocity-horizontal position for each time in the time period The curve performs cross-correlation judgment to obtain the propagation distance of the shear wave in the time period.
进一步的, 当根据所述回波数据进行数据计算以获取测量计算结果时, 获 取不同目标区域之间的杨氏模量比值, 当获取所述杨氏模量比值时, 满足以下 公式, Further, when data calculation is performed according to the echo data to obtain a measurement calculation result, a Young's modulus ratio between different target regions is obtained, and when the Young's modulus ratio is obtained, the following formula is satisfied.
Figure imgf000003_0001
dx 2 ld2 2
Figure imgf000003_0001
d x 2 ld 2 2
式中, di、 ί 7表示不同的目标区域在一段时间 7 ^内的所述剪切波的传播 距离, c表示所述剪切波的传播速度。 进一步的, 当显示所述测量计算结果时, 釆用波前电影图、 传播距离分布 图、 传播时间分布图、 传播距离曲线图、 传播距离示意图、 弹性比值示意图、 传播速度平方比示意图、 传播距离平方比示意图中的至少一种。 Where di, ί 7 represent the propagation distance of the shear wave for a different target region within a period of time 7 ^, and c represents the propagation velocity of the shear wave. Further, when the measurement calculation result is displayed, the wavefront film map, the propagation distance distribution map, the propagation time distribution map, the propagation distance curve diagram, the propagation distance diagram, the elastic ratio diagram, the squared ratio of the propagation velocity, and the propagation distance are used. At least one of the squared ratio diagrams.
进一步的, 当显示所述测量计算结果时, 釆用灰阶或彩色编码模式。 一种弹性测量检测系统, 包括超声探头、 控制模块、 信号处理模块、 计算 模块及显示装置, 所述超声探头设有收发模块, 所述超声探头的收发模块、 信 号处理模块、 计算模块、 显示装置依次连接, 所述控制模块连接于所述收发模 块,  Further, when the measurement calculation result is displayed, a gray scale or color coding mode is used. An elastic measurement and detection system includes an ultrasonic probe, a control module, a signal processing module, a calculation module and a display device, wherein the ultrasonic probe is provided with a transceiver module, a transceiver module, a signal processing module, a calculation module and a display device of the ultrasonic probe Connected in sequence, the control module is connected to the transceiver module,
所述收发模块用于发射推动脉冲、追踪脉冲及参考脉冲, 并接收所述追踪 脉冲与所述参考脉冲的回波数据;  The transceiver module is configured to transmit a push pulse, a tracking pulse and a reference pulse, and receive echo data of the tracking pulse and the reference pulse;
所述控制模块用于控制所述收发模块发射推动脉冲、 追踪脉冲及参考脉 冲;  The control module is configured to control the transceiver module to transmit a push pulse, a tracking pulse, and a reference pulse;
所述信号处理模块用于对回波数据进行信号预处理;  The signal processing module is configured to perform signal preprocessing on the echo data;
所述计算模块用于对波束合成输出的信号进行处理计算;  The calculation module is configured to process and calculate a signal output by the beam synthesis;
所述显示装置用于将所述计算模块生成的弹性图像进行显示。  The display device is configured to display an elastic image generated by the calculation module.
进一步的, 所述计算模块包括:  Further, the calculating module includes:
质点速度计算单元, 用于计算目标组织的质点速度;  a particle velocity calculation unit for calculating a particle velocity of the target tissue;
传播距离计算单元, 用于计算剪切波的传播距离;  a propagation distance calculation unit for calculating a propagation distance of the shear wave;
杨氏模量比值计算单元, 用于计算目标组织的杨氏模量比值杨氏模量比 值。  The Young's modulus ratio calculation unit is used to calculate the Young's modulus ratio Young's modulus ratio of the target tissue.
本发明的方法不需要操作者控制按压探头的力度即可对反映组织弹性的 相关参数进行计算获取, 使得测量的重复性和稳定性得到保证。  The method of the invention can calculate and acquire related parameters reflecting the elasticity of the tissue without the operator controlling the strength of pressing the probe, so that the repeatability and stability of the measurement are ensured.
附图说明 DRAWINGS
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施 例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地, 下面描述 中的附图仅仅是本发明的一些实施例, 对于本领域普通技术人员来讲, 在不付 出创造性劳动的前提下, 还可以根据这些附图获得其他的附图。  In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the embodiments or the description of the prior art will be briefly described below. Obviously, the drawings in the following description are only It is a certain embodiment of the present invention, and other drawings can be obtained from those skilled in the art without any creative work.
图 1是第一较佳实施方式提供的弹性测量检测方法流程示意图; 图 2至图 4是第一较佳实施方式提供的弹性测量检测方法中产生剪切波传 播的推动波发射示意图; 1 is a schematic flow chart of an elastic measurement detecting method provided by a first preferred embodiment; FIG. 2 to FIG. 4 are shear wave transmissions generated by the elastic measurement detecting method provided by the first preferred embodiment. a schematic diagram of the propagating wave emission of the broadcast;
图 5是第一较佳实施方式提供的弹性测量检测方法中的位移计算示意图; 图 6 是第一较佳实施方式提供的弹性测量检测方法中不同目标区域内的 传播巨离 -时间曲线示意图;  5 is a schematic diagram of displacement calculation in the elastic measurement detecting method provided by the first preferred embodiment; FIG. 6 is a schematic diagram of a propagation large off-time curve in different target regions in the elastic measurement detecting method provided by the first preferred embodiment;
图 7是第二较佳实施方式提供的弹性测量检测方法流程示意图;  7 is a schematic flow chart of an elastic measurement detecting method provided by a second preferred embodiment;
图 8是本发明提供的弹性测量检测系统的构成示意图。 具体实施方式  Fig. 8 is a schematic view showing the configuration of an elastic measurement detecting system provided by the present invention. detailed description
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清 楚、 完整地描述, 显然, 所描述的实施例仅仅是本发明一部分实施例, 而不是 全部的实施例。基于本发明中的实施例, 本领域普通技术人员在没有作出创造 性劳动前提下所获得的所有其他实施例, 都属于本发明保护的范围。  BRIEF DESCRIPTION OF THE DRAWINGS The technical solutions in the embodiments of the present invention will be described in detail below with reference to the accompanying drawings. All other embodiments obtained by a person of ordinary skill in the art based on the embodiments of the present invention without creative work are within the scope of the present invention.
请参阅图 1, 本发明第一较佳实施方式提供一种弹性测量检测方法, 包括 以下步骤:  Referring to FIG. 1, a first preferred embodiment of the present invention provides an elastic measurement and detection method, including the following steps:
步骤 S101, 于目标区域产生剪切波。 本实施例中, 目标区域可根据弹性 测量的需求确定, 目标区域后向目标区域发射推动脉冲 ( push脉冲) 以产生 剪切波。 所述目标区域可通过任意适用方式选择确定, 如可釆用常规二维 B 模式成像、常规弹性成像 E模式等各类适用成像检测方式初步检测后确定,也 可根据检测需求选定。  Step S101, generating a shear wave in the target area. In this embodiment, the target area can be determined according to the demand of the elastic measurement, and the target area is rearward to transmit a push pulse (push pulse) to the target area to generate a shear wave. The target area can be selected and determined by any suitable method, such as a conventional two-dimensional B mode imaging, a conventional elastography E mode, and the like, and various types of applicable imaging detection methods are used for preliminary detection, and may also be selected according to detection requirements.
所述目标区域的数量可以为一个或者多个。 当目标区域的数量为多个时, 多个目标区域的各自纵向深度或横向位置可以不同。可以理解的是, 当同时对 多个目标区域进行测量时,可通过后续步骤获取多个目标区域之间的平均距离 比值, 从而反映所述多个目标区域之间的弹性差异。  The number of target areas may be one or more. When the number of target areas is plural, the respective longitudinal depths or lateral positions of the plurality of target areas may be different. It can be understood that when measuring a plurality of target regions at the same time, an average distance ratio between the plurality of target regions can be obtained by a subsequent step, thereby reflecting a difference in elasticity between the plurality of target regions.
确定目标区域后, 即可通过探头向目标区域附近发射特定脉冲, 即推动脉 冲。所述推动脉冲发射时长比常规超声发射脉冲长,约为几十到几百微秒量级。 所述推动脉冲可以在发射位置产生剪切波源, 自剪切波源发射剪切波并沿着与 所述推动脉冲发射方向不同的方向传播。  Once the target area has been determined, a specific pulse can be transmitted through the probe near the target area, ie, the artery is pushed. The push pulse emission duration is longer than a conventional ultrasonic emission pulse, on the order of tens to hundreds of microseconds. The push pulse may generate a shear wave source at the emission position, the shear wave is emitted from the shear wave source and propagates in a direction different from the direction in which the push pulse is emitted.
如图 2至图 4所示,在本实施例中,发射所述推动脉冲可以釆用聚焦发射 或无聚焦发射, 所述推动脉冲也可釆用发射一次或者连续发射多次,还可以调 整连续多次发射的聚焦深度或 /和横向位置, 以形成特定的剪切波传播方向和 拓宽剪切波的传播范围。 As shown in FIG. 2 to FIG. 4, in the embodiment, the driving pulse may be transmitted by focusing emission or non-focusing, and the driving pulse may also be used for transmitting once or continuously for multiple times, and may also be adjusted. The depth of focus or/and the lateral position of the multiple successive shots are formed to form a specific shear wave propagation direction and broaden the propagation range of the shear wave.
当确认并设置多个目标区域时,各个目标区域可以发射具有相同或不同的 发射聚焦方式、 连续发射次数、 聚焦位置的多组推动脉冲, 每组推动脉冲针对 各自目标区域发射并在该目标区域附近产生剪切波,以便对相应剪切波的回波 数据进行分析计算。 由于推动脉冲发射需要较大的声场能量以提高剪切波强 度, 多组推动脉冲发射之间可能需要间隔一定时间以保证能量安全性,每组推 动脉冲的发射的强度和各组推动脉冲的发射的时间间隔由系统预先确定。  When a plurality of target areas are confirmed and set, each of the target areas may emit a plurality of sets of push pulses having the same or different emission focus modes, consecutive number of shots, and focus positions, each set of push pulses being transmitted for the respective target areas and in the target area Shear waves are generated nearby to analyze and calculate the echo data of the corresponding shear wave. Since the push pulse emission requires a large sound field energy to increase the shear wave intensity, it may take a certain time between multiple sets of push pulse emission to ensure energy safety, the intensity of each group of push pulses and the emission of each group of push pulses. The time interval is predetermined by the system.
多个目标区域也可以共用同一组推动脉冲发射, 即一组推动脉冲发射后, 产生剪切波传播, 剪切波传播会经过多个目标区域, 然后获取上述经过的多个 目标区域的回波数据进行分析计算。系统可以根据目标区域深度和横向位置的 远近来决定是否需要共用推动脉冲发射。  Multiple target areas can also share the same set of push pulse transmissions, that is, after a set of push pulses are transmitted, shear wave propagation occurs, and shear wave propagation passes through multiple target areas, and then echoes of the plurality of target regions passing through are acquired. The data is analyzed and calculated. The system can determine whether a common push pulse emission is required depending on the depth of the target area and the distance of the lateral position.
当向目标区域内发射多组脉冲时,每组发射可包括多次连续发射,每次发 射的脉冲可具有不同的聚焦位置、 纵向深度、横向位置。 可通过调节各次发射 的位置和 /或发射间隔时间, 来控制其产生的剪切波的传播方向、 传播宽度。 多组不同方向或者不同位置的剪切波传播后的计算结果, 可以进行复合、加权 等处理, 以增加结果准确性。  When multiple sets of pulses are transmitted into the target area, each set of shots may include multiple consecutive shots, each of which may have a different focus position, longitudinal depth, and lateral position. The direction of propagation and the propagation width of the shear wave generated by each of the shots can be controlled by adjusting the position of each shot and/or the interval between shots. The calculation results of multiple sets of shear wave propagation in different directions or different positions can be combined and weighted to increase the accuracy of the results.
步骤 S102, 向目标区域发射追踪脉冲 ( tracking波 ) 并接收其回波信息, 并获取回波数据。  Step S102, transmitting a tracking pulse (tracking wave) to the target area and receiving its echo information, and acquiring echo data.
在本步骤中, 所述回波数据指追踪脉冲的追踪脉冲的回波数据。 向目标区 域发射推动脉冲并产生剪切波后,发射追踪脉冲并接收其回波信息,从而得到 目标区域内一段时间内一段传播范围内的追踪脉冲的回波数据。追踪脉冲的发 射间隔时间可预先确定。 所述追踪脉冲的回波数据记录了剪切波传播过程中, 上述一段传播范围内各位置上的组织信息。  In this step, the echo data refers to echo data of a tracking pulse of a tracking pulse. After the push pulse is generated to the target area and the shear wave is generated, the tracking pulse is transmitted and its echo information is received, thereby obtaining echo data of the tracking pulse in a range of propagation within a target area for a period of time. The transmission interval of the tracking pulse can be predetermined. The echo data of the tracking pulse records the tissue information at each position within the above-mentioned propagation range during the shear wave propagation.
步骤 S103, 根据所述回波数据进行数据计算, 获取测量计算结果。 在本 实施例中, 可通过追踪脉冲的回波数据进行数据计算并获取位移、 质点速度、 一段时间内各时刻剪切波的传播距离、一段距离内剪切波的传播时间等计算结 果, 并加以显示。  Step S103: Perform data calculation according to the echo data, and obtain a measurement calculation result. In this embodiment, the data can be calculated by tracking the echo data of the pulse, and the calculation results of the displacement, the particle velocity, the propagation distance of the shear wave at each moment in a period of time, and the propagation time of the shear wave within a certain distance can be obtained, and Show it.
如图 5所示, 剪切波经过目标区域时, 目标区域的组织的相应位置会产生 纵向位移, 即位置变化。 所述位移计算即指对目标区域的组织位移进行计算并 获取位移数据。 位移计算可进一步包括以下步骤: As shown in Figure 5, when the shear wave passes through the target area, the corresponding position of the tissue of the target area is generated. Longitudinal displacement, ie positional change. The displacement calculation refers to calculating the tissue displacement of the target area and acquiring displacement data. The displacement calculation can further include the following steps:
步骤 S 1031,取目标区域内任一位置的任一时刻的回波数据作为基准回波 数据。  In step S1031, echo data at any time in any position in the target area is taken as the reference echo data.
步骤 S1032 , 将基准回波数据分段, 每段设为基准回波核数据 (kernel ), 设置所述基准回波核数据的基准回波核数据中心。 可以理解的是,每段所述基 准回波核数据的长度与所述基准回波核数据中心之间的间隔可由系统预先确 定。  Step S1032: Segment the reference echo data, each segment is set as a reference echo core data (kernel), and set a reference echo core data center of the reference echo core data. It will be appreciated that the interval between the length of each of the reference echo core data and the reference echo core data center may be predetermined by the system.
步骤 S1033 , 获取各个时刻的回波数据上与所述基准回波核数据中心之间 互相关性最大的位置,取该互相关性最大的位置与所述基准回波核数据的中心 位置之差为各个时刻的位移数据。  Step S1033: Obtain a position where the cross-correlation between the echo data of each time and the reference echo core data center is the largest, and obtain the difference between the position with the largest cross-correlation and the center position of the reference echo core data. Displacement data for each moment.
本步骤中,所述各个时刻的回波数据上与所述基准回波核数据中心互相关 性最大的位置与所述基准回波核数据的中心位置之差即为该基准回波核数据 位置对应的位移。 本步骤中, 可取目标区域内某一横向位置上不同时刻的追踪 脉冲的回波数据, 分别与该基准回波数据进行互相关比较, 可以得到各个深度 位置上的组织在不同时刻的位移数据。该位移反映了剪切波传播导致的组织内 质点相对于上述传播时刻的相对纵向位移。  In this step, the difference between the position of the echo data at each time and the position of the reference echo core data center and the center position of the reference echo core data is the reference echo core data position. Corresponding displacement. In this step, the echo data of the tracking pulse at different times in a certain horizontal position in the target area may be compared with the reference echo data, and the displacement data of the tissue at different depth positions may be obtained. This displacement reflects the relative longitudinal displacement of the tissue endoplasmic point relative to the above propagation moment caused by shear wave propagation.
在步骤 S1033中,基准回波数据与不同时刻的推动脉冲的回波数据可釆用 块匹配(block-matching )的方式进行互相关比较, 以获取位置变化的位移值。 可以理解的是, 块匹配的算法选择及块匹配位移的匹配标准也可自行设定。 并 且由于剪切波传播产生的位移非常小,也可以将回波数据上相同位置的分段数 据进行相移 (phase shift )计算, 直接通过相移计算出位移值, 该位置值即为 位移数据。  In step S1033, the reference echo data and the echo data of the push pulse at different times may be cross-correlated by block-matching to obtain the displacement value of the position change. It can be understood that the matching criteria of block matching algorithm and block matching displacement can also be set by themselves. Moreover, since the displacement caused by the shear wave propagation is very small, the segment data of the same position on the echo data can be phase-shifted, and the displacement value is directly calculated by the phase shift, and the position value is the displacement data. .
对目标区域内的一个或多个位置,都可以进行上述位移计算过程。对所得 位移数据, 可以在空间或时间维度进行一定的滤波处理以降低噪声。  The above displacement calculation process can be performed for one or more locations within the target area. For the resulting displacement data, a certain filtering process can be performed in the spatial or temporal dimension to reduce noise.
步骤 S103中根据所述回波数据进行数据计算时, 可包括质点速度计算。 所述质点速度计包括以下步骤:  When the data calculation is performed based on the echo data in step S103, the particle velocity calculation may be included. The particle velocity meter includes the following steps:
步骤 S1034 , 获取目标区域的组织的位移数据后, 取同一个位置不同时刻 的位移数据,沿时间方向求梯度,即可得到该位置处不同时刻的质点速度曲线。 该质点速度曲线反映了剪切波传播过程中各个位置上质点在不同时刻的 运动瞬时速度。对目标区域内的一个或多个位置, 都可以进行上述质点速度计 算过程。对所得质点速度数据, 可以在空间或时间维度进行一定的滤波处理以 降低噪声。 Step S1034: After acquiring the displacement data of the tissue of the target area, taking the displacement data of the same position at different times and obtaining the gradient along the time direction, the particle velocity curve at different times at the position can be obtained. The particle velocity curve reflects the instantaneous velocity of the particle at various points in the shear wave propagation process at different times. The above particle velocity calculation process can be performed for one or more locations within the target area. For the obtained particle velocity data, a certain filtering process can be performed in the spatial or temporal dimension to reduce noise.
步骤 S103中根据所述回波数据进行数据计算时, 可包括传播距离计算。 所述传播距离计算包括以下步骤:  When the data calculation is performed based on the echo data in step S103, the propagation distance calculation may be included. The propagation distance calculation includes the following steps:
步骤 S1035 , 取目标区域中任一位置对应的各个时刻的位移数据, 形成位 移 -时间曲线, 找出该曲线上峰值对应的时刻。 该峰值对应的时刻反映了剪切 波经过该位置的时刻。可以理解的是, 剪切波经过该位置的时刻也为质点速度 -时间曲线的峰值对应的时刻, 也可以利用该位置的质点速度-时间曲线来寻找 对应的剪切波经过该位置的时刻。本步骤中, 可以直接取出多个特定时刻附近 一段时间 Δ Τ内的位移或速度数据, 进行上述传播距离计算。 区别仅在于计算 出来的距离难以准确判断位于目标区域内的哪个位置,但能判断位移剪切波传 播多长时间之后。  Step S1035: Take displacement data of each time corresponding to any position in the target area, form a displacement-time curve, and find a time corresponding to the peak on the curve. The time corresponding to the peak reflects the moment when the shear wave passes through the position. It can be understood that the time when the shear wave passes through the position is also the time corresponding to the peak of the particle velocity-time curve, and the particle velocity-time curve of the position can also be used to find the time when the corresponding shear wave passes the position. In this step, the displacement or velocity data within a certain period of time Δ 附近 near a specific time can be directly taken out, and the above propagation distance calculation is performed. The only difference is that the calculated distance is difficult to accurately determine which position is within the target area, but it can be judged after how long the displacement shear wave is transmitted.
步骤 S1036, 获取目标区域内的所述峰值对应的时刻前后一段时间段 Δ Τ 内的位移 -横向位置曲线或者质点速度-横向位置曲线, 对所述时间段 Δ Τ内各 个时刻的位移-横向位置曲线或者质点速度-横向位置曲线进行互相关判断, 获 取所述时间段内的剪切波横向位置的差异,所述差异反映了剪切波在上述一段 时间 Δ Τ内的传播距离。 对目标区域内的所有位置均进行上述计算处理, 可以 得到剪切波经过目标区域内各位置附近时在一段时间 Δ Τ内的传播距离。在本 步骤中, 所述时间段的长度 Δ Τ可自行确定, 该时间段可以为上述时刻之前至 之后、 上述时刻之前或上述时刻之后。  Step S1036: Obtain a displacement-transverse position curve or a particle velocity-transverse position curve in a period of time Δ Τ before and after the peak corresponding to the peak in the target region, and a displacement-transverse position at each moment in the period Δ Τ The curve or the particle velocity-transverse position curve performs a cross-correlation judgment to obtain a difference in the lateral position of the shear wave in the time period, the difference reflecting the propagation distance of the shear wave within the above-mentioned period Δ Τ. The above calculation processing is performed on all the positions in the target area, and the propagation distance of the shear wave within a period of time Δ 时 when the shear wave passes near each position in the target area can be obtained. In this step, the length Δ Τ of the time period may be determined by itself, and the time period may be before or after the above time, before the above time or after the above time.
可以理解的是,在步骤 S1035中,可直接取出多个特定时刻附近一段时间 △ T内的位移或速度数据, 进行上述传播距离计算。 对上述传播距离求均值, 可以得到目标区域内的一段时间 Δ Τ内的平均传播距离,或者剪切波传播过程 中一段时间 Δ Τ内的平均传播距离。  It can be understood that, in step S1035, the displacement or velocity data within a period of time ΔT near a certain specific time can be directly taken out, and the above-described propagation distance calculation is performed. By averaging the above propagation distances, the average propagation distance within a period of time Δ Τ in the target region or the average propagation distance within a period of time Δ 剪切 during shear wave propagation can be obtained.
将多个连续时刻对应的 Δ Τ内的传播距离取出, 如图 6所示, 可以得到传 播距离 -时间曲线, 也可以累积起来得到较长一段时间 7 ^内的总传播距离。 还可以取出两个不同位置对应的各个时刻的位移数据, 进行互相关比较, 找出 两者的时间差异, 则可以得到剪切波在这两个位置之间的传播时间。 步骤 S103中根据所述回波数据进行数据计算时, 可包括弹性比值计算, 所述弹性比值主要指系统有多个目标区域时,可以计算出各个目标区域的杨氏 模量比值, 以反映各目标区域间的硬度差异程度。 The propagation distances in Δ 对应 corresponding to a plurality of consecutive moments are taken out. As shown in FIG. 6, the propagation distance-time curve can be obtained, or the total propagation distance within 7^ for a long period of time can be accumulated. It is also possible to take out the displacement data of each time corresponding to two different positions, and perform cross-correlation comparison to find out The time difference between the two can obtain the propagation time of the shear wave between these two positions. When the data is calculated according to the echo data in step S103, the elastic ratio calculation may be included. The elastic ratio mainly refers to a ratio of Young's modulus of each target region when the system has multiple target regions, to reflect each The degree of hardness difference between the target areas.
在一定条件下, 剪切波的传播速度与组织硬度有近似固定的关系:  Under certain conditions, the propagation velocity of the shear wave has a nearly fixed relationship with the tissue hardness:
E = 3pc2 上式中, P表示组织密度, £表示组织的杨氏模量值, c表示剪切波的传 播速度。 一定条件下, 杨氏模量越大, 意味着组织硬度越大。 假设两个组织间 的杨氏模量不同, 分别为 和 , 则其比值满足: E = 3pc 2 In the above formula, P represents the tissue density, £ represents the Young's modulus value of the tissue, and c represents the propagation velocity of the shear wave. Under certain conditions, the larger the Young's modulus, the greater the tissue hardness. Assuming that the Young's modulus between the two tissues is different, and they are respectively, the ratio is satisfied:
E E2 = 3plcl 2 /(3p2c2 2) 假设组织密度相近, 即 Α Α, 贝' J : EE 2 = 3p l c l 2 /(3p 2 c 2 2 ) Assuming tissue densities are similar, ie Α Α, 贝' J :
E^ l E2 cl I c2 可以使用一段时间 7。^内的剪切波传播距离的平方比值来计算, 即:
Figure imgf000009_0001
dx 2 ld2 2
E^ l E 2 c l I c 2 can be used for a period of time 7. Calculate the squared ratio of the shear wave propagation distance within ^, ie:
Figure imgf000009_0001
d x 2 ld 2 2
式中, ί 7、 ί 7代表不同的目标区域在一段时间 Tlong内的剪切波传播距离。 当然,也可以直接计算出各个目标区域内的传播速度, 然后计算速度平方 之比来得到弹性比值。 Where ί 7 , ί 7 represent the shear wave propagation distances of different target regions over a period of time T long . Of course, it is also possible to directly calculate the propagation velocity in each target region, and then calculate the ratio of the squares of the velocity to obtain the elastic ratio.
步骤 S103中根据所述回波数据进行数据计算时, 可包括传播速度计算。 对各个目标区域中的任意深度, 取出其各横向位置的质点速度-时间曲线 或位移 -时间曲线, 利用互相关比较找到各曲线之间的时移, 该时移代表剪切 波经过上述参与互相关比较的两个横向位置之间的时间差。横向位置差与时移 的比值即代表相应位置处的传播速度。还比如可以利用波的传播方程直接推导 出近似计算公式如下:  When the data calculation is performed based on the echo data in step S103, the propagation speed calculation may be included. For any depth in each target area, the particle velocity-time curve or displacement-time curve of each lateral position is taken out, and the time shift between the curves is found by cross-correlation comparison, which represents the shear wave passing through the above-mentioned mutual participation. The time difference between the two lateral positions of the correlation comparison. The ratio of the lateral position difference to the time shift represents the propagation speed at the corresponding position. For example, the approximate calculation formula can be directly derived by using the wave propagation equation as follows:
2 d2u I dt2 2 d 2 u I dt 2
C ~  C ~
d u dx + d u dz 式中, c表示传播速度, 可以视为纵向位移数据或纵向速度数据, X代 表横向坐标, z代表 έ从向坐标。  d u dx + d u dz where c is the propagation velocity and can be regarded as longitudinal displacement data or longitudinal velocity data, X represents the lateral coordinate, and z represents the έ slave coordinate.
步骤 S104, 显示测量计算结果。 经过上述弹性测量计算后, 有多种方式 将结果显示。 In step S104, the measurement calculation result is displayed. After the above elastic measurement calculation, there are many ways The result is displayed.
如釆用波前电影图, 即将某时刻的位移数据或质点速度数据取出, 形成当 前时刻在目标区域各位置的分布图,该分布图可以反映当前时刻下剪切波的传 播位置。 当多个时刻的位移数据或质点速度数据分布图连续播放时, 即可形成 剪切波波前传播电影图, 直观表现剪切波的传播过程。  For example, if the wavefront film map is used, the displacement data or the particle velocity data at a certain time is taken out to form a distribution map of the current time at each position of the target region, and the distribution map can reflect the propagation position of the shear wave at the current time. When the displacement data or the particle velocity data distribution map of multiple moments is continuously played, the shear wave wavefront propagating film map is formed, and the propagation process of the shear wave is visually represented.
如釆用传播距离分布图, 即将目标区域内各个位置处的 Δ Τ时间内的传播 距离进行成像显示。  For example, the propagation distance distribution map is used to image the propagation distance of Δ Τ at various positions in the target area.
如釆用传播时间分布图,即将目标区域内各位置附近某固定距离之间的传 播时间进行成像显示。 当然也可以曲线形式显示传播时间。  For example, if the propagation time distribution map is used, the propagation time between a fixed distance near each position in the target area is imaged and displayed. Of course, the propagation time can also be displayed in a curved form.
如釆用传播距离曲线图,即将目标区域内从某个时刻开始往后一段时间内 的传播距离-时间曲线绘制出来, 多个目标区域时可显示多条曲线。  For example, if the propagation distance curve is used, the propagation distance-time curve in the target area from a certain moment to a later period is drawn, and multiple curves can be displayed in multiple target areas.
如直接显示目标区域内的各位置处平均传播距离示意图,或者直接显示剪 切波在各个不同时刻时某小段时间内的平均传播距离示意图。  For example, the schematic diagram of the average propagation distance at each position in the target area is directly displayed, or the average propagation distance of the shear wave at a certain time in a certain period of time is directly displayed.
如当有多个目标区域时, 显示多个目标区域内的弹性比值示意图、或传播 速度平方比示意图、 或传播距离平方比示意图。  For example, when there are multiple target regions, a schematic diagram of the elastic ratio in multiple target regions, or a schematic diagram of the square of the propagation velocity, or a square of the propagation distance is displayed.
上述显示过程中, 可以使用灰阶或彩色编码, 或与其他模式的图像叠加、 融合等方式加强显示效果。  In the above display process, grayscale or color coding may be used, or the display effect may be enhanced by superimposing and blending with other modes of images.
如图 7所示, 本发明的第二较佳实施例提供一种弹性测量检测方法, 本实 施例的弹性测量检测方法与第一较佳实施例大致相同, 包括以下步骤:  As shown in FIG. 7, a second preferred embodiment of the present invention provides an elastic measurement and detection method. The elastic measurement and detection method of the present embodiment is substantially the same as the first preferred embodiment, and includes the following steps:
步骤 S201, 向目标区域发射参考脉冲( reference波), 并接收所述参考脉 冲于所述目标区域内的回波信息, 获取参考脉冲的回波数据,取参考脉冲的回 波数据为基准回波数据。在本实施例中,在每组推动脉冲发射前发射参考脉冲 并接受其在目标区域内的回波信息,获取参考脉冲的回波数据并设其为基准回 波数据。所述参考脉冲的回波信息记录了在产生剪切波并开始传播前的目标区 域的组织信息。  Step S201, transmitting a reference pulse (reference wave) to the target area, and receiving echo information of the reference pulse in the target area, acquiring echo data of the reference pulse, and taking echo data of the reference pulse as a reference echo. data. In the present embodiment, the reference pulse is transmitted and the echo information in the target region is received before each group of push pulses is transmitted, and the echo data of the reference pulse is acquired and set as the reference echo data. The echo information of the reference pulse records the tissue information of the target area before the shear wave is generated and starts to propagate.
步骤 S202, 于目标区域产生剪切波。  Step S202, generating a shear wave in the target area.
步骤 S203, 向目标区域发射追踪脉冲并接收其回波信息, 并获取回波数 据。在本实施例的本步骤中,所述回波数据指追踪脉冲的追踪脉冲的回波数据。  Step S203, transmitting a tracking pulse to the target area and receiving its echo information, and acquiring echo data. In the present step of the embodiment, the echo data refers to echo data of a tracking pulse of a tracking pulse.
步骤 S204, 根据所述回波数据进行数据计算, 获取测量计算结果。 在本 实施例中, 将参考脉冲的回波数据取为基准回波数据,通过追踪脉冲的回波数 据与作为基准回波数据的参考脉冲的回波数据进行数据计算并获取位移、质点 速度、一段时间内各时刻剪切波的传播距离、一段距离内剪切波的传播时间等 计算结果。 Step S204, performing data calculation according to the echo data, and obtaining a measurement calculation result. In this In the embodiment, the echo data of the reference pulse is taken as the reference echo data, and the data is calculated by the echo data of the tracking pulse and the echo data of the reference pulse as the reference echo data, and the displacement, the particle velocity, and the time are obtained. The calculation results of the propagation distance of the shear wave at each moment and the propagation time of the shear wave within a certain distance.
取目标区域内某一横向位置上不同时刻的追踪脉冲的回波数据,分别与该 横向位置上的基准回波数据进行互相关比较,可以得到各个深度位置上的组织 在不同时刻的位移数据。该位移主要反映了剪切波传播导致的组织内质点相对 于剪切波传播前的组织内质点的纵向位移。质点速度、一段时间内各时刻剪切 波的传播距离、一段距离内剪切波的传播时间等计算结果等计算过程于第一较 佳实施大致相同, 在此不再赘述。  The echo data of the tracking pulses at different times in a certain horizontal position in the target region are compared with the reference echo data at the lateral position, and the displacement data of the tissue at different depth positions can be obtained. This displacement mainly reflects the longitudinal displacement of the intragranular mass in the tissue caused by shear wave propagation relative to the intragranular mass before the shear wave propagates. The calculation process such as the particle velocity, the propagation distance of the shear wave at each time in a period of time, and the propagation time of the shear wave within a certain distance are substantially the same in the first preferred embodiment, and will not be described herein.
步骤 S205, 显示测量计算结果。  Step S205, displaying the measurement calculation result.
本实施例中, 在于目标区域产生剪切波之前先行向目标区域发射参考脉 冲, 并获取该参考脉冲的回波数据作为基准回波数据,从而获取剪切波传播导 致的组织内质点相对于剪切波传播前的组织内质点的纵向位移、 相应质点速 度、一段时间内各时刻剪切波的传播距离、一段距离内剪切波的传播时间等计 算结果等计算结果。 本实施例中的其他各个步骤与第一较佳实施例大致相同, 在此不再赘述。  In this embodiment, before the shearing wave is generated in the target region, the reference pulse is first transmitted to the target region, and the echo data of the reference pulse is acquired as the reference echo data, thereby acquiring the intragranular particle relative to the shear caused by the shear wave propagation. The calculation results such as the longitudinal displacement of the mass point in the tissue before the wave propagation, the corresponding particle velocity, the propagation distance of the shear wave at each moment in a period of time, and the propagation time of the shear wave in a distance. The other steps in this embodiment are substantially the same as those in the first preferred embodiment, and are not described herein again.
如图 8所示, 本发明还提供一种超声弹性测量系统, 包括超声探头 11、 控制模块 12、 信号处理模块 13、 计算模块 15及显示装置 17, 所述超声探头 11设有收发模块 110, 所述超声探头 11的收发模块 110、 信号处理模块 13、 计算模块 15、 显示装置 17依次连接, 所述控制模块 12连接于所述收发模块 110。 其中:  As shown in FIG. 8 , the present invention further provides an ultrasonic elastic measuring system, which includes an ultrasonic probe 11 , a control module 12 , a signal processing module 13 , a calculation module 15 , and a display device 17 . The ultrasonic probe 11 is provided with a transceiver module 110 . The transceiver module 110, the signal processing module 13, the calculation module 15, and the display device 17 of the ultrasound probe 11 are sequentially connected, and the control module 12 is connected to the transceiver module 110. among them:
所述收发模块 110用于发射推动脉冲、 追踪脉冲及参考脉冲, 并接收所述 追踪脉冲与所述参考脉冲的回波数据。  The transceiver module 110 is configured to transmit a push pulse, a tracking pulse, and a reference pulse, and receive echo data of the tracking pulse and the reference pulse.
所述控制模块 12用于控制所述收发模块 110发射推动脉冲、 追踪脉冲及 参考脉冲。 实际使用中, 所述控制模块 12按照预设时间间隔发射由推动脉冲、 参考脉冲、 追踪脉冲构成的特定的超声序列,从而便于在目标区域产生剪切波 并供所述超声探头 11的收发模块 110接收相应的回波数据。  The control module 12 is configured to control the transceiver module 110 to transmit a push pulse, a tracking pulse, and a reference pulse. In actual use, the control module 12 transmits a specific ultrasound sequence consisting of a push pulse, a reference pulse, and a tracking pulse according to a preset time interval, thereby facilitating generation of a shear wave in the target region and providing the transceiver module of the ultrasound probe 11 . 110 receives the corresponding echo data.
所述信号处理模块 13用于对回波数据进行信号预处理, 从而便于计算模 块 15进行后续计算, 所述信号预处理可包括波束合成处理, 还可包括如信号 放大、 模数转换、 正交分解等。 The signal processing module 13 is configured to perform signal preprocessing on the echo data, thereby facilitating computational mode Block 15 performs subsequent calculations, which may include beamforming processing, and may also include, for example, signal amplification, analog to digital conversion, orthogonal decomposition, and the like.
计算模块 15用于对波束合成输出的信号进行处理计算。 所述计算模块 15 可用于计算剪切波的各类传播参数, 具体的, 所述计算模块 15包括:  The calculation module 15 is configured to perform processing calculation on the signal of the beam synthesis output. The calculation module 15 can be used to calculate various types of propagation parameters of the shear wave. Specifically, the calculation module 15 includes:
质点速度计算单元 151, 用于计算目标组织的质点速度;  a particle velocity calculation unit 151, configured to calculate a particle velocity of the target tissue;
传播距离计算单元 152, 用于计算剪切波的传播距离;  a propagation distance calculation unit 152, configured to calculate a propagation distance of the shear wave;
杨氏模量比值计算单元 153, 用于计算目标组织的杨氏模量比值杨氏模量 比值。  The Young's modulus ratio calculation unit 153 is configured to calculate the Young's modulus ratio Young's modulus ratio of the target tissue.
计算模块 15通过计算某个位置开始各个时刻的剪切波传播距离, 或某段 时间内的平均距离, 或不同位置开始某段时间内的平均距离比值, 或几个目标 区域内的平均传播速度比值等,从而反映目标区域的物理特性, 并生成相应的 弹性图像。  The calculation module 15 calculates the shear wave propagation distance at each moment from a certain position, or the average distance within a certain period of time, or the average distance ratio of a certain period of time, or the average propagation speed in several target regions. The ratio, etc., reflects the physical characteristics of the target area and generates a corresponding elastic image.
显示装置 17用于将所述计算模块 15生成的弹性图像进行显示。  The display device 17 is for displaying the elasticity image generated by the calculation module 15.
本发明的弹性测量检测方法及系统,通过向待检测组织内部的目标区域内 发射一系列特定脉冲产生剪切波的传播,然后在剪切波传播路径范围内不断发 射接收回波信号, 从该回波信号中提取出一段时间内剪切波的传播距离曲线, 且计算出平均距离, 最终显示。 如果相同发射条件下, 在同一时间内, 剪切波 传播的越远, 则说明组织的硬度越高, 否则说明组织越软。  The elastic measurement detecting method and system of the present invention generates a shear wave propagation by transmitting a series of specific pulses in a target region inside a tissue to be detected, and then continuously transmits and receives an echo signal within a range of the shear wave propagation path, from The echo distance of the shear wave is extracted from the echo signal for a period of time, and the average distance is calculated and finally displayed. If the shear wave propagates at the same time under the same emission conditions, the higher the hardness of the tissue, otherwise the softer the tissue.
本发明的方法不需要操作者控制按压探头的力度即可对反映组织弹性的 相关参数进行计算获取, 使得测量的重复性和稳定性得到保证。 另一方面, 本 发明可以得到定量的参数值, 特别适用于对均勾组织的弹性测量。 同时, 本发 明还可提供不同组织间的弹性差异参数,用以反映不同组织之间的弹性差异程 度。  The method of the invention can calculate and acquire related parameters reflecting the elasticity of the tissue without the operator controlling the strength of pressing the probe, so that the repeatability and stability of the measurement are ensured. On the other hand, the present invention can obtain quantitative parameter values, and is particularly suitable for elastic measurement of uniform hook structures. At the same time, the present invention can also provide elastic difference parameters between different tissues to reflect the degree of elasticity difference between different tissues.
以上所揭露的仅为本发明一种较佳实施例而已,当然不能以此来限定本发 明之权利范围,本领域普通技术人员可以理解实现上述实施例的全部或部分流 程, 并依本发明权利要求所作的等同变化, 仍属于发明所涵盖的范围。  The above is only a preferred embodiment of the present invention, and of course, the scope of the present invention is not limited thereto, and those skilled in the art can understand all or part of the process of implementing the above embodiments, and according to the present invention. The equivalent changes required are still within the scope of the invention.

Claims

权 利 要 求 Rights request
1、 一种弹性测量检测方法, 其特征在于: 包括以下步骤, 1. An elastic measurement and detection method, characterized in that: the following steps are included,
于目标区域产生剪切波;  Producing shear waves in the target area;
向所述目标区域发射追踪脉冲并接收所述追踪脉冲的回波数据; 根据所述回波数据进行数据计算以获取测量计算结果;  Transmitting a tracking pulse to the target area and receiving echo data of the tracking pulse; performing data calculation according to the echo data to obtain a measurement calculation result;
显示所述测量计算结果。  The measurement calculation result is displayed.
2、 如权利要求 1所述的弹性测量检测方法, 其特征在于:  2. The elastic measurement detecting method according to claim 1, wherein:
所述弹性测量检测方法于目标区域产生剪切波前还包括以下步骤, 向目标区域发射参考脉冲,接收所述参考脉冲于所述目标区域内的回波信 息, 获取所述参考脉冲的回波数据。  The elastic measurement detecting method further includes the following steps of generating a shear wavefront in the target region, transmitting a reference pulse to the target region, receiving echo information of the reference pulse in the target region, and acquiring an echo of the reference pulse data.
3、 如权利要求 2所述的弹性测量检测方法, 其特征在于: 当根据所述回 波数据进行数据计算以获取测量计算结果时,通过数据计算获取位移数据, 当 获取位移数据时包括以下步骤,  3. The elastic measurement detecting method according to claim 2, wherein: when the data is calculated according to the echo data to obtain the measurement calculation result, the displacement data is acquired by the data calculation, and the following steps are included when the displacement data is acquired. ,
取目标区域内任一位置的任一时刻的回波数据作为基准回波数据; 将基准回波数据分为多段基准回波核数据,设置所述基准回波核数据的基 准回波核数据中心;  Taking echo data at any time in any position in the target area as reference echo data; dividing the reference echo data into multi-segment reference echo core data, and setting a reference echo core data center of the reference echo core data ;
获取各个时刻的回波数据上与所述基准回波核数据中心之间互相关性最 大的位置,取该互相关性最大的位置与所述基准回波核数据的中心位置之差为 各个时刻的位移数据。  Obtaining a position where the cross-correlation between the echo data of each time and the reference echo core data center is the largest, and taking the difference between the position with the largest cross-correlation and the center position of the reference echo core data for each moment. Displacement data.
4、 如权利要求 3所述的弹性测量检测方法, 其特征在于: 当取目标区域 内任一位置的任一时刻的回波数据作为基准回波数据时,所述基准回波数据釆 用任一时刻的追踪脉冲的回波数据或任一时刻的参考脉冲的回波数据。  4. The elastic measurement detecting method according to claim 3, wherein: when the echo data at any one of the positions in the target region is taken as the reference echo data, the reference echo data is used as the reference echo data. The echo data of the tracking pulse at one time or the echo data of the reference pulse at any time.
5、 如权利要求 3或 4所述的弹性测量检测方法, 其特征在于: 当根据所 述回波数据进行数据计算以获取测量计算结果时,通过数据计算获取不同时刻 的质点速度, 当获取质点速度时包括以下步骤,  The elastic measurement detecting method according to claim 3 or 4, wherein when the data calculation is performed according to the echo data to obtain the measurement calculation result, the particle velocity at different times is obtained by the data calculation, and the particle is acquired. Speed includes the following steps,
获取目标区域的组织的位移数据后, 取同一个位置不同时刻的位移数据; 沿时间方向求梯度以获取该位置处不同时刻的质点速度曲线。  After obtaining the displacement data of the tissue of the target area, the displacement data of the same position is taken at different times; the gradient is obtained along the time direction to obtain the particle velocity curve at different times at the position.
6、 如权利要求 3所述的弹性测量检测方法, 其特征在于: 当根据所述回 波数据进行数据计算以获取测量计算结果时,通过数据计算获取一段时间段内 的剪切波的传播距离, 当获取所述传播距离进行计算时包括以下步骤, 6. The elastic measurement detecting method according to claim 3, wherein: When the wave data is used for data calculation to obtain the measurement calculation result, the propagation distance of the shear wave in a period of time is obtained by data calculation, and the following steps are included when the propagation distance is obtained for calculation.
取目标区域中任一位置对应的各个时刻的位移数据,形成位移 -时间曲线, 找出该曲线上峰值对应的时刻;  Taking displacement data of each moment corresponding to any position in the target area, forming a displacement-time curve, and finding a time corresponding to the peak on the curve;
获取目标区域内的所述峰值对应的时刻前后一段时间段内的位移-横向位 置曲线或者质点速度-横向位置曲线, 对所述时间段内各个时刻的位移 -横向位 置曲线或者质点速度 -横向位置曲线进行互相关判断, 获取所述时间段内的剪 切波的传播距离。  Obtaining a displacement-transverse position curve or a particle velocity-transverse position curve in a period of time before and after the time corresponding to the peak in the target area, a displacement-lateral position curve or a particle velocity-horizontal position for each time in the time period The curve performs cross-correlation judgment to obtain the propagation distance of the shear wave in the time period.
7、 如权利要求 6所述的弹性测量检测方法, 其特征在于: 当根据所述回 波数据进行数据计算以获取测量计算结果时,获取不同目标区域之间的杨氏模 量比值, 当获取所述杨氏模量比值时, 满足以下公式, 7. The elastic measurement detecting method according to claim 6, wherein: when performing data calculation based on the echo data to obtain a measurement calculation result, obtaining a Young's modulus ratio between different target regions, when acquiring When the Young's modulus ratio is satisfied, the following formula is satisfied.
Figure imgf000014_0001
dx 2 ld2 2
Figure imgf000014_0001
d x 2 ld 2 2
式中, di、 ί 7表示不同的目标区域在一段时间 7 ^内的所述剪切波的传播 距离, c表示所述剪切波的传播速度。 Where di, ί 7 represent the propagation distance of the shear wave for a different target region within a period of time 7 ^, and c represents the propagation velocity of the shear wave.
8、 如权利要求 1所述的弹性测量检测方法, 其特征在于: 当显示所述测 量计算结果时, 釆用波前电影图、 传播距离分布图、 传播时间分布图、 传播距 离曲线图、 传播距离示意图、 弹性比值示意图、 传播速度平方比示意图、 传播 距离平方比示意图中的至少一种。  8. The elastic measurement detecting method according to claim 1, wherein: when displaying the measurement calculation result, using a wavefront film map, a propagation distance distribution map, a propagation time distribution map, a propagation distance graph, and a propagation The distance map, the elastic ratio diagram, the square of the propagation velocity, and the square of the propagation distance are at least one of the schematic diagrams.
9、 如权利要求 8所述的弹性测量检测方法, 其特征在于: 当显示所述测 量计算结果时, 釆用灰阶或彩色编码模式。  9. The elastic measurement detecting method according to claim 8, wherein: when the measurement calculation result is displayed, a gray scale or color coding mode is used.
10、 一种弹性测量检测系统, 其特征在于: 包括超声探头、 控制模块、 信 号处理模块、 计算模块及显示装置, 所述超声探头设有收发模块, 所述超声探 头的收发模块、 信号处理模块、 计算模块、 显示装置依次连接, 所述控制模块 连接于所述收发模块,  10 . An elastic measurement and detection system, comprising: an ultrasound probe, a control module, a signal processing module, a calculation module and a display device, wherein the ultrasound probe is provided with a transceiver module, the transceiver module of the ultrasound probe, and a signal processing module The computing module and the display device are sequentially connected, and the control module is connected to the transceiver module.
所述收发模块用于发射推动脉冲、追踪脉冲及参考脉冲, 并接收所述追踪 脉冲与所述参考脉冲的回波数据;  The transceiver module is configured to transmit a push pulse, a tracking pulse and a reference pulse, and receive echo data of the tracking pulse and the reference pulse;
所述控制模块用于控制所述收发模块发射推动脉冲、 追踪脉冲及参考脉 冲;  The control module is configured to control the transceiver module to transmit a push pulse, a tracking pulse, and a reference pulse;
所述信号处理模块用于对回波数据进行信号预处理; 所述计算模块用于对波束合成输出的信号进行处理计算; The signal processing module is configured to perform signal preprocessing on the echo data; The calculation module is configured to process and calculate a signal output by the beam synthesis;
所述显示装置用于将所述计算模块生成的弹性图像进行显示。  The display device is configured to display an elastic image generated by the calculation module.
11、 如权利要求 10所述的弹性测量检测系统, 其特征在于: 所述计算模 块包括:  11. The elastic measurement and detection system of claim 10, wherein: the calculation module comprises:
质点速度计算单元, 用于计算目标组织的质点速度;  a particle velocity calculation unit for calculating a particle velocity of the target tissue;
传播距离计算单元, 用于计算剪切波的传播距离;  a propagation distance calculation unit for calculating a propagation distance of the shear wave;
杨氏模量比值计算单元, 用于计算目标组织的杨氏模量比值杨氏模量比  Young's modulus ratio calculation unit for calculating Young's modulus ratio Young's modulus ratio of target tissue
PCT/CN2014/085821 2014-09-03 2014-09-03 Elasticity measurement and detection method and system WO2016033752A1 (en)

Priority Applications (4)

Application Number Priority Date Filing Date Title
PCT/CN2014/085821 WO2016033752A1 (en) 2014-09-03 2014-09-03 Elasticity measurement and detection method and system
CN201910741683.9A CN110432926B (en) 2014-09-03 2014-09-03 Elasticity measurement detection method and system
CN201480081082.5A CN106572838B (en) 2014-09-03 2014-09-03 Elasticity measurement detection method and system
CN201910741685.8A CN110368031B (en) 2014-09-03 2014-09-03 Elasticity measurement detection method and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/CN2014/085821 WO2016033752A1 (en) 2014-09-03 2014-09-03 Elasticity measurement and detection method and system

Publications (1)

Publication Number Publication Date
WO2016033752A1 true WO2016033752A1 (en) 2016-03-10

Family

ID=55439001

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2014/085821 WO2016033752A1 (en) 2014-09-03 2014-09-03 Elasticity measurement and detection method and system

Country Status (2)

Country Link
CN (3) CN110368031B (en)
WO (1) WO2016033752A1 (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108158610A (en) * 2018-01-16 2018-06-15 苏州国科昂卓医疗科技有限公司 A kind of elastograph imaging method, device, equipment and ultrasound imaging probe
CN109893172A (en) * 2019-02-22 2019-06-18 清华大学 Based on determination method and device, the computer equipment of the mechanics parameter of elastogram
WO2019196033A1 (en) * 2018-04-11 2019-10-17 深圳迈瑞生物医疗电子股份有限公司 Ultrasound elastography method and system
CN111399039A (en) * 2019-01-02 2020-07-10 无锡海斯凯尔医学技术有限公司 Slope parameter extraction method, device and storage medium
CN111399038A (en) * 2019-01-02 2020-07-10 无锡海斯凯尔医学技术有限公司 Slope parameter extraction method and device and computer readable storage medium
CN111401108A (en) * 2019-01-02 2020-07-10 无锡海斯凯尔医学技术有限公司 Slope parameter extraction method and device and computer readable storage medium

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6996035B2 (en) * 2017-11-02 2022-01-17 富士フイルムヘルスケア株式会社 Ultrasonic diagnostic device and method for evaluating physical properties of living tissue

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040167403A1 (en) * 2000-04-05 2004-08-26 Nightingale Kathryn R. Methods, systems, and computer program products for ultrasound measurements using receive mode parallel processing
CN102018533A (en) * 2009-09-11 2011-04-20 迈瑞控股(香港)有限公司 Method for shear wave field formation and ultrasonic image system
CN102469980A (en) * 2009-07-17 2012-05-23 皇家飞利浦电子股份有限公司 Spatially-fine shear wave dispersion ultrasound vibrometry sampling
CN102481137A (en) * 2009-06-30 2012-05-30 皇家飞利浦电子股份有限公司 Push/tracking sequences for shear wave dispersion vibrometry
CN103269639A (en) * 2010-12-22 2013-08-28 皇家飞利浦电子股份有限公司 Shear wave velocity estimation using center of mass
CN103300890A (en) * 2012-03-16 2013-09-18 通用电气公司 System and method for measuring tissue mechanical characteristics
CN103908289A (en) * 2012-12-31 2014-07-09 Ge医疗系统环球技术有限公司 Method for removing background noise in shear wave and corresponding ultrasonic imaging system

Family Cites Families (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6702744B2 (en) * 2001-06-20 2004-03-09 Advanced Cardiovascular Systems, Inc. Agents that stimulate therapeutic angiogenesis and techniques and devices that enable their delivery
EP1800603B1 (en) * 2004-08-05 2013-10-23 Hitachi Medical Corporation Method and apparatus for displaying elastic image and ultrasonograph
US8118744B2 (en) * 2007-02-09 2012-02-21 Duke University Methods, systems and computer program products for ultrasound shear wave velocity estimation and shear modulus reconstruction
US8602994B2 (en) * 2009-03-09 2013-12-10 Mayo Foundation For Medical Education And Research Method for ultrasound vibrometry using orthogonal basis functions
JP5559788B2 (en) * 2009-07-07 2014-07-23 株式会社日立メディコ Ultrasonic diagnostic equipment
CN101699280B (en) * 2009-10-15 2011-08-17 北京索瑞特医学技术有限公司 Method and device for ultrasonic and nondestructive detection of elasticity of viscoelastic medium
CN102667522B (en) * 2009-11-25 2014-10-08 皇家飞利浦电子股份有限公司 Ultrasonic shear wave imaging with focused scanline beamforming
US9678236B2 (en) * 2010-04-27 2017-06-13 Halliburton Energy Services, Inc. Fracture characterization by interferometric drillbit imaging, time reversal imaging of fractures using drill bit seismics, and monitoring of fracture generation via time reversed acoustics and electroseismics
CN101869485B (en) * 2010-06-23 2012-07-04 深圳大学 Ultrasonic imaging method and device
US10041910B2 (en) * 2010-11-14 2018-08-07 University Of Utah Research Foundation Method and system with oppositely-facing ultrasonic transducers for determining tissue pathology
CN103492855B (en) * 2011-02-25 2016-03-30 梅约医学教育与研究基金会 Use the ultrasonic vibration measuring of non-focused ultrasound
US20120253194A1 (en) * 2011-03-30 2012-10-04 Tadashi Tamura Methods and apparatus for ultrasound imaging
CN102824194B (en) * 2011-06-14 2016-08-03 深圳迈瑞生物医疗电子股份有限公司 Displacement detecting method in a kind of elastogram and device
CN102525568A (en) * 2012-01-17 2012-07-04 北京索瑞特医学技术有限公司 Subtraction elastography method
US10357226B2 (en) * 2012-03-12 2019-07-23 Mayo Foundation For Medical Education And Research System and method for model-independent quantification of tissue viscoelastic properties using ultrasound
US9220479B2 (en) * 2012-03-30 2015-12-29 Hitachi Aloka Medical, Ltd. Methods and apparatus for ultrasound imaging
JP5733835B2 (en) * 2012-08-27 2015-06-10 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Measuring device and ultrasonic diagnostic device
CN103845074B (en) * 2012-11-28 2017-12-12 深圳迈瑞生物医疗电子股份有限公司 A kind of ultrasonic elastograph imaging system and method
CN103462643B (en) * 2013-09-29 2015-02-11 深圳市开立科技有限公司 Shear wave speed measurement method, device and system
CN103720489B (en) * 2013-12-30 2015-10-28 中国科学院深圳先进技术研究院 Pathological tissues growth monitoring method and system

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040167403A1 (en) * 2000-04-05 2004-08-26 Nightingale Kathryn R. Methods, systems, and computer program products for ultrasound measurements using receive mode parallel processing
CN102481137A (en) * 2009-06-30 2012-05-30 皇家飞利浦电子股份有限公司 Push/tracking sequences for shear wave dispersion vibrometry
CN102469980A (en) * 2009-07-17 2012-05-23 皇家飞利浦电子股份有限公司 Spatially-fine shear wave dispersion ultrasound vibrometry sampling
CN102018533A (en) * 2009-09-11 2011-04-20 迈瑞控股(香港)有限公司 Method for shear wave field formation and ultrasonic image system
CN103269639A (en) * 2010-12-22 2013-08-28 皇家飞利浦电子股份有限公司 Shear wave velocity estimation using center of mass
CN103300890A (en) * 2012-03-16 2013-09-18 通用电气公司 System and method for measuring tissue mechanical characteristics
CN103908289A (en) * 2012-12-31 2014-07-09 Ge医疗系统环球技术有限公司 Method for removing background noise in shear wave and corresponding ultrasonic imaging system

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108158610A (en) * 2018-01-16 2018-06-15 苏州国科昂卓医疗科技有限公司 A kind of elastograph imaging method, device, equipment and ultrasound imaging probe
CN108158610B (en) * 2018-01-16 2024-02-06 苏州国科昂卓医疗科技有限公司 Elastic imaging method, device, equipment and ultrasonic imaging probe
WO2019196033A1 (en) * 2018-04-11 2019-10-17 深圳迈瑞生物医疗电子股份有限公司 Ultrasound elastography method and system
CN111399039A (en) * 2019-01-02 2020-07-10 无锡海斯凯尔医学技术有限公司 Slope parameter extraction method, device and storage medium
CN111399038A (en) * 2019-01-02 2020-07-10 无锡海斯凯尔医学技术有限公司 Slope parameter extraction method and device and computer readable storage medium
CN111401108A (en) * 2019-01-02 2020-07-10 无锡海斯凯尔医学技术有限公司 Slope parameter extraction method and device and computer readable storage medium
CN111401108B (en) * 2019-01-02 2021-04-06 无锡海斯凯尔医学技术有限公司 Slope parameter extraction method and device and computer readable storage medium
CN109893172A (en) * 2019-02-22 2019-06-18 清华大学 Based on determination method and device, the computer equipment of the mechanics parameter of elastogram
CN109893172B (en) * 2019-02-22 2020-06-19 清华大学 Method and device for determining mechanical parameters based on elastography and computer equipment

Also Published As

Publication number Publication date
CN110368031A (en) 2019-10-25
CN110432926B (en) 2022-06-07
CN106572838B (en) 2019-09-06
CN110368031B (en) 2022-10-04
CN106572838A (en) 2017-04-19
CN110432926A (en) 2019-11-12

Similar Documents

Publication Publication Date Title
WO2016033752A1 (en) Elasticity measurement and detection method and system
EP2926739B1 (en) Acquisition control for elasticity ultrasound imaging
JP5773833B2 (en) System for determining shear wave information in medical ultrasound imaging
KR102332464B1 (en) Shear wave estimation from analytic data
JP6129744B2 (en) Adjusting the measurement of the acoustic radiation force effect on the background motion effect
JP4795675B2 (en) Medical ultrasound system
US9907539B2 (en) Sparse tracking in acoustic radiation force impulse imaging
US11154277B2 (en) Tissue viscoelastic estimation from shear velocity in ultrasound medical imaging
CN107049361B (en) Sound velocity imaging using shear waves
JP2012081269A5 (en)
CN101273903A (en) Ultrasonic imaging apparatus and ultrasonic velocity optimization method
WO2015128656A1 (en) Methods and systems for measuring properties with ultrasound
CN111735526B (en) Ultrasonic elastography device and shear wave velocity measurement method for elasticity measurement
CN110507360B (en) Shear wave imaging method and system
WO2014157510A1 (en) Ultrasonic diagnostic device and ultrasonic diagnostic method
CN106955125A (en) Motion independence in acoustic radiation power Pulse Imageing
JP2013198635A (en) Ultrasonic diagnostic apparatus
JP2015128554A (en) Ultrasonic diagnostic equipment
KR20190132264A (en) Shear wave imaging based on ultrasound with increased pulse repetition interval
CN110946620B (en) Shear wave imaging method, apparatus and storage medium
JP5650430B2 (en) Ultrasound system and method for providing color M-mode video and luminance M-mode video
JP6133130B2 (en) Thickness measuring method and thickness measuring apparatus
WO2019218141A1 (en) Shear wave elasticity measurement method and shear wave elastography system
JP6165024B2 (en) Ultrasonic diagnostic equipment
JP2013244222A (en) Ultrasonic diagnostic equipment, method for determining reliability of set sound velocity, and program

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: 14901247

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

32PN Ep: public notification in the ep bulletin as address of the adressee cannot be established

Free format text: NOTING OF LOSS OF RIGHTS PURSUANT TO RULE 112(1) EPC (EPO FORM 1205A DATED 08.09.2017)

122 Ep: pct application non-entry in european phase

Ref document number: 14901247

Country of ref document: EP

Kind code of ref document: A1