CN106037815A - Ultrasonic echo statistical parameter imaging system and method for thermal coagulation monitoring - Google Patents
Ultrasonic echo statistical parameter imaging system and method for thermal coagulation monitoring Download PDFInfo
- Publication number
- CN106037815A CN106037815A CN201610326626.0A CN201610326626A CN106037815A CN 106037815 A CN106037815 A CN 106037815A CN 201610326626 A CN201610326626 A CN 201610326626A CN 106037815 A CN106037815 A CN 106037815A
- Authority
- CN
- China
- Prior art keywords
- ultrasonic echo
- statistical parameter
- ultrasonic
- imaging
- thermocoagulation
- 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.)
- Pending
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
-
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
- A61B8/5238—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image
- A61B8/5246—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image combining images from the same or different imaging techniques, e.g. color Doppler and B-mode
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Medical Informatics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biophysics (AREA)
- Physics & Mathematics (AREA)
- Heart & Thoracic Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
The invention discloses an ultrasonic echo statistical parameter imaging system and method for thermal coagulation monitoring. The imaging system comprises main control equipment, a power amplifier, a noninvasive power transducer, an ultrasonic detection and imaging transducer, full-digital ultrasonic equipment, data collection equipment, a sample, image collection equipment, an indicator lamp, a minimally invasive power transducer, a light-transmitting water tank and three-dimensional mobile equipment. According to the invention, the ultrasonic echo statistical parameter imaging method, which overcomes micro-bubble coverage and upper strong echo region sheltering, is adopted to monitor a process of forming thermal coagulation, so that early detection of the initiation of the thermal coagulation is achieved, lower part display of the thermal coagulation is enhanced and a thermal coagulation micro-bubble ratio in a thermal ablation process is improved so as to enhance a contrast ratio of the thermal coagulation monitoring.
Description
[technical field]
The invention belongs to ultrasound detection and Imaging for Monitoring technology, relate to biomedicine and message area, be specifically related to one
Ultrasonic echo statistical parameter imaging system and method for thermocoagulation monitoring.
[background technology]
Raise the thermal ablation techniques causing albumen non reversibility thermocoagulation to be mechanism with temperature and there is noinvasive or Wicresoft, fast
The advantages such as speed, side effect is little, use the heating ablation skill of the realizations of " green " physical energy such as radio frequency, microwave, ultrasonic, laser both at home and abroad
Art achieves challenging achievement in fields such as biology, medical science, chemical industries, uses high-strength focused super the most at home
Acoustic and thermal ablation techniques treatment tumor aspect achieves the achievement of initiative.The detection of thermocoagulation process and Imaging for Monitoring are to realize height
Imitate the basis with safe heating ablation and guarantee.Ultrasonic and nuclear magnetic resonance is the most common Imaging for Monitoring mode.Ultrasonic due to
Its low cost, real-time good, applied widely and be prone to many thermal ablation techniques are combined etc. advantage become currently used relatively
Wide detection and Imaging for Monitoring mode.Existing patent both domestic and external uses conventional ultrasound B-mode image to carry out thermal ablation process
Imaging for Monitoring, such as United States Patent (USP) US6425867B1 and Chinese patent CN1565671A etc..Except causing in thermal ablation process
Outside albumen non reversibility thermocoagulation, there is also cavitation and vaporize the microvesicle caused.What B ultrasonic image actually obtained is to organize to dissipate dorsad
Penetrating information, change and little before the backscatter signals in the thermocoagulation region that heating ablation causes and heating ablation, B ultrasonic image observation arrives
Information major part be the scattered information of bubble that cavitation vaporization causes.Therefore, B ultrasonic image for thermocoagulation process detection with
Imaging for Monitoring has certain limitation.Many scholars have carried out other ultrasonic imaging techniques at the detection of thermocoagulation process and prison
The research of control imaging, estimates and imaging, Ultrasonic tissue characterization parameter Imaging for Monitoring, elastogram etc. including ultrasonic temperature.
Owing to being affected by ablation procedure cavitation vaporization microvesicle and strong echo area, in thermal ablation process solidification initial with
Detection bottom solidification and Imaging for Monitoring have certain difficulty.Ultrasonic echo statistical parameter formation method depends on backscattering
The statistical distribution of echo envelope, is affected less by echo amplitude.
[summary of the invention]
For above-mentioned detection and the deficiency in monitoring, the present invention relates to a kind of ultrasonic echo statistics for thermocoagulation monitoring
Parametric imaging System and method for.Ultrasonic backscattering signal in this invention synchronous monitoring thermal ablation process and ultrasonic statistics ginseng
Amount；And utilize ultrasonic statistical parameter image and B-mode ultrasonoscopy that target area is carried out synchronous dynamic image supervisory control.
In order to realize above-mentioned task, the present invention takes following technical scheme:
A kind of ultrasonic echo statistical parameter imaging system for thermocoagulation monitoring, including main control device；Described master control
Power amplifier, data acquisition equipment, image capture device and three-dimensional mobile device is connected on equipment；Wherein, described merit
The outfan of rate amplifier connects noinvasive power changing device and Wicresoft's power changing device；The input of described data acquisition equipment
End connection has totally digitilized ultrasonic device, and the input of total digitalization ultrasonic device connects ultrasound detection and imaging transducer；
Described three-dimensional mobile device is for regulating sample locus in printing opacity tank；Described image capture device is used for adopting
Collection sample heating ablation dynamic image；
Described main control device controls power amplifier and drives noinvasive power changing device and/or Wicresoft's power changing device to sample
Product emitted energy；Total digitalization ultrasonic device gathers ultrasonic echo rf data by ultrasound detection with imaging transducer, and leads to
Cross data acquisition equipment to be input to main control device and carry out storing and processing.
The described high-definition camera that image capture device resolution is 1280 × 1024, it is arranged at sample side.
Also include the sequencing contro by main control device, in the dynamic image that image capture device collects, characterize the time
Display lamp.
Described noinvasive power changing device is for focusing on ultrasound thermal ablation irradiator；Described Wicresoft's power changing device is microwave
Heating ablation irradiator, radiofrequency irradiator or LASER HEAT melt irradiator；
Described printing opacity tank is transparent glass or lucite tank；Described three-dimensional mobile device is that micron order is swept
Retouch precision and programmable total digitalization 3-D scanning mobile device.
The formation method of a kind of ultrasonic echo statistical parameter imaging system for thermocoagulation monitoring, comprises the following steps:
1) main control device controls the locus of three-dimensional mobile device regulation sample, and sample is positioned over noinvasive power changing
In the Net long wave radiation district of device and/or Wicresoft's power changing device；
2) main control device controls power amplifier driving noinvasive power changing device and/or Wicresoft's power changing device emission energy
Amount, acts on sample；
3) ultrasound detection is directed at sample with imaging transducer；
4) main control device control total digitalization ultrasonic device is penetrated by ultrasound detection and imaging transducer collection ultrasonic echo
Frequency evidence；
5) the ultrasonic echo rf data that ultrasound detection and imaging transducer collect is input to by data acquisition equipment
Main control device is also stored on main control device；
6) the ultrasonic echo rf data on main control device obtains surpassing in sample according to ultrasonic echo statistical parameter algorithm
Sound echo statistical parameter and dynamic changing process thereof.
The imaging of a kind of ultrasonic echo statistical parameter imaging system for thermocoagulation monitoring according to claim 7
Method, it is characterised in that also include:
7) by step 5) in the ultrasonic echo rf data that obtains and step 6) in the ultrasonic echo statistical parameter number that obtains
According to synchronizing to build ultrasonic echo statistical parameter image and the B-mode ultrasonoscopy of dynamically change；
8) by step 5) in the ultrasonic echo rf data that obtains and step 6) in the ultrasonic echo statistical parameter number that obtains
According to, extract thermocoagulation region and the echo amplitude in microvesicle region and ultrasonic echo statistical parameter, according to the meter of thermocoagulation microvesicle ratio
Calculation method obtains thermocoagulation microvesicle ratio and dynamic changing process thereof.
Improve further as side of the present invention, step 6) it is specially and penetrates based on Nakagami statistical model statistics ultrasonic echo
Frequency evidence: according to predetermined, rf data being divided into multiple target area, the ultrasonic echo chosen in tissue in single region is penetrated
Frequency is according to the Nakagami parameter of the ultrasonic echo statistical analysis being calculated in this region according to the algorithm of Nakagami parameter；
Rf data individual processing to each region, ultrasonic echo statistical parameter and dynamic change thereof to multiple regions are carried out simultaneously
Detection.
Improve further as side of the present invention, use ultrasonic echo statistical parameter image based on Nakagami parameter to mesh
Mark region carries out imaging: first ultrasonic echo statistical parameter imaging determines the method for estimation chosen and sliding window size, makes slip
Window is carried out longitudinal direction with single pixel for step-length at imaging plane and is slid laterally, often slide a pixel, utilizes and chooses
Method of estimation calculates sliding window by the statistical parameter value of echo signal envelope, and this value is assigned to the central pixel point of sliding window,
Travel through whole imaging region to obtain ultrasonic echo statistical parameter image.Improve further as side of the present invention, build ultrasonic time
First ultrasonic echo radiofrequency signal is carried out denoising during ripple statistical parameter image, join based on maximal possibility estimation Nakagami
The denoising method of amount imaging is:
First, the ultrasonic echo rf data collected is separately added into twice random white noise, obtains noisy letter
Number S1 and S2, the white noise amplitude of addition is equal to system white noise amplitude, and system white noise amplitude is by no echo area ultrasonic signal
Estimate gained；
Secondly, using pulse width as the partial auto-correlation of sample window wide calculating S1 and S2, correlation coefficient square is obtained
Battle array；
Then, using 80% as threshold value, the value less than 80% in correlation matrix is found out, by whole for the value of they correspondences
It is set to the minima of echo radiofrequency signal, obtains the radiofrequency signal after system white noise；
Finally, will be deemed as the minima that the point of noise is set in echo-signal, y is less than the statistical parameter of 0.25 correspondence
Value is set to 0, the y corresponding statistical parameter value more than 0.2 is solved according to maximum-likelihood estimation.
The present invention compared with prior art, has the advantage that
The monitoring system of the present invention is set by main control device and the power amplifier, the data acquisition that are connected on main control device
Standby, image capture device and three-dimensional mobile device are constituted, and simple in construction, perfect in shape and function, the present invention detects the ultrasonic system of multizone
Meter parameter, it is achieved that ultrasonic statistical parameter synchronous detecting in different modes thermal ablation process, thermocoagulation microvesicle than synchronous detecting with
Damage dynamically changes synchronous detecting.
With existing thermocoagulation detection compared with Imaging for Monitoring technology, the method for the present invention to ultrasonic echo rf data according to
Ultrasonic echo statistical parameter algorithm obtains ultrasonic echo statistical parameter, by analyze ultrasonic echo statistical parameter dynamically change and
Ultrasonic statistical parameter imaging, can overcome the screening of the strong echo area of covering and top of thermal ablation process cavitation vaporization microvesicle
Gear, is used in different modes thermal ablation process detection and the Imaging for Monitoring of solidification, it is achieved heat by ultrasonic statistical parameter formation method
The detection ahead of time that solidification is initial, strengthens solidification bottom and shows, optimizes detection and the Imaging for Monitoring effect of thermocoagulation process.
Further, with existing thermocoagulation detection compared with Imaging for Monitoring technology, invention also improves thermocoagulation microvesicle ratio
Enhance the contrast of image.
[accompanying drawing explanation]
Fig. 1 is the block diagram of the ultrasonic echo statistical parameter monitoring system for thermocoagulation monitoring that the present invention proposes.
Fig. 2 is ultrasonic echo statistical parameter and echo amplitude synchronous detecting and imaging system and method in thermal ablation process institute
The flow chart controlled.
Fig. 3 is the Nakagami parametric image denoising based on maximal possibility estimation and imaging algorithm flow chart improved.
Fig. 4 is that thermal ablation process thermocoagulation initiates ultrasonic echo statistical parameter image, the knot of B-mode image synchronous monitoring
Really.
Fig. 5 is to occur ultrasonic echo statistical parameter image and B-mode when the strong echo area in top blocks in thermal ablation process
Image synchronization monitoring and statistical parameter and the dynamic result of variations of echo amplitude.
Fig. 6 is that in thermal ablation process, ultrasonic echo statistical parameter and B-mode synchronous monitoring are dynamic with thermal image solidification microvesicle ratio
The result of state change.
Wherein, 1, for main control device；2, power amplifier；3, noinvasive power changing device；4, ultrasound detection and imaging transducing
Device；5, total digitalization ultrasonic device；6, data acquisition equipment；7, sample；8, image capture device；9, display lamp；10, Wicresoft's merit
Rate transducer；11, printing opacity tank；12, three-dimensional mobile device.
[detailed description of the invention]
Understand the present invention in order to clearer, be described in further detail below in conjunction with drawings and the embodiments:
A kind of ultrasonic echo statistical parameter imaging system for thermocoagulation monitoring.
Seeing Fig. 1, a kind of ultrasonic echo statistical parameter imaging system for thermocoagulation monitoring of the present invention includes: master control
Equipment 1, power amplifier 2, noinvasive power changing device 3, ultrasound detection and imaging transducer 4, total digitalization ultrasonic device 5, number
Move according to collecting device 6, sample 7, image capture device 8, display lamp 9, Wicresoft's power changing device 10, printing opacity tank 11 and three-dimensional
Equipment 12；Main control device 1 connects power amplifier 2, data acquisition equipment 6, image capture device 8, display lamp 9 and three
Dimension mobile device 12；Wherein, the outfan of main control device 1 is connected with the control end of power amplifier 2, power amplifier 2 defeated
Go out end to be connected with the input of noinvasive power changing device 3 and Wicresoft's power changing device 10, noinvasive power changing device 3 and Wicresoft's power
Transducer 10 is arranged on the circumference of printing opacity tank 11；The outfan of main control device 1 is connected with the control end of data acquisition equipment 6,
The input of data acquisition equipment 6 is connected with the outfan of total digitalization ultrasonic device 5, the input of total digitalization ultrasonic device 5
End is connected with the outfan of ultrasound detection with imaging transducer 4；The outfan of main control device 1 and the control of three-dimensional mobile device 12
End is connected, and three-dimensional mobile device 12 is for regulating the sample 7 locus at printing opacity tank 11.
This system controls power amplifier 2 by main control device 1 and drives noinvasive power changing device 3 and/or Wicresoft's power to change
Can device 10 emitted energy；This system controls total digitalization ultrasonic device 5 by ultrasound detection and imaging transducing by main control device 1
Device 4 gathers ultrasonic echo rf data, and is input to main control device 1 by data acquisition equipment 6 and is stored in main control device 1
On；Image capture device 6 is acquired, and display lamp 9 is controlled by main control device 1, characterizes the time in the dynamic image collected.
Wherein, main control device 1 is time series and energy output two ore control and the many merits possessing data-storing and process
Energy main control device, wherein, main control device 1 can synchronize space-time to the time of power amplifier 2 energy output and amplitude and control, and
To total digitalization ultrasonic device 5 by the ultrasonic echo rf data that ultrasound detection and imaging transducer 4 gather carry out storage and
Process.
Described noinvasive power changing device 3 is for focusing on ultrasound thermal ablation irradiator；Described Wicresoft's power changing device 10 is
Microwave thermal melts irradiator；Described Wicresoft's power changing device 10 is radiofrequency irradiator；Described Wicresoft's power changing
Device 10 melts irradiator for LASER HEAT；
Described printing opacity tank 11 is transparent glass or lucite tank；Described three-dimensional mobile device 12 is micron
Level scanning accuracy and the totally digitilized 3-D scanning mobile device that can program.
Described image capture device 8 resolution is the high-definition camera of 1280 × 1024, and it is arranged at sample 7 side.
A kind of ultrasonic echo statistical parameter detection method overcoming microvesicle to cover reflection thermocoagulation initial period, concrete steps
As follows:
The present invention detects and solidifies the ultrasonic echo statistical parameter in initial target region in thermal ablation process and dynamically change
Process, is shown in Fig. 2.Main control device 1 controls three-dimensional mobile device 12 and regulates sample 7 and be at noinvasive power changing device 3 and/or micro-
In the Net long wave radiation district of wound power changing device 10；Ultrasound detection is directed at imaging transducer 4 Net long wave radiation of transducer simultaneously
Target area in district；System by default heating ablation state modulator power amplifier 2 drive noinvasive power changing device 3 and/or
Wicresoft's power changing device 10 emitted energy, acts on sample；System controls total digitalization ultrasonic device and passes through ultrasound detection and become
As transducer 4 gathers the ultrasonic echo rf data that target area heating ablation solidification is initial.
System uses total digitalization ultrasonic device 5 to gather ultrasonic echo rf data, is divided by rf data according to predetermined
For multiple target areas, the rf data individual processing to each region, simultaneously the ultrasonic echo statistical parameter to multiple regions
And dynamically change detects.The present invention uses Nakagami statistical model statistics ultrasonic backscattering signal.Choose group
Knit ultrasonic echo rf data in single region according to the algorithm of Nakagami parameter be calculated in this region ultrasonic
The Nakagami parameter of echo statistical analysis.Have a lot for obtaining the method for the estimator of Nakgami parameter, substantially can divide
For method for parameter estimation based on moments estimation and method for parameter estimation based on maximal possibility estimation.
A kind of ultrasonic echo statistical parameter formation method overcoming microvesicle covering detection thermocoagulation ahead of time initial, concrete steps
As follows:
The present invention uses B-mode image and ultrasonic echo statistical parameter image to detect heating ablation solidification is initial simultaneously
With Imaging for Monitoring.Main control device 1 controls power amplifier 2 by parameter preset and drives noinvasive power changing device 3 and/or Wicresoft's merit
Rate transducer 10 emitted energy, acts on sample 7；Control total digitalization ultrasonic device 5 to be changed with imaging by ultrasound detection simultaneously
The ultrasonic echo rf data that target area heating ablation solidification is initial can be gathered by device 4.Use the ultrasonic echo radio frequency number collected
According to synchronizing to build ultrasonic echo statistical parameter image and the B-mode ultrasonoscopy of dynamically change.
The ultrasonic backscattering information of the target area of conventional B mode image reaction, has reacted organization internal acoustic impedance difference
Difference, be at present the most ripe ultrasonic imaging method.Can produce cavitation vaporization microvesicle in thermal ablation process, cavitation vaporization is micro-
Bubble and damage all show as high echo in B-mode image, thus cause the initial detection of heating ablation solidification to have certain difficulty.
In Nakagmi statistical model, the Nakagami parameter that different scattering son distributions obtain is different, and Nakagmi parametric image has point
Distinguish the different ability scattering sub-concentration, and do not affected by echo amplitude.
For the impact overcoming heating ablation solidification initial cavitation vaporization microvesicle to block, use based on Nakagami parameter super
Sound echo statistical parameter image carries out imaging to target area.First ultrasonic echo statistical parameter imaging determines the estimation side chosen
Method and sliding window size, make sliding window be carried out longitudinal direction with single pixel for step-length at imaging plane and slide laterally, often slide
One pixel, the method for estimation calculating sliding window that utilization is chosen is by the statistical parameter value of echo signal envelope, and this value is composed
To the central pixel point of sliding window, travel through whole imaging region to obtain ultrasonic echo statistical parameter image.Due to thermocoagulation group
Knit concentration sub-with the scattering of microvesicle different, therefore can be by the detection initial to thermocoagulation of ultrasonic echo statistical parameter image and prison
Control imaging.
A kind of ultrasonic echo statistical parameter detection method overcoming the strong echo area in top to block reflection thermocoagulation bottom, tool
Body step is as follows:
The present invention monitors and solidifies the ultrasonic echo statistical parameter in bottom target region in thermal ablation process and dynamically change
Process.Main control device 1 controls three-dimensional mobile device 12 and regulates sample 7 and be at noinvasive power changing device 3 and/or Wicresoft's power
In the Net long wave radiation district of transducer 10；Ultrasound detection is directed in the Net long wave radiation district of transducer with imaging transducer 4 simultaneously
Target area；System drives noinvasive power changing device 3 and/or Wicresoft's power by default heating ablation state modulator power amplifier
Transducer 7 emitted energy, acts on sample；System is controlled total digitalization ultrasonic device and is adopted with imaging transducer by ultrasound detection
Ultrasonic echo rf data when being blocked bottom the heating ablation solidification of collection target area.
Along with the increase of heating ablation time, due to damage upper part strong scattering coefficient tissue and the impact of microvesicle, bottom surpasses
Sound deep fades so that damage is difficult to observe, and ultrasonic echo statistical parameter image can effectively suppress in backscattering echo
Top strong echo area occlusion effect, the weaker signal that B-mode image can not describe can be identified.White due to no echo area
Noise can affect ultrasonic echo statistical parameter image detection thermocoagulation region, needs to take in the case of echo scattering waveform is constant
Go out white noise.NCA denoising method obtains better effects in Nakagami parametric imaging based on moments estimation.NCA denoising method
Directly noise spot being set to 0, in ultrasonic echo statistical parameter imaging algorithm based on moments estimation, ultrasonic echo statistical parameter passes through
Utilize the k rank square μ of echo rf data envelope RkCalculate and obtain, make E () represent and expect, then,
μk=E (Rk)
Therefore noise spot is set to zero by the ultrasonic echo statistical parameter imaging algorithm of moments estimation can realize preferably going
Make an uproar effect, but in statistical parameter imaging algorithm based on maximal possibility estimation
From above formula, according to traditional NCA denoising method, noise spot is set to 0, can cause the m value calculating other points
Time, as long as sliding window comprises this point, may result in G=0, also result in y and trend towards infinity, cause the m value calculated to be not intended to
Justice, the present invention proposes the Nakagami image de-noising method based on maximal possibility estimation of a kind of improvement.See Fig. 3, first,
The ultrasonic echo rf data collected is separately added into twice random white noise, obtains noisy signal S1 and S2, add
White noise amplitude equal to system white noise amplitude, no echo area ultrasonic signal estimate gained；Secondly, using pulse width as
Sampling window width calculates the partial auto-correlation of S1 and S2, obtains correlation matrix；Then, using 80% as threshold value, phase is found out
Close coefficient matrix medium and small in the value of 80%, the value of they correspondences is all entered as the minima of echo radiofrequency signal, is gone
Radiofrequency signal after system white noise.Will be deemed as the minima that the point of noise is set in echo-signal, so can cause calculating
Y value out is far smaller than its normal value, therefore less than the statistical parameter value of 0.2 correspondence, y is set to 0, by right more than 0.2 of y
The statistical parameter value answered solves according to maximum-likelihood estimation.Preferable denoising effect can be reached.
A kind of overcome the strong echo area in top to block to strengthen the ultrasonic echo statistical parameter imaging side that shows, thermocoagulation bottom
Method, specifically comprises the following steps that
The present invention uses ultrasonic echo statistical parameter image and B-mode image to occurring in thermal ablation process that top is strong simultaneously
The thermocoagulation that echo area blocks carries out detecting and Imaging for Monitoring.Main control device controls power amplifier by parameter preset and drives
Noinvasive power changing device 3 and/or Wicresoft's power changing device 10 emitted energy, act on sample 7；Control total digitalization ultrasonic simultaneously
Equipment 5 gathers, by ultrasound detection and imaging transducer 4, the ultrasonic echo rf data that target area heating ablation solidification is initial.
The ultrasonic echo rf data collected is used to synchronize to build ultrasonic echo statistical parameter image and the B of dynamically change
Model ultrasonic image.Build ultrasonic echo statistical parameter image time first ultrasonic echo radiofrequency signal is carried out denoising.
After melting end, gather tangent plane along ultrasound detection and imaging transducer 4 and cut sample, and acquisition optics that tangent plane is taken pictures open
Image.Optical imagery, B-mode image and ultrasonic echo statistical parameter image are analyzed, when research is along with heating ablation
Between increase optical imagery when the strong echo area in top blocks, B-mode image and ultrasonic echo statistical parameter image and damage occur
Corresponding relation.
A kind of thermocoagulation microvesicle that improves monitors the ultrasonic echo statistical parameter formation method of contrast, bag than enhancing thermocoagulation
Include following steps:
The present invention detect target area in thermal ablation process thermocoagulation microvesicle ratio (lesion-to-bubble ratio,
And dynamic changing process LBR).System controls three-dimensional mobile device regulation sample by main control device and is at noinvasive power
In the Net long wave radiation district of transducer 3 and/or Wicresoft's power changing device 10；Ultrasound detection is directed at transducing with imaging transducing 4 simultaneously
Target area in the Net long wave radiation district of device；System drives noinvasive power to change by default heating ablation state modulator power amplifier
Energy device and/or Wicresoft's power changing device emitted energy, act on sample 7；System controls total digitalization ultrasonic device by ultrasonic
Detection and imaging transducer gather ultrasonic echo rf data when being blocked bottom the heating ablation solidification of target area.
LBR refers to the ratio of solidification and microbubbles scatter ability, is defined as:
Wherein, IlesionAnd IbubbleIt is respectively the average amplitude of solidification and microvesicle target area.Utilize the echo collected
Rf data and be calculated ultrasonic echo statistical parameter data.Extract area-of-interest (region of interest, ROI)
Echo amplitude and ultrasonic echo statistical parameter, be calculated LBR value and the dynamic changing process thereof of image according to formula.
The following is specific embodiment, but do not limit to and these embodiments.
Embodiment 1
See Fig. 4, the initial detection of solidification and ultrasonography monitoring imaging in thermal ablation process.At thermal ablation therapy in early days, use
Ultrasonic echo statistical parameter formation method detects and Imaging for Monitoring thermocoagulation is initial.
In the present embodiment, use focusing is ultrasonic carries out thermal ablation therapy.The acoustical power used is 85W, and dutycycle is 70%,
Total time is 4s, in thermal ablation process, uses image capture device collection to focus on transparent artificial in ultrasound thermal ablation therapeutic process
The dynamic changing process of damage in biological tissue's body mould, if shown in 4 (a).It is with digital super in thermal ablation process simultaneously
The ultrasonic echo backscattering data of target area during acoustic equipment synchronous acquisition thermal ablation therapy.These are used to collect
Ultrasonic echo rf data synchronizes to build ultrasonic B-mode image Fig. 4 (b) and ultrasonic echo statistical parameter image 4 (c), heat is disappeared
The ultrasonic echo statistical parameter image, B-mode image and the optical imagery that melt solidification initial are analyzed,
Embodiment 2
See Fig. 5, thermal ablation process occurs ultrasonic echo statistics in thermocoagulation bottom when the strong echo area in top blocks heat
Parameter dynamically changes and Imaging for Monitoring.In thermal ablation process, the ultrasonic echo amplitude in regions multiple to sample and ultrasonic echo
Statistical parameter carries out synchronous detecting.Synchronize to build ultrasonic echo statistical parameter image and ultrasonic B-mode image simultaneously.
The present embodiment use microwave therapy apparatus carry out thermal ablation therapy.Microwave thermal ablation parameter is: ablation power is
40W, the time of melting is 5min, in thermal ablation process, uses total digitalization ultrasonic device to control ultrasound detection and imaging transducing
Device gathers the ultrasonic echo rf data of heating ablation target area.Synchronize to build B-mode image Fig. 5 (a) and ultrasonic echo is added up
Parametric image Fig. 5 (b), in choosing area-of-interest monitored area respectively, ultrasonic echo amplitude 5 (c) is added up with ultrasonic echo simultaneously
The dynamic change of parameter 5 (d).
Embodiment 3
See Fig. 6, thermocoagulation Imaging for Monitoring based on the ultrasonic echo statistical parameter improving thermocoagulation microvesicle ratio.Disappear in heat
During melting, extract ultrasonic echo data and the ultrasonic echo statistical parameter data in resolidified region and microvesicle region respectively, calculate
The average amplitude of ROI region, is calculated LBR value according to formula and dynamically changes, and synchronizes to build ultrasonic echo statistics simultaneously
Parametric image and B-mode ultrasonoscopy.
Using in the present embodiment and focus on the ultrasonic thermal ablation therapy that carries out, use acoustical power is 140W, continuous radiation 32s,
During thermal ablation therapy, use target area ultrasonic echo radio frequency in total digitalization ultrasonic device synchronous acquisition thermal ablation process
Data.These rf datas are used to synchronize to build B-mode image Fig. 6 (a) and ultrasonic echo statistical parameter figure in thermal ablation process
As Fig. 6 (b).Utilize the rf data that collects to obtain dynamically the changing of LBR of B-mode image and statistical parameter image simultaneously
Journey Fig. 6 (c), concrete method, for being that to choose resolidified region and microvesicle region respectively be ROI region, calculates echo in ROI region
Amplitude and the meansigma methods of ultrasonic echo statistical parameter value.LBR value and dynamic changing process is obtained according to the computational methods of LBR value.
Above content is to combine concrete preferred implementation further description made for the present invention, it is impossible to assert
The detailed description of the invention of the present invention is only limitted to this, for general technical staff of the technical field of the invention, is not taking off
On the premise of present inventive concept, it is also possible to make some simple deduction or replace, all should be considered as belonging to the present invention by institute
The claims submitted to determine scope of patent protection.
Claims (10)
1. the ultrasonic echo statistical parameter imaging system for thermocoagulation monitoring, it is characterised in that: include main control device
(1)；Connect on described main control device (1) power amplifier (2), data acquisition equipment (6), image capture device (8) and
Three-dimensional mobile device (12)；Wherein, the outfan of described power amplifier (2) connects noinvasive power changing device (3) and micro-
Wound power changing device (10)；The input of described data acquisition equipment (6) connects has totally digitilized ultrasonic device (5), totally
The input of word ultrasonic device (5) connects ultrasound detection and imaging transducer (4)；Described three-dimensional mobile device (12) is used
In regulation sample (7) locus in printing opacity tank (11)；Described image capture device (8) is used for gathering sample (7)
Thermocoagulation dynamic image；
Described main control device (1) controls power amplifier (2) and drives noinvasive power changing device (3) and/or Wicresoft's power changing
Device (10) is to sample (7) emittance；Total digitalization ultrasonic device (5) gathers super by ultrasound detection and imaging transducer (4)
Sound echo rf data, and be input to main control device (1) carry out storing and processing by data acquisition equipment (6).
A kind of ultrasonic echo statistical parameter imaging system for thermocoagulation monitoring the most according to claim 1, its feature
It is: described image capture device (8) resolution is the high-definition camera of 1280 × 1024, and it is arranged at sample (7) side.
A kind of ultrasonic echo statistical parameter imaging system for thermocoagulation monitoring the most according to claim 1, its feature
Be: also include the sequencing contro by main control device (1), in the dynamic image that image capture device (8) collects characterize time
Between display lamp (9).
A kind of ultrasonic echo statistical parameter imaging system for thermocoagulation monitoring the most according to claim 1, its feature
It is: described noinvasive power changing device (3) is for focusing on ultrasound thermal ablation irradiator；Described Wicresoft's power changing device (10) is
Microwave thermal melts irradiator, radiofrequency irradiator or LASER HEAT and melts irradiator.
A kind of ultrasonic echo statistical parameter imaging system for thermocoagulation monitoring the most according to claim 1, its feature
It is: described printing opacity tank (11) is transparent glass or lucite tank；Described three-dimensional mobile device (12) is micro-
Meter level scanning accuracy and programmable total digitalization 3-D scanning mobile device.
6. imaging side based on a kind of ultrasonic echo statistical parameter imaging system for thermocoagulation monitoring described in claim 1
Method, it is characterised in that comprise the following steps:
1) main control device (1) controls the locus of three-dimensional mobile device (12) regulation sample (7), and sample (7) is positioned over nothing
In the Net long wave radiation district of wound power changing device (3) and/or Wicresoft's power changing device (10)；
2) main control device (1) controls power amplifier (2) driving noinvasive power changing device (3) and/or Wicresoft's power changing device
(10) emitted energy, acts on sample (7)；
3) ultrasound detection is directed at sample (7) with imaging transducer (4)；
4) main control device (1) control total digitalization ultrasonic device (5) gathers ultrasonic time by ultrasound detection and imaging transducer (4)
Ripple rf data；
5) the ultrasonic echo rf data that ultrasound detection and imaging transducer (4) collect is inputted by data acquisition equipment (6)
To main control device (1) and be stored on main control device (1)；
6) the ultrasonic echo rf data on main control device (1) obtains in sample (7) according to ultrasonic echo statistical parameter algorithm
Ultrasonic echo statistical parameter and dynamic changing process thereof.
The imaging side of a kind of ultrasonic echo statistical parameter imaging system for thermocoagulation monitoring the most according to claim 6
Method, it is characterised in that also include:
7) by step 5) in the ultrasonic echo rf data that obtains and step 6) in the ultrasonic echo statistical parameter data that obtain with
Step builds ultrasonic echo statistical parameter image and the B-mode ultrasonoscopy of dynamically change；
8) by step 5) in the ultrasonic echo rf data that obtains and step 6) in the ultrasonic echo statistical parameter data that obtain, carry
Take thermocoagulation region and the echo amplitude in microvesicle region and ultrasonic echo statistical parameter, according to the computational methods of thermocoagulation microvesicle ratio
Obtain thermocoagulation microvesicle ratio and dynamic changing process thereof.
8. according to the one-tenth of a kind of ultrasonic echo statistical parameter imaging system for thermocoagulation monitoring described in claim 6 or 7
Image space method, it is characterised in that step 6) it is specially and adds up ultrasonic echo rf data based on Nakagami statistical model: according in advance
Fixed rf data is divided into multiple target area, the ultrasonic echo rf data chosen in tissue in single region according to
The Nakagami parameter of the ultrasonic echo statistical analysis that the algorithm of Nakagami parameter is calculated in this region；To each region
Rf data individual processing, simultaneously to the ultrasonic echo statistical parameter in multiple regions and dynamically change detect.
The imaging side of a kind of ultrasonic echo statistical parameter imaging system for thermocoagulation monitoring the most according to claim 8
Method, it is characterised in that use ultrasonic echo statistical parameter image based on Nakagami parameter that target area is carried out imaging: super
First sound echo statistical parameter imaging determines the method for estimation chosen and sliding window size, makes sliding window with single pixel as step-length
Being carried out longitudinal direction at imaging plane and slide laterally, often slide a pixel, utilizes the method for estimation chosen to calculate sliding window
By the statistical parameter value of echo signal envelope, and this value is assigned to the central pixel point of sliding window, travel through whole imaging region with
Obtain ultrasonic echo statistical parameter image.
The imaging of a kind of ultrasonic echo statistical parameter imaging system for thermocoagulation monitoring the most according to claim 9
Method, it is characterised in that first ultrasonic echo radiofrequency signal is carried out denoising when building ultrasonic echo statistical parameter image,
Denoising method based on maximal possibility estimation Nakagami parametric imaging is:
First, the ultrasonic echo rf data collected is separately added into twice random white noise, obtains noisy signal S1
And S2, the white noise amplitude of addition is equal to system white noise amplitude, and system white noise amplitude is estimated by no echo area ultrasonic signal
Gained；
Secondly, calculate the partial auto-correlation of S1 and S2 using pulse width as sampling window width, obtain correlation matrix；
Then, using 80% as threshold value, find out the value less than 80% in correlation matrix, the value of they correspondences is all arranged
For the minima of echo radiofrequency signal, obtain the radiofrequency signal after system white noise；
Finally, will be deemed as the minima that the point of noise is set in echo-signal, y is put less than the statistical parameter value of 0.2 correspondence
It is 0, the y corresponding statistical parameter value more than 0.2 is solved according to maximum-likelihood estimation.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610326626.0A CN106037815A (en) | 2016-05-17 | 2016-05-17 | Ultrasonic echo statistical parameter imaging system and method for thermal coagulation monitoring |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610326626.0A CN106037815A (en) | 2016-05-17 | 2016-05-17 | Ultrasonic echo statistical parameter imaging system and method for thermal coagulation monitoring |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106037815A true CN106037815A (en) | 2016-10-26 |
Family
ID=57177025
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610326626.0A Pending CN106037815A (en) | 2016-05-17 | 2016-05-17 | Ultrasonic echo statistical parameter imaging system and method for thermal coagulation monitoring |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106037815A (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108553763A (en) * | 2018-01-19 | 2018-09-21 | 北京工业大学 | A kind of microwave heat therapeutic monitoring method based on ultrasonic echo decorrelation imaging technique |
| CN109171998A (en) * | 2018-10-22 | 2019-01-11 | 西安交通大学 | Heating ablation region recognition monitoring imaging method and system based on ultrasonic deep learning |
| CN112566694A (en) * | 2018-05-25 | 2021-03-26 | 卡尔迪亚韦弗公司 | Ultrasound treatment device comprising means for imaging cavitation bubbles |
| CN112966212A (en) * | 2021-02-07 | 2021-06-15 | 西安交通大学 | Multi-parameter real-time monitoring imaging system based on ultrasonic echo backscattering energy change |
| CN113117258A (en) * | 2019-12-30 | 2021-07-16 | 重庆融海超声医学工程研究中心有限公司 | Device and method for detecting tissue coagulation necrosis |
| CN113907792A (en) * | 2021-08-31 | 2022-01-11 | 西安交通大学 | Nonlinear decorrelation fusion ultrasonic multi-parameter thermal coagulation detection and imaging method |
| CN114305668A (en) * | 2021-12-22 | 2022-04-12 | 聚融医疗科技(杭州)有限公司 | Ultrasonic thermal ablation multi-parameter monitoring method and system based on demodulation domain parametric imaging |
| CN114587565A (en) * | 2022-03-01 | 2022-06-07 | 河南中医药大学 | Temperature control method and system in radio frequency ablation |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1847824A (en) * | 2006-02-27 | 2006-10-18 | 西安交通大学 | Analysis system and method of microbubble behavior in ultrasonic field based on superhigh speed photograph technology |
| US20080114242A1 (en) * | 2006-11-10 | 2008-05-15 | National Taiwan University | Ultrasonic imaging technique for differentiating the distribution of scatterers within a tissue |
| CN101623203A (en) * | 2009-08-07 | 2010-01-13 | 西安交通大学 | Multi-mode multi-parameter synchronous detection imaging monitoring system in transient physical process and monitoring method |
| CN103948402A (en) * | 2014-05-13 | 2014-07-30 | 中国科学院深圳先进技术研究院 | Tumor ultrasonic imaging feature extraction method and system |
| CN104398271A (en) * | 2014-11-14 | 2015-03-11 | 西安交通大学 | Method using three-dimensional mechanics and tissue specific imaging of blood vessels and plaques for detection |
-
2016
- 2016-05-17 CN CN201610326626.0A patent/CN106037815A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1847824A (en) * | 2006-02-27 | 2006-10-18 | 西安交通大学 | Analysis system and method of microbubble behavior in ultrasonic field based on superhigh speed photograph technology |
| US20080114242A1 (en) * | 2006-11-10 | 2008-05-15 | National Taiwan University | Ultrasonic imaging technique for differentiating the distribution of scatterers within a tissue |
| CN101623203A (en) * | 2009-08-07 | 2010-01-13 | 西安交通大学 | Multi-mode multi-parameter synchronous detection imaging monitoring system in transient physical process and monitoring method |
| CN103948402A (en) * | 2014-05-13 | 2014-07-30 | 中国科学院深圳先进技术研究院 | Tumor ultrasonic imaging feature extraction method and system |
| CN104398271A (en) * | 2014-11-14 | 2015-03-11 | 西安交通大学 | Method using three-dimensional mechanics and tissue specific imaging of blood vessels and plaques for detection |
Non-Patent Citations (4)
| Title |
|---|
| SIYUAN ZHANG 等: "Enhanced Lesion-to-Bubble Ratio on Ultrasonic Nakagami Imaging for Monitoring of High-Intensity Focused Ultrasound", 《J ULTRASOUND MED 2014》 * |
| SIYUAN ZHANG 等: "Feasibility of using Nakagami distribution in evaluating the formation of ultrasound-induced thermal lesions", 《2012 ACOUSTICAL SOCIETY OF AMERICA》 * |
| 宋晓阳: "《医学超声图像的轮廓波方法研究及其在相控HIFU治疗系统中的应用》", 《上海交通大学博士学位论文》 * |
| 张思远等: ""高强度聚焦超声热消融中参量监控成像研究"", 《中国超声医学工程学会超声诊疗、生物效应、仪器工程、重庆超声医学工程学会学术会议论文集》 * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108553763A (en) * | 2018-01-19 | 2018-09-21 | 北京工业大学 | A kind of microwave heat therapeutic monitoring method based on ultrasonic echo decorrelation imaging technique |
| CN112566694A (en) * | 2018-05-25 | 2021-03-26 | 卡尔迪亚韦弗公司 | Ultrasound treatment device comprising means for imaging cavitation bubbles |
| CN109171998A (en) * | 2018-10-22 | 2019-01-11 | 西安交通大学 | Heating ablation region recognition monitoring imaging method and system based on ultrasonic deep learning |
| US11058473B2 (en) | 2018-10-22 | 2021-07-13 | Xi'an Jiaotong University | Recognition, monitoring, and imaging method and system for thermal ablation region based on ultrasonic deep learning |
| CN113117258B (en) * | 2019-12-30 | 2023-04-28 | 重庆融海超声医学工程研究中心有限公司 | Detection device for tissue coagulation necrosis |
| CN113117258A (en) * | 2019-12-30 | 2021-07-16 | 重庆融海超声医学工程研究中心有限公司 | Device and method for detecting tissue coagulation necrosis |
| CN112966212A (en) * | 2021-02-07 | 2021-06-15 | 西安交通大学 | Multi-parameter real-time monitoring imaging system based on ultrasonic echo backscattering energy change |
| CN113907792A (en) * | 2021-08-31 | 2022-01-11 | 西安交通大学 | Nonlinear decorrelation fusion ultrasonic multi-parameter thermal coagulation detection and imaging method |
| CN113907792B (en) * | 2021-08-31 | 2022-08-05 | 西安交通大学 | Nonlinear decorrelation fusion ultrasonic multi-parameter thermal coagulation detection and imaging method |
| CN114305668A (en) * | 2021-12-22 | 2022-04-12 | 聚融医疗科技(杭州)有限公司 | Ultrasonic thermal ablation multi-parameter monitoring method and system based on demodulation domain parametric imaging |
| CN114305668B (en) * | 2021-12-22 | 2023-10-31 | 聚融医疗科技(杭州)有限公司 | Ultrasonic thermal ablation multi-parameter monitoring method and system based on demodulation domain parameter imaging |
| CN114587565A (en) * | 2022-03-01 | 2022-06-07 | 河南中医药大学 | Temperature control method and system in radio frequency ablation |
| CN114587565B (en) * | 2022-03-01 | 2023-09-29 | 河南中医药大学 | Temperature control method and system in radio frequency ablation |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106037815A (en) | Ultrasonic echo statistical parameter imaging system and method for thermal coagulation monitoring | |
| US11058399B2 (en) | Bubble-induced color doppler feedback during histotripsy | |
| US20070049827A1 (en) | Clinical feedback of ablation efficacy during ablation procedure | |
| DE69634714T2 (en) | Therapeutic ultrasound device | |
| US10226236B2 (en) | Ultrasound imaging system and method with automatic adjustment and/or multiple sample volumes | |
| JP5150267B2 (en) | Ultrasound diagnostic imaging system to detect liver damage | |
| JP4682149B2 (en) | Ultrasound imaging system and method for simultaneously displaying blood flow and perfusion parameters | |
| US20110208061A1 (en) | Ultrasonic lesion identification using temporal parametric contrast images | |
| CN101623203B (en) | Multi-mode multi-parameter synchronous detection imaging monitoring system in transient physical process and monitoring method | |
| CN107260217B (en) | The passive imaging method of three-dimensional and system for brain focused ultrasonic cavitation real time monitoring | |
| US10332250B2 (en) | Three-dimensional cavitation quantitative imaging method for microsecond-resolution cavitation spatial-temporal distribution | |
| EP3769690B1 (en) | Method for tissue characterization by ultrasound wave attenuation measurements and ultrasound system for tissue characterization | |
| Mast et al. | Ultrasound monitoring of in vitro radio frequency ablation by echo decorrelation imaging | |
| CN103202713B (en) | Image optimization method by blending of ultrasound fundamental wave and harmonic wave | |
| CN106264722A (en) | Window modulation combines Gauss polynomial matching monitoring radio frequency ablation device and method | |
| US20100125206A1 (en) | Methods for diagnosis and treatment of vessel occlusion | |
| JPS6355334B2 (en) | ||
| Gu et al. | Parametric perfusion imaging based on low-cost ultrasound platform | |
| Chen et al. | In vivo 3D power Doppler imaging using continuous translation and ultrafast ultrasound | |
| JP4909132B2 (en) | Optical tomography equipment | |
| CN112930145A (en) | Ultrasonic control unit | |
| Rahimian et al. | An acoustic backscatter-based method for estimating attenuation towards monitoring lesion formation in high intensity focused ultrasound | |
| Macé et al. | High sensitivity brain angiography using Ultrafast Doppler | |
| JP3534667B2 (en) | Ultrasonic measuring device | |
| US20210330226A1 (en) | Imaging method and imaging system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20161026 |
|
| RJ01 | Rejection of invention patent application after publication |