CN107518918A - Ultrasonic elastograph imaging method and its system - Google Patents
Ultrasonic elastograph imaging method and its system Download PDFInfo
- Publication number
- CN107518918A CN107518918A CN201710954146.3A CN201710954146A CN107518918A CN 107518918 A CN107518918 A CN 107518918A CN 201710954146 A CN201710954146 A CN 201710954146A CN 107518918 A CN107518918 A CN 107518918A
- Authority
- CN
- China
- Prior art keywords
- mrow
- module
- data
- adaptability
- present frame
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 29
- 238000004364 calculation method Methods 0.000 claims abstract description 33
- 230000004927 fusion Effects 0.000 claims abstract description 23
- 238000012545 processing Methods 0.000 claims abstract description 20
- 238000000034 method Methods 0.000 claims abstract description 16
- 239000000523 sample Substances 0.000 claims abstract description 10
- 238000001914 filtration Methods 0.000 claims description 29
- 238000006073 displacement reaction Methods 0.000 claims description 8
- 230000014509 gene expression Effects 0.000 claims description 3
- 230000037433 frameshift Effects 0.000 claims 1
- 230000002123 temporal effect Effects 0.000 claims 1
- 238000012285 ultrasound imaging Methods 0.000 abstract description 3
- 238000005070 sampling Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 3
- 241001269238 Data Species 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 238000002559 palpation Methods 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000010349 pulsation Effects 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 210000002966 serum Anatomy 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
- 230000002792 vascular Effects 0.000 description 1
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/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Animal Behavior & Ethology (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
The present invention relates to a kind of ultrasonic elastograph imaging method and its system, including：Probe, signal acquisition module, signal processing module and multiframe Fusion Module, signal acquisition module are used for the rf data for gathering target area object；Signal processing module handles the shift value for obtaining present frame according to the rf data of present frame with the rf data of the default consecutive frame of caching by computing cross-correlation, and the shift value of present frame is counted as to the adaptability to changes of present frame by gradient calculation turn；Multiframe Fusion Module reads some adaptability to changes frames adjacent with current adaptability to changes frame, by multiframe fusion treatment, obtains relative resilient information of the target area to surrounding background area.The present invention is had the characteristics of direction of motion uniformity using being organized in elastic ultrasound imaging procedure in continuous multiple frames image, is improved reliability and the degree of accuracy of adaptability to changes estimation using the ultrasonic elastograph imaging mode with multiframe fusion treatment.
Description
Technical field
It is especially a kind of super with multiframe fusion treatment the present invention relates to a kind of ultrasonic elastograph imaging method and its system
Acoustic elasticity imaging method and its system, belong to ultrasonic imaging technique field.
Background technology
The elasticity or hardness of biological tissue and the biological characteristics of focus are closely related, the hardness of biological tissue or elasticity
Change often implies the generation of lesion, and this is largely relevant with combining form with the molecule composition of biological tissue, for
The diagnosis of disease has important reference value.Clinically use the method for the consistency and elasticity of more detection biological tissue just
It is palpation, is directly used in by clinician and detected part is extruded, the hardening of tissue is judged by the sense of touch of finger
Degree is detected and diagnosed, and is easily invaded patient privacy and is easily caused mistaken diagnosis.
At present, the Ultrasonic Elasticity Imaging of the consistency and elasticity based on detection biological tissue is played the part of in terms of disease detection
Important role.Its principle is by mild compression destination organization or breathing by human body itself, the vascular pulsation of popping one's head in
Etc. process to organizing the formation of certain pressure, obtain and compress two forward and backward frame ultrasound echo signals.Believed according to two frame ultrasonic echos
Number calculate shift value of the test serum after being squeezed using Displacement Estimation algorithm and be converted into strain information and transmit to display
Interface is simultaneously intuitively shown with image format, to aid in diagnosis or treatment.However, prior art mainly uses base
Elastogram is not made full use of in the mode of single frames processing, or the mode based on simple multiframe linear filtering, these schemes
During the kinetic characteristic that shows is organized in multiple image.
The content of the invention
The purpose of the present invention is to overcome the deficiencies in the prior art, there is provided a kind of ultrasonic elastograph imaging method and its is
System is continuous more using being organized in elastic ultrasound imaging procedure using the ultrasonic elastograph imaging mode with multiframe fusion treatment
There is the characteristics of direction of motion uniformity in two field picture, improve reliability and the degree of accuracy of adaptability to changes estimation.
According to technical scheme provided by the invention, the ultrasonic elastograph imaging system, including：
Probe；
Signal acquisition module, for gathering the rf data of target area object reflection by the probe；
Signal processing module, for calculating the adaptability to changes of present frame according to the rf data of collection；And
Multiframe Fusion Module, the multiframe Fusion Module are connected with signal processing module, for reading and current adaptability to changes
The adjacent some adaptability to changes frames of frame, merged by multiframe and obtain elastic information of the target area with respect to surrounding background area.
Further, the signal processing module includes：
Rf data filtration module；
Rf data cache module, for caching the rf data of the rf data filtration module filtering process；
Displacement Estimation module, for the present frame rf data according to the rf data filtration module filtering process with delaying
The rf data that the consecutive frame in rf data cache module be present handles the shift value for obtaining present frame by computing cross-correlation；
Strain calculation module, the shift value for the present frame according to acquisition turn to be counted as answering for present frame by gradient calculation
Variable force；And
Elastic information cache module, it is connected with the strain calculation module.
Further, the signal processing module also includes beam-controller, and the beam-controller is adopted with the signal
Collect module connection, for control signal acquisition module to target area transmitting focusing ultrasonic wave and control signal acquisition module pair
Target area receives echo rf data.
Further, the signal acquisition module includes transmitting beamformer and receives Beam-former, launching beam
Shaper launches target area the focusing ultrasonic wave of certain carrier wave, receives Beam-former and obtains returning from target area reflection
Ripple rf data.
Further, in addition to elastic information display module, the elastic information display module connect with multiframe Fusion Module
Connect, for showing that the multiframe Fusion Module is transmitted to the elastic information of the elastic information display module.
The ultrasonic elastograph imaging method, comprises the following steps：
S1, the rf data of collection target area object reflection；
S2, the adaptability to changes according to the rf data of collection calculating present frame；
S3, some adaptability to changes frames adjacent with current adaptability to changes frame are read, it is relative to merge acquisition target area by multiframe
The elastic information of surrounding background area.
Further, the step S2 is concretely comprised the following steps：
S21, the rf data of caching rf data filtration module processing；
S22, according to the present frame rf data of rf data filtration module filtering process and it is buffered in rf data caching
The rf data of consecutive frame in unit module handles the shift value for obtaining present frame by computing cross-correlation；
S23, turned by gradient calculation according to the shift value of the present frame of acquisition to be counted as the adaptability to changes of present frame.
Further, the computing cross-correlation in the step S22 is：
Qab(t)=max-m/2≤τ≤m/2Rab(t, τ),
Wherein, n is that cross-correlation calculation accumulates siding-to-siding block length, and span is positive integer；M searches for for cross-correlation calculation result
Siding-to-siding block length, span are positive integer；τ is the searching position value in the cross-correlation calculation result region of search, and τ spans are
Positive integer；Rab(t, τ) is certain locus t, searching position τ cross-correlation calculation result, Qab(t) it is the displacement of current location
Value, t is positive integer.
Further, the gradient calculation in the step S23 is：
Strain (x, y, z)=(Qab(t-s/2)-Qab(t+s/2))/s,
Wherein, Strain (x, y, z) is the adaptability to changes of current location；S is gradient calculation gap length, and value is just whole
Number；X is image transverse axis coordinate, span 1 to X；Y is image ordinate of orthogonal axes, span 1 to Y；Z sits for image timeline
Mark, span 1 to Z；Strain (x, y, z) is signed number, and positive number represents current pixel and compares previous frame toward close to probe
Apparent motion, negative represent current pixel and compare previous frame toward away from detecting head surface motion.
Further, the step S3 is specially：
Obtain and include the continuous adjacent N frames of current pixel (x0, y0) adaptability to changes frame information, given threshold threshold_N,
When pixel strain force information is more than threshold value more than the 0 or statistics number statistics_N less than 0 in continuous adjacent N frames
During threshold_N, mark current pixel is valid pixel, calculates the effective strain force information of the pixel：
Wherein, z (1) expressions and the most short valid frame of present frame time gap, Q (*) is the anti-letter related to time gap
Number.
Ultrasonic elastograph imaging method of the present invention and its system, using the ultrasonic elastograph imaging side of multiframe fusion treatment
Formula, there is the characteristics of direction of motion uniformity using being organized in elastic ultrasound imaging procedure in continuous multiple frames image, improving should
The reliability of variable force estimation and the degree of accuracy.
Brief description of the drawings
Fig. 1 is the schematic diagram of ultrasonic elastograph imaging system of the present invention.
Fig. 2 is ultrasonic elastograph imaging system direction of displacement schematic diagram of the present invention.
Description of reference numerals：1- probes, 2- signal acquisition modules, 3- signal processing modules, 21- transmitting beamformer,
22- beam-controllers, 23- receive Beam-former, 31- rf datas filtration module, 32- rf datas cache module, 33- positions
Estimation module, 34- strain calculations module, 35- elastic informations frame buffer module, 4- multiframes Fusion Module, 5- elastic informations is moved to show
Show module.
Embodiment
With reference to specific accompanying drawing, the invention will be further described.
As shown in figure 1, ultrasonic elastograph imaging system of the present invention includes probe 1, signal acquisition module 2, signal transacting
Module 3, multiframe Fusion Module 4 and elastic information display module 5.Signal acquisition module 2 is used for by the collection of probe 1 target
The rf data of regional aim thing reflection.Signal processing module 3 is used to calculate present frame according to the rf data of collection
Adaptability to changes.Multiframe Fusion Module 4, the multiframe Fusion Module 4 are connected with signal processing module 3, for reading and current strain
The adjacent some adaptability to changes frames of power frame, merged by multiframe and obtain elastic information of the target area with respect to surrounding background area, bullet
Property information is mainly relative hardness of the target area to surrounding background area, and the strain ratio of harder target area is larger, softer
The strain ratio of target area is smaller.
The signal processing module 3 of the present invention is estimated including rf data filtration module 31, rf data cache module 32, displacement
Count module 33, strain calculation module 34 and elastic information cache module 35.Rf data cache module 32 is used to cache described penetrate
Frequency according to the filtering process of filtration module 31 the rf data.For example, penetrating when rf data cache module 32 is cached with T0
Frequency data message, rf data information when being now, in rf data filtration module 31 T1, then during T1 rf data letter
Breath is cached in rf data cache module 32, instead of rf data information during T0.Displacement Estimation module 33 is used for according to institute
State the present frame T of the filtering process of rf data filtration module 31n+1Rf data faWith being buffered in rf data cache module 32
In consecutive frame TnRf data fb, the shift value for obtaining present frame is handled by computing cross-correlation.As shown in Fig. 2 shift value
For signed number, positive number represents current pixel and compares the past close detecting head surface motion of previous frame, and negative represents current pixel and compared
Previous frame moves toward away from detecting head surface, and pixel is the concept of spatial sampling point, and adaptability to changes is the characteristic of spatial sampling point, frame bag
Many space pixels are contained.Strain calculation module 34 is used to turn to be counted as by gradient calculation according to the shift value of the present frame of acquisition
The adaptability to changes of present frame.Elastic information cache module 35 is connected with the strain calculation module 34.
The signal processing module 3 of the present invention also includes beam-controller 22, and the beam-controller 22 is adopted with the signal
Collection module 2 connects, for control signal acquisition module 2 to target area transmitting focusing ultrasonic wave and control signal acquisition module
2 pairs of target areas receive echo rf data.Signal acquisition module 2 includes transmitting beamformer 21 and receives Beam-former
23, transmitting beamformer 21 launches target area the focusing ultrasonic wave of certain carrier wave, receive Beam-former 23 obtain from
The echo rf data of target area reflection.Present invention additionally comprises elastic information display module 5, the elastic information display module
5 are connected with multiframe Fusion Module 4, for showing that the multiframe Fusion Module 4 is transmitted to the bullet of the elastic information display module 5
Property information.
The ultrasonic elastograph imaging method, comprises the following steps：
S1, the rf data of collection target area object reflection；
S2, the adaptability to changes according to the rf data of collection calculating present frame；
S3, some adaptability to changes frames adjacent with current adaptability to changes frame are read, it is relative to merge acquisition target area by multiframe
The elastic information of surrounding background area.
Further, the step S2 is concretely comprised the following steps：
S21, the rf data of caching rf data filtration module processing；
S22, according to the present frame rf data of rf data filtration module filtering process and it is buffered in rf data caching
The rf data of consecutive frame in unit module handles the shift value for obtaining present frame by computing cross-correlation；According to rf data
The present frame T of the filtering process of filtration module 31n+1Rf data faWith the consecutive frame being buffered in rf data cache module 32
TnRf data fb, the shift value for obtaining present frame is handled by computing cross-correlation.The shift value of present frame is signed number,
Positive number, which represents current pixel, compares previous frame and is moved toward close to detecting head surface, and negative represents current pixel and compares previous frame toward remote
Detecting head surface moves, and pixel is the concept of spatial sampling point, and adaptability to changes is the characteristic of spatial sampling point, and frame contains many spaces
Pixel；
S23, turned by gradient calculation according to the shift value of the present frame of acquisition to be counted as the adaptability to changes of present frame.
Further, the computing cross-correlation in the step S22 is：
Qab(t)=max-m/2≤τ≤m/2Rab(t, τ),
Wherein, n is that cross-correlation calculation accumulates siding-to-siding block length, and span is positive integer；M searches for for cross-correlation calculation result
Siding-to-siding block length, span are positive integer；τ is the searching position value in the cross-correlation calculation result region of search, and τ spans are
Positive integer；Rab(t, τ) is certain locus t, searching position τ cross-correlation calculation result, Qab(t) it is the displacement of current location
Value, t is positive integer.
Further, the gradient calculation in the step S23 is：
Strain (x, y, z)=(Qab(t-s/2)-Qab(t+s/2))/s,
Wherein, strain (x, y, z) is the adaptability to changes of current location；S is gradient calculation gap length, and value is just whole
Number；X is image transverse axis coordinate, span 1 to X；Y is image ordinate of orthogonal axes, span 1 to Y；Z sits for image timeline
Mark, span 1 to Z；Strain (x, y, z) is signed number, and positive number represents current pixel and compares previous frame toward close to probe
Apparent motion, negative represent current pixel and compare previous frame toward away from detecting head surface motion.
The step S3 is specially：
Obtain and include the continuous adjacent N frames of current pixel (x0, y0) adaptability to changes frame information, given threshold threshold_N,
When pixel strain force information is more than threshold value more than the 0 or statistics number statistics_N less than 0 in continuous adjacent N frames
During threshold_N, mark current pixel is that (purpose of this step is to filter out valid pixel to valid pixel, only as past N
The pixel that consistent directivity characteristic is shown in frame is only valid pixel), calculate the effective strain force information of the pixel：
Wherein, z (1) expressions and the most short valid frame of present frame time gap, Q (*) is the anti-letter related to time gap
Number, such as Q (dixtant)=e-distant/theta, e (*) is exponential function, and theta is at control multiframe adaptability to changes information smoothing
The parameter of reason, span are positive integer.This multiframe fusion treatment make use of quasistatic ultrasonic elastograph imaging principle, by more
The smoothing processing of sample statistics method and time-axis direction, suppress strain force information wrong in some frames, calculated so as to improve
As a result robustness.The strain force information calculated by multiframe fusion treatment is signed number, by asking for absolute value, is changed into
Unsigned number, further handled to obtain elastic information by log-compressed, scanning conversion, GTG mapping etc..
It is the description of specific embodiment of the present invention above, it should be understood that do not depart from this hair into those skilled in the art
Under bright true spirit and scope, substitution modification is carried out by other various simple changes and equivalent, reaches of the present invention
Purpose, such modification is all by scope of the following claims.
Claims (10)
1. a kind of ultrasonic elastograph imaging system, including：
Pop one's head in (1)；
Signal acquisition module (2), for gathering the rf data of target area object reflection by the probe (1)；
Signal processing module (3), for calculating the adaptability to changes of present frame according to the rf data of collection；And
Multiframe Fusion Module (4), the multiframe Fusion Module (4) are connected with signal processing module (3), for reading with currently should
The adjacent some adaptability to changes frames of variable force frame, merged by multiframe and obtain elastic information of the target area with respect to surrounding background area.
2. ultrasonic elastograph imaging system as claimed in claim 1, it is characterized in that：The signal processing module (3) includes：
Rf data filtration module (31)；
Rf data cache module (32), for caching the radio frequency number of rf data filtration module (31) filtering process
According to；
Displacement Estimation module (33), for the present frame rf data according to rf data filtration module (31) filtering process
The rf data of consecutive frame with being buffered in rf data cache module (32) is handled by computing cross-correlation obtains present frame
Shift value；
Strain calculation module (34), the shift value for the present frame according to acquisition turn to be counted as answering for present frame by gradient calculation
Variable force；And
Elastic information cache module (35), it is connected with the strain calculation module (34).
3. ultrasonic elastograph imaging system as claimed in claim 1, it is characterized in that：The signal processing module (3) also includes ripple
Beam controller (22), the beam-controller (22) is connected with the signal acquisition module (2), for control signal acquisition module
(2) echo rf data is received to target area to target area transmitting focusing ultrasonic wave and control signal acquisition module (2).
4. ultrasonic elastograph imaging system as claimed in claim 1, it is characterized in that：The signal acquisition module (2) includes transmitting
Beam-former (21) and reception Beam-former (23), transmitting beamformer (21) launch certain carrier wave to target area
Ultrasonic wave is focused on, Beam-former (23) is received and obtains the echo rf data reflected from target area.
5. ultrasonic elastograph imaging system as claimed in claim 1, it is characterized in that：Also include elastic information display module (5), institute
Elastic information display module (5) is stated to be connected with multiframe Fusion Module (4), for show the multiframe Fusion Module (4) transmit to
The elastic information of the elastic information display module (5).
6. a kind of ultrasonic elastograph imaging method, it is characterized in that, comprise the following steps：
S1, the rf data of collection target area object reflection；
S2, the adaptability to changes according to the rf data of collection calculating present frame；
S3, some adaptability to changes frames adjacent with current adaptability to changes frame are read, merged by multiframe and obtain target area with respect to surrounding
The elastic information of background area.
7. ultrasonic elastograph imaging method as claimed in claim 6, it is characterized in that：The step S2's concretely comprises the following steps：
S21, the rf data of caching rf data filtration module (31) processing；
S22, according to the present frame rf data of rf data filtration module (31) filtering process and it is buffered in rf data caching
The rf data of consecutive frame in unit module (32) handles the shift value for obtaining present frame by computing cross-correlation；
S23, turned by gradient calculation according to the shift value of the present frame of acquisition to be counted as the adaptability to changes of present frame.
8. ultrasonic elastograph imaging method as claimed in claim 7, it is characterized in that：Computing cross-correlation in the step S22 is：
<mrow>
<msub>
<mi>R</mi>
<mi>ab</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mo>&Integral;</mo>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mi>n</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mi>n</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msubsup>
<msub>
<mi>f</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>f</mi>
<mi>b</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>+</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mi>dx</mi>
<mo>,</mo>
</mrow>
Qab(t)=max-m/2≤τ≤m/2Rab(t, τ),
Wherein, n is that cross-correlation calculation accumulates siding-to-siding block length, and span is positive integer；M is the cross-correlation calculation result region of search
Length, span are positive integer；τ is the searching position value in the cross-correlation calculation result region of search, and τ spans are just whole
Number；Rab(t, τ) is certain locus t, searching position τ cross-correlation calculation result, Qab(t) it is the shift value of current location, t
For positive integer.
9. ultrasonic elastograph imaging method as claimed in claim 8, it is characterized in that：Gradient calculation in the step S23 is：
Strain (x, y, z)=(Qab(t-s/2)-Qab(t+s/2))/s,
Wherein, Strain (x, y, z) is the adaptability to changes of current location；S is gradient calculation gap length, and value is positive integer；X is
Image transverse axis coordinate, span 1 to X；Y is image ordinate of orthogonal axes, span 1 to Y；Z is image temporal axial coordinate, value
Scope 1 is to Z；Strain (x, y, z) is signed number.
10. ultrasonic elastograph imaging method as claimed in claim 6, it is characterized in that：The step S3 is specially：
Obtain the continuous adjacent N frames for including current pixel (x0, y0) adaptability to changes frame information, given threshold threshold_N, when even
Continue pixel strain force information in adjacent N frames and be more than threshold value more than the 0 or statistics number statistics_N less than 0
During threshold_N, mark current pixel is valid pixel, calculates the effective strain force information of the pixel：
<mrow>
<mi>S</mi>
<mi>t</mi>
<mi>r</mi>
<mi>a</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mn>0</mn>
<mo>,</mo>
<mi>y</mi>
<mn>0</mn>
<mo>,</mo>
<mi>z</mi>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>s</mi>
<mi>t</mi>
<mi>a</mi>
<mi>t</mi>
<mi>i</mi>
<mi>s</mi>
<mi>t</mi>
<mi>i</mi>
<mi>c</mi>
<mi>s</mi>
<mo>_</mo>
<mi>N</mi>
</mrow>
</munderover>
<mrow>
<mo>(</mo>
<mi>Q</mi>
<mo>(</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mi>z</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>z</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>*</mo>
<mi>s</mi>
<mi>t</mi>
<mi>r</mi>
<mi>a</mi>
<mi>i</mi>
<mi>n</mi>
<mo>(</mo>
<mrow>
<mi>x</mi>
<mn>0</mn>
<mo>,</mo>
<mi>y</mi>
<mn>0</mn>
<mo>,</mo>
<mi>z</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>s</mi>
<mi>t</mi>
<mi>a</mi>
<mi>t</mi>
<mi>i</mi>
<mi>s</mi>
<mi>t</mi>
<mi>i</mi>
<mi>c</mi>
<mi>s</mi>
<mo>_</mo>
<mi>N</mi>
</mrow>
</munderover>
<mrow>
<mo>(</mo>
<mi>Q</mi>
<mo>(</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mi>z</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>z</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
Wherein, z (1) expressions and the most short valid frame of present frame time gap, Q (*) is the inverse function related to time gap.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710954146.3A CN107518918B (en) | 2017-10-13 | 2017-10-13 | Ultrasonic elastography method and system thereof |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710954146.3A CN107518918B (en) | 2017-10-13 | 2017-10-13 | Ultrasonic elastography method and system thereof |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107518918A true CN107518918A (en) | 2017-12-29 |
| CN107518918B CN107518918B (en) | 2020-06-26 |
Family
ID=60684039
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710954146.3A Active CN107518918B (en) | 2017-10-13 | 2017-10-13 | Ultrasonic elastography method and system thereof |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107518918B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112826533A (en) * | 2021-01-11 | 2021-05-25 | 深圳华声医疗技术股份有限公司 | Ultrasonic imaging space compounding method and device, ultrasonic diagnostic apparatus and storage medium |
| CN112887022A (en) * | 2019-11-29 | 2021-06-01 | 无锡祥生医疗科技股份有限公司 | Ultrasonic image data transmission method and terminal ultrasonic equipment |
| CN113397588A (en) * | 2020-03-16 | 2021-09-17 | 深圳市理邦精密仪器股份有限公司 | Elastography method and device and medical equipment |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090292205A1 (en) * | 2006-07-18 | 2009-11-26 | Takashi Osaka | Ultrasonic diagnostic apparatus |
| JP2010119630A (en) * | 2008-11-20 | 2010-06-03 | Ge Medical Systems Global Technology Co Llc | Ultrasonograph |
| CN102423264A (en) * | 2011-09-01 | 2012-04-25 | 中国科学院深圳先进技术研究院 | Image-based biological tissue elasticity measuring method and device |
| CN104739451A (en) * | 2013-12-27 | 2015-07-01 | 深圳迈瑞生物医疗电子股份有限公司 | Elastic image imaging method and device and ultrasonic imaging device |
| CN104783835A (en) * | 2014-01-16 | 2015-07-22 | 富士胶片株式会社 | Ultrasound diagnostic apparatus and ultrasound image producing method |
| CN106373166A (en) * | 2016-10-13 | 2017-02-01 | 重庆理工大学 | Multi-shape plane shear wave composite imaging method based on Mach cone effects |
| CN106923864A (en) * | 2015-11-03 | 2017-07-07 | 东芝医疗系统株式会社 | Diagnostic ultrasound equipment, image processing apparatus and image processing program |
-
2017
- 2017-10-13 CN CN201710954146.3A patent/CN107518918B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090292205A1 (en) * | 2006-07-18 | 2009-11-26 | Takashi Osaka | Ultrasonic diagnostic apparatus |
| JP2010119630A (en) * | 2008-11-20 | 2010-06-03 | Ge Medical Systems Global Technology Co Llc | Ultrasonograph |
| CN102423264A (en) * | 2011-09-01 | 2012-04-25 | 中国科学院深圳先进技术研究院 | Image-based biological tissue elasticity measuring method and device |
| CN104739451A (en) * | 2013-12-27 | 2015-07-01 | 深圳迈瑞生物医疗电子股份有限公司 | Elastic image imaging method and device and ultrasonic imaging device |
| CN104783835A (en) * | 2014-01-16 | 2015-07-22 | 富士胶片株式会社 | Ultrasound diagnostic apparatus and ultrasound image producing method |
| CN106923864A (en) * | 2015-11-03 | 2017-07-07 | 东芝医疗系统株式会社 | Diagnostic ultrasound equipment, image processing apparatus and image processing program |
| CN106373166A (en) * | 2016-10-13 | 2017-02-01 | 重庆理工大学 | Multi-shape plane shear wave composite imaging method based on Mach cone effects |
Non-Patent Citations (1)
| Title |
|---|
| JINGFENG JIANG ; T.J. HALL ; A.M. SOMMER: "A novel performance descriptor for ultrasonic strain imaging: a preliminary study", 《IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112887022A (en) * | 2019-11-29 | 2021-06-01 | 无锡祥生医疗科技股份有限公司 | Ultrasonic image data transmission method and terminal ultrasonic equipment |
| CN112887022B (en) * | 2019-11-29 | 2022-11-25 | 无锡祥生医疗科技股份有限公司 | Ultrasonic image data transmission method and terminal ultrasonic equipment |
| CN113397588A (en) * | 2020-03-16 | 2021-09-17 | 深圳市理邦精密仪器股份有限公司 | Elastography method and device and medical equipment |
| CN112826533A (en) * | 2021-01-11 | 2021-05-25 | 深圳华声医疗技术股份有限公司 | Ultrasonic imaging space compounding method and device, ultrasonic diagnostic apparatus and storage medium |
| CN112826533B (en) * | 2021-01-11 | 2021-08-17 | 深圳华声医疗技术股份有限公司 | Ultrasonic imaging space compounding method and device, ultrasonic diagnostic apparatus and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107518918B (en) | 2020-06-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4898809B2 (en) | Ultrasonic diagnostic equipment | |
| JP5385533B2 (en) | Elastic imaging method and apparatus | |
| US8133179B2 (en) | Method of displaying an elastic image | |
| JP5721311B2 (en) | Contrast-enhanced ultrasound diagnostic imaging method, computer-readable storage medium, and contrast-enhanced ultrasound diagnostic imaging apparatus | |
| JP4795675B2 (en) | Medical ultrasound system | |
| US20080214934A1 (en) | Inter-frame processing for contrast agent enhanced medical diagnostic ultrasound imaging | |
| US20070093716A1 (en) | Method and apparatus for elasticity imaging | |
| US20050283076A1 (en) | Non-invasive diagnosis of breast cancer using real-time ultrasound strain imaging | |
| CN107518918A (en) | Ultrasonic elastograph imaging method and its system | |
| JP2007282932A (en) | Method of generating elastic image and ultrasonograph | |
| CN103169500A (en) | Ultrasonic diagnostic apparatus, medical image diagnostic apparatus, and medical image processing method | |
| JP2003070786A (en) | Ultrasonic imaging unit | |
| KR20120013585A (en) | Method for Estimating Acoustic Velocity of Ultrasonic Image and Ultrasonic Diagnosis Apparatus using the same | |
| WO2015099836A1 (en) | Ultrasound transducer and ultrasound imaging system with a variable thickness dematching layer | |
| US20070167772A1 (en) | Apparatus and method for optimized search for displacement estimation in elasticity imaging | |
| EP2138102B1 (en) | Formation of an enhanced elastic image in an ultrasound system | |
| KR101629541B1 (en) | Ultrasonic diagnostic apparatus and control program thereof | |
| JP5415669B2 (en) | Ultrasonic diagnostic equipment | |
| KR101014558B1 (en) | Elastic image formation method with the normallization of strain, and supersonic wave system for the same | |
| US9125618B2 (en) | Providing an elastic image in an ultrasound system | |
| KR20110088830A (en) | Ultrasonic diagnosis apparatus and method for processing ultrasonic image | |
| EP0934724A1 (en) | Ultrasonic diagnostic apparatus | |
| US20230061869A1 (en) | System and methods for beamforming sound speed selection | |
| Matsumura et al. | Diagnostic results for breast disease by real-time elasticity imaging system | |
| JP4005856B2 (en) | Ultrasonic diagnostic equipment |
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 |