WO2017144716A1 - Filtres anti-fouillis pour une imagerie élastographique et autre imagerie de déformation d'ultrasons - Google Patents

Filtres anti-fouillis pour une imagerie élastographique et autre imagerie de déformation d'ultrasons Download PDF

Info

Publication number
WO2017144716A1
WO2017144716A1 PCT/EP2017/054436 EP2017054436W WO2017144716A1 WO 2017144716 A1 WO2017144716 A1 WO 2017144716A1 EP 2017054436 W EP2017054436 W EP 2017054436W WO 2017144716 A1 WO2017144716 A1 WO 2017144716A1
Authority
WO
WIPO (PCT)
Prior art keywords
displacement
image
strain
tissue
images
Prior art date
Application number
PCT/EP2017/054436
Other languages
English (en)
Inventor
Abhay Vijay PATIL
Francois Guy Gerard Marie Vignon
Sheng-Wen Huang
Scott William DIANIS
Karl Erhard Thiele
Original Assignee
Koninklijke Philips N.V.
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Koninklijke Philips N.V. filed Critical Koninklijke Philips N.V.
Priority to US16/077,558 priority Critical patent/US20190029651A1/en
Publication of WO2017144716A1 publication Critical patent/WO2017144716A1/fr

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/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
    • A61B8/5276—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts due to motion
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0883—Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of the heart
    • 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/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
    • 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
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5238—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image
    • A61B8/5246—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image combining images from the same or different imaging techniques, e.g. color Doppler and B-mode
    • 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/8909—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration
    • G01S15/8915—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array
    • 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/52036—Details of receivers using analysis of echo signal for target characterisation
    • G01S7/52042—Details of receivers using analysis of echo signal for target characterisation determining elastic properties of the propagation medium or of the reflective target
    • 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/52077—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 with means for elimination of unwanted signals, e.g. noise or interference
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0825—Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of the breast, e.g. mammography
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/13—Tomography
    • A61B8/14—Echo-tomography
    • 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/467—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient characterised by special input means
    • 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
    • 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

Definitions

  • This invention relates to medical diagnostic ultrasound systems and, in particular, to ultrasound systems which assess tissue stiffness and function through the assessment of tissue deformation.
  • Elastography is the assessment of the elastic properties of tissue in the body. It has been found that the stiffness of tissue in the body can give an indication of whether the tissue may be malignant or benign.
  • the female breast for instance, can contain a variety of different lumps, cysts, and other growths, some of which may be malignant and some of which may be benign.
  • ultrasound is frequently used to assess tissue characteristics to determine whether to biopsy suspect tissue.
  • Elastography can be performed to determine whether the breast contains softer or harder (stiffer) regions. Since stiffer tissue correlates more greatly with malignant masses, the identification of regions of stiffer tissue can indicate a need to make a definitive diagnosis by biopsy.
  • a problem posed by elastography is the need to measure quantifiable properties of tissue noninvasively within the body. This means that the properties of the target tissue cannot be measured directly at the site of the tissue, but only through measurements made at the surface of the body through intervening tissues. Accordingly it is desirable to simplify the problem and make certain approximations and assumptions that will lead to valid data and analyses.
  • One set of assumptions that is frequently made is that the tissue being examined is homogeneous and isotropic. These assumptions enable certain property of materials equations to be applied to the problem, Poisson's ratio and Young's modulus.
  • Poisson's ratio is the ratio, when a sample is stretched or compressed in a given direction, of the expansion or contraction (strain) normal to the stretching or compressing force, to the expansion or contraction axially in the direction of the force.
  • a related measure is Young's modulus, which is a measure of stiffness, and is defined as the ratio of the uniaxial stress (pressure) applied to a sample over the resulting uniaxial strain (deformation).
  • the stress component at target tissue is generally unknown and difficult to measure noninvasively.
  • strain tissue deformation
  • characteristics of the myocardium such as the speckle characteristic of localized tissue resulting from the coherent nature of ultrasonic imaging.
  • reflections and reverberations from the chest wall, ribs and lungs can give rise to clutter in the echo data used for imaging.
  • clutter can manifest itself as static or slow-moving features that are superimposed over and interfere with the desired echo signals from the myocardial tissue. Clutter negatively affects image quality and hinders segmentation and automatic tracking algorithms used for the imaging of cardiac deformation.
  • the clutter or wall filter for blood flow imaging is a high pass filter which passes the higher velocity blood flow signals while filtering out the generally much stronger signals returned from slowly moving tissue such as the walls of the heart chambers and blood vessels.
  • the clutter filters used to image blood flow cannot be implemented for cardiac tissue deformation imaging for two reasons. First, the ensemble length for cardiac tissue Doppler is very low due to the low motional velocities of tissue, and there are insufficient signal samples to adequately operate a multi-tap filter such as those used in colorflow.
  • tissue of interest in cardiac deformation imaging is generally static or slow-moving for a significant portion of the cardiac cycle, which results in the cancellation by a wall filter of the signal of interest for the imaging of tissue deformation. Elimination of the signal of interest renders such a filter inapposite for the imaging or analysis of tissue motion.
  • a clutter filtering technique and apparatus which can remove clutter from ultrasonic signals returned from slowly moving tissue.
  • the present invention provides an ultrasonic diagnostic imaging system which reduces clutter in images of moving tissue.
  • the system includes a source of successive frames of ultrasonic image data of moving tissue, a clutter filter configured to filter the successive frames to reduce clutter, thereby generating clutter-reduced successive frames of the moving tissue; and a tissue displacement estimates selector.
  • the tissue displacement estimates selector is configured to
  • the inventive technique and apparatus receive successive frames of ultrasonic image data of the target anatomy, preferably at a relatively high frame rate of acquisition.
  • the data can be operated upon as r.f. data or as detected (e.g. , envelope detected) data.
  • Successive frames of data are clutter filtered in slow-time, that is, as a function of the acquisition frame rate, to eliminate signal data from stationary objects. Since the subject signals are from slowly moving tissue, shorter length filters (fewer taps) are generally employed.
  • Frame to frame tissue displacement is then estimated for both the filtered and unfiltered data.
  • the displacement estimate judged to be of the highest quality, e.g. , the best frame to frame correlation, is then selected as the displacement estimate to be used for imaging.
  • discontinuities in the image data due to the use of both filtered and unfiltered data can be removed by a spatial or spatial/temporal filter.
  • the moving tissue data thus processed is then used for analysis such as for the production of strain images as described below.
  • FIGURE 1 illustrates in block diagram form an ultrasonic diagnostic imaging system constructed in accordance with the principles of the present invention.
  • FIGURE 2 illustrates a corner-turning data table which can be used in an implementation of the present invention.
  • FIGURE 3 illustrates a two-tap FIR filter suitable for use in an implementation of the present invention.
  • FIGURE 4 is a flowchart of an implementation of the clutter filtering technique of the present invention.
  • An ultrasound probe 10 has an array transducer 12 which transmits ultrasound waves to and receives echoes from a region of the body.
  • the array transducer can be a one-dimensional array of transducer elements or a two-dimensional array of transducer elements for scanning a two dimensional image field or a three dimensional image field in the body.
  • the elements of the array transducer are driven by a transmit beamformer 16 which controls the steering, focusing and penetration of transmit beams from the array.
  • a receive beamformer 18 receives echoes from the transducer elements and combines them to form coherent echo signals from points in the image field.
  • the transmit and receive beamformers are coupled to the transducer array elements by transmit/receive switches 14 which protect sensitive receive circuitry during transmission.
  • a beamformer controller 20 synchronizes and controls the operation of the beamformers.
  • the received echo signals are demodulated into quadrature (I and Q) samples by a quadrature bandpass (QBP) filter 22.
  • the QBP filter can also provide band limiting and bandpass filtering of the received signals.
  • the received signals may then undergo further signal processing such as harmonic separation and frequency compounding by a signal processor 24.
  • the processed echo signals are stored in a frame memory 30 where the echo data is stored as a corner turning memory for use in a clutter filter of the present invention as discussed below.
  • the echo signals are also applied to a B mode processor 26 which performs amplitude detection of the echo signals by the equation ( I 2 + Q 2 ) 1 2 for formation of a B mode image, and to a Doppler processor 28 for Doppler shift detection of flow motion at points in the image field.
  • the outputs of the B mode processor 26 and the Doppler processor 28 are coupled to an image processor 42.
  • the image processor processes the image data from these sources, e.g. , by scan conversion to a desired image format, image overlay, etc. , and produces an image for display on a display 50.
  • the frame memory 30 stores successive time- separated samplings of the image field on a spatial basis.
  • the frame rate of the image acquisition is preferably at least 150 frames per second to avoid aliasing as per the Nyquist criterion.
  • a higher frame rate of at least 800Hz can also be used for aliasing phase shift (aka Doppler) estimators to estimate the motion.
  • Doppler phase shift estimators
  • For imaging other anatomy with different motional characteristics other frame rates may be used.
  • the image data of these frames is filtered by clutter filter 42 which filters the data as a function of the frame acquisition rate: samples (pixels) of the image field acquired from the same spatial location in different images are high-pass filtered by the clutter filter 42.
  • the clutter filter has a cutoff frequency around DC so as to reject static tissue motion.
  • Short-length FIR filters have been found effective for this filtering of tissue motion, such as a two-tap filter with tap weights of [+1] and [-1], or a three-tap filter with weights of [+1], [-2] and [+1] .
  • the number of taps used depends on the expected range of motion of the tissue, with faster motion requiring shorter filters, and the bandwidth around DC that one wishes to reject. In some embodiments, more than one clutter filter can also be used.
  • the filtered frame data and the unfiltered frame data from the frame memory 30 are applied to a displacement estimates selector 44.
  • displacement estimates selector estimates the tissue displacement indicated for the frame-to-frame data for both the clutter-filtered and the unfiltered image data. There can be one displacement estimate for the unfiltered data and one displacement estimate per clutter filter used to filter the data. As described further herein, only one clutter filter is described for simplicity, but the extension to using several clutter filters would be readily understood by one of ordinary skill in the art.
  • the inputs to the displacement estimates selector in the drawing show the input of the data of one image 1; and a subsequent image I i+i for the unfiltered images and the input of the corresponding filtered images I'i and I' i+i .
  • the algorithm used to estimate the frame-to-frame displacements also produces a quality metric for the displacement estimates.
  • This metric may take the form of the local normalized cross correlation coefficient from a correlation of the successive pixel data. Another possibility is to calculate the local normalized cross-correlation coefficient between frames aligned according to the displacement estimates. Other metrics which can be used include the non- normalized correlation coefficient, sum of absolute differences, or sum of square differences.
  • the displacement estimates selector 44 selects a displacement estimate to be used for the final image, either that from the clutter- filtered data, that from the unfiltered data, or a combination thereof. The selection can be made based on the quality metric, the more favorable metric indicating the displacement to be chosen. This selection is independently made for each pixel of the final image. For example, if at one pixel location the filtered data yields a higher cross correlation than that of the non-filtered data, the final displacement estimate is the one that was estimated from the filtered data.
  • a spatial filter 46 which may also be implemented as a spatial/temporal filter. Four-point or eight-point spatial filtering may be employed, as well as other spatial filtering or averaging techniques.
  • the images produced by the spatial filter 46 are used for the calculation of strain by a strain estimator 32, operating on the frame-to-frame displacement of particles in the image field; strain is calculated as a spatial derivative of displacement.
  • Strain may be calculated from radio frequency (r.f.) or baseband I and Q data, and may also be calculated from amplitude-detected (B mode) or tissue Doppler data. Any form of strain calculation such as strain, the ratio of lateral to axial strain, and strain velocity estimation may be employed. For instance, the echoes received at a common point in consecutive frames may be correlated to estimate displacement at the point. If no motion is present at the point, the echoes from consecutive frames will be the same. If motion is present, the echoes will be different and the motion vector indicates the displacement.
  • US Pat. 6,558,324 (Von Behren et al.) describes both amplitude and phase sensitive techniques for estimating strain and employs speckle tracking for strain estimation through block matching and correlation.
  • US Pat. 6,099,471 (Torp et al.) describes the estimation of strain velocity calculated as a gradient of tissue velocity.
  • US Pat. 5,800,356 (Criton et al.) describes the use of the Doppler vector to select points for strain estimation in the direction of the tissue motion.
  • a phase-sensitive technique is used since, as recognized by Von Behren et al., r.f. data will typically yield the most accurate estimates of strain.
  • Another reason for the preference of strain estimation with phase-sensitive techniques is that the slight motion produced by physiological activity and even from the small, virtually imperceptible motion occurring while holding a probe against the body can be sensed and used to estimate strain by the strain estimator 32.
  • the strain estimator 32 produces an estimated strain value at each point
  • strain estimator (pixel location) in the image field, and these values are spatially organized by the strain estimator as a strain image of the image field.
  • the strain image is coupled to the image processor 42, as are the outputs of the B mode processor 26 and the Doppler processor 28, and the strain image is displayed on display 50.
  • the strain image is also used to produce a strain ratio image at 36.
  • the strain ratio image is produced by dividing each strain value of the strain image by a strain value for normal tissue. This value may be provided automatically as by averaging or taking the mean or median value of a plurality of strain values, such as the strain values in a region in a corner of an image (on the assumption that the user will position suspect target anatomy in the center of the image.)
  • the strain value for normal tissue is taken from an indication of normal tissue in an image which is indicated by a user.
  • the strain ratio image is responsive to a reference cursor from a control panel 40 manipulated by the user to indicate a point or region of normal tissue in an image.
  • Each strain value in the strain image is divided by a strain value of normal tissue to produce a strain ratio image at 36.
  • the user may manipulate a control of the control panel 40 to set a threshold or range of values against which the strain ratio values are compared. Strain ratio values which exceed the threshold or the range of values are uniquely highlighted in the strain ratio image as by displaying them with unique colors or brightnesses. For example the user can set the threshold at 5, and all points in the strain ratio display with a value of 5 or greater, indicative of high stiffness, can be displayed in a bright red color. The user can quickly spot suspect regions in the image from the distinctive bright red color in the strain ratio image.
  • the user can manipulate a fade control on the control panel 40 which is coupled to the image processor 42.
  • the fade control enables the user to adjust the relative transparency of the B mode and strain ratio images.
  • the user can fade the strain ratio image to be completely transparent to see the corresponding structural image alone, or can fade the B mode tissue image to see only the strain ratio image, or an intermediate transparency setting for the two overlaid images.
  • FIGURE 2 illustrates a useful way of organizing the image data in frame memory 30, which in this example is that of a corner turning memory.
  • a corner turning memory is typically used to store echo data used for Doppler processing as it organizes the data in rows of common spatial locations, across which temporally different samples are stored.
  • each row of data comprises an ensemble of data for a spatial location (pixel) which may be Fourier-processed to yield a Doppler estimate.
  • the vertical direction (depth z) is the "fast time" direction in which echo signals received by sampling from along one transmit- receive A-line are stored.
  • the fast time sampling rate in a typical ultrasound system is on the order of many megaHertz.
  • each I,Q value is the fast time subscript, indicating the sequence in which the I,Q sample values are received along the line.
  • the horizontal axis of the data is the "slow time" axis indicating the pulses P n which generate the echoes from transmitting at successive times along the same spatial A-line, which is typically denominated in Hertz (frames per second.)
  • the second subscript of each I,Q value indicates the slow time sequence of each I,Q sample of a row of data, each row containing the data from temporal sampling of the same pixel location in the image field.
  • the data table of FIGURE 2 illustrates the quadrature data acquired from eight samples in depth along an A-line of an image field (an actual ultrasound system may sample a line scores or hundreds of times,) and the A-line has been sampled eight times in succession by transmitting eight temporally successive pulses P 1 -P 8 along the A-line.
  • the frame memory 30 will contain other data tables from the acquisition of echo signal data from all A- lines of an image.
  • a typical two-dimensional image may be formed by 128 adjacent A-lines, for example, and thus the frame memory for such an image will contain 128 tables of data like the one shown in FIGURE 2.
  • the frame memory 30 is generally an addressable digital random access memory, for instance.
  • FIGURE 3 illustrates a two-tap FIR filter 60 suitable for use in clutter filter 42.
  • a signal from a spatial location in one image pixel value, S(t,n) is applied to one tap of the filter and the pixel value of the same spatial location in a temporally different image, S(t,n-k) is applied to the other tap of the filter.
  • a signal S is denominated by a fast time acquisition time t and a slow time of acquisition n or n-k, respectively.
  • the pixel values at a common spatial location of two successively acquired images are applied to the taps of the filter and weighted by weights of +1 and - 1, respectively, as stated for one of the example filters described above.
  • Multipliers 62 and 64 multiply the image values by the weights and the weighted image values are summed by a summer 70 to produce the filtered output value of the clutter filter.
  • different length filters can be used; a three-tap
  • implementation would have a third tap with a third multiplier for a third temporally different pixel value, which is summed with the others by the summer 70.
  • FIGURE 4 is a flowchart illustration of the operation of frame memory 30, clutter filter 42, displacement estimates selector 44, and spatial filter 46 of FIGURE 1 in accordance with the principles of the present invention.
  • Blocks 1; and I' i+1 represent the frame data stored in frame memory 30 from the acquisition of successive images of the same image field of a body.
  • blocks I'i and I' i+1 representing the same images after clutter filtering.
  • block D i representing the pixel-by- pixel displacement estimates caused by a motional difference between the original two images I; and Ii + i.
  • a quality metric block G containing the cross-correlation coefficients of the cross- correlation of the two images.
  • the displacement estimates selector selects for output the displacement data with the best quality metric C, on a pixel-by- pixel basis, which is shown as block D"i.
  • the displacement data image is then smoothed by spatial filtering to produce the final image D"'i for strain estimation processing.
  • the operation of the clutter filter 42, the displacement estimates selector 44, and the spatial filter are all performed by digital software running on a microprocessor and programmed to execute algorithms that effect the desired processing of the acquired echo signal data stored in the frame memory 30.
  • Typical algorithms which may be employed are the following.
  • the local signal of one A-line for multiple frames can be described as
  • Tissue displacements from frame to frame can be locally estimated by computing the normalized cross-correlation function of multiple lags:
  • the high-pass-filtered version of is first calculated (in slow-time) according to
  • coefficients are those of a high-pass filter.
  • the objective is to eliminate stationary signal components from the input signal.
  • the displacement that maximizes is another estimate of the interframe displacement (the subscript stands for "high-pass").
  • cases may be a better estimate of the real displacement (i.e., when there is little inter-frame displacement and the high-pass filtering operation destroys all useful signals); in other cases may be a better estimate (i.e. , when there is a presence of stationary tissue clutter that can be eliminated by the high-pass filtering operation).
  • the objective is to find, for each pixel, which of these two estimates is the best. This is simply done by comparing the corresponding values of the cross-correlation functions and picking the estimate that gives the highest value, i.e. :
  • the final estimated displacement can be a blend of the two estimates according to where/ and g are blending functions of the correlation coefficients.
  • a blend or average of the two values may be the desired result.
  • the process is repeated for all image pixels and frames at which an estimate of the inter-frame displacement is desired.
  • the final output is a 2D or 3D image of displacement.
  • the map of displacement data is finally smoothed by applying a low-pass spatial (and possibly temporal) filter:
  • the sum can also be extended over slow-time and over the other spatial dimensions for added robustness, e.g. ,
  • the ultrasound system components and processors of FIGURE 1 may be implemented utilizing any combination of dedicated hardware boards, DSPs, microprocessors, etc. and software programs stored on a system disk drive or solid state memory.
  • the functionality and components of the system may be implemented utilizing an off-the-shelf PC with a single microprocessor or multiple microprocessors, with the functional operations distributed between the processors.
  • the functionality of FIGURE 1 may be implemented utilizing a hybrid configuration in which certain modular functions are performed utilizing dedicated hardware, while the remaining modular functions are performed utilizing an off-the shelf PC, software and the like.
  • FIGURE 1 An example of this configuration is a probe containing the transducer array and microbeamformer to produce beamformed echo signals, which are then further processed to produce images entirely by software programs of a tablet computer, on which the final images are displayed.
  • the Philips Healthcare Visiq ultrasound system is an example of such a system implementation, in which all of the ultrasound system functionality after beamforming is performed by software executed by the tablet microprocessor.
  • the various functions of the blocks shown in FIGURE 1 also may be implemented as software modules within a processing unit.
  • the various embodiments and/or components for example, the modules, or components and controllers therein, also may be implemented as part of one or more computers or microprocessors.
  • the computer or processor may include a computing device, an input device, a display unit and an interface, for example, for accessing the Internet.
  • the computer or processor may include a
  • the microprocessor may be connected to a communication bus, for example, to access a PACS system.
  • the computer or processor may also include a memory.
  • the memory may include Random Access Memory (RAM) and Read Only Memory (ROM).
  • the computer or processor further may include a storage device, which may be a hard disk drive or a removable storage drive such as a floppy disk drive, optical disk drive, solid-state thumb drive, and the like.
  • the storage device may also be other similar means for loading computer programs or other instructions into the computer or processor.
  • the term "computer” or “module” may include any processor-based or microprocessor-based system including systems using microcontrollers, reduced instruction set computers (RISC), ASICs, logic circuits, and any other circuit or processor capable of executing the functions described herein.
  • RISC reduced instruction set computers
  • ASIC application specific integrated circuit
  • logic circuits logic circuits, and any other circuit or processor capable of executing the functions described herein.
  • the above examples are exemplary only, and are thus not intended to limit in any way the definition and/or meaning of the term
  • the computer or processor executes a set of instructions that are stored in one or more storage elements, in order to process input data.
  • the storage elements may also store data or other information as desired or needed.
  • the storage element may be in the form of an information source or a physical memory element within a processing machine.
  • the set of instructions may include various commands that instruct the computer or processor as a processing machine to perform specific operations such as the methods and processes of the various embodiments of the invention.
  • the set of instructions may be in the form of a software program.
  • the software may be in various forms such as system software or application software and which may be embodied as a tangible and non-transitory computer readable medium.
  • the software may be in the form of a collection of separate programs or modules, a program module within a larger program or a portion of a program module.
  • the software also may include modular programming in the form of object-oriented programming.
  • the processing of input data by the processing machine may be in response to operator commands, or in response to results of previous processing, or in response to a request made by another processing machine.

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Medical Informatics (AREA)
  • Veterinary Medicine (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Surgery (AREA)
  • Molecular Biology (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Cardiology (AREA)
  • Physiology (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

La présente invention concerne un système et un procédé d'imagerie de diagnostic par ultrasons qui filtrent les données d'images ultrasonores successives de tissus mobiles. Au moins deux images successives sont filtrées contre les fouillis pour éliminer les fouillis dans les données d'images d'objets stationnaires. Les images successives et les images successives filtrées sont utilisées pour produire des images de déplacement d'un mouvement de tissus, image par image, et une mesure de qualité de la qualité du déplacement estimé est obtenue. L'image de déplacement présentant la meilleure mesure de qualité est sélectionnée en vue d'un traitement ultérieur tel que la production d'images élastographiques. L'image de déplacement sélectionnée peut en outre être filtrée spatialement pour réduire les artefacts aux bords.
PCT/EP2017/054436 2016-02-26 2017-02-27 Filtres anti-fouillis pour une imagerie élastographique et autre imagerie de déformation d'ultrasons WO2017144716A1 (fr)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US16/077,558 US20190029651A1 (en) 2016-02-26 2017-02-27 Clutter filters for strain and other ultrasonic deformation imaging

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US201662300233P 2016-02-26 2016-02-26
US62/300233 2016-02-26

Publications (1)

Publication Number Publication Date
WO2017144716A1 true WO2017144716A1 (fr) 2017-08-31

Family

ID=58162611

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/EP2017/054436 WO2017144716A1 (fr) 2016-02-26 2017-02-27 Filtres anti-fouillis pour une imagerie élastographique et autre imagerie de déformation d'ultrasons

Country Status (2)

Country Link
US (1) US20190029651A1 (fr)
WO (1) WO2017144716A1 (fr)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20190108071A (ko) * 2018-03-13 2019-09-23 지멘스 메디컬 솔루션즈 유에스에이, 인크. 음향 방사력-기반 초음파 이미징에서의 적응식 클러터 필터링

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10456115B2 (en) * 2013-04-22 2019-10-29 Samsung Electronics Co., Ltd. Ultrasound system and clutter filtering method thereof
US11006831B2 (en) * 2018-03-05 2021-05-18 Mothership Medical, Inc. Wireless biological monitoring
US11294052B2 (en) * 2018-06-18 2022-04-05 The Board Of Trustees Of The University Of Illinois Ultrasonic imaging with clutter filtering for perfusion
JP6924236B2 (ja) * 2019-09-26 2021-08-25 ゼネラル・エレクトリック・カンパニイ 超音波診断装置及びその制御プログラム
US11250564B2 (en) * 2019-12-19 2022-02-15 GE Precision Healthcare LLC Methods and systems for automatic measurement of strains and strain-ratio calculation for sonoelastography
WO2021216723A1 (fr) * 2020-04-22 2021-10-28 Bfly Operations, Inc. Procédés et appareils de formation de faisceau dans des systèmes à ultrasons
CN114376606B (zh) * 2022-01-18 2023-05-09 武汉联影医疗科技有限公司 一种超声成像的滤波方法和系统

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5623929A (en) * 1995-06-30 1997-04-29 Siemens Medical Systems, Inc. Ultrasonic doppler flow imaging method for eliminating motion artifacts
US5800356A (en) 1997-05-29 1998-09-01 Advanced Technology Laboratories, Inc. Ultrasonic diagnostic imaging system with doppler assisted tracking of tissue motion
WO1999017660A1 (fr) * 1997-10-07 1999-04-15 Ge Vingmed Ultrasound As Procede et appareil de calcul et d'affichage en temps reel de la deformation en imagerie echographique
US6099471A (en) 1997-10-07 2000-08-08 General Electric Company Method and apparatus for real-time calculation and display of strain in ultrasound imaging
US6309357B1 (en) * 1999-08-09 2001-10-30 Acuson Corporation Medical diagnostic ultrasound system and method for improved flow or movement detection with multiple clutter filters
US6558324B1 (en) 2000-11-22 2003-05-06 Siemens Medical Solutions, Inc., Usa System and method for strain image display
JP2006304965A (ja) * 2005-04-27 2006-11-09 Matsushita Electric Ind Co Ltd 超音波診断装置
WO2007030016A1 (fr) * 2005-09-08 2007-03-15 Angelsen Bjoern A J Imagerie acoustique par manipulation basse frequence non lineaire de proprietes de diffusion et de propagation haute frequence
US20090304246A1 (en) * 2008-05-15 2009-12-10 Walker William F Reduction of echo decorrelation facilitating motion estimation

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5623929A (en) * 1995-06-30 1997-04-29 Siemens Medical Systems, Inc. Ultrasonic doppler flow imaging method for eliminating motion artifacts
US5800356A (en) 1997-05-29 1998-09-01 Advanced Technology Laboratories, Inc. Ultrasonic diagnostic imaging system with doppler assisted tracking of tissue motion
WO1999017660A1 (fr) * 1997-10-07 1999-04-15 Ge Vingmed Ultrasound As Procede et appareil de calcul et d'affichage en temps reel de la deformation en imagerie echographique
US6099471A (en) 1997-10-07 2000-08-08 General Electric Company Method and apparatus for real-time calculation and display of strain in ultrasound imaging
US6309357B1 (en) * 1999-08-09 2001-10-30 Acuson Corporation Medical diagnostic ultrasound system and method for improved flow or movement detection with multiple clutter filters
US6558324B1 (en) 2000-11-22 2003-05-06 Siemens Medical Solutions, Inc., Usa System and method for strain image display
JP2006304965A (ja) * 2005-04-27 2006-11-09 Matsushita Electric Ind Co Ltd 超音波診断装置
WO2007030016A1 (fr) * 2005-09-08 2007-03-15 Angelsen Bjoern A J Imagerie acoustique par manipulation basse frequence non lineaire de proprietes de diffusion et de propagation haute frequence
US20090304246A1 (en) * 2008-05-15 2009-12-10 Walker William F Reduction of echo decorrelation facilitating motion estimation

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
LOUPAS T ET AL: "Experimental evaluation of velocity and power estimation for ultrasound blood flow imaging, by means of a two-dimensional autocorrelation approach", IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS AND FREQUENCY CONTROL, IEEE, US, vol. 42, no. 4, 1 July 1995 (1995-07-01), pages 689 - 699, XP011439272, ISSN: 0885-3010, DOI: 10.1109/58.393111 *
MAULDIN F W ET AL: "The Singular Value Filter: A General Filter Design Strategy for PCA-Based Signal Separation in Medical Ultrasound Imaging", IEEE TRANSACTIONS ON MEDICAL IMAGING, IEEE SERVICE CENTER, PISCATAWAY, NJ, US, vol. 30, no. 11, 1 November 2011 (2011-11-01), pages 1951 - 1964, XP011379749, ISSN: 0278-0062, DOI: 10.1109/TMI.2011.2160075 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20190108071A (ko) * 2018-03-13 2019-09-23 지멘스 메디컬 솔루션즈 유에스에이, 인크. 음향 방사력-기반 초음파 이미징에서의 적응식 클러터 필터링
CN110269643A (zh) * 2018-03-13 2019-09-24 美国西门子医疗解决公司 在基于声学辐射力的超声成像中的自适应杂波滤波
KR102245671B1 (ko) 2018-03-13 2021-04-27 지멘스 메디컬 솔루션즈 유에스에이, 인크. 음향 방사력-기반 초음파 이미징에서의 적응식 클러터 필터링
US11464495B2 (en) 2018-03-13 2022-10-11 Siemens Medical Solutions Usa, Inc. Adaptive clutter filtering in acoustic radiation force-based ultrasound imaging

Also Published As

Publication number Publication date
US20190029651A1 (en) 2019-01-31

Similar Documents

Publication Publication Date Title
US20190029651A1 (en) Clutter filters for strain and other ultrasonic deformation imaging
US8523776B2 (en) Ultrasonic doppler imaging apparatus and method with blood velocity waveform processing
JP7462816B2 (ja) ベクトルフローデータを使用する擾乱した血流の自動検出及び視覚化のためのシステム及び方法
KR102025328B1 (ko) 평면파 합성을 이용한 초음파 벡터 도플러 영상의 생성 장치 및 방법
JP2009530010A (ja) 心機能異常の分析のための心臓エコー検査装置及び方法
US9028413B2 (en) Prediction-based flow estimation for ultrasound diagnostic imaging
JP2006505321A (ja) 運動部分の３次元撮像用フェーズドアレイ音響システム
JP7232195B2 (ja) 血管内の壁せん断応力の同時視覚化及び定量化のためのシステム及び方法
CN106529561B (zh) 超声彩色流中的闪光伪像检测
JP4763588B2 (ja) 超音波診断装置
KR20120044267A (ko) 초음파 진단 장치 및 조직 움직임 추적 방법
JP2007518512A (ja) 心筋灌流を表示するための画像分割
US20190365344A1 (en) System and method for characterizing liver perfusion of contrast agent flow
EP3108456B1 (fr) Visualisation adaptée au mouvement dans l'imagerie médicale 4d
EP3600058B1 (fr) Système et procédé de visualisation et de quantification simultanées d'écoulement sanguin à l'aide d'une imagerie d'écoulement de vecteur ultrasonore
US10624608B2 (en) Ultrasonic diagnostic apparatus
JP2007222533A (ja) 超音波診断装置及び超音波画像処理方法
JP2003061958A (ja) 超音波診断装置
US7261695B2 (en) Trigger extraction from ultrasound doppler signals
WO2015104607A1 (fr) Modes d'imagerie ultrasonore pour quantification et analyse en temps réel automatisées
JP7366221B2 (ja) 解剖学的、機能的及び血行動態イメージングのためのトリプルモード超音波イメージング
US20100234730A1 (en) Image processing device, ultrasonic imaging apparatus including the same, and image processing method
US11497473B2 (en) Ultrasound cardiac processing
CN107920802B (zh) 利用自适应采样窗口尺寸的谱多普勒处理
JP2020014595A (ja) 超音波診断装置

Legal Events

Date Code Title Description
NENP Non-entry into the national phase

Ref country code: DE

121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 17707313

Country of ref document: EP

Kind code of ref document: A1

122 Ep: pct application non-entry in european phase

Ref document number: 17707313

Country of ref document: EP

Kind code of ref document: A1