CN104825195A - Shear wave viscoelasticity imaging method and system - Google Patents
Shear wave viscoelasticity imaging method and system Download PDFInfo
- Publication number
- CN104825195A CN104825195A CN201510272479.9A CN201510272479A CN104825195A CN 104825195 A CN104825195 A CN 104825195A CN 201510272479 A CN201510272479 A CN 201510272479A CN 104825195 A CN104825195 A CN 104825195A
- Authority
- CN
- China
- Prior art keywords
- tissue
- value
- density
- shearing wave
- viscoelasticity
- 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
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computerised tomographs
-
- 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/0858—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving measuring tissue layers, e.g. skin, interfaces
-
- 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/5261—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 different diagnostic modalities, e.g. ultrasound and X-ray
Abstract
The invention provides a shear wave viscoelasticity imaging method and system. The method includes the steps that a CT value of tissue is obtained; according to a tissue mark, body model CT value and body model density value relation information corresponding to the tissue mark is obtained; a tissue density value corresponding to the CT value of the tissue is determined; shear waves are stimulated inside the tissue, ultrasonic waves are generated in the tissue through an ultrasonic transducer, and an ultrasonic wave return signal is received; according to the ultrasonic wave return signal, a propagation characteristic parameter of the shear waves in the tissue is obtained, and according to the propagation characteristic parameter and the tissue density value, the viscoelasticity parameter of the tissue is obtained; according to the viscoelasticity parameter, a shear wave viscoelasticity image of the tissue is formed. In this way, the effect of viscoelasticity parameter detection is improved, and more reliable references are provided for pathology judgment on the tissue by a doctor according to the shear wave viscoelasticity image of the tissue.
Description
Technical field
The invention belongs to ultrasonic therapy technical field, specifically relate to a kind of shearing wave viscoelasticity formation method and system.
Background technology
Shearing wave viscoelasticity imaging technique is detected by shearing wave to measure and the technology of imaging the viscoelastic parameters of tissue.
Elastic Modulus E and shearing wave spread speed V is organized according to viscoelasticity in prior art
srelation obtain the elastic modelling quantity of this viscoelasticity tissue:
wherein, ρ is tissue density, and then obtains the elastic information of this viscoelasticity tissue.But prior art exists a drawback: because tissue density values not easily obtains, prior art is all that (density value of such as liver is 1010kg/m to the fixing density value of viscoelasticity setting tissue one
3), but the tissue density of Different Individual there are differences (such as fatty liver, hepatic fibrosis), therefore this set constant density value method meeting and tissue actual density between there is deviation, thus affect the precision of tissue viscoelasticity's parameter measurement, after carrying out affecting viscoelastic parameters imaging, doctor is to the judgement of histopathology, delay treatment.
Summary of the invention
In order to solve at least one problem mentioned in background technology, the invention provides a kind of shearing wave viscoelasticity formation method and system, by obtaining the actual density value of individual tissue, obtain measurement result tissue viscoelasticity's parameter more accurately, overcome prior art and obtain tissue viscoelasticity's parameter by organizing constant density value, the technological deficiency that accuracy is not high, thus improve the effect of viscoelastic parameters detection, judge to provide reference more reliably for doctor carries out pathology according to the shearing wave viscoelasticity image of tissue to tissue.
The invention provides a kind of shearing wave viscoelasticity formation method, comprising:
Obtain the computer tomography CT value of tissue; According to organization identification, obtain the body mould CT value corresponding with described organization identification and body density of modes value relation information; Determine the density value of the described tissue corresponding to CT value of described tissue;
Shearing wave is excited at described organization internal;
In described tissue, produce ultrasound wave by ultrasonic transducer, and receive ultrasonic echo signal; Obtain the propagation characteristic parameter of described shearing wave in described tissue according to described ultrasonic echo signal, according to the density value of described propagation characteristic parameter and described tissue, obtain the viscoelastic parameters of described tissue; According to described viscoelastic parameters, form the shearing wave viscoelasticity image of described tissue.
The present invention also provides a kind of shearing wave viscoelasticity imaging system, comprising:
Computer tomography CT value detection device, control processor, shearing wave excitation apparatus, ultrasonic transducer, imaging device;
Described computer tomography CT value detection device is connected with described control processor, for obtain tissue computer tomography CT value and send the organization identification of described tissue, described tissue CT value to described control processor;
Described control processor, for according to described organization identification, obtains the body mould CT value corresponding with described organization identification and body density of modes value relation information; Determine the density value of the described tissue corresponding to CT value of described tissue;
Described shearing wave excitation apparatus is connected with described control processor, for exciting shearing wave at described organization internal; Described ultrasonic transducer is connected with described control processor, for producing ultrasound wave and receive ultrasonic echo signal in described tissue;
Accordingly, described control processor, also for obtaining the propagation characteristic parameter of described shearing wave in described tissue according to described ultrasonic echo signal, according to the density value of described propagation characteristic parameter and described tissue, obtains the viscoelastic parameters of described tissue;
Described imaging device is connected with described control processor, for according to described viscoelastic parameters, forms the shearing wave viscoelasticity image of described tissue.
Shearing wave viscoelasticity formation method provided by the invention and system, by obtaining the CT value of tissue, and obtaining the body mould CT value corresponding with organization identification and body density of modes value relation information, determining the density value of tissue; Again by inspiring shearing wave at organization internal; Sent by ultrasonic transducer and receive ultrasonic echo signal to obtain shearing wave propagation characteristic parameter in the tissue, thus the density value of conjunctive tissue, obtain the viscoelastic parameters of tissue, form the shearing wave viscoelasticity image of this tissue.By obtaining the actual density value of individual tissue, obtain measurement result tissue viscoelasticity's parameter more accurately, overcome prior art and obtain tissue viscoelasticity's parameter by organizing constant density value, the technological deficiency that accuracy is not high, thus improve the effect of viscoelastic parameters detection, judge to provide reference more reliably for doctor carries out pathology according to the shearing wave viscoelasticity image of tissue to tissue.
Accompanying drawing explanation
Fig. 1 is the flow chart of shearing wave viscoelasticity formation method embodiment one of the present invention;
Fig. 2 is the flow chart of shearing wave viscoelasticity formation method embodiment two of the present invention;
Fig. 3 is the flow chart of shearing wave viscoelasticity formation method embodiment three of the present invention;
Fig. 4 is the structural representation of shearing wave viscoelasticity imaging system embodiment one of the present invention.
Detailed description of the invention
Fig. 1 is the flow chart of shearing wave viscoelasticity formation method embodiment one of the present invention, and as shown in Figure 1, the method comprises:
(Computed Tomography is called for short: CT) value in the computed tomography of step 101, acquisition tissue.
The Flexible change of biological tissue is usually closely related with pathological phenomenon, and pathological tissues and normal structure there are differences in elasticity or hardness, the elasticity of the measurement biological tissue medically using elastic modelling quantity quantitative.And determine that an important parameter of biological tissue elasticity modulus is exactly the density value of biological tissue self, but Medical Technology can't carry out biopsy directly to obtain actual density value to biological tissue at present, therefore, be a fixed value for calculating the density value of biological tissue elasticity modulus in prior art, this density value just causing tissue is not for individual variation, and then the elastic modelling quantity precision utilizing this constant density value to calculate is not high, thus be unfavorable for that doctor makes correct pathology and judges.The present invention adopts the method indirectly obtaining biological tissue's actual density value, first by CT value detection device, CT image scan is carried out to tissue to be detected, obtain the CT value of tissue, again this CT value is sent in the control processor of shearing wave viscoelasticity imaging system and processes, to obtain the actual density value of this tissue.
Step 102, according to organization identification, obtain the body mould CT value corresponding with organization identification and body density of modes value relation information; Determine the density value of the tissue corresponding to CT value of tissue.
Store the relation information of CT value that not consubstantiality module knits and density value in control processor, body mould refers to the model simulated according to different biological tissue, and this model can the physical characteristic of the actual biological tissue of Complete Characterization.Such as: the body mould of liver, the body mould etc. of mammary gland.In addition, linear ratio relation is there is between the CT value of biological tissue and density value, that is, the density of the CT value reflection tissue of tissue, namely the CT value organized is higher, and to be equivalent to the density of tissue higher, then can be calculated the actual density value of tissue to be detected by the CT value of the known reference phantom of density and body mould.Suppose that the density of nominal volume mould is ρ
1, its CT value is H
1; The CT value of tissue to be detected is H, then the actual density ρ of tissue to be detected calculates by following formula: ρ=(H/H
1) × ρ
1, thus obtain the actual density value of tissue.Wherein, organization identification is for identifying the type of tissue to be detected, and the type of tissue can according to the type of histoorgan, as liver, mammary gland; Further, can also be subdivided into male organ, woman organ or the organ type etc. in all ages and classes stage in organ type, the body mould type corresponding to organization identification specifically can by technical staff's sets itself as required, and the application is not restricted this.
Step 103, excite shearing wave at organization internal.
Excite the mode of shearing wave a lot of at organization internal, such as: adopt oscillation generator to produce according to the driving excitation parameters arranged and drive pumping signal, this driving excitation parameters comprises: the waveform, dutycycle, amplitude, frequency etc. that drive pumping signal; Usual signal frequency is 20 ~ 5000Hz, preferably, and low frequency signal frequency 50Hz; Under the effect of this driving pumping signal, oscillation generator vibrates to tissue surface, shearing wave is inspired to make organization internal, be organized in this vibration direction of wave travel, namely go deep into axially being pushed of organization internal from tissue surface, make the shearing wave that inspires in tissue vertically, laterally (perpendicular to axial direction) direction propagate.Or launch ultrasound wave by ultrasonic transducer, produce acoustic radiation force with concentration ultrasonic at organization internal.Wherein, the hardness that shearing velocity of wave propagation can be organized or elasticity effect, hardness height is such as organized then to shear velocity of wave propagation fast, such as: the canceration of tissue can increase the elasticity (hardness) of tissue, therefore can by determining that tissue shears the pathological state that velocity of wave propagation comes diagnostic organization.
Step 104, produce ultrasound wave in the tissue by ultrasonic transducer, and receive ultrasonic echo signal.
Wave propagation information is sheared for getting tissue, adopt ultrasonic transducer to tissue emissions ultrasound wave, and the ultrasound echo signal that the ultrasound wave that collection is launched is reflected back from the tested tissue having shearing wave to propagate, shear velocity of wave propagation information to obtain from this echo-signal.
Step 105, obtain shearing wave propagation characteristic parameter in the tissue according to ultrasonic echo signal, according to the density value of propagation characteristic parameter and tissue, obtain the viscoelastic parameters of tissue.
The information that propagation characteristic parameter can be propagated based on the shearing wave got in echo-signal by empirical equation in tissue calculates, then calculates the viscoelastic parameters of tissue according to propagation characteristic parameter and the actual density value of organizing obtained before.
Step 106, according to viscoelastic parameters, the shearing wave viscoelasticity image of formative tissue.
According to the viscoelastic parameters got in the echo-signal that never shearing wave is propagated in the same time, then by color map, the viscoelastic parameters value of tissue is mapped to gray scale or coloured image, forms reflection and organize the viscoelastic image of shearing wave, thus location pathological changes.
The shearing wave viscoelasticity formation method of the present embodiment, by obtaining the CT value of tissue, and obtaining the body mould CT value corresponding with organization identification and body density of modes value relation information, determining the density value of tissue; Again by inspiring shearing wave at organization internal; Sent by ultrasonic transducer and receive ultrasonic echo signal to obtain shearing wave propagation characteristic parameter in the tissue, thus the density value of conjunctive tissue, obtain the viscoelastic parameters of tissue, form the shearing wave viscoelasticity image of this tissue.By obtaining the actual density value of individual tissue, obtain measurement result tissue viscoelasticity's parameter more accurately, overcome prior art and obtain tissue viscoelasticity's parameter by organizing constant density value, the technological deficiency that accuracy is not high, thus improve the effect of viscoelastic parameters detection, judge to provide reference more reliably for doctor carries out pathology according to the shearing wave viscoelasticity image of tissue to tissue.
On the basis of above-described embodiment, further, body mould CT value and body density of modes relation information are: the linear scale factor between body mould CT value and body density of modes value; Or be: body mould CT value and the corresponding numerical tabular between body density of modes value.Accordingly, step 102, according to organization identification, obtain the body mould CT value corresponding with organization identification and body density of modes value relation information; Determine the density value of the tissue corresponding to CT value of tissue, specifically comprise:
According to organization identification, obtain the linear scale factor between the body mould CT value corresponding with organization identification and body density of modes value, according to the CT value of tissue, calculate the density value of tissue.
Or, comprising:
According to organization identification, obtain the corresponding numerical tabular between the body mould CT value corresponding with organization identification with body density of modes value, according to the CT value of tissue, in numerical tabular, search the density value of tissue.
Concrete, linear ratio relation is there is between the CT value of biological tissue and density value, that is, the density of the CT value reflection tissue of tissue, namely the CT value organized is higher, and to be equivalent to the density of tissue higher, then can by the density of the CT Zhi Lai computation organization of the known reference phantom of density and body mould: suppose that the density of nominal volume mould is ρ
1, its CT value is H
1; The CT value of tissue to be detected is H, then the actual density ρ of tissue to be detected calculates by following formula: ρ=(H/H
1) × ρ
1, thus obtain the actual density value of tissue.Therefore, first according to organization identification, the body mould corresponding with organization identification can be got, such as: organization identification is women's liver mark, then correspond to the body mould information of women's liver, get the linear scale factor C between the CT value of this body mould and density value, i.e. C=ρ
1/ H
1; Again according to the CT value H of tissue actual measurement to be detected, calculate the actual density value ρ=C × H of this tissue.
Can also set up various body mould CT value and the corresponding numerical tabular of body density of modes value, by look-up table, obtain organizing and organize actual density value corresponding to CT value, look-up table can save the process that control processor carries out computing, speed up processing.In addition, can also based on the manual input density value of the empirical value of medical personnel etc.
Further, step 103, excite shearing wave at organization internal, specifically comprise:
By vibrator organizing the vibration of outer surface applying low frequency instant, excite shearing wave to make organization internal.
Concrete, vibrator can comprise: low-frequency pulse producer, amplifier; Low-frequency pulse producer is connected with control processor, receives the control signal that control processor sends; One end of amplifier is connected with low-frequency pulse producer, for amplifying the low frequency electrical number that low-frequency pulse producer produces; The other end of amplifier is connected with vibrator contact, vibrator contact is used for and organizes exterior surface, the vibration signal that low-frequency pulse producer sends amplifies through amplifier, then sends the wave of oscillation via vibrator contact and be transmitted into organization internal, makes organization internal oscillatory excitation go out shearing wave.Wherein, low-frequency pulse producer can be made up of sinusoidal signal generator circuit, also can be any circuit board that can send sinusoidal low-frequency pumping signal, as miniature sound card, the concrete formation of low-frequency pulse producer is by those skilled in the art according to existing generator designed, designed, and the application is not construed as limiting this.The low frequency signal frequency that low-frequency pulse producer produces and amplitude are controlled by control processor, and usual signal frequency is 20 ~ 5000Hz, preferably, and low frequency signal frequency 50Hz; The low frequency signal that amplifier can send low-frequency pulse producer carries out power amplification, makes it can inspire oscillatory shear ripple in tissue.Or vibrator can also comprise: low-frequency vibration driver element, motor; Motor is connected with low-frequency vibration driver element; Low-frequency vibration driver element is connected with control processor, for receiving the drive singal of control processor to make motor pulsation mechanical movement.
Or tissue surface being vibrated, to make organization internal excite shearing wave, specifically comprises: by microphone organizing outer surface generation sound wave, exciting shearing wave to make organization internal.
Concrete, by microphone organizing outer surface generation acoustic vibration, produce shearing wave with excite tissue internal motion.
Or produce acoustic radiation force with concentration ultrasonic at organization internal, to make organization internal excite shearing wave, specifically comprise: launch ultrasound wave by ultrasonic transducer, focusing ultrasonic wave is at organization internal, produce acoustic radiation force, acoustic radiation force promotes histokinesis, forms shearing wave; Wherein, this ultrasonic transducer and " step 104, producing ultrasound wave in the tissue by ultrasonic transducer, and receiving ultrasonic echo signal before." in ultrasonic transducer can be same or different ultrasonic transducers.
Concrete, the ultrasonic transducer producing acoustic radiation force both can use the ultrasonic transducer detecting shearing wave, also can use ultrasonic transducers different with it.Acoustic radiation force is the radiation power produced by the ultrasonic pulse of some strength, and it is launched into organization internal, produces displacement under making to be organized in the impetus of acoustic radiation force, forms shearing wave.
Fig. 2 is the flow chart of shearing wave viscoelasticity formation method embodiment two of the present invention, on the basis of above-described embodiment, further, as shown in Figure 2, step 105, obtains shearing wave propagation characteristic parameter in the tissue specifically comprise according to ultrasonic echo signal:
Step 1051, signal processing is carried out to ultrasonic echo signal, determine the kinematic parameter of tissue.
Concrete, kinematic parameter comprises: be organized in the displacement information under shearing wave effect and/or strain information.Wherein, signal processing comprises: from time domain cross-correlation, spectrum cross-correlation, square error and, trackings of SPECKLE TRACKING, scale invariant feature point, dynamic programming, zero crossing tracking, at least one in peak value searching.Meanwhile, echo-signal also needs to carry out the process such as filtering, amplification, digitized and beam shaping, to ensure the precision of kinematic parameter.
Step 1052, obtain the propagation characteristic parameter of shearing wave according to beginning parameter transform model.
Concrete, propagation characteristic parameter comprises: shear velocity of wave propagation and/or propagation attenuation coefficient.According to the displacement produced under being organized in shearing wave effect or strain (deformation), rule of thumb formulae discovery obtains the propagation characteristic parameter of shearing wave, velocity of wave propagation and/or propagation attenuation coefficient substitution empirical equation will be sheared again, obtain the viscoelastic parameters of tissue, thus obtain the pathological information of tissue.Wherein, displacement obtains by the integration of Tissue velocity and time, and strain can obtain by carrying out differential to displacement of tissue.Such as: the shearing wave spread speed obtaining tissue is V
s, and be ρ according to the density that CT value obtains tissue, then according to Young's modulus
obtain the viscoelastic parameters E of this viscoelasticity tissue.
Fig. 3 is the flow chart of shearing wave viscoelasticity formation method embodiment three of the present invention, and on the basis of above-described embodiment, further, as shown in Figure 3, step 106, according to viscoelastic parameters, the shearing wave viscoelasticity image of formative tissue, specifically comprises:
Step 1061, basis preset the mapping relations of viscoelastic parameters and color, search the color coding that viscoelastic parameters is corresponding.
Concrete, viscoelastic parameters comprises: modulus of shearing, Young's modulus, shearing elasticity, shear viscosity, mechanical impedance, mechanical relaxation time, anisotropy etc.Set up mapping relations one by one by between default viscoelastic parameters and color, thus find the color-coded information corresponding to it according to different viscoelastic parameters.
The shearing wave viscoelasticity image of step 1062, formative tissue.
Concrete, mapping comprises: grey scale mapping, color mapped.Utilize color map, the viscoelastic parameters of tissue is mapped to the color coding of gray scale or coloured image, thus the elastic information at each focus place of clear display organization in the form of images.
Fig. 4 is the structural representation of shearing wave viscoelasticity imaging system embodiment one of the present invention, and as shown in Figure 4, this system comprises:
Computer tomography CT value detection device 1, control processor 2, shearing wave excitation apparatus 3, ultrasonic transducer 4, imaging device 5;
Wherein, computer tomography CT value detection device 1 is connected with control processor 2, for obtain tissue computer tomography CT value and send tissue organization identification, tissue CT value to control processor 2;
Control processor 2, for according to organization identification, obtains the body mould CT value corresponding with organization identification and body density of modes value relation information; Determine the density value of the tissue corresponding to CT value of tissue;
Shearing wave excitation apparatus 3 is connected with control processor 2, for exciting shearing wave at organization internal;
Ultrasonic transducer 4 is connected with control processor 2, for producing ultrasound wave in the tissue and receiving ultrasonic echo signal;
Accordingly, control processor 2, also for obtaining shearing wave propagation characteristic parameter in the tissue according to ultrasonic echo signal, according to the density value of propagation characteristic parameter and tissue, obtains the viscoelastic parameters of tissue;
Imaging device 5 is connected with control processor 2, for according to viscoelastic parameters, and the shearing wave viscoelasticity image of formative tissue.
The shearing wave viscoelasticity imaging system of the present embodiment may be used for the technical scheme performing embodiment of the method shown in Fig. 1, and it realizes principle and technique effect is similar, repeats no more herein.
On the basis of above-described embodiment, further, body mould CT value and body density of modes relation information are: the linear scale factor between body mould CT value and body density of modes value; Or be: body mould CT value and the corresponding numerical tabular between body density of modes value.
Control processor 2, specifically for according to organization identification, obtains the linear scale factor between the body mould CT value corresponding with organization identification and body density of modes value, according to the CT value of tissue, calculates the density value of tissue; Or specifically for according to organization identification, obtain the corresponding numerical tabular between the body mould CT value corresponding with organization identification with body density of modes value, according to the CT value of tissue, in numerical tabular, search the density value of tissue.
Further, shearing wave excitation apparatus 3 is specially vibrator, for organizing outer surface to apply low frequency instant vibration, excites shearing wave to make organization internal; Or shearing wave excitation apparatus 3 is specially microphone, for organizing outer surface to produce sound wave, excite shearing wave to make organization internal; Or shearing wave excitation apparatus 3 is specially ultrasonic transducer, for launching ultrasound wave, focusing ultrasonic wave, at organization internal, produces acoustic radiation force, and acoustic radiation force promotes histokinesis, forms shearing wave; Wherein, this ultrasonic transducer and ultrasonic transducer before 4 can be same or different ultrasonic transducers.
Control processor 2, specifically for carrying out signal processing to ultrasonic echo signal, determine the kinematic parameter of tissue, wherein, kinematic parameter comprises: be organized in the displacement information under shearing wave effect and/or strain information; The propagation characteristic parameter of shearing wave is obtained according to beginning parameter transform model; Propagation characteristic parameter comprises: shear velocity of wave propagation and/or propagation attenuation coefficient; Signal processing comprises: from time domain cross-correlation, spectrum cross-correlation, square error and, trackings of SPECKLE TRACKING, scale invariant feature point, dynamic programming, zero crossing tracking, at least one in peak value searching.
Imaging device 5, specifically for according to the mapping relations presetting viscoelastic parameters and color, searches the color coding that viscoelastic parameters is corresponding, the shearing wave viscoelasticity image of formative tissue; Wherein, mapping comprises: grey scale mapping, color mapped; Viscoelastic parameters comprises: modulus of shearing, Young's modulus, shearing elasticity, shear viscosity, mechanical impedance, mechanical relaxation time, anisotropy.
The shearing wave viscoelasticity imaging system of the present embodiment may be used for the technical scheme performing embodiment of the method shown in Fig. 1 ~ 3, and it realizes principle and technique effect is similar, repeats no more herein.
One of ordinary skill in the art will appreciate that: all or part of step realizing said method embodiment can have been come by the hardware that programmed instruction is relevant, aforesaid program can be stored in a computer read/write memory medium, this program, when performing, performs the step comprising said method embodiment; And aforesaid storage medium comprises: ROM, RAM, magnetic disc or CD etc. various can be program code stored medium.
Last it is noted that above each embodiment is only in order to illustrate technical scheme of the present invention, be not intended to limit; Although with reference to foregoing embodiments to invention has been detailed description, those of ordinary skill in the art is to be understood that: it still can be modified to the technical scheme described in foregoing embodiments, or carries out equivalent replacement to wherein some or all of technical characteristic; And these amendments or replacement, do not make the essence of appropriate technical solution depart from the scope of various embodiments of the present invention technical scheme.
Claims (10)
1. a shearing wave viscoelasticity formation method, is characterized in that, comprising:
Obtain the computer tomography CT value of tissue; According to organization identification, obtain the body mould CT value corresponding with described organization identification and body density of modes value relation information; Determine the density value of the described tissue corresponding to CT value of described tissue;
Shearing wave is excited at described organization internal;
In described tissue, produce ultrasound wave by ultrasonic transducer, and receive ultrasonic echo signal; Obtain the propagation characteristic parameter of described shearing wave in described tissue according to described ultrasonic echo signal, according to the density value of described propagation characteristic parameter and described tissue, obtain the viscoelastic parameters of described tissue; According to described viscoelastic parameters, form the shearing wave viscoelasticity image of described tissue.
2. method according to claim 1, is characterized in that, described relation information is described body mould CT value and the linear scale factor between described body density of modes value or the corresponding numerical tabular between described body mould CT value with described body density of modes value.
3. method according to claim 2, is characterized in that, described according to organization identification, obtains the body mould CT value corresponding with described organization identification and body density of modes value relation information; Determine that the density value of the described tissue corresponding to CT value of described tissue comprises:
According to described organization identification, obtain the linear scale factor between the body mould CT value corresponding with described organization identification and body density of modes value, according to the CT value of described tissue, calculate the density value of described tissue; Or, according to described organization identification, obtain the corresponding numerical tabular between the body mould CT value corresponding with described organization identification with body density of modes value, according to the CT value of described tissue, in described numerical tabular, search the density value of described tissue.
4. method according to claim 1, is characterized in that, describedly obtains the propagation characteristic parameter of described shearing wave in described tissue according to described ultrasonic echo signal and comprises:
Carry out signal processing to described ultrasonic echo signal, determine the kinematic parameter of described tissue, described kinematic parameter comprises: described in be organized in displacement information under described shearing wave effect and/or strain information; The propagation characteristic parameter of described shearing wave is obtained according to described beginning parameter transform model; Described propagation characteristic parameter comprises: described shearing velocity of wave propagation and/or propagation attenuation coefficient;
Described signal processing comprises: from time domain cross-correlation, spectrum cross-correlation, square error and, trackings of SPECKLE TRACKING, scale invariant feature point, dynamic programming, zero crossing tracking, at least one in peak value searching.
5. method according to claim 1, is characterized in that, described according to described viscoelastic parameters, and the shearing wave viscoelasticity image forming described tissue comprises:
According to the mapping relations of default viscoelastic parameters and color, search the color coding that described viscoelastic parameters is corresponding, form the shearing wave viscoelasticity image of described tissue; Described mapping comprises: grey scale mapping, color mapped;
Described viscoelastic parameters comprises: modulus of shearing, Young's modulus, shearing elasticity, shear viscosity, mechanical impedance, mechanical relaxation time, anisotropy.
6. a shearing wave viscoelasticity imaging system, is characterized in that, comprising:
Computer tomography CT value detection device, control processor, shearing wave excitation apparatus, ultrasonic transducer, imaging device;
Described computer tomography CT value detection device is connected with described control processor, for obtain tissue computer tomography CT value and send the organization identification of described tissue, described tissue CT value to described control processor;
Described control processor, for according to described organization identification, obtains the body mould CT value corresponding with described organization identification and body density of modes value relation information; Determine the density value of the described tissue corresponding to CT value of described tissue;
Described shearing wave excitation apparatus is connected with described control processor, for exciting shearing wave at described organization internal; Described ultrasonic transducer is connected with described control processor, for producing ultrasound wave and receive ultrasonic echo signal in described tissue;
Accordingly, described control processor, also for obtaining the propagation characteristic parameter of described shearing wave in described tissue according to described ultrasonic echo signal, according to the density value of described propagation characteristic parameter and described tissue, obtains the viscoelastic parameters of described tissue;
Described imaging device is connected with described control processor, for according to described viscoelastic parameters, forms the shearing wave viscoelasticity image of described tissue.
7. system according to claim 6, is characterized in that,
Described relation information is described body mould CT value and the linear scale factor between described body density of modes value or the corresponding numerical tabular between described body mould CT value with described body density of modes value.
8. system according to claim 6, is characterized in that,
Described control processor, specifically for according to described organization identification, obtains the linear scale factor between the body mould CT value corresponding with described organization identification and body density of modes value, according to the CT value of described tissue, calculates the density value of described tissue; Or specifically for according to described organization identification, obtain the corresponding numerical tabular between the body mould CT value corresponding with described organization identification with body density of modes value, according to the CT value of described tissue, in described numerical tabular, search the density value of described tissue.
9. system according to claim 6, is characterized in that,
Described control processor, specifically for carrying out signal processing to described ultrasonic echo signal, determine the kinematic parameter of described tissue, described kinematic parameter comprises: described in be organized in displacement information under described shearing wave effect and/or strain information; The propagation characteristic parameter of described shearing wave is obtained according to described beginning parameter transform model; Described propagation characteristic parameter comprises: described shearing velocity of wave propagation and/or propagation attenuation coefficient; Described signal processing comprises: from time domain cross-correlation, spectrum cross-correlation, square error and, trackings of SPECKLE TRACKING, scale invariant feature point, dynamic programming, zero crossing tracking, at least one in peak value searching.
10. system according to claim 6, is characterized in that,
Described imaging device, specifically for according to the mapping relations presetting viscoelastic parameters and color, searches the color coding that described viscoelastic parameters is corresponding, forms the shearing wave viscoelasticity image of described tissue; Described mapping comprises: grey scale mapping, color mapped; Described viscoelastic parameters comprises: modulus of shearing, Young's modulus, shearing elasticity, shear viscosity, mechanical impedance, mechanical relaxation time, anisotropy.
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510272479.9A CN104825195A (en) | 2015-05-25 | 2015-05-25 | Shear wave viscoelasticity imaging method and system |
| HK16101282.7A HK1213452A1 (en) | 2015-05-25 | 2016-02-03 | Method and system for shear wave viscoelastic imaging |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510272479.9A CN104825195A (en) | 2015-05-25 | 2015-05-25 | Shear wave viscoelasticity imaging method and system |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN104825195A true CN104825195A (en) | 2015-08-12 |
Family
ID=53803678
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510272479.9A Pending CN104825195A (en) | 2015-05-25 | 2015-05-25 | Shear wave viscoelasticity imaging method and system |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN104825195A (en) |
| HK (1) | HK1213452A1 (en) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105212968A (en) * | 2015-10-29 | 2016-01-06 | 无锡海斯凯尔医学技术有限公司 | Elastomeric check method and apparatus |
| CN105455851A (en) * | 2015-12-24 | 2016-04-06 | 无锡海斯凯尔医学技术有限公司 | Method and equipment for detecting viscoelasticity parameters of viscoelastic media |
| CN105910956A (en) * | 2016-04-15 | 2016-08-31 | 中国科学院上海硅酸盐研究所 | Industrial CT analysis method of object density and density distribution |
| WO2017121170A1 (en) * | 2016-01-11 | 2017-07-20 | 无锡海斯凯尔医学技术有限公司 | Tissue parameter detection method and system |
| CN107095692A (en) * | 2016-02-19 | 2017-08-29 | 乐普（北京）医疗器械股份有限公司 | Ultrasonic imaging system, method for ultrasonic imaging and one dimension displacement scan method |
| CN107510474A (en) * | 2017-09-21 | 2017-12-26 | 深圳开立生物医疗科技股份有限公司 | Shearing wave elastograph imaging method and system |
| WO2019015396A1 (en) * | 2017-07-21 | 2019-01-24 | 无锡海斯凯尔医学技术有限公司 | Method and device for measuring viscoelasticity of medium |
| CN110292399A (en) * | 2018-05-04 | 2019-10-01 | 深圳迈瑞生物医疗电子股份有限公司 | A kind of method and system of shearing wave elasticity measurement |
| CN111728642A (en) * | 2017-06-06 | 2020-10-02 | 深圳迈瑞生物医疗电子股份有限公司 | Method, device and system for imaging in ultrasonic scanning |
| CN112135567A (en) * | 2018-05-03 | 2020-12-25 | 皇家飞利浦有限公司 | Shear wave amplitude reconstruction for tissue elasticity monitoring and display |
-
2015
- 2015-05-25 CN CN201510272479.9A patent/CN104825195A/en active Pending
-
2016
- 2016-02-03 HK HK16101282.7A patent/HK1213452A1/en unknown
Cited By (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2688294C1 (en) * | 2015-10-29 | 2019-05-21 | Вуси Хиски Медикал Текнолоджис Ко., Лтд. | Method and device for determining elasticity |
| WO2017071605A1 (en) * | 2015-10-29 | 2017-05-04 | 无锡海斯凯尔医学技术有限公司 | Elasticity detection method and device |
| US11426144B2 (en) | 2015-10-29 | 2022-08-30 | Wuxi Hisky Medical Technologies Co., Ltd. | Method and device for elasticity detection |
| CN105212968A (en) * | 2015-10-29 | 2016-01-06 | 无锡海斯凯尔医学技术有限公司 | Elastomeric check method and apparatus |
| CN105455851A (en) * | 2015-12-24 | 2016-04-06 | 无锡海斯凯尔医学技术有限公司 | Method and equipment for detecting viscoelasticity parameters of viscoelastic media |
| WO2017107660A1 (en) * | 2015-12-24 | 2017-06-29 | 无锡海斯凯尔医学技术有限公司 | Method and apparatus for measuring viscoelastic parameter of viscoelastic medium |
| KR101914021B1 (en) | 2015-12-24 | 2018-10-31 | 우시 히스키 메디칼 테크놀로지스 컴퍼니., 리미티드. | Method and Apparatus for Detecting Viscoelastic Parameters of Viscoelastic Media |
| US10420528B2 (en) | 2015-12-24 | 2019-09-24 | Wuxi Hisky Medical Technologies Co., Ltd. | Method and device for detecting viscoelastic parameter of viscoelastic medium |
| WO2017121170A1 (en) * | 2016-01-11 | 2017-07-20 | 无锡海斯凯尔医学技术有限公司 | Tissue parameter detection method and system |
| CN107095692B (en) * | 2016-02-19 | 2024-02-09 | 乐普（北京）医疗器械股份有限公司 | Ultrasonic imaging system, ultrasonic imaging method and one-dimensional displacement scanning method |
| CN107095692A (en) * | 2016-02-19 | 2017-08-29 | 乐普（北京）医疗器械股份有限公司 | Ultrasonic imaging system, method for ultrasonic imaging and one dimension displacement scan method |
| CN105910956B (en) * | 2016-04-15 | 2019-02-26 | 中国科学院上海硅酸盐研究所 | A kind of method of industry CT object analysis density and Density Distribution |
| CN105910956A (en) * | 2016-04-15 | 2016-08-31 | 中国科学院上海硅酸盐研究所 | Industrial CT analysis method of object density and density distribution |
| CN111728642A (en) * | 2017-06-06 | 2020-10-02 | 深圳迈瑞生物医疗电子股份有限公司 | Method, device and system for imaging in ultrasonic scanning |
| CN111728642B (en) * | 2017-06-06 | 2021-03-19 | 深圳迈瑞生物医疗电子股份有限公司 | Method, device and system for imaging in ultrasonic scanning |
| US11717270B2 (en) | 2017-06-06 | 2023-08-08 | Shenzhen Mindray Bio-Medical Electronics Co., Ltd. | Method, apparatus and system for imaging in ultrasonic scanning |
| WO2019015396A1 (en) * | 2017-07-21 | 2019-01-24 | 无锡海斯凯尔医学技术有限公司 | Method and device for measuring viscoelasticity of medium |
| CN107510474B (en) * | 2017-09-21 | 2020-07-10 | 深圳开立生物医疗科技股份有限公司 | Shear wave elastic imaging method and system |
| CN107510474A (en) * | 2017-09-21 | 2017-12-26 | 深圳开立生物医疗科技股份有限公司 | Shearing wave elastograph imaging method and system |
| CN112135567A (en) * | 2018-05-03 | 2020-12-25 | 皇家飞利浦有限公司 | Shear wave amplitude reconstruction for tissue elasticity monitoring and display |
| CN110292399A (en) * | 2018-05-04 | 2019-10-01 | 深圳迈瑞生物医疗电子股份有限公司 | A kind of method and system of shearing wave elasticity measurement |
| CN110292399B (en) * | 2018-05-04 | 2022-03-08 | 深圳迈瑞生物医疗电子股份有限公司 | Method and system for measuring shear wave elasticity |
Also Published As
| Publication number | Publication date |
|---|---|
| HK1213452A1 (en) | 2016-07-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104825195A (en) | Shear wave viscoelasticity imaging method and system | |
| CN106175831B (en) | A kind of tissue viscoelasticity's measurement method based on shearing wave amplitude and phase-detection | |
| AU2016347143B2 (en) | Elasticity detection method and device | |
| Zhang et al. | Noninvasive method for estimation of complex elastic modulus of arterial vessels | |
| Wells et al. | Medical ultrasound: imaging of soft tissue strain and elasticity | |
| CN102641137B (en) | Measurement of viscoelastic property using amplitude-phase modulation | |
| JP4990963B2 (en) | Method and apparatus for imaging a viscoelastic medium | |
| JP5485508B2 (en) | Method and apparatus for improved ultrasonic distortion measurement of soft tissue | |
| EP3231369A1 (en) | Ultrasound diagnostic device and elasticity evaluation method | |
| CN105662473A (en) | Tissue parameter detection method and system | |
| CN105266851A (en) | Method and device for detecting viscosity parameter of viscoelastic medium | |
| US11439367B2 (en) | Hybrid elastography Method, probe and device for hybrid elastography | |
| CN104771192A (en) | Method for processing form and elasticity information of tissue and elasticity detection apparatus | |
| Collin et al. | Quantitative observations of cavitation activity in a viscoelastic medium | |
| CN105147327A (en) | Generation method for multi-frequency shear wave | |
| CN204600529U (en) | Elastomeric check equipment | |
| TW201907871A (en) | Method and device for measuring viscoelasticity of medium | |
| JP4381118B2 (en) | Ultrasonic diagnostic equipment | |
| US11808898B2 (en) | Method for measuring an ultrasonic attenuation parameter guided by harmonic elastography, probe and device for the implementation of the method | |
| RU2794039C2 (en) | Hybrid elastography method, probe and device for hybrid elastography | |
| RU2786332C2 (en) | Method for measurement of parameter of ultrasound attenuation by harmonic elastography, probe and device for implementation of method | |
| CN114869329A (en) | Elastic imaging method and system | |
| Nenadic et al. | Measuring bladder viscoelasticity using ultrasound | |
| CN113749687A (en) | Test apparatus, test method, and test phantom | |
| Zhang et al. | Generation of ring resonant mode for measuring artery elastic modulus |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| EXSB | Decision made by sipo to initiate substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| REG | Reference to a national code |
Ref country code: HK Ref legal event code: DE Ref document number: 1213452 Country of ref document: HK |
|
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20150812 |
|
| REG | Reference to a national code |
Ref country code: HK Ref legal event code: WD Ref document number: 1213452 Country of ref document: HK |