CN108553763B - Microwave thermotherapy monitoring method based on ultrasonic echo decorrelation imaging technology - Google Patents
Microwave thermotherapy monitoring method based on ultrasonic echo decorrelation imaging technology Download PDFInfo
- Publication number
- CN108553763B CN108553763B CN201810052311.0A CN201810052311A CN108553763B CN 108553763 B CN108553763 B CN 108553763B CN 201810052311 A CN201810052311 A CN 201810052311A CN 108553763 B CN108553763 B CN 108553763B
- Authority
- CN
- China
- Prior art keywords
- decorrelation
- ultrasonic
- data
- cumdecol
- imaging
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/02—Radiation therapy using microwaves
- A61N5/022—Apparatus adapted for a specific treatment
- A61N5/025—Warming the body, e.g. hyperthermia treatment
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0833—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Animal Behavior & Ethology (AREA)
- Pathology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Physics & Mathematics (AREA)
- Biophysics (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
The invention aims to provide a microwave thermotherapy solidification region detection method based on an ultrasonic decorrelation imaging technology, which is characterized by continuously collecting ultrasonic echo signals under different moments of microwave thermotherapy, respectively carrying out instantaneous decorrelation imaging (decorrelation imaging of every two adjacent frames of ultrasonic signals) and maximum accumulated decorrelation imaging (maximum value of each pixel point in all instantaneous decorrelation images), setting a threshold value, drawing a solidification region in the accumulated decorrelation images, and combining a polynomial fitting technology to realize quantitative detection of the microwave thermotherapy solidification region. The invention utilizes all ultrasonic signals collected in the thermal therapy process, can reduce the influence of signal mutation caused by cavitation and increase the accuracy of detection of a coagulation area.
Description
Technical Field
The invention belongs to the field of signal processing, and particularly relates to a microwave thermotherapy monitoring method based on an ultrasonic echo decorrelation imaging technology.
Background
Liver cancer is one of the common malignant tumors, and surgical resection is considered as the first choice method for treating liver cancer at present. However, due to the influence of the tumor location, size, presence or absence of blood vessels around the tumor, and the like, only a few patients are suitable for surgical operations. In recent years, thermoablation has become an effective new therapy for treating liver tumors. Microwave thermotherapy is widely concerned with its characteristics of no influence from current conduction, little influence from carbonization and blood perfusion, fast temperature rise, large ablation range, etc. Quantitative real-time monitoring of the coagulation zone plays a crucial role in order to ensure maximum killing of tumor cells and protection of normal tissue cells during hyperthermia. Currently, the current practice is. The tumor thermotherapy monitoring technology mainly comprises: magnetic Resonance Imaging (MRI) technology, ultrasonic imaging (USI) technology, and Computed Tomography (CT) technology. MRI has high precision for measuring the tissue temperature in the thermal therapy process, but the real-time performance and the imaging resolution are not high, and the cost is high; CT has certain radiation to human body, and is not suitable for long-term monitoring; USI has become the subject of research by many scholars in terms of its non-destructive, real-time, and low cost. The traditional ultrasonic imaging is easily affected by air bubbles generated in the thermotherapy process, artifacts are generated below the air bubbles, and a coagulation area cannot be accurately depicted. The ultrasonic decorrelation imaging technology is used for tracking tissue degeneration caused by hyperthermia by calculating a decorrelation coefficient of ultrasonic echo data of adjacent frames and imaging. Compared with other quantitative ultrasonic technologies, the decorrelation imaging technology is easier to realize real-time monitoring of tumor hyperthermia.
The invention aims to provide a microwave thermotherapy solidification region detection method based on an ultrasonic decorrelation imaging technology, which is characterized by continuously collecting ultrasonic echo signals under different moments of microwave thermotherapy, respectively carrying out instantaneous decorrelation imaging (decorrelation imaging of every two adjacent frames of ultrasonic signals) and maximum accumulated decorrelation imaging (maximum value of each pixel point in all instantaneous decorrelation images), setting a threshold value, drawing a solidification region in the accumulated decorrelation images, and combining a polynomial fitting technology to realize quantitative detection of the microwave thermotherapy solidification region.
Disclosure of Invention
The invention provides microwave thermotherapy non-invasive monitoring based on an ultrasonic echo decorrelation imaging technology. The ultrasonic echo signals under different moments of microwave thermotherapy are continuously collected, instantaneous decorrelation imaging (decorrelation imaging of every two adjacent frames of ultrasonic signals) and maximum accumulated decorrelation imaging (maximum value of each pixel point in all instantaneous decorrelation images) are respectively carried out, a threshold value is set, a solidification region is depicted in the accumulated decorrelation images by combining a polynomial fitting technology, and quantitative detection of the solidification region of the microwave thermotherapy is realized.
A microwave thermotherapy monitoring method based on ultrasonic echo decorrelation imaging technology comprises the following steps:
step 2, respectively calculating instantaneous decorrelation of every two adjacent frames of ultrasonic data to obtain Insdecol1 and Insdecol K … Insdecol N-1;
step 3, taking the maximum value of each pixel point of all instantaneous decorrelation images Insdecol1 and Insdecol2 … Insdecol N-1 to form a maximum accumulated decorrelation image Cumdecol;
step 4, carrying out logarithmic transformation on Cumdecol for direct comparison with the B-type ultrasonic image to obtain log _ Cumdecol;
step 5, selecting a threshold value according to a coagulation area (gold standard) in the liver tissue section after the thermal therapy is finished, and distinguishing normal liver tissue and heat injury tissue in the maximum range of log _ Cumdecol;
and 6, describing a solidification region by combining a polynomial fitting technology.
Preferably, step 2 specifically comprises:
step 2.1, performing Hilbert transform on two adjacent frames of ultrasonic data _ K and data _ K +1 to obtain complex analytic signals I _ K and I _ K + 1;
step 2.2, performing Gaussian convolution filtering on the I _ K · conj (I _ K +1), the abs (I _ K) · 2 and the abs (I _ K +1) · 2 respectively to obtain R01, R00 and R11;
step 2.3, calculating the instantaneous decorrelation instecolk of the two adjacent frames of ultrasonic data according to the formula 2 (R00. times. R11-abs (R01. times. 2)/(R00. times. R11+ mean (mean (R00. times. R11)));
step 2.4, traversing K to 1: n-1, and repeating the steps 2.1-2.3.
Preferably, step 6 specifically comprises:
step 6.1, carrying out polynomial fitting on each column of Cumdecol1, and then carrying out polynomial fitting on each column to obtain a polynomial fitting decorrelation image poly _ cumecol;
and 6.2, selecting a parameter equipotential line, and drawing a solidification region in the poly _ cumdecol.
Compared with the prior art, the invention has the following advantages and beneficial effects:
1. the invention utilizes all ultrasonic echo signals collected in the thermal therapy process, and can reduce the phenomenon of tissue cavitation
The effect of the mutation causing the signal is used to increase the accuracy of detection of the coagulated region.
Drawings
FIG. 1: a flow chart of the method of the invention;
FIG. 2: flow chart of instantaneous decorrelation algorithm in method
FIG. 3: maximum cumulative decorrelation flow chart in the method of the invention
FIG. 4: method for representing a frozen region based on cumulative decorrelation images
Detailed Description
Fig. 1 shows a flow chart of the method of the present invention, which mainly includes maximum cumulative decorrelation imaging based on ultrasound echo signals and coagulated region delineation based on the cumulative decorrelation images.
As shown in fig. 2, the ultrasound echo transient decorrelation imaging specifically includes:
step 2, performing Hilbert transform on two adjacent frames of ultrasonic data _ K and data _ K +1 to obtain complex analysis signals I _ K and I _ K + 1;
step 3, performing Gaussian convolution filtering on the I _ K · conj (I _ K +1), the abs (I _ K) · 2 and the abs (I _ K +1) · 2 respectively to obtain R01, R00 and R11;
step 4, calculating instantaneous decorrelation Insdeclok of two adjacent frames of ultrasonic data according to a formula 2 (R00. R11-abs (R01. Lambda 2)/(R00. R11+ mean (R00. R11)));
and 5, traversing K to 1: n-1, repeating the steps 2-4;
and 6, taking the maximum value of each pixel point of all instantaneous decorrelation images Insdecol1 and Insdecol2 … Insdecol N-1 to form a maximum accumulated decorrelation image Cumdecol.
As shown in fig. 3, the maximum cumulative decorrelation imaging of the ultrasound echo specifically includes:
Step 2, carrying out logarithmic compression on Cumdecol to obtain log _ Cumdecol
As shown in fig. 4, the depiction of the coagulated region based on the accumulated decorrelated image specifically includes:
Step 2, carrying out polynomial fitting on each column of Cumdecol1, and then carrying out polynomial fitting on each column to obtain a polynomial fitting decorrelation image poly _ cumecol
And 3, selecting a parameter equipotential line, and drawing a solidification region in the poly _ cumdecol.
Claims (3)
1. A microwave thermotherapy monitoring method based on ultrasonic echo decorrelation imaging technology is characterized by comprising the following steps:
step 1, continuously collecting ultrasonic echo signals data _1 and data _2 … data _ N at different moments of microwave thermotherapy;
step 2, respectively calculating instantaneous decorrelation of every two adjacent frames of ultrasonic data to obtain Insdecol1 and Insdecol K … Insdecol N-1;
step 3, taking the maximum value of each pixel point of all instantaneous decorrelation images Insdecol1 and Insdecol2 … Insdecol N-1 to form a maximum accumulated decorrelation image Cumdecol;
step 4, carrying out logarithmic transformation on Cumdecol for direct comparison with the B-type ultrasonic image to obtain log _ Cumdecol;
step 5, selecting a threshold value according to a coagulation area (gold standard) in the liver tissue section after the thermal therapy is finished, and distinguishing normal liver tissue and heat injury tissue in the maximum range of log _ Cumdecol;
and 6, describing a solidification region by combining a polynomial fitting technology.
2. The microwave hyperthermia monitoring method based on ultrasound echo decorrelation imaging technology according to claim 1, wherein the step 2 is specifically:
step 2.1, performing Hilbert transform on two adjacent frames of ultrasonic data _ K and data _ K +1 to obtain complex analytic signals I _ K and I _ K + 1;
step 2.2, performing Gaussian convolution filtering on the I _ K · conj (I _ K +1), the abs (I _ K) · 2 and the abs (I _ K +1) · 2 respectively to obtain R01, R00 and R11;
step 2.3, calculating the instantaneous decorrelation instecolk of the two adjacent frames of ultrasonic data according to the formula 2 (R00. times. R11-abs (R01. times. 2)/(R00. times. R11+ mean (mean (R00. times. R11)));
step 2.4, traversing K to 1: n-1, and repeating the steps 2.1-2.3.
3. The microwave hyperthermia monitoring method based on ultrasound echo decorrelation imaging technique according to claim 1, wherein step 6 is specifically:
step 6.1, carrying out polynomial fitting on each column of Cumdecol1, and then carrying out polynomial fitting on each column to obtain a polynomial fitting decorrelation image poly _ cumecol;
and 6.2, selecting a parameter equipotential line, and drawing a solidification region in the poly _ cumdecol.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810052311.0A CN108553763B (en) | 2018-01-19 | 2018-01-19 | Microwave thermotherapy monitoring method based on ultrasonic echo decorrelation imaging technology |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810052311.0A CN108553763B (en) | 2018-01-19 | 2018-01-19 | Microwave thermotherapy monitoring method based on ultrasonic echo decorrelation imaging technology |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108553763A CN108553763A (en) | 2018-09-21 |
| CN108553763B true CN108553763B (en) | 2020-02-18 |
Family
ID=63530963
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810052311.0A Active CN108553763B (en) | 2018-01-19 | 2018-01-19 | Microwave thermotherapy monitoring method based on ultrasonic echo decorrelation imaging technology |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108553763B (en) |
Family Cites Families (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6193657B1 (en) * | 1998-12-31 | 2001-02-27 | Ge Medical Systems Global Technology Company, Llc | Image based probe position and orientation detection |
| CN100515345C (en) * | 2003-06-12 | 2009-07-22 | 伯拉考开发股份有限公司 | Blood flow estimate through replenishment curve fitting in ultrasound contrast imaging |
| CN1894594B (en) * | 2003-12-16 | 2010-05-26 | 皇家飞利浦电子股份有限公司 | Ultrasonic diagnostic imaging system with automatic control of penetration, resolution and frame rate |
| CN100386057C (en) * | 2005-10-31 | 2008-05-07 | 西安交通大学 | System and method of perfusion imaging and ultrasonic control releasing based-on capsule microbubble |
| CN101854537B (en) * | 2009-03-31 | 2014-03-26 | 深圳迈瑞生物医疗电子股份有限公司 | Ultrasound image data optimization and radiography quantitative analysis method and system |
| US20120143042A1 (en) * | 2010-12-06 | 2012-06-07 | Palmeri Mark L | Ultrasound Methods, Systems and Computer Program Products for Imaging Fluids Using Acoustic Radiation Force |
| CN102499711B (en) * | 2011-09-28 | 2013-07-10 | 无锡祥生医学影像有限责任公司 | Three-dimensional or four-dimensional automatic ultrasound image optimization and adjustment method |
| CN102764139B (en) * | 2012-07-12 | 2014-10-29 | 复旦大学 | Medical ultrasonic beam forming method based on feature space analysis and region identification |
| US20140066759A1 (en) * | 2012-09-04 | 2014-03-06 | General Electric Company | Systems and methods for parametric imaging |
| CN103519848A (en) * | 2013-10-25 | 2014-01-22 | 中国科学院深圳先进技术研究院 | Tissue displacement estimation method and system based on ultrasonic echo radio frequency signals |
| RU2716430C2 (en) * | 2015-12-31 | 2020-03-11 | Конинклейке Филипс Н.В. | Automated blood pool identification system and method of controlling said system |
| CN106037815A (en) * | 2016-05-17 | 2016-10-26 | 西安交通大学 | Ultrasonic echo statistical parameter imaging system and method for thermal coagulation monitoring |
-
2018
- 2018-01-19 CN CN201810052311.0A patent/CN108553763B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN108553763A (en) | 2018-09-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Chen et al. | Ultrasound super-resolution imaging provides a noninvasive assessment of renal microvasculature changes during mouse acute kidney injury | |
| Ireland et al. | Microwave head imaging for stroke detection | |
| US10475192B2 (en) | Motion tracking during non-invasive therapy | |
| CN107569256B (en) | Ultrasonic method for measuring temperature change of biological tissue based on thermal expansion and gating algorithm | |
| Mahoney et al. | Volumetric contrast‐enhanced ultrasound imaging of renal perfusion | |
| Wang et al. | Monitoring radiofrequency ablation with ultrasound Nakagami imaging | |
| Yang et al. | Correlations between B‐mode ultrasonic image texture features and tissue temperature in microwave ablation | |
| US20190351261A1 (en) | Selective resampling during non-invasive therapy | |
| WO2012120498A1 (en) | Non-invasive thermal treatment monitoring | |
| Abbass et al. | In vivo ultrasound thermal ablation control using echo decorrelation imaging in rabbit liver and VX2 tumor | |
| CN106264722A (en) | Window modulation combines Gauss polynomial matching monitoring radio frequency ablation device and method | |
| Li et al. | Ultrasound entropy imaging for detection and monitoring of thermal lesion during microwave ablation of liver | |
| Yang et al. | A systemic study on the performance of different quantitative ultrasound imaging techniques for microwave ablation monitoring of liver | |
| CN108553763B (en) | Microwave thermotherapy monitoring method based on ultrasonic echo decorrelation imaging technology | |
| McCall et al. | Longitudinal 3-D visualization of microvascular disruption and perfusion changes in mice during the evolution of glioblastoma using super-resolution ultrasound | |
| Perlman et al. | Target visualisation and microwave hyperthermia monitoring using nanoparticle-enhanced transmission ultrasound (NETUS) | |
| Wang et al. | High contrast ultrasonic method with multi-spatiotemporal compounding for monitoring catheter-based ultrasound thermal therapy: Development and ex vivo evaluations | |
| Takeuchi et al. | Statistical analysis of ultrasonic scattered echoes enables the non-invasive measurement of temperature elevations inside tumor tissue during oncological hyperthermia | |
| Wang et al. | In vivo Nakagami‐m parametric imaging of microbubble‐enhanced ultrasound regulated by RF and VF processing techniques | |
| Bottiglieri et al. | Monitoring microwave thermal ablation using electrical impedance tomography: an experimental feasibility study | |
| Wang et al. | Ultrasound single-phase CBE imaging for monitoring radiofrequency ablation of the liver tumor: A preliminary clinical validation | |
| CN102525568A (en) | Subtraction elastography method | |
| Yin et al. | Separated respiratory phases for in vivo ultrasonic thermal strain imaging | |
| Wang et al. | Magnetic resonance imaging‑navigated argon‑helium cryoablation therapy against a rabbit VX2 brain tumor | |
| Yuxin et al. | Microwave Ablation Monitoring based on Automatic Segmentation of Ultrasonic Harmonic Envelope Nakagami Parametric Imaging |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |