WO1995020357A1 - Method for reducing speckle in vessel imaging - Google Patents
Method for reducing speckle in vessel imaging Download PDFInfo
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 21
- 238000003384 imaging method Methods 0.000 title abstract description 8
- 210000004369 blood Anatomy 0.000 claims abstract description 31
- 239000008280 blood Substances 0.000 claims abstract description 31
- 238000002592 echocardiography Methods 0.000 claims abstract description 12
- 210000004204 blood vessel Anatomy 0.000 claims abstract description 9
- 238000005070 sampling Methods 0.000 claims description 7
- 230000002708 enhancing effect Effects 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims 3
- 238000010408 sweeping Methods 0.000 claims 1
- 230000001131 transforming effect Effects 0.000 claims 1
- 238000002604 ultrasonography Methods 0.000 abstract description 13
- 238000012545 processing Methods 0.000 abstract description 12
- 238000012285 ultrasound imaging Methods 0.000 abstract description 5
- 239000000523 sample Substances 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 3
- 238000003708 edge detection Methods 0.000 description 3
- 239000013598 vector Substances 0.000 description 3
- 210000000601 blood cell Anatomy 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000009877 rendering Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 244000208734 Pisonia aculeata Species 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000007429 general method Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/06—Measuring blood flow
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/12—Diagnosis using ultrasonic, sonic or infrasonic waves in body cavities or body tracts, e.g. by using catheters
-
- 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
-
- Y—GENERAL 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
- Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
- Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10S128/00—Surgery
- Y10S128/916—Ultrasound 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
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.
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)
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 |
EP1105044A1 (en) | 1998-07-21 | 2001-06-13 | 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 |
US7359554B2 (en) * | 2002-08-26 | 2008-04-15 | Cleveland Clinic Foundation | System and method for identifying a vascular border |
US7074188B2 (en) * | 2002-08-26 | 2006-07-11 | The Cleveland Clinic Foundation | System and method of characterizing vascular tissue |
US7927275B2 (en) * | 2002-08-26 | 2011-04-19 | The Cleveland Clinic Foundation | System and method of aquiring blood-vessel data |
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 |
CN102512200A (en) * | 2003-11-26 | 2012-06-27 | 艾玛克公司 | Transesophageal ultrasound using a narrow 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 |
CA2763794C (en) | 2009-06-04 | 2017-01-17 | Rotation Medical, Inc. | Apparatus having bowstring-like staple delivery to a target tissue |
US8483488B2 (en) * | 2009-08-07 | 2013-07-09 | Medinol Ltd. | Method and system for stabilizing a series of intravascular ultrasound images and extracting vessel lumen from the images |
US10123796B2 (en) | 2014-11-04 | 2018-11-13 | Rotation Medical, Inc. | Medical implant delivery system and related methods |
AU2015343273B2 (en) | 2014-11-04 | 2017-12-14 | Rotation Medical, Inc. | Medical implant delivery system and related methods |
EP3215025B1 (en) | 2014-11-04 | 2020-12-23 | Rotation Medical, Inc. | Medical implant delivery system |
WO2017117437A1 (en) | 2015-12-31 | 2017-07-06 | Rotation Medical, Inc. | Fastener delivery system and related methods |
CN109561877A (en) | 2016-08-03 | 2019-04-02 | 皮埃-哈韦斯特控股有限公司 | A kind of includes the system and method for the internal pressure of intravascular blood pressure for non-invasive measurement |
EP3278735A1 (en) | 2016-08-03 | 2018-02-07 | PI-Harvest Holding AG | A system and method for non-invasive measurement of pressure inside a body including intravascular blood pressure |
CN110349256B (en) * | 2019-07-16 | 2023-05-23 | 深圳大学 | Vascular reconstruction method and device and computer terminal |
Citations (3)
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)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5000185A (en) * | 1986-02-28 | 1991-03-19 | Cardiovascular Imaging Systems, Inc. | Method for intravascular two-dimensional ultrasonography and recanalization |
US4794931A (en) * | 1986-02-28 | 1989-01-03 | Cardiovascular Imaging Systems, Inc. | Catheter apparatus, system and method for intravascular two-dimensional ultrasonography |
US5285788A (en) * | 1992-10-16 | 1994-02-15 | Acuson Corporation | Ultrasonic tissue imaging method and apparatus with doppler velocity and acceleration processing |
-
1994
- 1994-01-26 US US08/187,584 patent/US5363850A/en not_active Expired - Lifetime
- 1994-11-15 US US08/340,178 patent/US5520185A/en not_active Expired - Lifetime
- 1994-12-22 WO PCT/US1994/014917 patent/WO1995020357A1/en active Application Filing
Patent Citations (3)
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 | |
CN110325119B (en) | Ovarian follicle count and size determination | |
US5260871A (en) | Method and apparatus for diagnosis of breast tumors | |
US11064972B2 (en) | Systems and methods for detecting and displaying body lumen bifurcations | |
EP1101128B1 (en) | Method and apparatus for spatial and temporal filtering of intravascular ultrasonic image data | |
US6385332B1 (en) | Automated segmentation method for 3-dimensional ultrasound | |
US8202220B2 (en) | Ultrasound diagnostic apparatus | |
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 | |
CN111462077A (en) | Method for characterizing biological tissues by utilizing nonlinear information entropy | |
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 | |
Zatari et al. | In vivo liver differentiation by ultrasound using an artificial neural network |
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 |