The present invention relates to the technical field of medical instruments, in particularly to a method and a device for elasticity detection.
The elasticity of biological tissues is closely related to a characteristic of a lesion, thereby having important reference value for disease diagnosis. Therefore, the elasticity detection technique for non-destructive elasticity detection of viscoelastic medium such as human liver has been widely used in recent years.
Elasticity detection technique is used to quantitatively detect elasticity information of a tissue, i.e., elasticity modulus. The basic principle is that applying an internal (including its own) or an external dynamic/static/quasi-static excitation on a tissue, and under the action of elastic mechanics, biomechanics and other physical laws, the tissue will produce a response. Since different tissues (including normal and pathological tissues) have different elasticity coefficients (stress/strain), they have different strains after being subjected to an external force. The image transformation is performed on echo signals before and after tissue deformation, thereby visually displaying the elasticity of the tissue by the color of the image, and locating the lesion.
At present, in order to ensure that the elasticity detection results of the tissue are as accurate as possible, it is often necessary to repeat the tissue measurement several times to obtain a final elasticity detection result. That is, an operator places a probe in the corresponding position of the tissue, starts a first measurement, and controls the host for a related processing to obtain a first elasticity detection result; and then places the probe in the corresponding position of the tissue again and starts a second measurement to obtain a second elasticity detection result. The above process is repeated for times. A final elasticity detection result is obtained by considering multiple elasticity detection results obtained above.
In the above measurement, as the operator performs multiple measurements by using a hand-held probe, he needs to locate the measurement position for several times, so it is hard to ensure that the located measurement position for every time is the same position, leading a deviation of the measurement position, which affects accuracy of the measurement result.
In order to solve at least one of the problems mentioned in the background, the present invention provides a method and a device for elasticity detection to overcome the drawback of the conventional multi-measurement method in which the elasticity detection result is poor in accuracy.
The present invention provides a method for elasticity detection, comprising:
step 1: maintaining an excitation device and an ultrasonic transducer in contact with a surface of a tissue; controlling the excitation device to periodically excite N shear waves in the tissue at a preset time interval, and controlling the ultrasonic transducer to transmit an ultrasonic wave to the tissue, N being an integer greater than 1;
step 2: receiving, by the ultrasonic transducer, an ultrasonic echo signal corresponding to each of the shear waves;
step 3: obtaining a propagation characteristic parameter of each of the shear waves according to the ultrasonic echo signal corresponding to each of the shear waves;
step 4: calculating obtaining an elasticity parameter of the tissue, according to the propagation characteristic parameters of the N shear waves and a tissue density of the tissue.
The present invention provides an elasticity detection device, comprising:
a control host, an excitation device and an ultrasonic transducer; the excitation device and the ultrasound transducer are maintained in contact with a surface of a tissue.
The control host comprises:
a controlling module, configured to control the excitation device to periodically excite N shear waves in the tissue at a preset time interval and, control the ultrasonic transducer to transmit an ultrasonic wave to the tissue, where N is an integer greater than 1;
a receiving module, configured to receive by the ultrasonic transducer an ultrasonic echo signal corresponding to each of the shear waves;
an acquiring module, configured to acquire a propagation characteristic parameter of each of the shear waves, according to an ultrasonic echo signal corresponding to each of the shear waves;
a calculating module, configured to calculate an elasticity parameter of the tissue, according to propagation characteristic parameters of the N shear waves and a tissue density of the tissue.
The present invention provides a method and a device for elasticity detection. When an elasticity detection of a tissue is required, the excitation device and the ultrasonic transducer are maintained in contact with the surface of the tissue, and the control host controls the excitation device to continuously and periodically excite N shear waves in the tissue at a preset time interval without manually locating the excitation device several times on the surface of the tissue, thereby ensuring that the N shear waves are excited in respect to the same position, which is helpful to ensure the accuracy of the final elasticity parameter detection result. After transmitting an ultrasonic wave corresponding to each of the shear waves and receiving an ultrasonic echo signal corresponding to each of the shear waves, by the ultrasonic transducer, a propagation characteristic parameter of each of the shear waves are acquired according to the ultrasonic echo signal corresponding to each of the shear waves. Finally, the elasticity parameter of the tissue is calculated according to the propagation characteristic parameter of the N shear waves and the density of the tissue. Since the excitation device and the ultrasonic transducer are maintained in contact with the surface of the tissue, the receiving of the ultrasonic echo signal is continuous, and thereby the elasticity detection efficiency can be improved.
step 101: controlling an excitation device to periodically excite N shear waves in a tissue at a preset time interval, and controlling an ultrasound transducer to transmit an ultrasound wave to the tissue.
Wherein the excitation device and the ultrasonic transducer are maintained in contact with a surface of the tissue, and N is an integer greater than 1.
In the present embodiment, when an elasticity detection of the tissue is required, the excitation device and the ultrasonic transducer are placed on the surface of the tissue, and are maintained in contact with the surface of the tissue.
Further, the control host controls the excitation device to periodically excite shear waves that will propagate in the tissue, at a preset time interval. After excitation of each of the shear waves, the ultrasonic transducer is controlled to transmit ultrasonic waves corresponding to the shear waves and receive corresponding ultrasonic echo signals for analyzing the propagation characteristics of the shear waves in the tissue.
In the present embodiment, the excitation device of the shear waves includes any one of a vibrator, an ultrasonic transducer, and a loudspeaker.
Correspondingly, an excitation mode of the shear waves may be that the vibrator applies a low frequency instantaneous vibration to an outer surface of the tissue so as to produce a shear wave in the tissue; an ultrasonic wave transmitted by the ultrasonic transducer focuses inside the tissue to generate an acoustic radiation force so as to produce shear waves in the tissue; and the loudspeaker generates sound waves with a certain frequency on an outer surface of the tissue so as to produce shear waves in the tissue.
It should be noted that in the mode in which a shear wave is generated by the ultrasonic transducer, the ultrasonic transducer for generating the shear wave and the ultrasonic transducer for transmitting an ultrasonic wave and receiving an ultrasonic echo signal can be the same or different.
Step 102: receiving, by the ultrasonic transducer, an ultrasonic echo signal corresponding to each of the shear waves.
The ultrasonic transducer transmits the received ultrasonic echo signal corresponding to each of the shear waves to a main control host, thereby allowing the control host to perform subsequent processing on each of the ultrasonic echo signals.
Step 103: acquiring a propagation characteristic parameter of each of the shear waves, according to the ultrasonic echo signal corresponding to each of the shear waves.
Step 104: calculating an elasticity parameter of the tissue, according to propagation characteristic parameters of the N shear waves and a tissue density of the tissue.
In the present embodiment, after receiving the ultrasonic echo signal corresponding to each of the shear waves, the control host can perform analysis and processing on each ultrasonic echo signal to obtain a propagation characteristic parameter corresponding to each of the shear waves, and then elasticity parameters of the tissue are acquired according to the propagation characteristic parameter of each of the shear waves.
Specifically, the control host can acquire a tissue motion parameter corresponding to each of the shear waves by performing signal processing on the ultrasonic echo signal corresponding to each of the shear waves, and then the propagation characteristic parameter of each of the shear waves is acquired according to the tissue motion parameter corresponding to each of the shear waves.
Wherein the signal processing on each ultrasonic echo signal includes at least one of the following: time-domain cross-correlation, spectral cross-correlation, square error sum, spot tracking, scale invariant feature point tracking, dynamic programming, zero cross-tracking, and peak searching. The tissue motion parameter corresponding to each of the shear waves, such as displacement or strain, can be acquired through signal processing on each ultrasonic echo signal. Thus, the propagation characteristic parameter corresponding to each of the shear waves, such as propagation speed or propagation attenuation coefficient, can be acquired based on the tissue motion parameter.
After the acquirement of the propagation characteristic parameter of each of the shear waves, the elasticity parameter of the tissue can be calculated according to the propagation characteristic parameter of each of the shear waves and the tissue density of the tissue. The elasticity parameter includes at least one of the following parameters: shear modulus, Young's modulus, shear elasticity, shear viscosity, mechanical impedance, mechanical relaxation time and anisotropy, wherein in general the shear modulus or the Young's modulus is mainly acquired.
In the present embodiment, since the excitation device and the ultrasonic transducer are always maintained in contact with a surface of the tissue during the process of f multiple elasticity detections of the tissue, and start of multiple detection processes by using the excitation device and the ultrasonic transducer are controlled by the control host, there is no need for the operator to relocate the position of the excitation device and the ultrasonic transducer on the surface of the tissue when starting the elasticity detection process every time. Furthermore, the control host can continuously receive ultrasonic echo signals received by the ultrasonic transducer during each elasticity detection process, and after receiving the ultrasonic echo signal corresponding to the each of the shear waves, that is the multiple elasticity detection processes, the control host performs concentration process on the ultrasonic echo signals so as to acquire the propagation characteristic parameter corresponding to each of the shear waves. Finally the elasticity detection results of the tissue, i.e., the elasticity parameters, are acquired according to the propagation characteristic parameter corresponding to each of the shear waves and the tissue density.
Specifically, the control host can calculate the elasticity parameter of the tissue, according to the propagation characteristic parameters of the N shear waves and the tissue density of the tissue, in two alternative ways as shown in
step 201: calculating an average value of the propagation characteristic parameters of the N shear waves; and
step 202: calculating an elasticity parameter of the tissue according to the average value of the propagation characteristic parameters of the N shear waves and the tissue density.
In the present embodiment, it is possible to first calculate an average value of the propagation characteristic parameters, such as displacements, of the N shear waves, and then a correlation calculation is performed based on the average value and the tissue density, to obtain the elasticity parameters of the tissue, such as elasticity modulus and Young's modulus.
step 301: calculating the elasticity parameter of the tissue corresponding to each of the shear waves, according to the propagation characteristic parameter of each of the shear waves and the tissue density;
step 302: calculating an average value of elasticity parameters of the tissue respectively corresponding to the N shear waves, and the elasticity parameter of the tissue is an average value of the elasticity parameters of the tissue respectively corresponding to the N shear waves.
In the present embodiment, it is possible to first calculate the elasticity parameter of the tissue corresponding to each shear wave, that is every elasticity detection process, and then an average value of the N elasticity parameters may be calculated, and this average value may be used as the final elasticity parameter of the tissue.
In the above embodiment, when an elasticity detection of the tissue is required, the excitation device and the ultrasonic transducer are maintained in contact with the surface of the tissue, and the control host controls the excitation device to continuously and periodically excite N shear waves in the tissue at a preset time interval without manually locating the excitation device several times on the surface of the tissue, thereby ensuring that the N shear waves are excited in respect to the same position, which is helpful to ensure accuracy of the final elasticity parameter detection result. After transmitting the ultrasonic wave corresponding to each of the shear waves and receiving the ultrasonic echo signal corresponding to each of the shear waves by the ultrasonic transducer, the propagation characteristic parameter of each of the shear waves is acquired according to the ultrasonic echo signal corresponding to each of the shear waves. Finally, an elasticity parameter of the tissue is calculated according to propagation characteristic parameters of the N shear waves and the density of the tissue. Since the excitation device and the ultrasonic transducer are maintained in contact with the surface of the tissue, the receiving of the ultrasonic echo signals is continuous, and thereby elasticity detection efficiency can be improved.
step 401: controlling the excitation device to periodically excite N shear waves having one or M frequencies in the tissue at a preset time interval, and controlling the ultrasonic transducer to transmit ultrasonic waves to the tissue at the same or different angles and/or frequencies, wherein 2≤M≤N.
Wherein, 2≤M≤N.
In the present embodiment, in order to ensure a more accurate elasticity detection result of the tissue, frequencies of the N shear waves may be different, and further, the angles and frequencies of transmitting the ultrasonic waves and receiving the ultrasonic echo signals by the ultrasonic transducer can be different.
For example, in a mode in which a vibrator excites a shear wave, it is possible to control the vibrator to transmit mechanical vibrations with different frequencies, so that shear waves with different frequencies can be excited in the tissue. For another example, when exciting the shear waves by a loudspeaker, the loudspeaker can be controlled to transmit sound waves with different frequencies, thereby exciting the shear waves with different frequencies.
It should be noted that the control of each of the shear waves and the control of the ultrasound waves may include any one of the following combinations: exciting the shear waves with a single frequency, and transmitting the ultrasonic waves with the same frequency to the tissue at the same or different angles; exciting the shear waves with a single frequency, and transmitting the ultrasonic waves with different frequencies to the tissue at the same angle; exciting the shear waves with M frequencies, and transmitting the ultrasonic waves to the tissue at the same angle or frequency; and exciting the shear waves with M frequencies, and transmitting the ultrasonic waves to the tissue at different angles or frequencies or at different angles and frequencies, and so on.
In the present embodiment, by controlling the excitation device and the ultrasonic transducer to generate the shear waves with different frequencies and the ultrasonic waves with different angles in the tissue, elasticity detection dimension of the tissue is diversified, which is helpful to improve accuracy and reliability of the detection results.
Step 402: receiving, by the ultrasonic transducer, an ultrasonic echo signal corresponding to each of the shear waves.
Step 403: acquiring the propagation characteristic parameter of each of the shear waves according to the ultrasonic echo signal corresponding to each of the shear waves.
Step 404: calculating the elasticity parameter of the tissue according to propagation characteristic parameters of the N shear waves and the tissue density of the tissue.
Step 405: mapping the elasticity parameter of the tissue into a corresponding grayscale or color image by using a grayscale or color mapping.
In the present embodiment, in order to more intuitively understand the elasticity status of the tissue, the calculated elasticity parameter can be mapped in color to obtain a corresponding elasticity image.
a control host 1, an excitation device 2 and an ultrasonic transducer 3.
The control host 1 includes a controlling module 11, a receiving module 12, an acquiring module 13, and a calculating module 14.
The controlling module 11 is configured to control the excitation device to periodically excite N shear waves in a tissue at preset time interval and control the ultrasonic transducer to transmit an ultrasound wave to the tissue, wherein the excitation device and the ultrasonic transducer are maintained in contact with a surface of the tissue, and N is an integer greater than 1.
The receiving module 12 is configured to receive an ultrasonic echo signal corresponding to each of the shear waves by the ultrasonic transducer.
The acquiring module 13 is configured to acquire a propagation characteristic parameter of each of the shear waves according to the ultrasonic echo signal corresponding to each of the shear waves.
The calculating module 14 is configured to calculate an elasticity parameter of the tissue, according to propagation characteristic parameters of the N shear waves and the tissue density of the tissue.
The elasticity detection device of the present embodiment can be used in the technical solutions of the embodiments of the above-mentioned method, and has similar realization principle and technical effect thereto, and thus will not be repeated here.
It will be understood by those of ordinary skill in the art that all or part of the steps of the embodiments of the method described above may be accomplished by using program related hardware, and the program may be stored in a computer readable storage medium. When the program is executed, the steps included in the embodiments of the above-mentioned method will be executed. The above-mentioned storage medium includes various kinds of media that can store a program code, such as a ROM, a RAM, a magnetic disk, or an optical disk.
Finally, it should be noted that, the above embodiments are merely intended to illustrate rather than limit the technical solutions of the present invention; although the present invention has been described in detail with reference to the above embodiments, one with ordinary skill in the art shall understand that modifications can still be made to the technical solutions recorded in the foregoing embodiments, or equivalent substitutions can still be made to part or all of the technical features therein. Neither these modifications nor these substitutions shall make the essence of the corresponding technical solutions depart from the scope of the technical solutions in the above embodiments of the present application.
Number | Date | Country | Kind |
---|---|---|---|
201510718719.3 | Oct 2015 | CN | national |
This application is a continuation of International Application No. PCT/CN2016/103485, filed on Oct. 27, 2016, which claims the priority benefit of China Patent Application No. 201510718719.3, filed on Oct. 29, 2015. The contents of the above identified applications are incorporated herein by reference in their entireties.
Number | Name | Date | Kind |
---|---|---|---|
5606971 | Sarvazyan | Mar 1997 | A |
5955671 | Gilmore et al. | Sep 1999 | A |
20070049824 | Konofagou | Mar 2007 | A1 |
20090056453 | McAleavey | Mar 2009 | A1 |
20110040185 | Matsumura | Feb 2011 | A1 |
20110319756 | Zheng et al. | Dec 2011 | A1 |
20120108968 | Freiburger | May 2012 | A1 |
20120278005 | Sumi | Nov 2012 | A1 |
20130245442 | Hazard | Sep 2013 | A1 |
20130317361 | Tabaru | Nov 2013 | A1 |
20130317362 | Shi et al. | Nov 2013 | A1 |
20150173718 | Tabaru | Jun 2015 | A1 |
20150374338 | Sandrin | Dec 2015 | A1 |
Number | Date | Country |
---|---|---|
101699280 | Apr 2010 | CN |
104380134 | Feb 2015 | CN |
104825195 | Aug 2015 | CN |
105212968 | Jan 2016 | CN |
105662473 | Jun 2016 | CN |
2005-283265 | Oct 2005 | JP |
2005-283265 | Oct 2005 | JP |
2009201989 | Sep 2009 | JP |
2015-023913 | Feb 2015 | JP |
2016-506846 | Mar 2016 | JP |
2012 105 462 | Aug 2013 | RU |
WO2011004661 | Jan 2011 | WO |
2014128182 | Aug 2014 | WO |
2015029651 | Mar 2015 | WO |
WO 2015116893 | Aug 2015 | WO |
Entry |
---|
Publisher: IT'IS Foundation; Title: Tissue Properties Database V1.0; Release Date: Sep. 26, 2011; << https://itis.swiss/virtual-population/tissue-properties/database/density/ >> ; 3 pages total. (Year: 2011). |
Nevo et al., “Automated Tracking of the Mitral Valve Annulus Motion in Apical Echocardiographic Images Using Multidimensional Dynamic Programming”, Ultrasound in Med. & Biol., vol. 33, No. 9, 2007, pp. 1389-1399 (Year: 2007). |
Ning, Minghao. “SIFT (Scale-invariant feature transform)”. Towards Data Science, Apr. 9, 2019, pp. 1-23. << https://towardsdatascience.com/sift-scale-invariant-feature-transform-c7233dc60f37 >> (Year: 2019). |
The Chinese First Examination Report of corresponding Chinese application No. 201510718719.3, dated Aug. 30, 2017. |
The Chinese Second Examination Report of corresponding Chinese application No. 201510718719.3, dated Apr. 13, 2018. |
The extended European Search Report of corresponding European patent application No. 16859040.4-1124/3369382, dated Nov. 18, 2018. |
The Japanese Examination Report of corresponding Japanese application No. 2018-540201, dated Mar. 13, 2019. |
The Russian Examination Report of corresponding Russian Federation Application No. 2018119513/14(030693), dated Jan. 17, 2019. |
The Chinese First Examination Report of corresponding China patent application No. 201510718719.3, dated Aug. 30, 2017. |
International Search Report of corresponding International PCT Application No. PCT/CN2016/103485, dated Dec. 28, 2016. |
Office Action of the corresponding Japan application. |
Number | Date | Country | |
---|---|---|---|
20170311929 A1 | Nov 2017 | US |
Number | Date | Country | |
---|---|---|---|
Parent | PCT/CN2016/103485 | Oct 2016 | US |
Child | 15655844 | US |