CN110477948A - Elastograph imaging method and device, imaging device, storage medium - Google Patents
Elastograph imaging method and device, imaging device, storage medium Download PDFInfo
- Publication number
- CN110477948A CN110477948A CN201910773249.9A CN201910773249A CN110477948A CN 110477948 A CN110477948 A CN 110477948A CN 201910773249 A CN201910773249 A CN 201910773249A CN 110477948 A CN110477948 A CN 110477948A
- Authority
- CN
- China
- Prior art keywords
- data
- matrix
- cross
- correlation coefficient
- biological tissue
- 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.)
- Granted
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/48—Diagnostic techniques
- A61B8/485—Diagnostic techniques involving measuring strain or elastic properties
-
- 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
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 present invention relates to elastograph imaging method and device, electronic equipment, storage mediums.Wherein, elastograph imaging method is used in the imaging device of elastogram system, and elastogram system further includes probe, pops one's head in for emitting rf wave, elastograph imaging method includes: multiframe I/Q data corresponding to echo-signal of the biological tissue based on radio-wave reflection during acquisition biological tissue is extruded；Calculate the first cross-correlation coefficient of adjacent two frames I/Q data；A frame speed data are calculated according to the first cross-correlation coefficient；The VELOCITY DISTRIBUTION data of biological tissue are obtained according to multiframe I/Q data, VELOCITY DISTRIBUTION data include multiframe speed data；Diagram of strains is generated according to VELOCITY DISTRIBUTION data.Present invention combination Doppler effect and cross correlation algorithm calculate the VELOCITY DISTRIBUTION during biological tissue is extruded, generate the diagram of strains of biological tissue, realize elastogram, and precision is high, speed is fast, can improve healthcare givers to the diagnostic reliability and efficiency of disease.
Description
Technical field
The present invention relates to technical field of medical equipment more particularly to a kind of elastograph imaging method and device, imaging device, deposit
Storage media.
Background technique
Ultrasonic Elasticity Imaging has become a research hotspot of medical ultrasound image, is widely used in breast, preceding
The organellas such as column gland, thyroid gland imaging, the basic principle of elastogram: to tissue apply a dynamic either internally or externally or
Static/quasi-static excitation, under the effect of the physics laws such as Elasticity, biomethanics, tissue will generate a response.Simply
For, elasticity modulus larger (i.e. harder tissue) strain is smaller, or the amplitude vibrated is smaller or speed is larger, utilizes ultrasound
Imaging is estimated that the corresponding situation of organization internal, to reflect the difference of the mechanical attributes such as the elasticity modulus of organization internal.
Currently, elastogram is mainly based upon the strain pattern algorithm of Displacement Estimation, by calculating adjacent two frames RF data
Between displacement, and then calculate Strain Distribution image.But the precision of the algorithm is low, and it is affected by noise serious, cause into image quality
It measures poor.
Summary of the invention
The present invention provides a kind of elastograph imaging method and device, electronic equipment, storage medium, to solve the bullet of the prior art
Property the imaging algorithm precision that uses of imaging it is low, it is affected by noise serious, cause the defect that image quality is poor.
The present invention is to solve above-mentioned technical problem by following technical proposals:
According to the first aspect of the invention, a kind of elastograph imaging method, the elastograph imaging method are used in elastogram system
The imaging device of system, the elastogram system further include probe, and the probe is for emitting rf wave, the elastogram side
Method includes:
During acquisition biological tissue is extruded, echo-signal institute of the biological tissue based on the radio-wave reflection is right
The multiframe I/Q data answered；
Calculate the first cross-correlation coefficient of adjacent two frames I/Q data；
A frame speed data are calculated according to first cross-correlation coefficient；
The VELOCITY DISTRIBUTION data of the biological tissue are obtained according to the multiframe I/Q data, the VELOCITY DISTRIBUTION data include
Multiframe speed data；
Diagram of strains is generated according to the VELOCITY DISTRIBUTION data.
Optionally, the first cross-correlation coefficient of adjacent two frames I/Q data is calculated, comprising:
The first IQ matrix and the 2nd IQ matrix are generated according to the adjacent two frames I/Q data；
The conjugate matrices of the first IQ matrix and the 2nd IQ matrix are carried out respectively based on convolution algorithm window
Convolution algorithm；
Product calculation will be done by the corresponding element of the conjugate matrices of the first IQ matrix and the 2nd IQ matrix of convolution algorithm,
Obtain first cross-correlation coefficient.
Optionally, a frame speed data are calculated according to first cross-correlation coefficient, comprising:
The first phase difference of the adjacent two frames I/Q data is calculated according to first cross-correlation coefficient；
The speed data is calculated according to the first phase difference and the frequency of the rf wave.
Optionally, the speed data is calculated according to the first phase difference and the frequency of the rf wave, comprising:
The second cross-correlation coefficient of the submatrix of the different offset rows of the first IQ matrix is calculated, and calculates described second
The third cross-correlation coefficient of the submatrix of the different offset rows of IQ matrix；
Sum operation is carried out to second cross-correlation coefficient and the third cross-correlation coefficient, and according to sum operation
As a result it is poor to calculate second phase；
The frequency is updated according to the second phase difference；
The speed data is calculated according to the first phase difference and updated frequency.
Optionally, diagram of strains is generated according to the VELOCITY DISTRIBUTION data, comprising:
B mould image is generated according to the multiframe I/Q data；
Strain Distribution data are generated according to the VELOCITY DISTRIBUTION data, and the Strain Distribution data are carried out at imaging
Reason generates diagram of strains；
According to the dimension difference between the B mould image and the diagram of strains, scales the B mould image or described answer
Variation Butut, so that the size of the B mould image is identical as the size of the diagram of strains；
Merge the B mould image and diagram of strains Jing Guo size adjusting.
According to the second aspect of the invention, a kind of elastogram device is provided, the elastogram device is applied in elasticity
The imaging device of imaging system, the elastogram system further include probe, and the probe is for emitting rf wave, the elasticity
Imaging device includes:
First data acquisition module, for obtaining during biological tissue is extruded, the biological tissue is based on described penetrate
Multiframe I/Q data corresponding to the echo-signal of frequency wave reflection；
First computing module, for calculating the first cross-correlation coefficient of adjacent two frames I/Q data；
Second computing module, for calculating a frame speed data according to first cross-correlation coefficient；
Second data acquisition module, for obtaining the VELOCITY DISTRIBUTION number of the biological tissue according to the multiframe I/Q data
According to the VELOCITY DISTRIBUTION data include multiframe speed data；
Image generation module, for generating diagram of strains according to the VELOCITY DISTRIBUTION data.
Optionally, first computing module is specifically used for:
The first IQ matrix and the 2nd IQ matrix are generated according to the adjacent two frames I/Q data；
The conjugate matrices of the first IQ matrix and the 2nd IQ matrix are carried out respectively based on convolution algorithm window
Convolution algorithm；
Product calculation will be done by the corresponding element of the conjugate matrices of the first IQ matrix and the 2nd IQ matrix of convolution algorithm,
Obtain first cross-correlation coefficient.
Optionally, second computing module is specifically used for:
The first phase difference of adjacent two frames I/Q data is calculated according to first cross-correlation coefficient；
The speed data is calculated according to the first phase difference and the frequency of rf wave.
Optionally, the number of speed is calculated according to the first phase difference and the frequency of rf wave in second computing module
According to when, be specifically used for:
The second cross-correlation coefficient of the submatrix of the different offset rows of the first IQ matrix is calculated, and calculates described second
The third cross-correlation coefficient of the submatrix of the different offset rows of IQ matrix；
Sum operation is carried out to second cross-correlation coefficient and the third cross-correlation coefficient, and according to sum operation
As a result it is poor to calculate second phase；
The frequency is updated according to the second phase difference；
The speed data is calculated according to the first phase difference and updated frequency.
Optionally, described image generation module is specifically used for:
B mould image is generated according to the multiframe I/Q data；
Strain Distribution data are generated according to the VELOCITY DISTRIBUTION data, and the Strain Distribution data are carried out at imaging
Reason generates initial strain distribution map；
According to the dimension difference between the B mould image and the initial strain distribution map, the B mould image or institute are scaled
Initial strain distribution map is stated, so that the size of the B mould image is identical as the size of the initial strain distribution map；
Merge the B mould image and initial strain distribution map Jing Guo size adjusting, generates final diagram of strains.
According to the third aspect of the invention we, a kind of imaging device is provided, including memory, processor and is stored in memory
Computer program that is upper and can running on a processor, the processor realize any of the above-described when executing the computer program
The elastograph imaging method.
According to the fourth aspect of the invention, a kind of computer readable storage medium is stored thereon with computer program, described
The step of elastograph imaging method described in any of the above embodiments is realized when computer program is executed by processor.
The positive effect of the present invention is that: present invention combination Doppler effect and cross correlation algorithm calculate biological group
The VELOCITY DISTRIBUTION during being extruded to be knitted, the diagram of strains of biological tissue is generated, realizes elastogram, precision is high, speed is fast,
Healthcare givers can be improved to the diagnostic reliability and efficiency of disease.
It should be understood that above general description and following detailed description be only it is exemplary and explanatory, not
It can the limitation present invention.
Detailed description of the invention
The drawings herein are incorporated into the specification and forms part of this specification, and shows and meets implementation of the invention
Example, and be used to explain the principle of the present invention together with specification.
Fig. 1 is Doppler effect principle explanatory diagram.
Fig. 2 is the first pass figure of the elastograph imaging method of inventive embodiments 1.
Fig. 3 is the second flow chart of the elastograph imaging method of inventive embodiments 1.
Fig. 4 a is the first schematic diagram of implementation process of step 102 in Fig. 2.
Fig. 4 b is the second schematic diagram of implementation process of step 102 in Fig. 2.
Fig. 5 is the flow chart of the elastograph imaging method of inventive embodiments 2.
Fig. 6 is the module diagram of the elastogram device of inventive embodiments 3.
Fig. 7 is the structural schematic diagram of the imaging device of the embodiment of the present invention 5.
Specific embodiment
Example embodiments are described in detail here, and the example is illustrated in the accompanying drawings.Following description is related to
When attached drawing, unless otherwise indicated, the same numbers in different drawings indicate the same or similar elements.Following exemplary embodiment
Described in embodiment do not represent all embodiments consistented with the present invention.On the contrary, they be only with it is such as appended
The example of device and method being described in detail in claims, some aspects of the invention are consistent.
It is only to be not intended to limit the invention merely for for the purpose of describing particular embodiments in terminology used in the present invention.
It is also intended in the present invention and the "an" of singular used in the attached claims, " described " and "the" including majority
Form, unless the context clearly indicates other meaning.It is also understood that term "and/or" used herein refers to and wraps
It may be combined containing one or more associated any or all of project listed.
It will be appreciated that though various information, but this may be described using term first, second, third, etc. in the present invention
A little information should not necessarily be limited by these terms.These terms are only used to for same type of information being distinguished from each other out.For example, not departing from
In the case where the scope of the invention, the first information can also be referred to as the second information, and similarly, the second information can also be referred to as
One information.Depending on context, word as used in this " if " can be construed to " ... when " or " when ...
When " or " in response to determination ".
The present embodiment provides a kind of elastograph imaging methods, used in the imaging device of elastogram system, the elastogram system
System further includes probe, and the probe is for emitting rf wave.This method combination Doppler effect and cross correlation algorithm generate biological group
The diagram of strains knitted realizes elastogram, provides reference to the diagnosis of disease for healthcare givers.This implementation in order to facilitate understanding
The elastograph imaging method of example, is first briefly described its theoretical foundation.
About Doppler effect, referring to Fig. 1, the sensor in Fig. 1 includes transmitting transducer and reception energy converter.In Fig. 1,
Scatterer P is moved with speed v, namely with axial velocity va=vcos θ is moved to static sensor opposite direction.Assuming that transmitting
The frequency for the rf wave that energy converter issues is fC, the frequency (echo frequency) that the scatterer P moved reflects this pulse is
fR.Doppler frequency f is exactly the difference of the two frequencies, is expressed as follows:
Wherein, c indicates the velocity of sound.
When the movement velocity of scatterer P is far smaller than the velocity of sound, formula (1) can be converted are as follows:
In formula (2), due to radio frequency wave frequency rate fCIt is known quantity with velocity of sound c, therefore acquire Doppler frequency f to ask
The axial velocity that solution scatterer P is moved relative to sensor opposite direction.
The cross-correlation that frequency f, definition signal x (t) and signal x (t+ τ) how is estimated from two frame echo-signals is described below
Function R (τ) is expressed as follows:
Wherein, x*(t+ τ) is the conjugated signal of x (t+ τ)；Δ t indicates integrating range；D τ indicates the differential symbol of time variable
Number.
By the phase change of time domain in calculation formula (3) can approximation obtain frequency f variation:
Wherein, arg { R (τ) } indicates the phase difference of signal x (t) and x (t+ τ)；
To bring formula (4) into formula (2), can acquire:
Referring again to during elastogram, probe (ultrasonic probe or the probe that can issue other rf waves) squeezes life
Object tissue, this process may be considered the relative motion between probe and biological tissue, i.e., biological tissue is relative to probe
It is axially moved.By taking ultrasonic probe as an example, when ultrasonic probe squeezes biological tissue, ultrasonic probe is with TPRIt is issued for the period super
Sound rf wave receives the ultrasonic echo of tested biological tissue, to collect multiple frames of ultrasonic RF (radio frequency) in this process
Signal.The each frame ultrasonic radio frequency signal collected may include multi-strip scanning line, which can be perpendicular to tested life
The detection line on the surface of object tissue, every scan line may include multiple sampled points.Two frame ultrasonic radio frequency signals of arbitrary neighborhood
The moment is acquired, biological tissue's stress difference is tested, generates different deformation, the ultrasonic radio frequency signal collected can reflect the shape
Become.To be handled according to the analysis of ultrasonic radio frequency wave and echo to transmitting, the biography of rf wave in biological tissues can be obtained
Characteristic is broadcast, and then analyzes the soft or hard degree of biological tissue.The radio frequency wave frequency rate for the echo frequency and transmitting that receive of popping one's head in is not
With, the process is similar with Doppler effect principle, therefore, Doppler effect principle can be used, and calculates in conjunction with cross correlation algorithm
Biological tissue is extruded the VELOCITY DISTRIBUTION of direction (axial direction relative to probe).
When probe vertical squeezes biological tissue, the negative side of the scatterer P towards the sensor that are equivalent in Fig. 1 transports in the axial direction
Dynamic, θ=1 cos, the axial velocity acquired using formula (5) are the movement speed with biological tissue at this time.
Embodiment 1
As shown in Fig. 2, the elastograph imaging method of the present embodiment the following steps are included:
During step 101, acquisition squeeze biological tissue, corresponding to echo-signal of the biological tissue based on radio-wave reflection
Multiframe I/Q data.
Specifically, squeezing biological tissue using ultrasonic probe, ultrasonic probe is with T in the present embodimentPRFor period sending
Ultrasonic wave receives the ultrasonic echo of tested biological tissue in this process, to collect multiple frames of ultrasonic RF (radio frequency) letter
Number.RF signal is demodulated to obtain complex matrix, the real and imaginary parts of complex matrix respectively correspond I datum and Q data, respectively
Low-pass filtering is carried out to I datum and Q data, I/Q data can be obtained.In the present embodiment, for ease of calculation, according to I/Q data
The sequencing of acquisition carries out sequence number to it.
During elastogram, the I/Q data that samples is discrete signal, can be indicated with matrix, by formula (3),
(4) approximate description in discrete domain are as follows:
To,
Wherein, Num indicates the total quantity of I/Q data, and the similarity of adjacent two frames I/Q data is only calculated in the present embodiment, because
This Num=2；The offset row of h expression I/Q data matrix；0≤u≤H-h-1；0≤v≤W-1；0≤w≤Num-num-1；Num table
Show that the sequence of two frame I/Q datas is poor, num is the 0 same frame I/Q data of expression, and num is the adjacent two frames I/Q data of 1 expression.γ[h,
Num] expression sequence difference be num, deviate behavior h two frame I/Q datas similarity, arg { γ [h, num] } be its phase difference.
In the present embodiment, the sin cos functions transformation based on hardware end processing is can be used in demodulation, or is handled based on software end
Hilbert (Hilbert) transformation.It should be noted that the demodulation mode fast speed based on sin cos functions transformation, and
Demodulation mode precision based on Hilbert transformation is higher, and user can voluntarily select according to actual needs.Since the present embodiment uses
Elastogram is realized in conjunction with Doppler effect and cross correlation algorithm, and accuracy is higher, therefore realizes and demodulate even with hardware end, most
The precision being imaged afterwards there will not be too much influence, can also greatly improve the arithmetic speed of algorithm in this way.
Step 102, the first cross-correlation coefficient for calculating adjacent two frames I/Q data.
Below by taking adjacent the first I/Q data and the second I/Q data as an example, illustrate the calculating process of the first cross-correlation coefficient,
As shown in figure 3, step 102 specifically includes:
Step 102-1, the first IQ matrix and the 2nd IQ matrix are generated according to the first I/Q data and the second I/Q data respectively.
Step 102-1 namely the first I/Q data and the second I/Q data are indicated with matrix, obtain the first IQ matrix and the 2nd IQ
Matrix.Wherein, the first IQ matrix and the 2nd IQ matrix are plural two-dimensional matrixes.
Step 102-2, the conjugate matrices of the first IQ matrix and the 2nd IQ matrix are carried out respectively based on convolution algorithm window
Convolution algorithm.
In the present embodiment, use a behavior I, be classified as J, two-dimensional matrix that value is 1 is as convolution algorithm window.Step
102-2 is for noise reduction process.
It should be noted that due to the present embodiment primarily to obtaining biological tissue is extruded direction (relative to probe
Axial direction) movement speed, on the image namely longitudinal speed, therefore the relatively J that I can be arranged is more greatly for characterization, by examination
It tests and learns, when I is 20 times or more of J, conveniently.If also needing to obtain biological tissue relative to the shifting in probe horizontal direction
Dynamic speed, since the movement speed in the movement speed and axial direction in horizontal direction is mutually indepedent, convolution algorithm can
It is assigned on two threads and independently carries out, and then promote calculating speed.Certainly, carrying out convolution algorithm to different IQ matrixes can also lead to
The mode for crossing multi-threading parallel process can significantly improve calculation process speed while guaranteeing elastogram accuracy.
For the two-dimensional complex number IQ matrix shown in Fig. 4 a, Fig. 4 b, the detailed process of convolution algorithm is said
It is bright.Assuming that Fig. 4 a, Fig. 4 b show 12 × 8 complex matrixs, a stain in figure represents an element in matrix.It is black
Frame A indicates the convolution algorithm window having a size of I × J.Convolution algorithm is carried out to the first IQ matrix, also covers convolution algorithm window
On the first IQ matrix, and using n element as step-length, with arrow direction shown in scheming respectively from laterally, longitudinal two dimensions
Mobile convolution algorithm window, to traverse entire first IQ matrix.Convolution algorithm window often slides into some position, the IQ just covered
Submatrix and the operation window matrix multiple.Fig. 4 a is initial position of the convolution algorithm window in the first IQ matrix, it is covered at this time
The IQ submatrix of lid and the operation window matrix multiple.Fig. 4 b is that convolution algorithm windowsill vertically moves the position after 1 element step-length,
The IQ submatrix and the operation window matrix multiple that it is covered at this time.And so on, until traversing entire first IQ matrix, realize
To the convolution algorithm of the first IQ matrix.Data matrix advantage of this is that the influence that can reduce noise, after process of convolution
Size can become smaller.Wherein, n can self-setting according to actual needs, such as be set as 1 or 2 etc..To the conjugate matrices of the second matrix
The process for carrying out convolution algorithm is similar with the first matrix, and details are not described herein again.
Step 102-3, to the corresponding element of the conjugate matrices of the first IQ matrix and the 2nd IQ matrix by convolution algorithm
Product calculation is done, the first cross-correlation coefficient is obtained.
In step 102-3, so-called corresponding element does the element that product calculation also arranges the u row v of the first IQ matrix
The element multiplication arranged with the u row v of the conjugate matrices of the 2nd IQ matrix, the result of product calculation are the first I/Q data and the
As a result the cross-correlation coefficient of two I/Q datas and a complex matrix are expressed as γ [0,1], i.e., (sequence is poor for adjacent two IQ matrix
Num is the 1) related coefficient of error-free bit arithmetic (offset row h is 0).
Step 103, the VELOCITY DISTRIBUTION data that biological tissue is calculated based on the first cross-correlation coefficient.
In the present embodiment, step 103 is specifically included:
Step 103-1, the first phase difference per adjacent two frames I/Q data is calculated according to the first cross-correlation coefficient.
After related coefficient is calculated, using the equations phase difference in software tool, such as matlab (one is used
Kind of business mathematics software) the function angle for seeking plural number phase, acquire the phase difference of adjacent two frames I/Q data.
It should be noted that calculating phase-frequency characteristic, need to use arctan function, wherein provide, coordinate first and second as
Angle in limit is 0~π, and the angle of third and fourth quadrant is 0~-π.If the knot that an angle changes to 2 π, but actually obtain from 0
Fruit is 0~π, then by-π~0, phase can be jumped, and hopping amplitude is 2 π, this is just the winding of phase, be will affect final
The quality of imaging.In order to avoid the winding of phase, in step 103, solution winding processing is carried out to the phase difference being calculated, is made
Phase does not jump at π, to reflect true phase change, improves image quality.
Step 103-2, multiframe speed data is calculated according to the frequency of first phase difference and rf wave, according to multiframe speed
Degree is according to the VELOCITY DISTRIBUTION data for obtaining biological tissue.
After the phase difference for acquiring adjacent two frames I/Q data, formula (8) are carried it into, the frame speed of biological tissue can be acquired
Degree.Wherein, the frequency f for the rf wave that probe emits in formula (8)CNamely biological tissue receives the frequency of rf wave.Per adjacent
Two frame data require to calculate by formula (8), obtain multiframe speed data, obtain VELOCITY DISTRIBUTION number according to multiframe speed data
According to the VELOCITY DISTRIBUTION data are represented by a two-dimensional matrix.
Step 104, the Strain Distribution data that biological tissue is generated according to VELOCITY DISTRIBUTION data.
Specifically, differentiating in step 104 to VELOCITY DISTRIBUTION data, the Strain Distribution number of tissue can be obtained
According to.
The detailed process differentiated is described further below:
If the VELOCITY DISTRIBUTION data that step 103 obtains are indicated with one 1000 × 256 two-dimensional matrix.One is done now
Length is Δ m (it is assumed that Δ m is equal to 90) one-dimension array, wherein first element of one-dimension array is equal to the 1, the 90th position
Element is equal to -1, and other elements are 0, using the one-dimension array as operation window, traverses each of 1000 × 256 two-dimensional matrix
Convolution algorithm is arranged and carries out, this process is differential calculation, and the result differentiated is Strain Distribution data.This calculating
Method, which is simple and efficient, is not take up the excessive time.
Step 105 carries out imaging to Strain Distribution data, generates diagram of strains.
Specifically, step 105 includes:
Step 105-1, Strain Distribution data are converted into 8 gray level images to get preliminary Strain Distribution image is arrived.
Step 105-2, B mould image is obtained according to I/Q data.
B mould image can be obtained using existing method, details are not described herein again for specific implementation process.
In the implementation case, Strain Distribution image and B mould image are calculated due to using different filtering during processing
Method necessarily causes the picture size finally exported variant, therefore before merging Strain Distribution image and B mould image, needs
Adjust picture size.
Step 105-3, according to the dimension difference between B mould image and diagram of strains, B mould image or Strain Distribution are scaled
Figure, so that the size of B mould image is identical as the size of diagram of strains.
Wherein, zoomed image algorithm can be, but not limited to using bilinear interpolation or spline interpolation etc..Size it is identical namely
Guarantee that the position of biological tissue's organ and profile are completely coincident in image.
Preferably, size adjusting is carried out to diagram of strains in the present embodiment.By convolution algorithm, the ruler of diagram of strains
It is very little to become smaller, it needs to amplify diagram of strains processing: the size finally shown according to B mould image, image length and width
Size finds proportionate relationship, by Strain Distribution image proportionally relationship zoom into B mould picture size same size, make it
It can merge with B mould image.
Step 105-4, Pseudo Col ored Image is carried out to the diagram of strains after size adjusting, and is merged with B mould image.
Specifically, step 105-4 namely to diagram of strains adjusted carry out Pseudo Col ored Image, be taped against B mould figure
It is shown simultaneously as upper, by adjusting transparency, can be obtained and check diagram of strains.
In the present embodiment, in conjunction with Doppler effect and cross correlation algorithm, the speed during biological tissue is extruded is calculated
Distribution generates the diagram of strains of biological tissue, realizes elastogram, and precision is high, speed is fast, can improve healthcare givers to disease
Diagnostic reliability and efficiency.
In the present embodiment, during elastogram, the VELOCITY DISTRIBUTION data of generation and adjacent two frames I/Q data it is mutual
Relationship number can be used as the parameter of diagram of strains trust evaluation.The mean approximation pressure profile of VELOCITY DISTRIBUTION data, such as
The mean value fluctuation of the every frame speed distributed data of fruit is in (difference is smaller) under steady state, so that it may think the pressure that the process applies
Power is more steady, if mean value fluctuating change is more violent, so that it may think the pressure anomaly that the process applies, doctor can root
Force situation at that time is judged according to the fluctuation transformation of curve.Cross-correlation coefficient can be used to evaluate data acquisition quality, cross correlation
Number is high to illustrate that collected adjacent frame data do not have loss of data, otherwise explanation has loss of data, this numerical value can also be used as
Judge the parameter of current Strain Distribution image reliability.VELOCITY DISTRIBUTION data and cross-correlation coefficient are all in algorithm implementation procedure
Intermediate product, there is no additional time loss.
Embodiment 2
In the actual process, ripple attenuation, namely biology can occur in the communication process of biological tissue for ultrasonic radio frequency wave
Received ultrasonic radio frequency wave is organized gradually to decay with the increase fluctuation of depth, in order to further increase the accurate of elastogram
Degree needs during calculating speed distribution to frequency fCCarry out real-time update.
On the basis of embodiment 1, in the present embodiment, the frequency f based on real-time updateC, calculating speed distribution.
Relationship between the average frequency and phase of ultrasonic wave is as follows:
With the Δ θ in arg { γ [1,0] } replacement formula (10)i, updated frequency f can be obtainedC':
The speed data formula of adjacent two frames I/Q data is expressed as follows:
Wherein, 1/TsWhen indicating the data sampling frequency of probe, namely obtaining RF signal data, sampling radiofrequency signal scanning
The frequency of line；γ [0,1] indicates the cross-correlation coefficient of adjacent the first I/Q data and the second I/Q data of offset behavior 0, arg
{ γ [0,1] } is its phase difference；γ [1,0] indicates that the first I/Q data offset row h is the related coefficient and the 2nd IQ number of 1 matrix
The sum of the related coefficient of matrix for being 1 according to offset row h, arg { γ [1,0] } are its phase difference.
It can be seen that from formula (12), need to calculate phase difference arg { γ [1,0] } in the present embodiment.Or with adjacent the
For one I/Q data and the second I/Q data, illustrate the calculating process of [1,0] γ, as shown in figure 5, method further include:
Step 102-1 ', convolution fortune is carried out to the first IQ matrix and the 2nd IQ matrix respectively based on the first convolution algorithm window
It calculates, based on the second convolution algorithm window, the conjugate matrices of the conjugate matrices to the first IQ matrix and the 2nd IQ matrix are rolled up respectively
Product operation.
Wherein, the first convolution algorithm window is identical as the size of the second convolution algorithm window, moving step length is identical；The first volume
Product operation window the first IQ matrix initial position and the second convolution algorithm window the first IQ matrix conjugate matrices at the beginning of
Beginning position is different；First convolution algorithm window is in the initial position of the 2nd IQ matrix and the second convolution algorithm window in the 2nd IQ square
The initial position of the conjugate matrices of battle array is different.
Below in conjunction with Fig. 4 a, Fig. 4 b, the detailed process of the convolution algorithm of the first IQ matrix and its conjugate matrices is said
It is bright.Assuming that Fig. 4 a show the first IQ matrix and convolution algorithm during the first convolution algorithm window in the first IQ matrix
Initial position, by the IQ submatrix and the first convolution algorithm window matrix multiple of the covering of the first convolution algorithm window, operation result is made
The element arranged for the 1st row the 1st of new IQ matrix；Fig. 4 b shows the conjugate matrices and convolution algorithm process of the first IQ matrix
In initial position of the second convolution algorithm window in conjugate matrices, by the second convolution algorithm window covering IQ submatrix and this second
Convolution algorithm window matrix multiple, the element that operation result is arranged as the 1st row the 1st of new conjugate matrices.First convolution algorithm window with
Second convolution algorithm window using n element as step-length, with arrow direction shown in scheming respectively from laterally, longitudinal two dimensions move
Dynamic convolution algorithm window, to traverse entire first IQ matrix.The first convolution algorithm with n=1, for moving longitudinally, in Fig. 4 a
The upper side frame of window is moved to the 6th row from the first row of the first IQ matrix, terminates longitudinal traversal, not last to the first IQ matrix
Row data do operation；The second convolution algorithm window in Fig. 4 b is moved to the 7th row from the second row, terminates longitudinal traversal, not to first
The first row data of IQ matrix do operation；This is because the data of line direction are more in actual operation, even if losing a line number
Too big influence is not had according on imaging.Convolution algorithm is done in the same fashion for the 2nd IQ matrix and its conjugate matrices, this
Place repeats no more.
It can be seen that the convolution algorithm for the first IQ matrix and its conjugate matrices, the first convolution algorithm window from two figures
There is dislocation on initial position between the second convolution algorithm window, namely differs 1 element (h=1) on matrix line direction, this
It is the result for obtaining constant with its conjugate complex number matrix multiple in order to avoid complex matrix.It is multiplied, is multiplied by such dislocation
Matrix afterwards is still complex matrix.It should be noted that the dislocation of the first convolution algorithm window and the second convolution algorithm window
Distance can be one element distance of difference, can also differ multiple elements, but distance will lead to greatly very much algorithm failure, in turn
Final Strain Distribution picture quality is influenced, design parameter setting needs depending on specific input data.
Step 102-2 ', it is done to by the first IQ matrix of convolution algorithm and the conjugate matrices corresponding element of the first IQ matrix
Product calculation obtains the second cross-correlation coefficient；To the conjugate matrices of the 2nd IQ matrix and the 2nd IQ matrix by convolution algorithm
Corresponding element does product calculation, obtains third related coefficient.
Step 102-3 ', sum operation, the knot of sum operation are carried out to the second cross-correlation coefficient and third cross-correlation coefficient
Fruit is [1,0] γ, and poor according to γ [1,0] calculating second phase.
By step 102-3 ' be calculated as a result, bring into formula (11), realize to frequency fCReal-time update.
In step 103-2, then according to first phase difference and updated frequency fCCalculating speed data, referring to formula
(12)。
It should be noted that similar with first phase difference, it is also desirable to carry out solution winding processing to second phase difference.
In the present embodiment, when calculating the VELOCITY DISTRIBUTION during biological tissue is extruded, to frequency fCReal-time perfoming updates
It calculates, with the diagram of strains of this biological tissue generated, image quality is further increased.
Embodiment 3
Corresponding with the elastograph imaging method of embodiment 1, the present embodiment provides a kind of elastogram device, which is answered
Used in the imaging device of elastogram system, which further includes probe, is popped one's head in for emitting rf wave.Such as Fig. 6
Shown, the elastogram device of the present embodiment includes: the first data acquisition module 1, the first computing module 2, the second computing module
3, the second data acquisition module 4, image generation module 5, solution winding module 6.
For first data acquisition module 1 for obtaining during biological tissue is extruded, biological tissue is based on radio-wave reflection
Echo-signal corresponding to multiframe I/Q data.
First computing module 2 is used to calculate the first cross-correlation coefficient of adjacent two frames I/Q data.
In the present embodiment, the first computing module is specifically used for:
The first IQ matrix and the 2nd IQ matrix are generated according to adjacent two frames I/Q data；
Convolution fortune is carried out to the conjugate matrices of the first IQ matrix and the 2nd IQ matrix respectively based on convolution algorithm window
It calculates；
Product calculation will be done by the corresponding element of the conjugate matrices of the first IQ matrix and the 2nd IQ matrix of convolution algorithm,
Obtain the first cross-correlation coefficient.
Second computing module 3 is used to calculate a frame speed data according to the first cross-correlation coefficient.
In the present embodiment, the second computing module 3 is specifically used for:
The first phase difference of adjacent two frames I/Q data is calculated according to the first cross-correlation coefficient；
Speed data is calculated according to first phase difference and the frequency of rf wave.
In the present embodiment, in order to avoid the winding of phase, uncoiling is called after the phase difference that the second computing module 3 is calculated
It jumps phase not at π to carry out solution winding processing to phase difference around module 6, becomes to reflect true phase
Change, improves image quality.
Second data acquisition module 4 is used to obtain the VELOCITY DISTRIBUTION data of biological tissue, speed point according to multiframe I/Q data
Cloth data include multiframe speed data.
Image generation module 5 is used to generate diagram of strains according to VELOCITY DISTRIBUTION data.
In the present embodiment, image generation module 5 is specifically used for:
B mould image is generated according to multiframe I/Q data；
Strain Distribution data are generated according to VELOCITY DISTRIBUTION data, and imaging is carried out to Strain Distribution data, are generated just
Beginning diagram of strains；
According to the dimension difference between B mould image and initial strain distribution map, B mould image or initial strain distribution are scaled
Figure, so that the size of B mould image is identical as the size of initial strain distribution map；
Merge the B mould image and initial strain distribution map Jing Guo size adjusting, generates final diagram of strains.
In the present embodiment, in conjunction with Doppler effect and cross correlation algorithm, the speed during biological tissue is extruded is calculated
Distribution generates the diagram of strains of biological tissue, realizes elastogram, and precision is high, speed is fast, can improve healthcare givers to disease
Diagnostic reliability and efficiency.
Embodiment 4
Embodiment 4 is substantially the same manner as Example 3, the difference is that, the second computing module is based on real-time in the present embodiment
The frequency f of updateCCalculating speed distribution.
In the present embodiment, when the second computing module calculates speed data according to first phase difference and the frequency of rf wave,
It is specifically used for:
The second cross-correlation coefficient of the submatrix of the different offset rows of the first IQ matrix is calculated, and calculates the 2nd IQ matrix
The third cross-correlation coefficient of the submatrix of different offset rows；
Sum operation is carried out to the second cross-correlation coefficient and third cross-correlation coefficient, and is calculated according to the result of sum operation
Second phase is poor；
According to second phase difference renewal frequency；
Speed data is calculated according to first phase difference and updated frequency.
It should be noted that similar with first phase difference, it is also desirable to carry out solution winding processing to second phase difference.
In the present embodiment, when calculating the VELOCITY DISTRIBUTION during biological tissue is extruded, to frequency fCReal-time perfoming updates
It calculates, with the diagram of strains of this biological tissue generated, image quality is further increased.
Embodiment 5
Fig. 7 is a kind of structural schematic diagram of imaging device provided in an embodiment of the present invention, shows and is suitable for being used to realizing this
The block diagram of the exemplary imaging device 90 of invention embodiment.The imaging device 90 that Fig. 7 is shown is only an example, should not be right
The function and use scope of the embodiment of the present invention bring any restrictions.
As shown in fig. 7, imaging device 90 can be showed in the form of universal computing device, such as it can set for server
It is standby.The component of imaging device 90 can include but is not limited to: at least one above-mentioned processor 91, above-mentioned at least one processor
92, the bus 93 of different system components (including memory 92 and processor 91) is connected.
Bus 93 includes data/address bus, address bus and control bus.
Memory 92 may include volatile memory, such as random access memory (RAM) 921 and/or cache
Memory 922 can further include read-only memory (ROM) 923.
Memory 92 can also include (or the practical work of program means 925 with one group of (at least one) program module 924
Tool), such program module 924 includes but is not limited to: operating system, one or more application program, other program modules
And program data, it may include the realization of network environment in each of these examples or certain combination.
Processor 91 by the computer program that is stored in memory 92 of operation, thereby executing various function application and
Data processing, such as elastograph imaging method provided by the embodiment of the present invention 1 or 2.
Imaging device 90 can also be communicated with one or more external equipments 94 (such as keyboard, sensing equipment etc.).It is this
Communication can be carried out by input/output (I/O) interface 95.Also, the imaging device 90 that model generates can also pass through network
Adapter 96 and one or more network (such as local area network (LAN), wide area network (WAN) and/or public network, such as because of spy
Net) communication.As shown, the other modules for the imaging device 90 that network adapter 96 is generated by bus 93 and model communicate.
It should be understood that although not shown in the drawings, the imaging device 90 that can be generated with binding model uses other hardware and/or software mould
Block, including but not limited to: microcode, device driver, redundant processor, external disk drive array, RAID (disk array)
System, tape drive and data backup storage system etc..
It should be noted that although being referred to several units/modules or subelement/mould of imaging device in the above detailed description
Block, but it is this division be only exemplary it is not enforceable.In fact, embodiment according to the present invention, is retouched above
The feature and function for two or more units/modules stated can embody in a units/modules.Conversely, above description
A units/modules feature and function can with further division be embodied by multiple units/modules.
Embodiment 6
A kind of computer readable storage medium is present embodiments provided, computer program, described program quilt are stored thereon with
The step of elastograph imaging method provided by embodiment 1 or 2 is realized when processor executes.
Wherein, what readable storage medium storing program for executing can use more specifically can include but is not limited to: portable disc, hard disk, random
Access memory, read-only memory, erasable programmable read only memory, light storage device, magnetic memory device or above-mentioned times
The suitable combination of meaning.
In possible embodiment, the present invention is also implemented as a kind of form of program product comprising program generation
Code, when described program product is run on the terminal device, said program code is realized in fact for executing the terminal device
Apply the step in elastograph imaging method described in example 1 or 2.
Wherein it is possible to be write with any combination of one or more programming languages for executing program of the invention
Code, said program code can be executed fully on a user device, partly execute on a user device, is only as one
Vertical software package executes, part executes on a remote device or executes on a remote device completely on a user device for part.
Those skilled in the art are considering specification and after practicing disclosure disclosed herein, will readily occur to of the invention its
Its embodiment.The present invention is directed to cover any variations, uses, or adaptations of the invention, these modifications, purposes or
Person's adaptive change follows general principle of the invention and including the undocumented common knowledge in the art of the present invention
Or conventional techniques.The description and examples are only to be considered as illustrative, and true scope and spirit of the invention are by following
Claim is pointed out.
It should be understood that the present invention is not limited to the precise structure already described above and shown in the accompanying drawings, and
And various modifications and changes may be made without departing from the scope thereof.The scope of the present invention is limited only by the attached claims.
Claims (12)
1. a kind of elastograph imaging method, which is characterized in that the elastograph imaging method is applied to be set in the imaging of elastogram system
Standby, the elastogram system further includes probe, and for emitting rf wave, the elastograph imaging method includes: the probe
During acquisition biological tissue is extruded, corresponding to echo-signal of the biological tissue based on the radio-wave reflection
Multiframe I/Q data；
Calculate the first cross-correlation coefficient of adjacent two frames I/Q data；
A frame speed data are calculated according to first cross-correlation coefficient；
The VELOCITY DISTRIBUTION data of the biological tissue are obtained according to the multiframe I/Q data, the VELOCITY DISTRIBUTION data include multiframe
Speed data；
Diagram of strains is generated according to the VELOCITY DISTRIBUTION data.
2. elastograph imaging method as described in claim 1, which is characterized in that calculate the first cross-correlation of adjacent two frames I/Q data
Coefficient, comprising:
The first IQ matrix and the 2nd IQ matrix are generated according to the adjacent two frames I/Q data；
Convolution is carried out to the conjugate matrices of the first IQ matrix and the 2nd IQ matrix respectively based on convolution algorithm window
Operation；
Product calculation will be done by the corresponding element of the conjugate matrices of the first IQ matrix and the 2nd IQ matrix of convolution algorithm, obtained
First cross-correlation coefficient.
3. elastograph imaging method as claimed in claim 2, which is characterized in that calculate a frame according to first cross-correlation coefficient
Speed data, comprising:
The first phase difference of the adjacent two frames I/Q data is calculated according to first cross-correlation coefficient；
The speed data is calculated according to the first phase difference and the frequency of the rf wave.
4. elastograph imaging method as claimed in claim 3, which is characterized in that according to the first phase difference and the rf wave
Frequency calculate the speed data, comprising:
The second cross-correlation coefficient of the submatrix of the different offset rows of the first IQ matrix is calculated, and calculates the 2nd IQ square
The third cross-correlation coefficient of the submatrix of the different offset rows of battle array；
Sum operation is carried out to second cross-correlation coefficient and the third cross-correlation coefficient, and according to the result of sum operation
It is poor to calculate second phase；
The frequency is updated according to the second phase difference；
The speed data is calculated according to the first phase difference and updated frequency.
5. the elastograph imaging method as described in any one of claim 1-4, which is characterized in that according to the VELOCITY DISTRIBUTION number
According to generation diagram of strains, comprising:
B mould image is generated according to the multiframe I/Q data；
Strain Distribution data are generated according to the VELOCITY DISTRIBUTION data, and imaging is carried out to the Strain Distribution data, it is raw
At initial strain distribution map；
According to the dimension difference between the B mould image and the initial strain distribution map, the B mould image or described first is scaled
Beginning diagram of strains, so that the size of the B mould image is identical as the size of the initial strain distribution map；
Merge the B mould image and initial strain distribution map Jing Guo size adjusting, generates final diagram of strains.
6. a kind of elastogram device, which is characterized in that the elastogram device is applied to be set in the imaging of elastogram system
Standby, the elastogram system further includes probe, and for emitting rf wave, the elastogram device includes: the probe
First data acquisition module, for obtaining during biological tissue is extruded, the biological tissue is based on the rf wave
Multiframe I/Q data corresponding to the echo-signal of reflection；
First computing module, for calculating the first cross-correlation coefficient of adjacent two frames I/Q data；
Second computing module, for calculating a frame speed data according to first cross-correlation coefficient；
Second data acquisition module, for obtaining the VELOCITY DISTRIBUTION data of the biological tissue, institute according to the multiframe I/Q data
Stating VELOCITY DISTRIBUTION data includes multiframe speed data；
Image generation module, for generating diagram of strains according to the VELOCITY DISTRIBUTION data.
7. elastogram device as claimed in claim 6, which is characterized in that first computing module is specifically used for:
The first IQ matrix and the 2nd IQ matrix are generated according to the adjacent two frames I/Q data；
Convolution is carried out to the conjugate matrices of the first IQ matrix and the 2nd IQ matrix respectively based on convolution algorithm window
Operation；
Product calculation will be done by the corresponding element of the conjugate matrices of the first IQ matrix and the 2nd IQ matrix of convolution algorithm, obtained
First cross-correlation coefficient.
8. elastogram device as claimed in claim 7, which is characterized in that second computing module is specifically used for:
The first phase difference of adjacent two frames I/Q data is calculated according to first cross-correlation coefficient；
The speed data is calculated according to the first phase difference and the frequency of rf wave.
9. elastogram device as claimed in claim 8, which is characterized in that in second computing module according to described first
When the frequency of phase difference and rf wave calculates the speed data, it is specifically used for:
The second cross-correlation coefficient of the submatrix of the different offset rows of the first IQ matrix is calculated, and calculates the 2nd IQ square
The third cross-correlation coefficient of the submatrix of the different offset rows of battle array；
Sum operation is carried out to second cross-correlation coefficient and the third cross-correlation coefficient, and according to the result of sum operation
It is poor to calculate second phase；
The frequency is updated according to the second phase difference；
The speed data is calculated according to the first phase difference and updated frequency.
10. the elastogram device as described in any one of claim 6-9, which is characterized in that described image generation module
It is specifically used for:
B mould image is generated according to the multiframe I/Q data；
Strain Distribution data are generated according to the VELOCITY DISTRIBUTION data, and imaging is carried out to the Strain Distribution data, it is raw
At initial strain distribution map；
According to the dimension difference between the B mould image and the initial strain distribution map, the B mould image or described first is scaled
Beginning diagram of strains, so that the size of the B mould image is identical as the size of the initial strain distribution map；
Merge the B mould image and initial strain distribution map Jing Guo size adjusting, generates final diagram of strains.
11. a kind of imaging device including memory, processor and stores the calculating that can be run on a memory and on a processor
Machine program, which is characterized in that the processor is realized described in any one of claim 1 to 5 when executing the computer program
Elastograph imaging method.
12. a kind of computer readable storage medium, is stored thereon with computer program, which is characterized in that the computer program
The step of elastograph imaging method described in any one of claim 1 to 5 is realized when being executed by processor.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910773249.9A CN110477948B (en) | 2019-08-21 | 2019-08-21 | Elasticity imaging method and device, imaging equipment and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910773249.9A CN110477948B (en) | 2019-08-21 | 2019-08-21 | Elasticity imaging method and device, imaging equipment and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110477948A true CN110477948A (en) | 2019-11-22 |
| CN110477948B CN110477948B (en) | 2022-05-10 |
Family
ID=68552527
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910773249.9A Active CN110477948B (en) | 2019-08-21 | 2019-08-21 | Elasticity imaging method and device, imaging equipment and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110477948B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113397588A (en) * | 2020-03-16 | 2021-09-17 | 深圳市理邦精密仪器股份有限公司 | Elastography method and device and medical equipment |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030163044A1 (en) * | 2002-02-28 | 2003-08-28 | Andreas Heimdal | Method and apparatus for spectral strain rate visualization |
| CN102824194A (en) * | 2011-06-14 | 2012-12-19 | 深圳迈瑞生物医疗电子股份有限公司 | Displacement detecting method and device thereof in elasticity imaging |
| CN102860842A (en) * | 2012-09-26 | 2013-01-09 | 浙江大学 | Real-time quasi-static ultrasonic elastography method |
| CN103458800A (en) * | 2011-03-30 | 2013-12-18 | 日立阿洛卡医疗株式会社 | Methods and apparatus for ultrasound imaging |
| CN106456121A (en) * | 2014-05-02 | 2017-02-22 | B-K医疗公司 | Elastography visualization |
-
2019
- 2019-08-21 CN CN201910773249.9A patent/CN110477948B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030163044A1 (en) * | 2002-02-28 | 2003-08-28 | Andreas Heimdal | Method and apparatus for spectral strain rate visualization |
| CN103458800A (en) * | 2011-03-30 | 2013-12-18 | 日立阿洛卡医疗株式会社 | Methods and apparatus for ultrasound imaging |
| CN102824194A (en) * | 2011-06-14 | 2012-12-19 | 深圳迈瑞生物医疗电子股份有限公司 | Displacement detecting method and device thereof in elasticity imaging |
| CN102860842A (en) * | 2012-09-26 | 2013-01-09 | 浙江大学 | Real-time quasi-static ultrasonic elastography method |
| CN106456121A (en) * | 2014-05-02 | 2017-02-22 | B-K医疗公司 | Elastography visualization |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113397588A (en) * | 2020-03-16 | 2021-09-17 | 深圳市理邦精密仪器股份有限公司 | Elastography method and device and medical equipment |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110477948B (en) | 2022-05-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102641137B (en) | Measurement of viscoelastic property using amplitude-phase modulation | |
| CN101491450B (en) | Synchronized combining for contrast agent enhanced medical diagnostic ultrasound imaging | |
| CN104582582B (en) | Ultrasonic image-forming system memory architecture | |
| CN102640014B (en) | Image generating apparatus and image generating method | |
| Lu et al. | Unsupervised super-resolution framework for medical ultrasound images using dilated convolutional neural networks | |
| GB2502997A (en) | Suppression of reverberation and/or clutter in ultrasonic imaging | |
| CN106529561A (en) | Sparkle artifact detection in Ultrasound color flow | |
| CN109875606B (en) | Ultrasonic CT sound velocity imaging method based on prior reflection imaging | |
| CN104914440B (en) | A kind of compressed sensing ultrasonic beam synthetic method for merging dynamic aperture | |
| US20090062649A1 (en) | Scan conversion for ultrasonic imaging and apparatus using the same | |
| CN106680825A (en) | Acoustic array imaging system and method thereof | |
| US10772609B2 (en) | Ultrasonic imaging processing method and system based on RF data | |
| WO2019162485A1 (en) | Reconstruction system and method | |
| CN102631216A (en) | Ultrasound diagnostic apparatus and method of producing ultrasound image | |
| CN107260210A (en) | Diffraction source compensation in medical diagnostic ultrasonic imaging | |
| CN103945771A (en) | Ultrasonic diagnostic device and image processing device | |
| CN110301936A (en) | Frequency scanning for acoustic radiation force pulse | |
| CN110477948A (en) | Elastograph imaging method and device, imaging device, storage medium | |
| CN110536647A (en) | For handling system, method and computer-readable medium with composite ultraphonic image in the presence of motion | |
| US20130158403A1 (en) | Method for Obtaining a Three-Dimensional Velocity Measurement of a Tissue | |
| CN115666401A (en) | System and method for non-invasive pressure measurement | |
| Waag et al. | An eigenfunction method for reconstruction of large-scale and high-contrast objects | |
| CN111449681B (en) | Shear wave imaging method, device, equipment and readable storage medium | |
| CN105286917A (en) | Three-dimensional ultrasonic real-time imaging method and system based on area-array probes | |
| US20210088639A1 (en) | System and methods for flash suppression in ultrasound imaging |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |