US11096671B2 - Sparkle artifact detection in ultrasound color flow - Google Patents
Sparkle artifact detection in ultrasound color flow Download PDFInfo
- Publication number
- US11096671B2 US11096671B2 US14/850,707 US201514850707A US11096671B2 US 11096671 B2 US11096671 B2 US 11096671B2 US 201514850707 A US201514850707 A US 201514850707A US 11096671 B2 US11096671 B2 US 11096671B2
- Authority
- US
- United States
- Prior art keywords
- locations
- color flow
- sparkle
- similarity
- flow
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active, expires
Links
- 238000001514 detection method Methods 0.000 title claims description 10
- 238000002604 ultrasonography Methods 0.000 title description 11
- 230000033001 locomotion Effects 0.000 claims abstract description 79
- 208000000913 Kidney Calculi Diseases 0.000 claims abstract description 7
- 206010029148 Nephrolithiasis Diseases 0.000 claims abstract description 7
- 238000001914 filtration Methods 0.000 claims description 25
- 238000005070 sampling Methods 0.000 claims description 24
- 238000000034 method Methods 0.000 claims description 22
- 239000012530 fluid Substances 0.000 claims description 17
- 238000013507 mapping Methods 0.000 claims description 6
- 230000002708 enhancing effect Effects 0.000 claims description 3
- 208000001130 gallstones Diseases 0.000 claims description 3
- 210000003734 kidney Anatomy 0.000 claims 1
- 238000003384 imaging method Methods 0.000 abstract description 19
- 230000015654 memory Effects 0.000 description 23
- 230000006870 function Effects 0.000 description 15
- 238000012545 processing Methods 0.000 description 12
- 230000005540 biological transmission Effects 0.000 description 10
- 230000004044 response Effects 0.000 description 10
- 230000000875 corresponding effect Effects 0.000 description 5
- 230000009467 reduction Effects 0.000 description 5
- 239000000872 buffer Substances 0.000 description 4
- 230000001934 delay Effects 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 4
- 239000004575 stone Substances 0.000 description 4
- 238000012285 ultrasound imaging Methods 0.000 description 4
- 230000017531 blood circulation Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000002592 echocardiography Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000003321 amplification Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000000295 complement effect Effects 0.000 description 2
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 201000001883 cholelithiasis Diseases 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 239000012528 membrane Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000035479 physiological effects, processes and functions Effects 0.000 description 1
- 238000011524 similarity measure Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 230000003746 surface roughness Effects 0.000 description 1
- 230000026676 system process Effects 0.000 description 1
Images
Classifications
-
- 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/5269—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/22—Matching criteria, e.g. proximity measures
-
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/06—Measuring blood flow
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8979—Combined Doppler and pulse-echo imaging systems
- G01S15/8988—Colour Doppler imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52019—Details of transmitters
- G01S7/5202—Details of transmitters for pulse systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52023—Details of receivers
- G01S7/52025—Details of receivers for pulse systems
- G01S7/52026—Extracting wanted echo signals
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/5205—Means for monitoring or calibrating
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52053—Display arrangements
- G01S7/52057—Cathode ray tube displays
- G01S7/52071—Multicolour displays; using colour coding; Optimising colour or information content in displays, e.g. parametric imaging
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/20—Analysis of motion
- G06T7/254—Analysis of motion involving subtraction of images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/30—Noise filtering
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/46—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
- A61B8/461—Displaying means of special interest
- A61B8/463—Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
-
- 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/488—Diagnostic techniques involving Doppler signals
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
- A61B8/5223—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
-
- 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/10—Image acquisition modality
- G06T2207/10132—Ultrasound image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V2201/00—Indexing scheme relating to image or video recognition or understanding
- G06V2201/03—Recognition of patterns in medical or anatomical images
- G06V2201/031—Recognition of patterns in medical or anatomical images of internal organs
Definitions
- This present embodiments relate to sparkle detection in ultrasound color flow imaging.
- Color flow imaging is susceptible to artifacts when there is a loss of correlation in the backscattered echoes.
- Correlation loss may be indicative of underlying physiologies, such as kidney stones, due to surface roughness.
- Correlation loss may instead be due to system imperfections, such as phase noises or reverberation noises.
- gain or transmit power is reduced to reduce this sparkle artifact, but at the cost of losing sensitivity.
- Using spatial variance of the color flow may identify the sparkle, but may reduce aliased flow or turbulence as well as artifacts.
- Color flow data is estimated with different pulse repetition frequency (PRF).
- PRF pulse repetition frequency
- sparkle is identified.
- Color flow images may be filtered to reduce motion while maintaining the sparkle region (e.g., kidney stone imaging) or reduce the sparkle region while maintaining motion (e.g., remove sparkle as system noise).
- a method for sparkle artifact detection in color flow First color flow data representing locations in a patient is generated.
- the first color flow data is generated with a first pulse repetition frequency.
- Second color flow data representing the locations in the patient is generated.
- the second color flow data is generated with a second pulse repetition frequency.
- a degree of similarity between the first and second color flow data is determined for each of the locations.
- a color flow image is filtered where the filtering is based on the degree of similarity as a function of location. The filtered color flow image is displayed.
- a non-transitory computer readable storage medium has stored therein data representing instructions executable by a programmed processor for sparkle detection in Doppler imaging.
- the storage medium includes instructions for: scanning a plurality of locations with ultrasound, the scanning of each of the locations providing a flow sample count of return samples; estimating first Doppler values for the locations using the return samples of the flow sample count; estimating second Doppler values for the locations using a sub-sampling of the return samples of the flow sample count; multiplying the first Doppler values by a factor that is a function of the sub-sampling; correlating the second Doppler values with the multiplied first Doppler values; and detecting the sparkle from results of the correlating.
- a system for sparkle-based processing in flow images.
- a transducer and beamformer are provided for scanning a scan region.
- a Doppler estimator is configured to estimate, from the scanning, first motion values representing locations of the scan region and second motion values representing the locations of the scan region. The first motion values are estimated with a different number of samples from the scanning than the second motion values.
- a processor is configured to identify which of the locations have sparkle from a comparison of the first and second motion values.
- FIG. 1 is a flow chart of one embodiment of a method for sparkle detection in motion imaging
- FIGS. 2A-D are example graphs of weights for filtering as a function of the similarity
- FIGS. 3A and 3B are example velocity images generated with different PRF
- FIG. 3C is a binary filter or mask removing sparkle
- FIG. 3D is the flow image of FIG. 3A after applying the filter of FIG. 3C ;
- FIGS. 4A and 4B are example velocity images generated with different PRF
- FIG. 4C is a binary filter or mask
- FIG. 4D is the flow image of FIG. 4A after filtering to enhance sparkle regions
- FIG. 5 is a block diagram of one embodiment of a system for sparkle-based processing in flow images.
- Sparkling artifacts are detected in ultrasound imaging. Sparkle, sometimes called twinkle, is more random than fluid or tissue motion. Varying pulse repetition frequency (PRF) may produce different information for sparkle artifact. Velocity estimated from blood flow in general is independent of PRF. By generating two images at different PRFs and taking a normalized cross-correlation, a weighting matrix may be applied to produce images with reduced sparkle (e.g., flow only or cleaner flow image) or with reduced flow (e.g., artifact only or enhancing stone). By using PRF variation, sparkle artifacts are detected without affecting sensitivity.
- PRF pulse repetition frequency
- FIG. 1 shows one embodiment of a method for sparkle artifact detection in color flow.
- Color flow is used to indicate spatial motion imaging, such as fluid or tissue motion. “Color” is used to distinguish from spectral Doppler imaging, where the power spectrum for a range gate is estimated.
- the color “flow” data may not be of fluid (e.g., may be of tissue motion) and/or may not represent color (e.g., may be a scalar).
- the sparkle artifact is detected to enhance return from stones or to reduce system noise for color flow imaging.
- the method is performed by the ultrasound imaging system 10 of FIG. 5 , the processor 24 , or a different system, filter, and/or processor.
- the ultrasound imaging system 10 performs the acts.
- the processor 24 controls a beamformer for scanning in act 30 , causes generation of motion data by a Doppler estimator in acts 32 - 36 , and causes generation of the image by a scan converter, graphics memory, and/or display in act 46 , but itself determines similarity in act 40 , detects sparkle in act 42 , and filters in act 44 .
- a separate filter performs act 44 .
- FIG. 1 The acts of FIG. 1 are performed in the order shown or a different order.
- acts 34 and 36 are performed in the order shown or an opposite order.
- acts 30 and 32 for one PRF are repeated after performing acts 30 and 32 for another PRF.
- act 42 is not separately performed, but instead the filtering of act 44 is performed where the similarity indicates the sparkle locations as detection.
- an image is not generated in act 46 .
- Acts 34 and 36 are one embodiment for performing act 32 , but other embodiments may be provided.
- various locations within a scan region of the patient are scanned with ultrasound.
- a patient or region is scanned in real-time with the imaging.
- the scanned region is an interior of an object, such as the patient.
- the scan is of a volume, plane, or line region.
- Scanning a plane provides data representing different locations or samples of the plane.
- the data representing the region is formed from spatial sampling of the object.
- the spatial samples are for locations distributed in an acoustic sampling grid.
- the region for the color flow scan is a region of interest smaller than a field of view or for the entire field of view.
- the ultrasound system may scan the field of view using B-mode imaging.
- the color flow region is a sub-set of that field of view.
- the user or a processor determines the region of interest in which color flow scanning occurs. Alternatively, the color flow region is the full field of view.
- Spatial samples along one or more scan lines are received. Where the transmit beam insonifies just one receive scan line, then return or echo samples along that scan line are received. Where the transmit beam insonifies multiple scan lines, then samples along the multiple scan lines may be received. To generate the samples for different receive beams at a same time, parallel receive beamformation is performed. For example, a system may be capable of forming two or more, tens, or hundreds of receive beams in parallel. Alternatively, signals received from the elements are stored and sequentially processed. Spatial samples are acquired for a plurality of receive lines in the region of interest in response to one and/or in response to sequential transmit beams.
- the scanning may be performed a plurality of times to cover the region.
- the acts are repeated to scan different portions of the region of interest. Alternatively, performing once acquires the data for the entire region of interest.
- the complete region of interest is scanned multiple times in sequence. Scanning at different times in sequence acquires spatial samples associated with motion. Any now known or later developed pulse sequences may be used.
- a sequence of at least two (flow sample count) transmissions is provided along each scan line. For example, the flow sample count is 10-20, resulting in 10-20 samples for each location.
- Any pulse repetition frequency (i.e., rate of sampling for a location), flow sample count (i.e., number of samples for a location or used to estimate), and pulse repetition interval (i.e., time between each sample acquisition for a location) may be used.
- the echo responses to the transmissions of the sequence or return samples are used to estimate velocity, energy (power), and/or variance at a given time.
- the transmissions along one line(s) may be interleaved with transmissions along another line(s).
- the spatial samples for a given time are acquired using transmissions from different times.
- the estimates from different scan lines may be acquired sequentially, but rapidly enough to represent a same time from a user perspective. Multiple scans are performed to acquire estimates for different times.
- the samples may be acquired with scanning at one PRF or flow sample count. A sub-set of these samples may be used to estimate at the lower PRF.
- the scanning is repeated. For a first repetition, the return samples are acquired at one PRF. For subsequent repetition, the return samples are acquired at the other PRF.
- the return samples or estimates of flow are acquired by transfer over a network and/or loading from memory. Data previously acquired by scanning is acquired.
- an estimator or detector generates color flow data representing locations in the patient.
- Color flow data includes estimates of velocity, energy (e.g., power), and/or variance. In one embodiment, at least velocity is estimated.
- the color flow data may be for fluid or tissue. Estimates of velocity, energy, and/or variance of tissue motion may be generated. Any motion data, whether from fluid or tissue movement, may be acquired. Color flow data is used in examples below as motion of fluid, but may alternatively or additionally be tissue motion data.
- the received spatial samples may be clutter filtered.
- the clutter filter passes frequencies associated with fluid and not tissue motion or with tissue motion and not fluid.
- the clutter filtering is of signals in the pulse sequence for estimating motion at a given time (e.g., samples of the flow sample count).
- a given signal may be used for estimates representing different times, such as associated with a moving window for clutter filtering and estimation.
- Different filter outputs are used to estimate motion for a location at different times.
- Color flow data is generated from the spatial samples.
- Doppler processing such as autocorrelation, may be used. In other embodiments, temporal correlation may be used. Another process may be used to estimate the color flow data.
- Color Doppler parameter values e.g., velocity, energy, or variance values
- the change in frequency e.g., Doppler shift
- a sequence flow sample count of two or more samples may be used to estimate the color Doppler parameter values. Estimates are formed for different groupings of received signals, such as completely separate or independent groupings or overlapping groupings. The estimates for each grouping represent the spatial location at a given time.
- the estimation is performed for the different sampled spatial locations. For example, velocities for the different locations in a plane are estimated from echoes responsive to the scanning. Multiple frames of color flow data may be acquired to represent the region of interest at different times, respectively.
- the estimates may be thresholded. Thresholds are applied to the velocities and/or powers. For example, a low velocity threshold is applied. Velocities below the threshold are removed or set to another value, such as zero. As another example, where the energy is below a threshold, the velocity value for the same spatial location is removed or set to another value, such as zero. Alternatively, the estimated velocities are used without thresholding.
- the acquired motion or color flow data is a frame of data or image representing the patient at a given time, despite being estimated from received signals over the flow sample count.
- Other data may be generated, such as B-mode data.
- a B-mode image may be overlaid with or have an incorporated region of interest showing the color Doppler velocities. Within the region of interest, locations with no flow are shown as B-mode data.
- two or more frames of data representing the same locations but at different PRF are generated.
- the locations are the same.
- the same locations may not be exactly the same due to transducer and/or tissue motion, but are generally the same.
- Acts 34 and 36 show one embodiment where the two frames with different PRF are generated using a factor of 2.
- one frame of data is generated for one PRF.
- Doppler values are estimated using the return samples of one flow sample count (e.g., 16). Velocities or other motion values are estimated for each of the locations using the PRF and corresponding flow sample count.
- the true velocity v t of a flow pixel may be expressed as:
- z is the color ensemble samples
- v d is the Doppler velocity estimate
- i is an index of the samples
- n is an integer
- the PRF for one of the frames is the full or all of the flow samples, providing the maximum PRF given the acquired samples.
- sub-sampling is used (e.g., acquire 20 samples for the flow sample count, but use just 16).
- the estimated velocities are multiplied by a factor.
- the factor is related to or the same as the amount of sub-sampling for the other PRF or is a reflection of the difference in PRFs. For example, sub-sampling the flow sample count by a factor of 2 doubles the resulting velocity values. For more exact comparison, the estimates of velocity using the full or greater flow sample count (i.e., twice the PRF) are multiplied by 2. ⁇ 2v d ⁇ mod 2 ⁇ is simply the original image multiplied by 2 and then wrapped by 2 ⁇ .
- any multiplication may be used, such as multiplying the estimated velocities by the factor.
- the velocities are multiplied by 2 using finite precision arithmetic so that the wrapping is inherent during overflow (i.e., if the original velocities are 8-bit 2s complement, performing the multiply by 2 in 2s complement arithmetic automatically achieves the desired multiplication).
- the estimated velocities are not multiplied by the factor.
- the correlation for detecting sparkle relies on variation by location, so having similar magnitudes for each location may not be needed in some cases of no wrapping.
- motion data is generated with another PRF in act 36 .
- the Doppler estimator generates a second frame of color flow data representing the locations in the patient.
- the second frame of color flow data is generated with a second pulse repetition frequency.
- a different flow sample count is used.
- the same set of return samples or a different set of return samples are used.
- the different set has a different flow sample count and/or PRF.
- sub-sampling is used. The same return samples are sub-sampled to decrease the PRF and reduce the flow sample count.
- any reduction or sub-sampling may be used.
- the flow sample count is reduced by half or a factor of 2. Every other (e.g., samples 1-16 used for one frame and only even or odd number samples used for the other frame) or other grouping results in use of half of the return samples. Every third, or other grouping using fewer than all, fewer than half, and/or fewer than used for the other frame of velocity values may be used. Any integer or real value sub-sampling factor may be used.
- a subsampled image vel subsampled is created as:
- vel subsampled ⁇ arg ⁇ ⁇ i ⁇ ⁇ z ⁇ ( 2 ⁇ ⁇ i ) ⁇ z * ⁇ ( 2 ⁇ ⁇ i - 2 ) ⁇ ⁇ ⁇ mod ⁇ ⁇ 2 ⁇ ⁇ ⁇
- the subsampled image using a factor of 2 is expected to be twice the true velocity, as represented as:
- vel subsampledAVG ⁇ AVG( ⁇ v d ⁇ mod 2 ⁇ , ⁇ v d ′ ⁇ mod 2 ⁇ where term where ⁇ is the subsample factor.
- a processor determines a degree of similarity between the color flow data estimated with the different PRF and flow sample counts. Two frames of data, such as frames of velocity, are compared to determine similarity.
- vel times2 ⁇ 2v d ⁇ mod 2 ⁇
- vel subsampled and vel times2 should have high spatial correlation.
- artifacts that are random in nature, such as sparkle their spatial correlation is low.
- the similarity is determined as a function of location.
- the similarity between data from different frames is determined for each of the locations.
- the similarity may be difference.
- the similarity is based on kernels of any size (e.g., 5 axial and 3 lateral or the 2-D kernel size is 2A+1 samples axially ⁇ 2B+1 samples laterally) centered at each location.
- the similarity is determined as a function of location.
- a minimum sum of absolute differences is used.
- a correlation such as a cross-correlation
- the Doppler values from the sub-sampled estimation are correlated with the Doppler values from the full estimation multiplied by the sub-sampling factor.
- the cross-correlation or other similarity may be normalized.
- the degree of similarity may be quantified using normalized cross-correlation (NCC) between the two signals vel subsampledAVG and vel times ⁇ . 2D Normalized cross-correlation using segments of data along the axial direction and lateral direction at zero lag is performed as follows:
- the normalized cross-correlation coefficient ⁇ (i,j) at zero lag is calculated for every sample. Other correlation or similarity measures between the frames of data estimated with different PRF or flow sample count may be used.
- the processor detects the sparkle from the similarity. Locations with a high degree of correlation between data with different PRF are not sparkle, but are motion. Locations with a low degree of correlation are sparkle. The results of the correlation identify the locations with less correlation. Any threshold, such as an empirically developed threshold for a given imaging application, may be used to distinguish sparkle from motion.
- the processor does not specifically detect the sparkle or locations with sparkle.
- the similarity as a function of location may indicate the locations without specifically identifying a given location as sparkle.
- the detection is provided by creation of the similarity as a function of location.
- the processor or a filter filters a color flow image.
- the color flow image is velocity or other motion data to be used to generate an image (e.g., scalar estimate or Doppler values) or is RGB or other data for display.
- the color flow image is created from the estimates generated with the greater PRF, such as the velocities generated with the full sample count.
- the color flow image is created from the sub-sampling.
- the color flow image is created from samples from a different scan, so is not created from the data used for the determination of similarity.
- the filtering is based on the degree of similarity as a function of location. To suppress the sparkle, the filtering more heavily reduces values at locations of less similarity. To enhance the sparkle, the filtering more heavily reduces values at locations of greater similarity. The filtering avoids altering some Doppler values, but reduces others. Alternatively, the filtering results in alteration of the magnitude of all of the values.
- the reduction is moving the velocity closer to zero.
- the negative or positive velocities are reduced to lower negative or positive velocities while maintaining the sign.
- the estimates are not signed. The reduction is moving the estimates closer to zero.
- the motion data is multiplied by weights.
- the weights are fractional values between and/or including 0 and 1. Other weighting may be used, such as adding, subtracting, or dividing by weights.
- the color flow image or other motion data is weighted by the weights as a function of location.
- the weights are mapped from the degree of similarity.
- a pixel-by-pixel or location specific weighting matrix is created.
- the degree of similarity for each location is mapped to a weight for that location.
- the weighting is target-dependent.
- the weights vary as a function of the level of similarity, providing for reduction as a function of sparkle, such as to pass motion signals and remove or attenuate sparkle artifact.
- normalized cross-correlation coefficients may range from ⁇ 1 to 1 or from 0 to 1. Two signals are identical if the cross-correlation coefficient is 1 and they are considered uncorrelated if the coefficient is near or below zero.
- Any threshold within the range of normalized values may be used. In one embodiment, the threshold is used for binary weights. Above the threshold, the weight is 1 or 0. Below the threshold, the weight is the other of 0 or 1. Equal to the threshold may be mapped to either 0 or 1. The result is a binary mask as the filter. Values associated with sparkle are the only ones remaining or the only ones removed after filtering with the binary mask. The values associated with non-sparkle are removed or remain.
- the sample value is multiplied by 0. If the coefficient is less than the threshold value ⁇ , the sample value is multiplied by 1. This case may be used for enhancing kidney stones or gallstones where sparkle indicates desired information. The locations associated with greater similarity (i.e., with motion) are removed or reduced.
- the sample value will be multiplied by 1. If the coefficient is less than the threshold value ⁇ , the sample value is multiplied by 0. The motion data for locations associated with lesser similarity (i.e., with sparkle) are removed or reduced, while motion data for motion locations is maintained.
- weight mapping than a binary mask may be used.
- the weighting matrix may be generalized instead of using a mask of 0 and 1.
- FIGS. 2A-D Some examples of similarity to weight mapping are shown in FIGS. 2A-D .
- the normalized cross-correlation coefficient is mapped to weights.
- the maps of FIGS. 2A-D are for reducing velocity or other motion data at locations of sparkle. By applying the weights of FIGS. 2A-D , locations with greater similarity are weighted more strongly than locations with lesser similarity, reducing the Doppler values for the locations of the sparkle.
- Inverse maps may be used, weighting locations with lesser similarity more strongly than locations with greater similarity to reduce the Doppler values for the locations other than sparkle.
- Other linear or non-linear mapping of similarity to weights may be used.
- the filtered color flow image is displayed.
- the ultrasound system processes the filtered frame of data to create the image. Spatial filtering, temporal filtering, scan conversion, or other image processing is performed.
- the scalar values are mapped to display values, such as mapping to color values using a velocity scale.
- the resulting image is buffered for display. The display values are provided from the buffer to the display.
- Color flow e.g., Doppler energy or Doppler velocity
- Doppler tissue motion or other motion image
- the image may include other information.
- the image is an overlay of the color flow data on B-mode data.
- the color flow data e.g., velocities
- the B-mode data is used for tissue locations or low/no flow locations.
- the image includes motion values (e.g., velocities or energy) that have been altered to remove or reduce the sparkle artifact or motion information not associated with a rough surface.
- the color flow image such as a velocity image, is generated from the Doppler values after the filtering based on similarity. As a result, the color flow image highlights stones or other rough surface targets without other motion information or highlights flow with less sparkle artifact. For example, a kidney stone or gallstone image is displayed. As another example, a velocity of fluid image is displayed.
- FIGS. 3A-D and 4 A-C show two examples.
- FIGS. 3A-B and 4 A-B show velocity images in combination with B-mode. The velocity information is color mapped. The velocity portions of these images represent frames of motion data.
- the images use estimates from a full flow sample count or higher PRF.
- the images use estimates from a sub-sampling of the same return samples by a factor of 2 (e.g., estimating velocity using every other sample of the full flow sample count).
- FIGS. 3B and 4B are from motion data estimated with a lower PRF.
- FIG. 3C shows two white regions of no detected sparkle. Flow regions with correlation below a threshold are weighted by 0, resulting in the mask for passing clean flow of FIG. 3C .
- a velocity image of FIG. 3D results.
- the image has the sparkle removed or reduced.
- B-mode information may be provided for non-motion locations.
- FIG. 4A shows a color flow image.
- a binary mask e.g., FIG. 4C ) passing (e.g., weight of 1) locations associated with sparkle and removing flow information not associated with sparkle is created by correlating the motion information of FIGS. 4A and 4B .
- the color flow or motion information for locations associated with sparkle is maintained while removing motion information for other locations.
- the resulting FIG. 4D shows color information as sparkle representing a kidney stone without flow from blood or tissue motion.
- B-mode information may be provided for non-motion locations.
- FIG. 5 shows one embodiment of a system 10 for sparkle-based processing in flow images (i.e., in color flow, tissue motion, or other motion imaging).
- the system 10 implements the method of FIG. 1 or another method.
- the motion data may be compared to indicate sparkle.
- the results of the comparison are used for processing, filtering, or other imaging based on detected sparkle locations.
- the system 10 includes a transmit beamformer 12 , a transducer 14 , a receive beamformer 16 , a memory 18 , a filter 20 , a flow estimator 22 , another memory 28 , a processor 24 , and a display 27 . Additional, different or fewer components may be provided.
- the system includes a B-mode detector.
- the flow estimator 22 and processor 24 are provided without the front-end components, such as the transmit and receive beamformers 12 , 16 .
- the memories 18 and 28 are one component.
- the system 10 is a medical diagnostic ultrasound system.
- the system 10 is a computer or workstation.
- the flow estimator 22 is part of a medical diagnostic ultrasound system or other medical imaging system
- the processor 24 is part of a separate workstation or remote system, making of the ultrasound imaging system.
- the transducer 14 is an array of a plurality of elements.
- the elements are piezoelectric or capacitive membrane elements.
- the array is configured as a one-dimensional array, a two-dimensional array, a 1.5D array, a 1.25D array, a 1.75D array, an annular array, a multidimensional array, a wobbler array, combinations thereof, or any other now known or later developed array.
- the transducer elements transduce between acoustic and electric energies.
- the transducer 14 connects with the transmit beamformer 12 and the receive beamformer 16 through a transmit/receive switch, but separate connections may be used in other embodiments.
- the transmit and receive beamformers 12 , 16 are a beamformer for scanning a region of the patient with the transducer 14 .
- the transmit beamformer 12 using the transducer 14 , transmits one or more beams to scan a region.
- Vector®, sector, linear or other scan formats may be used.
- the receive lines and/or transmit beams are distributed in the scan region.
- the receive beamformer 16 samples the receive beams at different depths. Sampling the same locations at different times obtains a sequence for flow estimation.
- the transmit beamformer 12 is a processor, delay, filter, waveform generator, memory, phase rotator, digital-to-analog converter, amplifier, combinations thereof or any other now known or later developed transmit beamformer components.
- the transmit beamformer 12 digitally generates envelope samples. Using filtering, delays, phase rotation, digital-to-analog conversion, and amplification, the desired transmit waveform is generated.
- Other waveform generators may be used, such as switching pulsers or waveform memories.
- the transmit beamformer 12 is configured as a plurality of channels for generating electrical signals of a transmit waveform for each element of a transmit aperture on the transducer 14 .
- the waveforms are unipolar, bipolar, stepped, sinusoidal or other waveforms of a desired center frequency or frequency band with one, multiple, and/or fractional number of cycles.
- the waveforms have relative delay and/or phasing and amplitude for focusing the acoustic energy.
- the transmit beamformer 12 includes a controller for altering an aperture (e.g.
- an apodization profile e.g., type or center of mass
- a delay profile across the plurality of channels e.g., a delay profile across the plurality of channels
- a phase profile across the plurality of channels center frequency, frequency band, waveform shape, number of cycles and/or combinations thereof.
- a transmit beam focus is generated based on these beamforming parameters.
- the receive beamformer 16 is a preamplifier, filter, phase rotator, delay, summer, base band filter, processor, buffers, memory, combinations thereof or other now known or later developed receive beamformer components.
- the receive beamformer 16 is configured into a plurality of channels for receiving electrical signals representing echoes or acoustic energy impinging on the transducer 14 .
- a channel from each of the elements of the receive aperture within the transducer 14 connects to an amplifier and/or delay.
- An analog-to-digital converter digitizes the amplified echo signal.
- the digital radio frequency received data is demodulated to a base band frequency. Any receive delays, such as dynamic receive delays, and/or phase rotations are then applied by the amplifier and/or delay.
- a digital or analog summer combines data from different channels of the receive aperture to form return samples for one or a plurality of receive beams.
- the summer is a single summer or cascaded summer.
- the beamform summer is configured to sum in-phase and quadrature channel data in a complex manner such that phase information is maintained for the formed beam.
- the beamform summer sums data amplitudes or intensities without maintaining the phase information.
- the receive beamformer 16 is configured to form receive beams in response to the transmit beams. For example, the receive beamformer 16 receives one, two, or more receive beams in response to each transmit beam. The receive beams are collinear, parallel and offset or nonparallel with the corresponding transmit beams. The receive beamformer 16 outputs spatial samples representing different spatial locations of a scanned region. Once the channel data is beamformed or otherwise combined to represent spatial locations along the scan lines 11 , the data is converted from the channel domain to the image data domain. The phase rotators, delays, and/or summers may be repeated for parallel receive beamformation. One or more of the parallel receive beamformers may share parts of channels, such as sharing initial amplification.
- a velocity sample group corresponds to multiple transmissions for each of a plurality of scan lines 11 .
- the number of times a substantially same spatial location, such as a scan line 11 , is scanned within a velocity sample group is the velocity or flow sample count.
- the transmissions for different scan lines 11 , different velocity sample groupings or different types of imaging may be interleaved.
- the amount of time between transmissions to a substantially same scan line 11 within the velocity sample count is the pulse repetition interval.
- the pulse repetition interval establishes the pulse repetition frequency or vice versa.
- the memory 18 is video random access memory, random access memory, removable media (e.g. diskette or compact disc), hard drive, database, corner turning memory, or other memory device for storing data or video information.
- the memory 18 is a corner turning memory of a motion parameter estimation path.
- the memory 18 is configured to store signals responsive to multiple transmissions along a substantially same scan line.
- the memory 22 is configured to store ultrasound data formatted in an acoustic grid, a Cartesian grid, both a Cartesian coordinate grid and an acoustic grid, or ultrasound data representing a volume in a 3D grid. The return samples of the flow sample count for each of a plurality of locations are stored.
- the filter 20 is a clutter filter, finite impulse response filter, infinite impulse response filter, analog filter, digital filter, combinations thereof or other now known or later developed filter.
- the filter 20 includes a mixer to shift signals to baseband and a programmable low pass filter response for removing or minimizing information at frequencies away from the baseband.
- the filter 20 is a low pass, high pass, or band pass filter.
- the filter 20 reduces velocities from fluids or alternatively reduces the influence of data from tissue while maintaining velocity information from fluids.
- the filter 20 has a set response or may be programmed, such as altering operation as a function of signal feedback or other adaptive process.
- the memory 18 and/or the filter 20 are part of the flow estimator 22 .
- the Doppler or flow estimator 22 is a Doppler processor or cross-correlation processor for estimating the color flow data.
- another device now known or later developed for estimating velocity, power (e.g., energy), and/or variance from any or various input data may be provided.
- the flow estimator 22 receives a plurality of signals associated with a substantially same location at different times and estimates a Doppler shift frequency, based on a change or an average change in phase between consecutive signals from the same location. Velocity is calculated from the Doppler shift frequency. Alternatively, the Doppler shift frequency is used as a velocity. The power and variance may also be calculated.
- Color flow data (e.g., velocity, power, and/or variance) is estimated for spatial locations in the scan region from the beamformed scan samples.
- the flow data represents a plurality of different locations in a plane.
- the color flow data is motion data for tissue and/or fluid.
- the flow estimator 22 may apply one or more thresholds to identify sufficient motion information. For example, velocity and/or power thresholding for identifying velocities is used. In alternative embodiments, a separate processor or filter applies thresholds. In other embodiments, the thresholding is applied after any motions suppression, such as by the processor 24 .
- the flow estimator 22 outputs frames of data representing the scan region at different times.
- the beamformed samples for a given flow sample count are used to estimate for a time.
- a moving window with overlap of the data is used to estimate for other times. Velocities for each location at different times are output.
- Two or more frames are estimated from the scanning for the same locations.
- the frames are for a same or different time.
- the frames represent the locations at different times. Where the samples for the same flow sample count are used to estimate the frames, then the frames represent a same time.
- the two or more frames are estimated differently.
- the number of samples used in the estimation for each frame is different.
- the PRF or flow sample count in the estimation varies from frame-to-frame. For example, 16 samples are obtained and used for estimating velocity in one frame. Every other sample (e.g., 8 samples) is used to estimate the velocity in the other frame.
- the sub-sampling factor for the other frame is 2. A different number of samples are used to estimate the velocities of the two frames.
- the flow estimator 22 or the processor 24 is configured to normalize the motion data of the two or more frames. Due to sub-sampling, the velocities estimation may be different, such as the sub-sampling by the factor of 2 leading to the results of the autocorrelation being doubled relative to the frame estimated from the full flow sample count. The motion data from the full flow sample count is multiplied to normalize. Alternatively, multiplication or normalization is not performed, and the frames are used despite having different magnitude ranges.
- the processor 24 is a digital signal processor, a general processor, an application specific integrated circuit, field programmable gate array, control processor, digital circuitry, analog circuitry, graphics processing unit, filter, combinations thereof or other now known or later developed device for implementing calculations, algorithms, programming or other functions.
- the processor 24 operates pursuant to instruction provided in the memory 28 , or a different memory for sparkle-based image processing. Additional or multiple processors may be used.
- the processor 24 is configured by software, firmware, and/or hardware.
- the processor 24 receives color flow data from the flow estimator 22 , the memory 28 , and/or another source. Using the received motion data, the processor 24 is configured to identify which of the locations of the scan region correspond to sparkle in the estimates. The processor 24 is configured to compare the motion values from the different frames. For example, a spatial cross-correlation between the data of the frames is performed. The cross-correlation indicates a level of similarity between the different PRF frames. Similarity is determined for each of multiple locations represented in the frames. The locations with less correlation between motion values estimated from different PRF are locations of sparkle. The locations with more correlation are actual motion rather than sparkle.
- the processor 24 specifically identifies locations of sparkle. Alternatively or additionally, the processor 24 uses the indication of sparkle provided by the level of similarity for filtering or other image processing.
- the level of similarity may be mapped to filter weights.
- the weights are applied to a frame of motion values, such as to the frame estimated using the greater number of samples.
- the mapping and corresponding weights reduce locations not associated with sparkle or locations associated with sparkle. Other locations are maintained the same, increased, or not reduced as mush.
- the weight matrix and/or locations of sparkle may be filtered, such as with a low pass filter, to remove outlier locations or reduce outlier weights.
- the processor 24 or other component uses the filtered motion values to generate an image.
- the frame is scan converted and color mapped.
- the resulting color values are added to a B-mode image, such as an overlay, or used alone.
- the color values are placed in a display buffer to display an image on the display 27 .
- the display 27 is a CRT, LCD, plasma, projector, monitor, printer, touch screen, or other now known or later developed display device.
- the display 27 receives RGB, other color values, or other motion values and outputs an image.
- the image may be a gray scale or color image.
- the image represents the region of the patient scanned by the beamformer and transducer 14 , but with returns from sparkle locations removed, reduced, or maintained while reducing returns for other locations.
- the display 27 displays a Doppler or other color image from the motion values as filtered.
- the resulting image may represent rough surfaces, such as from stones, without or with less information from fluid flow or tissue motion. Alternatively, the resulting image may represent fluid flow or tissue motion without or with less information from sparkle artifact.
- the memory 28 is video random access memory, random access memory, removable media (e.g. diskette or compact disc), hard drive, database, or other memory device for storing color flow or other motion data.
- the stored data is in a polar or Cartesian coordinate format.
- the memory 28 is used by the processor 24 for the various filtering, detecting, identifying, correlation, calculating, or other acts.
- Non-transitory computer readable storage media include various types of volatile and nonvolatile storage media.
- the functions, acts or tasks illustrated in the figures or described herein are executed in response to one or more sets of instructions stored in or on computer readable storage media.
- the functions, acts or tasks are independent of the particular type of instructions set, storage media, processor or processing strategy and may be performed by software, hardware, integrated circuits, firmware, micro code and the like, operating alone or in combination.
- processing strategies may include multiprocessing, multitasking, parallel processing and the like.
- the instructions are stored on a removable media device for reading by local or remote systems.
- the instructions are stored in a remote location for transfer through a computer network or over telephone lines.
- the instructions are stored within a given computer, CPU, GPU or system.
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Acoustics & Sound (AREA)
- Multimedia (AREA)
- Hematology (AREA)
- Data Mining & Analysis (AREA)
- Bioinformatics & Computational Biology (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Physiology (AREA)
Priority Applications (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US14/850,707 US11096671B2 (en) | 2015-09-10 | 2015-09-10 | Sparkle artifact detection in ultrasound color flow |
JP2016172630A JP6385992B2 (ja) | 2015-09-10 | 2016-09-05 | 超音波カラーフローにおけるスパークルアーチファクトの検出 |
DE102016116658.8A DE102016116658A1 (de) | 2015-09-10 | 2016-09-06 | Sparkleartefakt-erkennung im ultraschall-farbfluss |
CN201610812815.9A CN106529561B (zh) | 2015-09-10 | 2016-09-09 | 超声彩色流中的闪光伪像检测 |
KR1020160116702A KR102041433B1 (ko) | 2015-09-10 | 2016-09-09 | 초음파 컬러 흐름에서의 스파클 아티팩트 검출 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US14/850,707 US11096671B2 (en) | 2015-09-10 | 2015-09-10 | Sparkle artifact detection in ultrasound color flow |
Publications (2)
Publication Number | Publication Date |
---|---|
US20170071577A1 US20170071577A1 (en) | 2017-03-16 |
US11096671B2 true US11096671B2 (en) | 2021-08-24 |
Family
ID=58160924
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US14/850,707 Active 2040-04-21 US11096671B2 (en) | 2015-09-10 | 2015-09-10 | Sparkle artifact detection in ultrasound color flow |
Country Status (5)
Country | Link |
---|---|
US (1) | US11096671B2 (ja) |
JP (1) | JP6385992B2 (ja) |
KR (1) | KR102041433B1 (ja) |
CN (1) | CN106529561B (ja) |
DE (1) | DE102016116658A1 (ja) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6604272B2 (ja) * | 2016-06-06 | 2019-11-13 | コニカミノルタ株式会社 | 超音波診断装置及び超音波信号処理方法 |
CA3039666C (en) | 2016-10-28 | 2022-08-23 | Ppg Industries Ohio, Inc. | Coatings for increasing near-infrared detection distances |
US11419584B2 (en) | 2017-06-28 | 2022-08-23 | Tokitae, LLC | Devices, systems, and methods for diagnosis of pulmonary conditions through detection of b-lines in lung sonography |
WO2019038210A1 (en) * | 2017-08-21 | 2019-02-28 | Koninklijke Philips N.V. | DETECTION, PRESENTATION AND REPORTING OF B LINES IN PULMONARY ULTRASONES |
JP7052313B2 (ja) * | 2017-11-20 | 2022-04-12 | コニカミノルタ株式会社 | 超音波信号処理装置、超音波診断装置、および、超音波信号処理方法 |
EP3727159B1 (en) * | 2017-12-20 | 2024-02-28 | Verathon INC. | Echo window artifact classification and visual indicators for an ultrasound system |
EP3881236A1 (en) | 2018-11-13 | 2021-09-22 | PPG Industries Ohio Inc. | Method of detecting a concealed pattern |
US11561329B2 (en) | 2019-01-07 | 2023-01-24 | Ppg Industries Ohio, Inc. | Near infrared control coating, articles formed therefrom, and methods of making the same |
EP4021306A4 (en) * | 2019-08-28 | 2023-08-23 | View Point Medical, Inc. | ULTRASONIC MARKER DETECTION, MARKER AND ASSOCIATED SYSTEMS, METHODS AND ARTICLES |
US20220362411A1 (en) * | 2019-09-20 | 2022-11-17 | Mayo Foundation For Medical Education And Research | Non-Metallic Ultrasound-Detectable Markers |
CN111436971B (zh) * | 2020-04-20 | 2022-06-10 | 汕头市超声仪器研究所股份有限公司 | 一种超声混响伪像抑制方法 |
CN112256474B (zh) * | 2020-09-11 | 2022-12-13 | 安徽中科新辰技术有限公司 | 一种针对软件实时同步的双机热备份系统 |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000308642A (ja) | 1998-12-31 | 2000-11-07 | General Electric Co <Ge> | セグメント化およびフロー・ダイナミクスを強化するためのカラー・フロー撮像 |
US20070078347A1 (en) * | 2005-09-30 | 2007-04-05 | Siemens Medical Solutions Usa, Inc. | Ultrasound color flow imaging at high frame rates |
US20080021326A1 (en) | 2003-12-15 | 2008-01-24 | Siemens Medical Solutions Usa, Inc. | Automatic optimization for ultrasound medical imaging |
US20090253986A1 (en) | 2005-11-10 | 2009-10-08 | Bracco Research Sa | Detection of Immobilized Contrast Agent in Medical Imaging Applications Based on Flow Dynamics Analysis |
US20100185085A1 (en) * | 2009-01-19 | 2010-07-22 | James Hamilton | Dynamic ultrasound processing using object motion calculation |
CN101897600A (zh) | 2009-05-26 | 2010-12-01 | 通用电气公司 | 用于自动超声图像优化的系统和方法 |
US20110164794A1 (en) * | 2010-01-05 | 2011-07-07 | Shenzhen Mindray Bio-Medical Electronics Co., Ltd. | Methods and systems for color flow dynamic frame persistence |
US20110263967A1 (en) | 2010-04-22 | 2011-10-27 | of higher education having a principal place of bussiness | Ultrasound based method and apparatus for stone detection and to facilitate clearance thereof |
US20120046548A1 (en) * | 2010-08-17 | 2012-02-23 | Imsonic Medical, Inc. | Handheld ultrasound color flow imaging system with mechanically scanned, mechanically focused multi-element transducers |
CN103123721A (zh) | 2011-11-17 | 2013-05-29 | 重庆海扶医疗科技股份有限公司 | 一种实时减少图像中伪影的方法以及装置 |
US20130336560A1 (en) | 2012-06-14 | 2013-12-19 | Siemens Medical Solutions Usa, Inc. | Motion Artifact Suppression in Ultrasound Diagnostic Imaging |
US20150230777A1 (en) | 2014-02-17 | 2015-08-20 | Siemens Medical Solutions Usa, Inc. | Twinkle artifact suppression in ultrasound color flow |
US20150320384A1 (en) * | 2014-05-06 | 2015-11-12 | University Of Washington | Methods and Systems for Detecting an Object in a Subject with Ultrasound |
US9318770B2 (en) | 2011-08-09 | 2016-04-19 | Toyota Jidosha Kabushiki Kaisha | Battery and method of manufacturing battery |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5961462A (en) | 1998-05-18 | 1999-10-05 | Atl Ultrasound | Ultrasonic doppler imaging at high frame rates of display |
JP2008171472A (ja) * | 2007-01-09 | 2008-07-24 | Fuji Electric Device Technology Co Ltd | 情報記録媒体用ガラス基板の製造方法、磁気記録ディスクの製造方法および磁気記録ディスク |
-
2015
- 2015-09-10 US US14/850,707 patent/US11096671B2/en active Active
-
2016
- 2016-09-05 JP JP2016172630A patent/JP6385992B2/ja active Active
- 2016-09-06 DE DE102016116658.8A patent/DE102016116658A1/de active Pending
- 2016-09-09 CN CN201610812815.9A patent/CN106529561B/zh active Active
- 2016-09-09 KR KR1020160116702A patent/KR102041433B1/ko active IP Right Grant
Patent Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000308642A (ja) | 1998-12-31 | 2000-11-07 | General Electric Co <Ge> | セグメント化およびフロー・ダイナミクスを強化するためのカラー・フロー撮像 |
US20080021326A1 (en) | 2003-12-15 | 2008-01-24 | Siemens Medical Solutions Usa, Inc. | Automatic optimization for ultrasound medical imaging |
US20070078347A1 (en) * | 2005-09-30 | 2007-04-05 | Siemens Medical Solutions Usa, Inc. | Ultrasound color flow imaging at high frame rates |
US20090253986A1 (en) | 2005-11-10 | 2009-10-08 | Bracco Research Sa | Detection of Immobilized Contrast Agent in Medical Imaging Applications Based on Flow Dynamics Analysis |
US20100185085A1 (en) * | 2009-01-19 | 2010-07-22 | James Hamilton | Dynamic ultrasound processing using object motion calculation |
CN101897600A (zh) | 2009-05-26 | 2010-12-01 | 通用电气公司 | 用于自动超声图像优化的系统和方法 |
US20110164794A1 (en) * | 2010-01-05 | 2011-07-07 | Shenzhen Mindray Bio-Medical Electronics Co., Ltd. | Methods and systems for color flow dynamic frame persistence |
US20110263967A1 (en) | 2010-04-22 | 2011-10-27 | of higher education having a principal place of bussiness | Ultrasound based method and apparatus for stone detection and to facilitate clearance thereof |
JP2013527782A (ja) | 2010-04-22 | 2013-07-04 | ザ ユニバーシティ オブ ワシントン スルー イッツ センター フォー コマーシャライゼーション | 結石を検出し、その除去を促進する超音波ベースの方法及び装置 |
US20120046548A1 (en) * | 2010-08-17 | 2012-02-23 | Imsonic Medical, Inc. | Handheld ultrasound color flow imaging system with mechanically scanned, mechanically focused multi-element transducers |
US9318770B2 (en) | 2011-08-09 | 2016-04-19 | Toyota Jidosha Kabushiki Kaisha | Battery and method of manufacturing battery |
CN103123721A (zh) | 2011-11-17 | 2013-05-29 | 重庆海扶医疗科技股份有限公司 | 一种实时减少图像中伪影的方法以及装置 |
US20130336560A1 (en) | 2012-06-14 | 2013-12-19 | Siemens Medical Solutions Usa, Inc. | Motion Artifact Suppression in Ultrasound Diagnostic Imaging |
US20150230777A1 (en) | 2014-02-17 | 2015-08-20 | Siemens Medical Solutions Usa, Inc. | Twinkle artifact suppression in ultrasound color flow |
US20150320384A1 (en) * | 2014-05-06 | 2015-11-12 | University Of Washington | Methods and Systems for Detecting an Object in a Subject with Ultrasound |
Non-Patent Citations (15)
Title |
---|
Chun, Doo-Man et al: "Web-based material database for material selction and its application programming interface (API) for CAD"; Key Engineering Materials vol. 345-346; pp. 1593-1596; Dec. 31, 2007; Trans Tech Publications; Switzerland. |
Cicek K et al: "Multiple attribute decision-making solution to material selection problem based on modified fuzzy axiomatic design-model selection interface algorithm", Materials and Design, London, GB, vol. 31, No. 4, pp. 2129-2133, XP026821195, ISSN: 0261-3069. |
CICEK, K. ; CELIK, M.: "Multiple attribute decision-making solution to material selection problem based on modified fuzzy axiomatic design-model selection interface algorithm", MATERIALS AND DESIGN, LONDON, GB, vol. 31, no. 4, 1 April 2010 (2010-04-01), GB, pages 2129 - 2133, XP026821195, ISSN: 0261-3069 |
Gao, Jing, et al. "Flow turbulence or twinkling artifact? A primary observation on the intrarenal color Doppler sonography." Clinical imaging 34.5, pp. 355-360, 2010. |
Hindi, Ammar et al.; "Artifacts in diagnostic ultrasound." In: Reports in Medical Imaging 6, 2013, S. 29-48. |
Jamzad, Amoon, and Seyed Kamaledin Setarehdan. "Simulation of the twinkling artifact in color flow doppler sonography: a phase noise hypothesis validation." Signal and Image Processing Applications (ICSIPA), Nov. 2011. |
Kim, Hyun Cheol, et al. "Color Doppler twinkling artifacts in various conditions during abdominal and pelvic sonography." Journal of Ultrasound in Medicine 29.4, pp. 621-632, 2010. |
Lan, Yuanpei, et al. "A web-based computer-aided material-selection system for aircraft design." 2010 WASE International Conference on Information Engineering. vol. 3. IEEE, 2010. |
Mohammad Hosein Fazel Zarandi et al: "A material selection methodology and expert system for sustainable product design"; The International Journal of Advanced Manufacturing Technology; Springer; Berlin, DE; vol. 57; No. 9-12; pp. 885-903; XP019972388; ISSN: 1433-3015; DOI: 10.I007/S00170-11I-3362-Y. |
MOHAMMAD HOSEIN FAZEL ZARANDI; SAEID MANSOUR; SEID ALI HOSSEINIJOU; MILAD AVAZBEIGI: "A material selection methodology and expert system for sustainable product design", THE INTERNATIONAL JOURNAL OF ADVANCED MANUFACTURING TECHNOLOGY, SPRINGER, BERLIN, DE, vol. 57, no. 9 - 12, 12 May 2011 (2011-05-12), Berlin, DE, pages 885 - 903, XP019972388, ISSN: 1433-3015, DOI: 10.1007/s00170-011-3362-y |
Office Action dated May 29, 2018 in corresponding Korean Patent Application No. 2016-0116702. |
Sakundarini Novita et al: "Optimal multi-material selection for lightweight design of automotive body assembly incorporating recyclability"; Materials and Design; London, GB; vol. 50; pp. 846-857; XP028552632; ISSN: 0261-3069; DOI: 10.1016/J.MATDES.2013.03.085. |
SAKUNDARINI NOVITA; TAHA ZAHARI; ABDUL-RASHID SALWA HANIM; GHAZILA RAJA ARIFFIN RAJA: "Optimal multi-material selection for lightweight design of automotive body assembly incorporating recyclability", MATERIALS AND DESIGN, LONDON, GB, vol. 50, 4 April 2013 (2013-04-04), GB, pages 846 - 857, XP028552632, ISSN: 0261-3069, DOI: 10.1016/j.matdes.2013.03.085 |
Tanabe, Masayuki et al.; "Effect of pulse repetition frequency on microcalcification detection in color flow imaging." In: Japanese Journal of Applied Physics 53, 07KF15, 2014, pp. 07KF15-1-07KF15-5. |
Tsao, Teng-Fu, et al. "Correlation study of the strength of the color Doppler twinkling artifact with the roughness of the reflecting surface and the Doppler angles." Journal of Medical Ultrasound 12.4, pp. 119-124, 2004. |
Also Published As
Publication number | Publication date |
---|---|
CN106529561A (zh) | 2017-03-22 |
KR102041433B1 (ko) | 2019-11-06 |
DE102016116658A1 (de) | 2017-03-16 |
US20170071577A1 (en) | 2017-03-16 |
JP6385992B2 (ja) | 2018-09-05 |
CN106529561B (zh) | 2020-03-17 |
JP2017051617A (ja) | 2017-03-16 |
KR20170035796A (ko) | 2017-03-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11096671B2 (en) | Sparkle artifact detection in ultrasound color flow | |
US9420997B2 (en) | Motion artifact suppression in ultrasound diagnostic imaging | |
CN107616809B (zh) | 医学诊断超声中的组织表征 | |
US9554770B2 (en) | High pulse repetition frequency for detection of tissue mechanical property with ultrasound | |
US9579083B2 (en) | Automatic doppler gate positioning in spectral doppler ultrasound imaging | |
US9468421B2 (en) | Visualization of associated information in ultrasound shear wave imaging | |
US9332962B2 (en) | Ultrasound ARFI displacement imaging using an adaptive time instance | |
US9011338B2 (en) | Gap filling for spectral doppler ultrasound | |
US20130258805A1 (en) | Methods and systems for producing compounded ultrasound images | |
US9028413B2 (en) | Prediction-based flow estimation for ultrasound diagnostic imaging | |
US10945701B2 (en) | Segment-based flash suppression in ultrasound color flow | |
US10799208B2 (en) | Compressional sound speed imaging using ultrasound | |
US11166699B2 (en) | Diffraction correction for attenuation estimation in medical diagnostic ultrasound | |
KR101939644B1 (ko) | 의료 진단 초음파 점탄성 이미징에서의 디프랙션 소스 보상 | |
US20230012638A1 (en) | Adaptive clutter filtering in acoustic radiation force-based ultrasound imaging | |
US8394027B2 (en) | Multi-plane/multi-slice processing for 2-D flow imaging in medical diagnostic ultrasound | |
US9320498B2 (en) | Twinkle artifact suppression in ultrasound color flow | |
US11717256B2 (en) | Motion independence in acoustic radiation force impulse imaging | |
US20140336510A1 (en) | Enhancement in Diagnostic Ultrasound Spectral Doppler Imaging | |
US11129598B2 (en) | Calibration for ARFI imaging | |
CN108338808B (zh) | 使用相干的剪切速度成像 | |
US11141138B2 (en) | Kalman filtering for flash artifact suppression in ultrasound imaging | |
US20180310918A1 (en) | Variable focus for shear wave imaging | |
US10624611B2 (en) | Dynamic power doppler imaging in ultrasound |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: SIEMENS MEDICAL SOLUTIONS USA, INC., PENNSYLVANIA Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:SEO, CHI HYUNG;WONG, KING YUEN;REEL/FRAME:036708/0222 Effective date: 20150902 |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: RESPONSE TO NON-FINAL OFFICE ACTION ENTERED AND FORWARDED TO EXAMINER |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: NON FINAL ACTION MAILED |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: RESPONSE TO NON-FINAL OFFICE ACTION ENTERED AND FORWARDED TO EXAMINER |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: FINAL REJECTION MAILED |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: FINAL REJECTION MAILED |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: RESPONSE AFTER FINAL ACTION FORWARDED TO EXAMINER |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: NOTICE OF ALLOWANCE MAILED -- APPLICATION RECEIVED IN OFFICE OF PUBLICATIONS |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: PUBLICATIONS -- ISSUE FEE PAYMENT RECEIVED |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |