WO1995020357A1 - Method for reducing speckle in vessel imaging - Google Patents

Method for reducing speckle in vessel imaging Download PDF

Info

Publication number
WO1995020357A1
WO1995020357A1 PCT/US1994/014917 US9414917W WO9520357A1 WO 1995020357 A1 WO1995020357 A1 WO 1995020357A1 US 9414917 W US9414917 W US 9414917W WO 9520357 A1 WO9520357 A1 WO 9520357A1
Authority
WO
WIPO (PCT)
Prior art keywords
image
pixels
pixel
blood
region
Prior art date
Application number
PCT/US1994/014917
Other languages
French (fr)
Inventor
Bobby Soni
Richard Didday
Kean Hurley
Original Assignee
Cardiovascular Imaging Systems, Inc.
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 Cardiovascular Imaging Systems, Inc. filed Critical Cardiovascular Imaging Systems, Inc.
Publication of WO1995020357A1 publication Critical patent/WO1995020357A1/en

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/06Measuring blood flow
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/12Diagnosis using ultrasonic, sonic or infrasonic waves in body cavities or body tracts, e.g. by using catheters
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52023Details of receivers
    • G01S7/52036Details of receivers using analysis of echo signal for target characterisation
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10STECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10S128/00Surgery
    • Y10S128/916Ultrasound 3-D imaging

Definitions

  • This invention relates to high resolution intravascular imaging and more particularly to ultrasound imaging and techniques for enhancing image quality.
  • echoes from blood molecules degrade the lumen-to-vessel wall contrast.
  • the production of high resolution images of vessel wall structures requires imaging at high ultrasound frequencies.
  • the back scatter from blood in such an image is a significant problem in high frequency ultrasound imaging, since the scattering of ultrasound from blood is proportional to the fourth power of the frequency such that the higher the ultrasound frequency the more pronounced is the back scatter from blood. It is desirable to provide imaging over a broad range of frequencies. Therefore, echoes in the ultrasound image due to blood cells must either be removed or suppressed to a level at which wall structures can be distinguished from blood. Since blood is typically in motion relative to the image of interest, an individual frame will contain speckle due to the interference of blood constituents.
  • One technique for enhancing image quality is to average successive image frames thereby to smooth out the impact of speckle (that is, the irregular pattern of backscatter from blood cells) in the ultrasound images.
  • Thresholding-based techniques have been suggested wherein blood-speckle induced smearing is removed based on a signal intensity compared with a threshold.
  • Such techniques are not effective in reducing the mean echo amplitude from the region of blood flow, and they cannot totally remove blood echoes from the image.
  • What is needed is a mechanism for identifying the signatures of blood- induced echoes and for subsequently suppression which allows a display of ultrasound images free of blood-induced echoes. Background information on the subject of intravascular ultrasonography is found in U.S. Patent No. 4,794,931 and related U.S. Patent No. 5,000,185.
  • ultrasound echoes representing blood vessel walls are distinguished from ultrasound echoes from blood.
  • a classifier which employs the mean and variance of raw data of gray intensities as acquired directly from an ultrasound scanner detector.
  • the statistical parameters of mean and variance are determined in the course of real-time processing and prior to attempted classification and then are used as parameters representing the true image as acquired.
  • Parameters are computed by first obtaining an image, then distinguishing between the blood region as opposed to the wall region by statistical analysis, namely, obtaining a mean and variance at candidate positions, multiplying the mean by the variance to obtain a "factor" and noting that the product of mean and variance of a blood region is substantially higher than the product of the mean and variance of a wall region.
  • a user can control the appearnce of the blood region of the image.
  • Figure l is a block diagram of an intravascular ultrasonic imaging system in accordance with the invention.
  • Figure 2 is a representative diagram of an array of a full density pixel image.
  • Figure 3 is a representative diagram of an array of a working image.
  • Figure 4 is a flow chart of the method according to the invention.
  • Figure 5 is a flow chart of a general method of processing in accordance with the invention.
  • a specialized signal processing device 10 is used with an ultrasonic imaging system 12 including a catheter probe 13 wherein ultrasonic beams 14 are emitted by an exciter 16, and spokes or vectors 18 of information are collected from a target 20 (the interior walls of a blood vessel) based on ultrasonic reflections at a transducer 22.
  • the reflections scale in amplitude over a range, and pixels representing points in a scanned (swept and pulled out) two- dimensional image are assigned a value on a gray scale between opaque to transparent.
  • the image is representative of a
  • slice of the structure of blood vessel 20 and includes wall structures and lumens of blood.
  • the information gathered from the sampling spoke 18 is used to determine a mean value and a variance value useful for calculation purposes as hereinafter explained.
  • a total of eight spokes radially directed from the catheter 13 is generally sufficient to obtain data to process the information.
  • the average mean and avera ⁇ e standard deviation (variance) are determined.
  • the raw image data is sampled for every point using a local neighborhood, where the size of the window is relatively dense, e.g., eleven by eleven pixels (Step A) .
  • a distance metric is evaluated for a two-class problem, the two classes being blood region and wall region (Step B) .
  • the distance metrics are "Fisher" distances.
  • Other distance metrics may be employed, such as a least squares distance metric.
  • the "Fisher” distance metric is preferred because it takes into account statistical variance.
  • the "Fisher" distances are represented by the following equations:
  • X L is the gray scale intensity of the (i,j)th pixel; n- j is the mean of class 1; var x is the variance of class 1; m 2 is the mean of class 2; and var 2 is the variance of class 2.
  • the decision rule for determining whether a pixel with Fisher distance d x belongs to class 1 or to class 2 is given by: If (d ⁇ is less than d 2 ) , then the (i,j)th belongs to class 1; otherwise the (i,j)th pixel belongs to class 2.
  • Steps A and B The entire image is thus processed with overlapping windows, segmenting each window into the two classes.
  • a local binary copy of the image is maintained in memory of the signal processing device 10 with each pixel belonging to either to the wall region or the blood region.
  • a pullback (longitudinal displacement of the tranducer along the blood vessel) is performed to produce a series of sliced images. Then the Fisher distances are computed for each slice, and the pixels of each image slice are classified so that a set of binary images is obtained.
  • the initially processed binary image (of one slice) 28 comprises an an array of pixels measuring 220 by 220 pixels (Fig. 2—not all pixels shown) .
  • the initial binary image 28 (Fig. 2) , which is typically too cumbersome to be subjected to a long and complex calculation) is sampled or abstracted at a low density offor example one pixel per eleven pixels of the initial image 28 to yield an abstracted image or working image 32 (Fig.
  • Step E 3) of for example twenty by twenty pixels, each pixel in the abstracted image 32 representing an eleven by eleven group of the full image 28 (Step E) .
  • This working image 32 is also stored in the image processing device 10 for manipulation. Then, based on knowledge of morphology, e.g., assuming that the lumen under examination is enclosed by the wall tissue 28 in the probed vessel structure or target 20, the working image 32 is modified, namely, a region therein is reconstructively "grown" from a seed point 34, if a growth path 36 exists of connected pixels 38 of the same class (Step F) .
  • a seed point 40 for the lumen a pixel in the center of the transducer 22 (Fig.
  • Step G The inventive process is employed to develop a viewable image of the other types of morphological elements (Step H, F, G) .
  • Wall seed point 34 is selected so as to "grow” to include all connected pixels 38 of the abstracted image which belong to the "wall" class.
  • intermediate reconstruction 44 then consists of interconnected interior pixels 42, interconnected exterior pixels 38 and untouched unclassified pixels 46.
  • the untouched pixels thereafter need to be classified. They are classified by checking the classification of their neighbors and then assigning to them the classification of their immediate neighbors (Step I) .
  • edge pixels 48 comprising the interface between the lumen and the wall must be processed.
  • the resolution of the abstracted image is generally insufficient to classify such points correctly.
  • the edge pixels must therefore be at least initially tagged as being edge pixels, that is, as being neither interior nor exterior
  • An edge detection processing technique is applied to the image to look for transitions from interior points to exterior points. Pixels found between interior points and exterior points are first tagged as edge pixels (Step J) .
  • the edge detection processing technique is first applied from the blood region looking out (Step K) .
  • the edge detection processing technique is applied from the wall region looking in (Step L) .
  • the classified abstracted image or reconstruction 44 is thereafter remapped into a full resolution gray scale image 50 (Step M) , with the following rules applied to assignment of gray scale values:
  • the gray scale value is set to zero or some other high contrast value which is a percentage of the original gray scale intensity.
  • the Fisher classification for the high resolution pixels in the full resolution image is used to identify the pixels as either in the blood region or the wall region, and the corresponding gray scale value is assigned to the resultant high resolution pixel.
  • Fig. 5 is a flow chart of program flow according to the invention.
  • Step AA the system is initialized (Step AB) , and it acquires data (Step AC) , or alternatively is calibrated (Step AD) or is terminated (Step AE) .
  • Step AF it either iteratively applies the algorithm (Steps AG, AH) and quits, or applies it and tests for the requirement of a rendering process (Step Al) .
  • the system either skips to the display step (Step AJ) or performs a rendering and a 3D animation (Steps AK, AL) , after which it displays the results (Step AJ) .
  • Steps AK, AL rendering and a 3D animation

Abstract

This invention is a method for using an ultrasound imaging system comrpising: a signal processing device (10) and an ultrasonic imaging system (12) including a catheter probe (13) wherein ultrasonic beams (14) are emitted by an exciter (16), and spokes (18) of information are collected from a target (20), the interior walls of a blood vessel based on ultrasonic reflections at a transducer (22). Ultrasound echoes representing vessel walls are distinguished from ultrasound echoes from blood by using a classifier which employs the mean and variance of the raw data of gray intensities as acquired directly from an ultrasound scanner detector. In a specific embodiment, the statistical parameters of mean and variance are determined in the course of real time processing and prior to attempted classification and then are used as parameters representing the true image as acquired.

Description

METHOD FOR REDUCING SPECKLE IN VESSEL IMAGING
BACKGROUND OF THE INVENTION This invention relates to high resolution intravascular imaging and more particularly to ultrasound imaging and techniques for enhancing image quality.
In ultrasound imaging, echoes from blood molecules degrade the lumen-to-vessel wall contrast. However, the production of high resolution images of vessel wall structures requires imaging at high ultrasound frequencies. The back scatter from blood in such an image is a significant problem in high frequency ultrasound imaging, since the scattering of ultrasound from blood is proportional to the fourth power of the frequency such that the higher the ultrasound frequency the more pronounced is the back scatter from blood. It is desirable to provide imaging over a broad range of frequencies. Therefore, echoes in the ultrasound image due to blood cells must either be removed or suppressed to a level at which wall structures can be distinguished from blood. Since blood is typically in motion relative to the image of interest, an individual frame will contain speckle due to the interference of blood constituents. One technique for enhancing image quality is to average successive image frames thereby to smooth out the impact of speckle (that is, the irregular pattern of backscatter from blood cells) in the ultrasound images. Thresholding-based techniques have been suggested wherein blood-speckle induced smearing is removed based on a signal intensity compared with a threshold. However, such techniques are not effective in reducing the mean echo amplitude from the region of blood flow, and they cannot totally remove blood echoes from the image. What is needed is a mechanism for identifying the signatures of blood- induced echoes and for subsequently suppression which allows a display of ultrasound images free of blood-induced echoes. Background information on the subject of intravascular ultrasonography is found in U.S. Patent No. 4,794,931 and related U.S. Patent No. 5,000,185.
SUMMARY OF THE INVENTION
According to the invention, in an ultrasound imaging system, ultrasound echoes representing blood vessel walls are distinguished from ultrasound echoes from blood. This is achieved by using a classifier which employs the mean and variance of raw data of gray intensities as acquired directly from an ultrasound scanner detector. In a specific embodiment, the statistical parameters of mean and variance are determined in the course of real-time processing and prior to attempted classification and then are used as parameters representing the true image as acquired.
Parameters are computed by first obtaining an image, then distinguishing between the blood region as opposed to the wall region by statistical analysis, namely, obtaining a mean and variance at candidate positions, multiplying the mean by the variance to obtain a "factor" and noting that the product of mean and variance of a blood region is substantially higher than the product of the mean and variance of a wall region. As a result, a user can control the appearnce of the blood region of the image. The invention will be better understood by reference to the following detailed description in connection with the accompanying drawings.
BRIEF DESCRIPTION OF THE DRAWINGS Figure l is a block diagram of an intravascular ultrasonic imaging system in accordance with the invention.
Figure 2 is a representative diagram of an array of a full density pixel image.
Figure 3 is a representative diagram of an array of a working image.
Figure 4 is a flow chart of the method according to the invention. Figure 5 is a flow chart of a general method of processing in accordance with the invention.
DESCRIPTION OF SPECIFIC EMBODIMENTS The invention may be understood by reference to the
Figures 1 through 5. A specialized signal processing device 10 is used with an ultrasonic imaging system 12 including a catheter probe 13 wherein ultrasonic beams 14 are emitted by an exciter 16, and spokes or vectors 18 of information are collected from a target 20 (the interior walls of a blood vessel) based on ultrasonic reflections at a transducer 22. The reflections scale in amplitude over a range, and pixels representing points in a scanned (swept and pulled out) two- dimensional image are assigned a value on a gray scale between opaque to transparent. The image is representative of a
"slice" of the structure of blood vessel 20 and includes wall structures and lumens of blood.
More specifically, information is gathered by projecting narrow ultrasonic sampling beams 14 from exciter 16 as it is rotated within catheter 13 within blood vessel 20, the reflections of which are recorded by transducer 22 as amplitude as a function of distance along the radius of each vector. "Factors" are computed at every unit radius in a sampling window (e.g., five by five pixels in size). Depending upon the absolute value of the factors along the entire radius of each reflected spoke (vector) 18, the end of the spoke 18 is identified as lying entirely within the blood region 24 or at the blood-wall interface 26 or within the vessel wall 28 or at some indeterminant position 30 beyond the vessel wall 28. So long as a reasonable length of a sampling spoke is determined to be within a known region, such as the blood region, the information gathered from the sampling spoke 18 is used to determine a mean value and a variance value useful for calculation purposes as hereinafter explained. A total of eight spokes radially directed from the catheter 13 is generally sufficient to obtain data to process the information. Using the accumulated information from as many as all eight spokes, the average mean and averaσe standard deviation (variance) are determined.
Referring to Fig. 4, as a first estimation of the (location of the) blood region and the wall region in a volume image, the raw image data is sampled for every point using a local neighborhood, where the size of the window is relatively dense, e.g., eleven by eleven pixels (Step A) . Then for each point a distance metric is evaluated for a two-class problem, the two classes being blood region and wall region (Step B) . In a specific embodiment, the distance metrics are "Fisher" distances. Other distance metrics may be employed, such as a least squares distance metric. The "Fisher" distance metric is preferred because it takes into account statistical variance. The "Fisher" distances are represented by the following equations:
d l = (χij " m1)2/var1 (Eq. 1) d2 = (xi;j - m2)2/var2 (Eq. 2)
where;
XL is the gray scale intensity of the (i,j)th pixel; n-j is the mean of class 1; varx is the variance of class 1; m2 is the mean of class 2; and var2 is the variance of class 2.
The decision rule for determining whether a pixel with Fisher distance dx belongs to class 1 or to class 2 is given by: If (dλ is less than d2) , then the (i,j)th belongs to class 1; otherwise the (i,j)th pixel belongs to class 2.
These Fisher distances are computed for every point within the sampling window. Votes are accumulated for the window for both the blood region class and the wall region class, depending whether the pixel is classified as belonging to "blood" or "wall." If the percentage of pixels classified as blood within the sampling window is more than the empirically-estimated threshold, then the center pixel of the window is classified as belonging to the blood region. Otherwise the center pixel is classified as belonging to the wall region. Thereafter the window is moved (by one pixel) (Steps C, D) , and the center pixel is again classified
(repeating Steps A and B) . The entire image is thus processed with overlapping windows, segmenting each window into the two classes.
A local binary copy of the image is maintained in memory of the signal processing device 10 with each pixel belonging to either to the wall region or the blood region. A pullback (longitudinal displacement of the tranducer along the blood vessel) is performed to produce a series of sliced images. Then the Fisher distances are computed for each slice, and the pixels of each image slice are classified so that a set of binary images is obtained.
After the first estimation is made, knowledge of blood vessel morphology is used to refine the initial binary classification of each pixel in each sliced image. (Up to 150 slices are processed, preferably within sixty seconds time, to obtain information which is sufficiently time invariant to be clinically significant.) The initially processed binary image (of one slice) 28 comprises an an array of pixels measuring 220 by 220 pixels (Fig. 2—not all pixels shown) . According to the invention, the initial binary image 28 (Fig. 2) , which is typically too cumbersome to be subjected to a long and complex calculation) is sampled or abstracted at a low density offor example one pixel per eleven pixels of the initial image 28 to yield an abstracted image or working image 32 (Fig. 3) of for example twenty by twenty pixels, each pixel in the abstracted image 32 representing an eleven by eleven group of the full image 28 (Step E) . This working image 32 is also stored in the image processing device 10 for manipulation. Then, based on knowledge of morphology, e.g., assuming that the lumen under examination is enclosed by the wall tissue 28 in the probed vessel structure or target 20, the working image 32 is modified, namely, a region therein is reconstructively "grown" from a seed point 34, if a growth path 36 exists of connected pixels 38 of the same class (Step F) . As a seed point 40 for the lumen, a pixel in the center of the transducer 22 (Fig. l) may be selected, thus assuring that the seed point is in the blood region 24 in the image 32. Pixels considered to be in the blood region are tagged as being "interior" pixels 42. These pixels are then assigned a distinct gray scale value in a reconstruction 44 of abstracted image 32 stored also in the signal processing device 10 (Step G). The inventive process is employed to develop a viewable image of the other types of morphological elements (Step H, F, G) . Wall seed point 34 is selected so as to "grow" to include all connected pixels 38 of the abstracted image which belong to the "wall" class. The growing procedure yields the connections, and the wall pixels 38 in the abstracted or working image 32 are so tagged as being "exterior" and assigned a distinct gray scale value different from the lumen-assigned gray scale value in the reconstruction 44 (Fig. 1) . An intermediately-processed abstracted image
(intermediate reconstruction 44) then consists of interconnected interior pixels 42, interconnected exterior pixels 38 and untouched unclassified pixels 46. The untouched pixels thereafter need to be classified. They are classified by checking the classification of their neighbors and then assigning to them the classification of their immediate neighbors (Step I) .
As a next step, edge pixels 48 comprising the interface between the lumen and the wall must be processed. The resolution of the abstracted image is generally insufficient to classify such points correctly. The edge pixels must therefore be at least initially tagged as being edge pixels, that is, as being neither interior nor exterior An edge detection processing technique is applied to the image to look for transitions from interior points to exterior points. Pixels found between interior points and exterior points are first tagged as edge pixels (Step J) . The edge detection processing technique is first applied from the blood region looking out (Step K) . Then the edge detection processing technique is applied from the wall region looking in (Step L) .
The classified abstracted image or reconstruction 44 is thereafter remapped into a full resolution gray scale image 50 (Step M) , with the following rules applied to assignment of gray scale values:
1) If the pixel in the reconstructed image is tagged as being an exterior point, then the corresponding (surrounding) pixels of the full resolution image are assigned the original gray scale value of the raw image.
2) If the pixel in the reconstructed image is tagged as being interior, then the corresponding (surrounding) pixels of the full resolution image are presumed to be blood and the gray scale value is set to zero or some other high contrast value which is a percentage of the original gray scale intensity.
3) If the pixel in the reconstructed image is tagged as an edge pixel, then the Fisher classification for the high resolution pixels in the full resolution image is used to identify the pixels as either in the blood region or the wall region, and the corresponding gray scale value is assigned to the resultant high resolution pixel.
In this manner, a high resolution image is obtained for display on a display device 52 with sharp edge identification, while minimizing the amount of required computation of Fisher distances.
Fig. 5 is a flow chart of program flow according to the invention. After start or following a "quit" state (Step AA) , the system is initialized (Step AB) , and it acquires data (Step AC) , or alternatively is calibrated (Step AD) or is terminated (Step AE) . After testing for the need for the Fisher Algorithm (Step AF) , it either iteratively applies the algorithm (Steps AG, AH) and quits, or applies it and tests for the requirement of a rendering process (Step Al) . The system either skips to the display step (Step AJ) or performs a rendering and a 3D animation (Steps AK, AL) , after which it displays the results (Step AJ) . There are several possible switches or options to the processing: To step through raw data or to display a three-dimensional movie; to save raw data or to display the three-dimensional movie; or to animate the raw data or display a three-movie.
While the invention has been particularly shown and described with reference to preferred embodiments thereof, it will be understood by those skilled in the art that the foregoing and other changes in the form and details may be made therein without departing from the spirit or scope of the invention. For example, variants on the distance metric may be employed. It is therefore not intended that this invention be limited, except as indicated by the appended claims.

Claims

WHAT IS CLAIMED IS:
1. A method for enhancing an intravascular ultrasonic blood vessel image comprising the steps of: a) illuminating an intravascular target from within a blood vessel with ultrasonic energy to generate ultrasonic echoes from said intravascular target; thereafter b) intercepting the ultrasonic echoes from the intravascular target; c) transforming the ultrasonic echoes into electrical signals; d) assigning a gray scale value to each increment (pixel) in an image based on amplitude of each said electrical signal, said image being constructed by illuminating and sweeping said intravascular target and circumferentially and radially sampling the electrical signal; e) computing a distance metric for each pixel using seed values for each said pixel; f) abstracting said image to obtain an abstracted image; g) growing a first interior classification for points in said abstracted image from a first seed pixel preselected to be entirely within a first region of a first known class; h) growing a second exterior classification for exterior points from a second seed pixel in a second region of a second known class; i) identifying and classifying other pixels in the abstracted image by searching for points between interior classified pixels and exterior classified pixels and tagging said other pixels as said edge pixels; j) remapping all interior pixels, exterior pixels and edge pixels so classified into a full resolution image 1) by mapping exterior pixels with first gray scale values, 2) by mapping interior pixels with second gray scale values and 3) by mapping edge pixels according to a high resolution Fisher classification for said edge pixels; and k) outputting a resultant reconstruction of said image for display.
2. The method according to claim 1 wherein said step a) comprises illuminating the intravascular target with a plurality of beams in radially-directed as spokes around a central point, and incrementally advancing the central point along an intravascular axis so as to generate information to construct a three-dimensional image.
3. The method according to claim 2 wherein said abstracting step comprises collapsing said image by a factor which yields said abstracted image with minimum detail sufficient to identify and classify between a blood region and a vessel wall region.
4. The method according to claim 1 wherein said abstracting step comprises collapsing said image by a factor which yields said abstracted image with minimum detail sufficient to identify and classify between a blood region and a vessel wall region.
5. The method according to claim 4 wherein said distance metric computing step comrises computing a Fisher distance for each pixel using said seed values for mean and variance to obtain a Fisher classification.
6. The method according to claim 1 wherein said distance metric computing step comrises computing a Fisher distance for each pixel using said seed values for mean and variance to obtain a Fisher classification.
PCT/US1994/014917 1994-01-26 1994-12-22 Method for reducing speckle in vessel imaging WO1995020357A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US08/187,584 1994-01-26
US08/187,584 US5363850A (en) 1994-01-26 1994-01-26 Method for recognition and reduction of blood speckle in blood vessel imaging system

Publications (1)

Publication Number Publication Date
WO1995020357A1 true WO1995020357A1 (en) 1995-08-03

Family

ID=22689574

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US1994/014917 WO1995020357A1 (en) 1994-01-26 1994-12-22 Method for reducing speckle in vessel imaging

Country Status (2)

Country Link
US (2) US5363850A (en)
WO (1) WO1995020357A1 (en)

Families Citing this family (40)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5453575A (en) * 1993-02-01 1995-09-26 Endosonics Corporation Apparatus and method for detecting blood flow in intravascular ultrasonic imaging
US5446800A (en) * 1994-06-13 1995-08-29 Diasonics Ultrasound, Inc. Method and apparatus for displaying angiographic data in a topographic format
JP3406106B2 (en) * 1995-02-06 2003-05-12 ジーイー横河メディカルシステム株式会社 Ultrasonic image display method and ultrasonic diagnostic apparatus
US5592941A (en) * 1995-10-20 1997-01-14 Diagnostic Ultrasound Corporation Apparatus and method for non-invasively detecting urinary tract infection using ultrasound
EP0839497A1 (en) * 1996-11-01 1998-05-06 EndoSonics Corporation A method for measuring volumetric fluid flow and its velocity profile in a lumen or other body cavity
US7789841B2 (en) 1997-02-06 2010-09-07 Exogen, Inc. Method and apparatus for connective tissue treatment
US6095976A (en) * 1997-06-19 2000-08-01 Medinol Ltd. Method for enhancing an image derived from reflected ultrasound signals produced by an ultrasound transmitter and detector inserted in a bodily lumen
US5957138A (en) * 1997-08-25 1999-09-28 Diasonics Ultrasound, Inc. Method and apparatus for three-dimensional flow lumen imaging
US5876343A (en) * 1997-09-23 1999-03-02 Scimed Life Systems, Inc. Methods and apparatus for blood speckle detection in an intravascular ultrasound imaging system
US5921934A (en) * 1997-11-25 1999-07-13 Scimed Life Systems, Inc. Methods and apparatus for non-uniform rotation distortion detection in an intravascular ultrasound imaging system
JP2889568B1 (en) * 1998-05-18 1999-05-10 正男 伊藤 Vascular thickness measurement device and arteriosclerosis diagnosis device
CA2338735A1 (en) 1998-07-21 2000-02-03 Acoustic Sciences Associates Synthetic structural imaging and volume estimation of biological tissue organs
US6181810B1 (en) * 1998-07-30 2001-01-30 Scimed Life Systems, Inc. Method and apparatus for spatial and temporal filtering of intravascular ultrasonic image data
DE19846241C2 (en) * 1998-10-08 2000-12-21 Inst Angewandte Informatik Im Method and device for moisture measurement on masonry materials with ultrasonic signal analysis
US6554774B1 (en) * 2000-03-23 2003-04-29 Tensys Medical, Inc. Method and apparatus for assessing hemodynamic properties within the circulatory system of a living subject
US6454715B2 (en) 2000-04-11 2002-09-24 Scimed Life Systems, Inc. Methods and apparatus for blood speckle detection in an intravascular ultrasound imaging system
NL1019612C2 (en) * 2001-12-19 2003-06-20 Gemeente Amsterdam Steam superheater.
US6817982B2 (en) * 2002-04-19 2004-11-16 Sonosite, Inc. Method, apparatus, and product for accurately determining the intima-media thickness of a blood vessel
US7927275B2 (en) * 2002-08-26 2011-04-19 The Cleveland Clinic Foundation System and method of aquiring blood-vessel data
US7074188B2 (en) * 2002-08-26 2006-07-11 The Cleveland Clinic Foundation System and method of characterizing vascular tissue
US7359554B2 (en) * 2002-08-26 2008-04-15 Cleveland Clinic Foundation System and method for identifying a vascular border
US6835177B2 (en) * 2002-11-06 2004-12-28 Sonosite, Inc. Ultrasonic blood vessel measurement apparatus and method
GB2396406A (en) * 2002-12-17 2004-06-23 Qinetiq Ltd Image analysis
US7727153B2 (en) * 2003-04-07 2010-06-01 Sonosite, Inc. Ultrasonic blood vessel measurement apparatus and method
JP4913601B2 (en) * 2003-11-26 2012-04-11 イマコー・インコーポレーテッド Transesophageal ultrasound using a thin probe
US7874990B2 (en) * 2004-01-14 2011-01-25 The Cleveland Clinic Foundation System and method for determining a transfer function
US7215802B2 (en) * 2004-03-04 2007-05-08 The Cleveland Clinic Foundation System and method for vascular border detection
US20060030777A1 (en) * 2004-07-30 2006-02-09 Liang David H T-statistic method for suppressing artifacts in blood vessel ultrasonic imaging
US7460716B2 (en) * 2004-09-13 2008-12-02 Boston Scientific Scimed, Inc. Systems and methods for producing a dynamic classified image
US7680307B2 (en) * 2005-04-05 2010-03-16 Scimed Life Systems, Inc. Systems and methods for image segmentation with a multi-stage classifier
US7831081B2 (en) * 2005-08-15 2010-11-09 Boston Scientific Scimed, Inc. Border detection in medical image analysis
AU2010256415B2 (en) 2009-06-04 2015-04-02 Rotation Medical, Inc. Apparatus having bowstring-like staple delivery to a target tissue
WO2011015952A1 (en) * 2009-08-07 2011-02-10 Medinol Ltd. Method and system for stabilizing a series of intravascular ultrasound images and extracting vessel lumen from the images
EP3215025B1 (en) 2014-11-04 2020-12-23 Rotation Medical, Inc. Medical implant delivery system
AU2015343273B2 (en) 2014-11-04 2017-12-14 Rotation Medical, Inc. Medical implant delivery system and related methods
US10675019B2 (en) 2014-11-04 2020-06-09 Rotation Medical, Inc. Medical implant delivery system and related methods
EP3397175B1 (en) 2015-12-31 2021-11-24 Rotation Medical, Inc. Fastener delivery system
EP3478184A1 (en) 2016-08-03 2019-05-08 PI-Harvest Holding AG A system and method for non-invasive measurement of pressure inside a body including intravascular blood pressure
EP3524164A1 (en) 2016-08-03 2019-08-14 PI-Harvest Holding AG A system, method and software for non-invasive measurement of intravascular, in particular intracardiac blood pressure
CN110349256B (en) * 2019-07-16 2023-05-23 深圳大学 Vascular reconstruction method and device and computer terminal

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5040225A (en) * 1987-12-07 1991-08-13 Gdp, Inc. Image analysis method
US5224483A (en) * 1992-05-19 1993-07-06 Hewlett-Packard Company Adaptive contrast enhancement for scanned ultrasonic image
US5243988A (en) * 1991-03-13 1993-09-14 Scimed Life Systems, Inc. Intravascular imaging apparatus and methods for use and manufacture

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4794931A (en) * 1986-02-28 1989-01-03 Cardiovascular Imaging Systems, Inc. Catheter apparatus, system and method for intravascular two-dimensional ultrasonography
US5000185A (en) * 1986-02-28 1991-03-19 Cardiovascular Imaging Systems, Inc. Method for intravascular two-dimensional ultrasonography and recanalization
US5285788A (en) * 1992-10-16 1994-02-15 Acuson Corporation Ultrasonic tissue imaging method and apparatus with doppler velocity and acceleration processing

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5040225A (en) * 1987-12-07 1991-08-13 Gdp, Inc. Image analysis method
US5243988A (en) * 1991-03-13 1993-09-14 Scimed Life Systems, Inc. Intravascular imaging apparatus and methods for use and manufacture
US5224483A (en) * 1992-05-19 1993-07-06 Hewlett-Packard Company Adaptive contrast enhancement for scanned ultrasonic image

Also Published As

Publication number Publication date
US5363850A (en) 1994-11-15
US5520185A (en) 1996-05-28

Similar Documents

Publication Publication Date Title
US5363850A (en) Method for recognition and reduction of blood speckle in blood vessel imaging system
US5260871A (en) Method and apparatus for diagnosis of breast tumors
US11064972B2 (en) Systems and methods for detecting and displaying body lumen bifurcations
US6385332B1 (en) Automated segmentation method for 3-dimensional ultrasound
US8202220B2 (en) Ultrasound diagnostic apparatus
EP1101128B1 (en) Method and apparatus for spatial and temporal filtering of intravascular ultrasonic image data
US6454715B2 (en) Methods and apparatus for blood speckle detection in an intravascular ultrasound imaging system
US6309353B1 (en) Methods and apparatus for tumor diagnosis
Bleck et al. Random field models in the textural analysis of ultrasonic images of the liver
Zahalka et al. An automated segmentation method for three-dimensional carotid ultrasound images
US20110182495A1 (en) System and method for automatic defect recognition of an inspection image
JPH11313822A (en) Three-dimensional imaging system, and method and device for following movement of scanning plane
WO2011015952A1 (en) Method and system for stabilizing a series of intravascular ultrasound images and extracting vessel lumen from the images
US9039620B2 (en) Ultrasound diagnostic apparatus
US6744911B1 (en) Tomographic segmentation
JP2000126178A (en) Method of quantifying stereo surface shape and automatic identification method of malignant tumor
JP2000126181A (en) Method for extracting and treating tumor
CN113570594A (en) Method and device for monitoring target tissue in ultrasonic image and storage medium
JPH08293025A (en) Image sorting device
Schlaps et al. Ultrasonic tissue characterization using a diagnostic expert system
Choy et al. Morphological image analysis of left-ventricular endocardial borders in 2D echocardiograms
US6358206B1 (en) Ultrasound process for the determination of the location of a parietal surface in a tissue and of the absolute radius of an artery, and ultrasound apparatus for carrying out such process
CN111462077A (en) Method for characterizing biological tissues by utilizing nonlinear information entropy
Zatari et al. In vivo liver differentiation by ultrasound using an artificial neural network
JPH07108043A (en) Ultrasonic blood flow observing device

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): CA JP

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): AT BE CH DE DK ES FR GB GR IE IT LU MC NL PT SE

DFPE Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed before 20040101)
121 Ep: the epo has been informed by wipo that ep was designated in this application
122 Ep: pct application non-entry in european phase
NENP Non-entry into the national phase

Ref country code: CA