CA2239330A1 - Ultrasonic diagnostic imaging system with doppler assisted tracking of tissue motion - Google Patents
Ultrasonic diagnostic imaging system with doppler assisted tracking of tissue motion Download PDFInfo
- Publication number
- CA2239330A1 CA2239330A1 CA002239330A CA2239330A CA2239330A1 CA 2239330 A1 CA2239330 A1 CA 2239330A1 CA 002239330 A CA002239330 A CA 002239330A CA 2239330 A CA2239330 A CA 2239330A CA 2239330 A1 CA2239330 A1 CA 2239330A1
- Authority
- CA
- Canada
- Prior art keywords
- lines
- velocity
- ultrasonic
- doppler
- image
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Abandoned
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/48—Diagnostic techniques
- A61B8/486—Diagnostic techniques involving arbitrary m-mode
-
- 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
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8979—Combined Doppler and pulse-echo imaging systems
-
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/48—Other medical applications
- A61B5/4887—Locating particular structures in or on the body
- A61B5/489—Blood vessels
-
- 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/52073—Production of cursor lines, markers or indicia by electronic means
Abstract
A method and apparatus are presented for tracing the border of tissue through temporally acquired scanlines comprising the steps of reducing noise in the scanlines, producing a map of tissue edges from the scanlines, denoting a tissue border to be traced, and using velocity information corresponding to tissue edges to trace the denoted border.
Description
CA 02239330 1998-0~-29 ULTRASONIC DIAGNOSTIC IMAGING SYSTEM WITH
DOPPLER ASSISTED TRACKING OF TISSUE MOTION
This invention relates to ultrasonic diagnostic imaging systems and, in particular, to ultrasonic diagnostic imaging systems in which motion depicted in ultrasonic images is tracked through the assistance of ultrasonic velocity information.
One of the primary uses of medical diagnostic ultrasound is the diagnosis of cardiac function and disease. Echocardiography revolves around the functioning of the heart and the diagnosis of abnormal conditions of cardiac performance and pathology. Of significant interest to the cardiologist in making these diagnoses is clear delineation of the heart structure itself, and in particular the motion of the well-defined heart muscle as it cyclically pumps blood. But many patients are what are known as "difficult patients,"
meaning that the patient's physiology impedes the creation of sharply defined images from strong, clear echo returns. In those and other, more normal situations, the echocardiographer is limited to diagnosing from echocardiograms of poorly defined cardiac structures. To aid in diagnoses in these situations, ultrasound systems facilitate the tracing of the outlines of cardiac structures on the displayed echocardiogram by the echocardiographer, generally with a joystick or trackball. Quantified measures and analyses can then be performed based upon the dimensions and measurement of the outlined structures. But such tracing can be a painstaking task, since the heart is constantly contracting and relaxing and valves are continually opening and closing. For a complete and accurate diagnosis, it CA 02239330 1998-0~-29 may be necessary to define the structures of the heart in every echocardiogram throughout the cardiac cycle so that the heart muscle is continually defined while it is in motion. Realtime cardiac images are produced at a nominal rate of 30 images per second or more, presenting the echocardiographer with a daunting number of images with image structures to trace. It would be desirable, of course, for the ultrasound system to outline the structures of the heart automatically, and to maintain such outlines by tracking these structures automatically as they are engaged in the dynamics of contraction and relaxation, opening and closing.
In accordance with the principles of the present invention, an ultrasonic diagnostic imaging system is provided in which moving structures in a realtime ultrasonic image sequence are automatically tracked by use of velocity information. In a constructed embodiment the velocity information is derived from ultrasonic Doppler signals. Such automatic tracking enables the heart wall to be outlined or traced continually as it moves with the motion of the beating heart cycle. A preferred method for automatically tracing a moving structure in a realtime image sequence comprises the steps of detecting the edges of structures in the images, identifying an edge which is to be tracked or traced, and tracking or tracing the identified edge during its motional sequence with the use of velocity information concerning the identified edge.
In an illustrated embodiment, the borders of the heart wall are identified and traced in an M-mode sequence by use of the velocity characteristics of the beating heart. This embodiment leads to a number of useful diagnostic tools, such as analysis of the CA 02239330 1998-0~-29 mean velocity of the heart wall or measurement of the velocity gradient across the heart muscle.
In the drawings: .
FIGURE 1 is a block diagram of an ultrasonic diagnostic imaging system constructed in accordance with the principles of the present invention;
FIGURE 2 is a two dimensional image of the left ventricle of the heart with an M-line cursor across the heart chamber;
FIGURE 3 is a greyscale M-mode display acquired from the M-line location in FIGURE 2;
FIGURE 3a represents one M-line of the M-mode display of FIGURE 3;
FIGURE 4 is a block diagram illustrating a preferred edge detection technique for the ultrasound system of FIGURE l;
FIGURES 5a-5c are plots illustrating the operation of the preferred edge detection technique of FIGURE 4;
FIGURE 6 is a color M-mode display correlated to the greyscale M-mode display of FIGURE 3;
FIGURE 7 illustrates the tracking of moving tissue structure in accordance with the principles of the present invention;
FIGURE 8 illustrates the tracing of the heart wall in a greyscale M-mode display in accordance with the present invention;
FIGURE 9 is a graph of the mean heart wall velocity produced from the color M-mode data of FIGURE 6 in accordance with the principles of the present invention; and FIGURE 10 illustrates the application of the principles of the present invention to two dimensional greyscale images.
Referring first to FIGURE 1, an ultrasonic CA 02239330 1998-0~-29 diagnostic imaging system constructed in accordance with the principles of the present invention is shown in block diagram form. An ultrasonic array probe 10 is operated under control of a beamformer 12 to transmit ultrasonic waves into the body of a patient and receives echoes returning from the body in response to those waves. The beamformer 12 is controlled by a beamformer controller 14 which controls the time of activation of the probe transducer array elements to steer and focus a transmitted beam. The beamformer is also controlled to appropriately delay and sum received echo signals to form received beams of coherent echoes along each desired scanline.
The echoes from the received scanlines may then be processed in a variety of ways. One way is by means of a B mode processor 18, which detects the echo amplitudes. The echo amplitudes may then be displayed spatially as a greyscale or B mode image of the body. Another processing technique is by means of a Doppler processor, which processes ensembles of scanlines from a common location to detect the Doppler shift of moving blood and other tissues along the scanlines. The Doppler signals can be continually received from a sample volume in the body and displayed in a spectral Doppler display. A
second means of display is a colorflow display, in which the Doppler signals are spatially overlaid for display with a greyscale image, thereby depicting flow parameters within the structures of the greyscale image.
In accordance with one aspect of the present invention, the detected and processed echo signals are displayed in a time-motion or M-mode display format. In an M-mode display a one dimensional CA 02239330 1998-0~-29 "_ image, a single scanline termed an A-line, is repetitively scanned and displayed over time. The resultant display reveals changes occurring along the scanline location with time. M-mode displays can be formed using either Doppler or greyscale scanlines.
In the system of FIGURE 1, Doppler scanlines produced by the Doppler processor 16 are stored as color M-lines in storage device 20, and greyscale scanlines produced by the B mode processor 18 are stored as greyscale M-lines in the storage device 20. The M-mode display is most useful for displaying structures in the body which are moving, such as the beating heart, as this movement can be seen in the time-displaced sequence of scanlines of the display.
The color M-lines stored in the storage device 20 may be read from the device and converted into video signals by a video processor 22, then displayed as a color M-mode display on a display 26. The greyscale M-lines may also be applied to the video processor 22 for display as a greyscale M-mode display.
In accordance with the principles of the present invention, the edges or borders of structures in the greyscale scanlines are traced by a motion tracker 30. The motion tracker 30 uses the Doppler signals to track the movement of moving structures such as blood cells and other tissue in the greyscale scanlines. The motion tracker first removes extraneous artifacts and noise from the greyscale scanlines by means of a noise eliminator 32 and, optionally, a noise eliminator 34 eliminates noise from the Doppler scanlines. The greyscale scanlines then undergo edge detection by an edge detector 36, which identifies the locations of edges or boundaries of structures along the scanlines. Finally, the edges or boundaries depicted in a sequence of scanlines are traced with an edge tracer 38. The edge tracing is initiated by user initialization, whereby a user points to one point along the edge which is to be traced. The edge tracer then traces the edge of which the identified point is a part. It does this through use of the velocity data of the Doppler information. The Doppler information provides data as to the velocity at locations along the greyscale scanline, including the boundary point identified by the user. By knowing the time interval which lapsed between the acquisition of each greyscale scanline, and the velocity and direction of motion of an identified boundary point, the system can compute the location to which the boundary point has moved by the time of the subsequently acquired scanline. The system then looks for the boundary point in the vicinity of the computed position and produces a trace connecting the points of the two scanlines. In this way the motion tracker 30 can continuously track moving structures or structure boundaries from one scanline to the next, and can trace or outline them in the M-mode display.
The trace defined by the edge tracer 38 is applied to the video processor in concert with the greyscale M-lines, where the trace is overlaid over the M-mode display to define the traced edge or border or other characteristic in the display. The trace can also be overlaid over the temporally corresponding color M-mode display.
Once the border of a structure such as the heart wall has been defined by edge tracing, a number of quantified measurements can be made of the now precisely delineated structure. Such quantified measures are performed by a motion analysis module 24 CA 02239330 1998-0~-29 , ''-as described below, in which the computation of measures such as mean velocities or velocity gradients are discussed.
An understanding of the operation of the system of FIGURE 1 may be had with reference to the following drawings. FIGURE 2 illustrates a two dimensional ultrasound image of the heart. Typically a sector image of the heart is produced with a phased array probe, although curved and linear array probes may also be used. At the top of the image is a portion of the right ventricle, marked RV in the drawing. A septum 40 separates the right ventricle from the left ventricle LV. The chamber of the left ventricle is bounded by the posterior wall 42. At the right side of the left ventricle is the mitral valve 44.
Extending from the apex of the image and across the right and left ventricles is an M-mode cursor 46.
The M-mode cursor marks the scanline of the two dimensional (2D) image which is repetitively scanned and reproduced to form an M-mode image as shown in FIGURE 3. An M-mode image is formed by repetitively scanning the same spatial scanline and reproducing each new scan of the scanline alongside the previous scanlines. Thus, the M-mode image is a time sequence of one dimensional images of a single A-line. The benefit provided by such a display is that the image reveals the motion of tissues and tissue boundaries intersected by the single scanline. In the present example, as the heart beats the left ventricle will cyclically contract, then relax. As the left ventricle contract, the septum 40 and the posterior wall will move closer to each other. And as the left ventricle relaxes these two boundaries of the left ventricle will move apart again. This cyclical CA 02239330 1998-0~-29 , ~._ motion is revealed in the M-mode display of FIGURE 3, where the septum 40 and the posterior wall are seen to periodically approach each other and move apart again. Since tissue boundaries tend to return relatively strong echoes, the boundaries of the septum 40 and the posterior wall 42 will be defined by strong echo returns 50, 52 in each scanline.
FIGURE 3a, for instance, illustrates the strong echo returns by dots on a single scanline 48, one of many such scanlines that are returned from along the location of the M-mode cursor 46. In FIGURE 3a the smaller dots between the strong echo returns 50, 52 represent relatively weaker echo returns from the blood pool within the left ventricle (between the dots 50, 52), and from the right ventricle above the septum 40 (50) and the abdominal tissue below the posterior wall 42 (52).
In accordance with the principles of the present invention, the ultrasound system of FIGURE 1 automatically traces one or more of the borders in the M-mode display of FIGURE 3. First, the preferred embodiment performs the optional step of noise elimination on the greyscale data of the M-mode display. Any of a number of noise reduction techniques may be employed, such as filtering the greyscale signals or threshold elimination of low level noise or averaging a number of consecutively received M-lines to form each M-line of the M-mode display. However, the present inventor have found a preferable noise reduction technique to be the mathematical morphologic operations of dilation and erosion. This technique is described in the article Haralick et al., "Image Analysis Using Mathematical Morphology" in IEEE Trans. Pattern Anal. and Machine Intell., vol. PAMI-9, (July 1987), which provides a CA 02239330 1998-0~-29 , ~
_g _ number of examples of mathematical operands which may be use in erosion and dilation. Dilation and erosion are performed along each M-line from the minimum depth of the scanline to the maximum depth of the scanline using the dilation equation A(i,j) = maXi-~kSi+~(A(k'i)) followed by the erosion equation of A(i,j) = mini_~k~i+~(A(k,j)) for each point i along each of j M-lines. Dilation and erosion will remove much of the unwanted noise from the greyscale data such as speckle artifacts and noise in the blood pool of the heart. Each M-line is thus smoother with decreased noise, but without loss of the sharpness of the edges of tissue defined by i5 the greyscale data.
The second step is to develop an "edge map" of the edges of tissue boundaries represented by the greyscale M-line data. A variety of techniques may be used to detect tissue edges, such as detecting transitions between successive data points above a certain level, or differentiating the M-line data. A
preferred technique is to use an algorithm which detects edges in the greyscale data by computing the variance of the mean value in a sliding window applied to each M-line from the minimum depth to the maximum depth of each scanline. The preferred algorithm is illustrated in FIGURE 4, which is a succession of summation functions ~ (60) and squaring functions ( )2 (62), ending in a comparison (subtraction) computation 64. In the preferred algorithm each summation function is of the form k=i+~
~ (A(k,j)) (2*~+1) CA 02239330 1998-0~-29 and each squaring function is of the form ( ) = (A(i,j)) The resolution of the edge map is a function of the size of the window 2~ of successive data points upon S which the algorithm is performed, with a smaller window producing a relatively large number of sharply defined edges and a larger window producing a relatively smaller number of more broadly defined edges. For example, suppose-that the data values of the data points of one M-line are as shown in FIGURE
5a, which represents the amplitudes of the data points proc~e~;ng from the minimum depth to the maximum depth. In this plot the first peak 72 represents data points returned from the septum 40 of FIGURE 2 or 3, the second peak 76 represents data points returned from the posterior wall 42 of these drawings, and the region 74 between these two peaks represents data points returned from the blood pool of the left ventricle. When a window of three successive data points is used in the algorithm of FIGURE 4 to produce the edge map, an edge map for the M-line such as that shown in FIGURE 5b is produced.
In this plot the peaks 71 and 73 are the leading and trailing edges of the septum 70 of FIGURE 5a, and the peaks 75 and 77 are the leading and trailing edges of the posterior wall 76 of FIGURE 5a. The smaller peaks 79 represent smaller transition in the data values from the septum and blood pool.
When a wider window such as five or six data point is used for the window, a more broadly defined group of edge peaks are produced as illustrated by FIGURE 5c, in which peaks 71' and 73' are the edges of the septum 72 and peaks 75' and 77' are the edges of the posterior wall. Experimentation with a number of different window sizes and M-mode data sets has CA 02239330 1998-0~-29 , ~
led the present inventors to prefer a smaller window of three data points for the edge detection algorithm of FIGURE 4.
The next step is to identify a point on an edge which is to be traced by the edge tracker. This is shown in FIGURE 8, where an arrow cursor 80 is manipulated over the M-mode display until the arrow is pointing at one point on an edge. In this example the arrow cursor 80 is pointing at the inner border of the septum 40. When the arrow cursor is positioned as desired, the user presses a control on the system to actuate the edge tracer 38, which traces the edge selected by the user.
In accordance with the principles of the present invention, the edge indicated by the user in the M-mode display is traced using a corresponding map of ultrasonic velocity data. An example of a velocity data map is shown in FIGURE 6, which represents a color M-mode display. The color M-mode display of FIGURE 6 is produced by the ultrasound system in an interleaved fashion with the M-lines of the M-mode display of FIGURE 8. For instance, a pulse is transmitted along the M-mode cursor 46 and echoes are returned from the scanline and their amplitudes detected and displayed as one M-line in the greyscale M-mode display of FIGURE 8. A succession of further scanlines are repetitively acquired from along the M-mode cursor over time as a Doppler ensemble of scanlines. The ensemble of scanlines is Doppler processed to detect the velocity of motion at each point along the M-mode cursor at substantially the same time as the time of acquisition of the greyscale M-line. The echoes used to produce the greyscale M-line can of course also be Doppler processed as one of the ensemble of Doppler lines. The ensemble of CA 02239330 1998-0~-29 Doppler lines can also be acquired both before and after the time of acquisition of the greyscale M-line. The succession of lines of velocity data lines can be mapped in the same manner as the M-mode display as indicated in FIGURE 6, in which each new line of velocity data is displayed next to the previously acquired line of velocity data, and with color coding of the direction of the motion, hence the term color M-mode display. When the M-mode cursor intersects a chamber of the heart as does the cursor 46, the color M-mode display will vary in color with the direction of the motion of blood and the heart walls. When blue is used for the color of motion away from the transducer and red for the color of motion toward the transducer, the color M-mode display will appear as alternating bands of blue and red as indicated by the alternating B and R bands indicated in the drawing. Since opposite walls of the heart are constantly moving in opposite directions relative to the transducer, i.e., both simultaneously toward the center of the chamber or both simultaneously away from the center of the chamber, the opposite walls will exhibit opposite colors on the same vertical scanline in the color M-mode display. The color M-mode display of FIGURE 6 comprises a velocity map corresponding to the greyscale M-mode display of FIGURE 8.
The velocity information of the color M-mode display is used to track the location of edges in the greyscale M-mode display as the edges appear at different depths along the M-line from one M-line to another as illustrated by FIGURE 7. In this drawing scanlines 88 and 98 represent successively acquired M-lines of a greyscale M-mode display. The large dots 81 and 83 represent detected outer and inner ATL-lS5 CA 02239330 l998-0~-29 edges of the septum 40 of FIGURE 8, and the large dots 85 and 87 represent detected inner and outer edges of the posterior wall 42 of FIGURE 8. The second greyscale M-line 98 is acquired after a time interval t~ from the time of acquisition of M-line 88, and its edge map identifies outer and inner edges 91,93 of the septum and inner and outer edges 95,97 of the posterior wall. Between the times of acquisition of these greyscale M-lines, a color Doppler M-line 100 is acquired and processed. Since the heart is contracting during these times the opposite walls of the heart are moving toward each other as depicted by the arrows on the color M-line 100. The arrows are larger in the vicinity of the heart walls at the top and bottom of the color M-line, depicting greater velocities at these locations. In the center of the heart there is little or no motion along the color M-line, as indicated by the dots without arrows in the center of the color M-line 100.
Suppose that the user has pointed the arrow cursor 80 at the inner edge of the septum as shown in FIGURE 8. When the greyscale M-line (or edge map line) 88 is so designated, the arrow cursor 80 is pointing at septum inner edge dot 83. The edge tracker 38 now produces a trace to the location of the septum inner edge on the adjacent greyscale line.
To do this the edge tracker must identify the inner edge of the septum on the adjacent greyscale M-line.
It does this by observing the velocity at the location of inner edge dot 83, which is the velocity value of velocity arrow 102 on the color M-line 100.
As the velocity arrow 102 indicates, the velocity of the inner edge of the septum is moving toward the center of the heart chamber (down in the drawing) at CA 02239330 l998-0~-29 , this time. The quantified value of the velocity arrow will provide both speed and direction. This velocity is multiplied by the time interval t~
between greyscale lines, yielding a displacement value. The system now has computed that the inner septum edge on line 98 should be displaced downward by the calculated displacement from the corresponding location of the inner edge 83 on line 88. Moving down M-line 98 by this displacement, the system finds inner edge dot 93 on the next M-line 98.
In a constructed embodiment the edge tracker 30 follows an edge from one M-line to the next by use of this technique, but with the added benefit of a window function. This window function is indicated by bracket W1 for M-line 88, which is centered around the edge dot 83 which the system is tracking. By calculating the velocity times time displacement value shown above, the system relocates the window function to the displaced bracket location W2 as shown for M-line 98. The system then identifies an edge of the edge map which is within the window W2 as the edge which is being tracked.
The system of the present invention proceeds in the above manner to track the indicated edge from one greyscale M-line to the next, using the velocity data to compute the new location of the edge in each successive M-line. The technique of the present invention may also be used to track the interior region of tissue such as inside the myocardial wall, enabling the motion of the heart muscle to be defined by a tracing. In a constructed embodiment the edge is tracked both forward and backward in time through the M-mode display. For example, if the arrow cursor 80 is pointing at edge 93, the location of edge 83 in the previous line can be computed using the interval CA 02239330 l998-0~-29 t~, and the velocity of velocity arrow 102 but with its sign (direction) reversed. The finished edge tracing will appear as a highlighted line 110 drawn over the M-mode display in FIGURE 8.
The ability to track an edge either forward or backward allows erroneous tracings to be corrected.
The greater the distance over which an edge is traced, the greater the possibility that the tracing will move to an adjacent image border. This raises the possibility of the endpoint of a long tracing being on a border other than the desired border. To overcome this condition, the user can point to the desired border at one end of the image and execute a tracing which proceeds forward in time through the image. Then the user can point to the desired border at the other end of the image and execute a tracing which proceeds backward in time through the image.
Since both traces are likely to be accurate near their starting points and less so at their end points, the two tracings can be combined on a weighted basis, with the weighting being a function of the proximity of each trace to its starting point.
The weighted average of the two traces will generally be an accurate trace of the desired border.
The edge tracing can be displayed over the color M-mode display as shown in FIGURE 6, since the color M-mode display corresponds to the greyscale M-mode display. In FIGURE 6 the user has used the ultrasound system of the present invention to trace four edges, the inner and outer borders of both the septum 40 and the posterior heart wall 42. With the borders of the heart walls thus delineated, several quantified measures of cardiac performance can be readily calculated. For instance, the color pixels of a heart wall on each vertical color Doppler line CA 02239330 l998-0~-29 , ~-between the tracings can be summed or averaged to compute the instantaneous or mean heart wall velocity. An exemplary plot 120 of the mean velocity of the heart wall which is computed in this manner is shown in FIGURE 9 in correspondence with the color M-mode display of FIGURE 6. Alternatively, the velocity gradient along each vertical color Doppler line between the tracings can be computed to produce a quantified measure of the velocity gradient across the heart wall. Other quantified measures of heart performance will be readily apparent to those skilled in the art.
It will be appreciated that tissue boundaries can also be identified in color M-mode images by transitions in color, hue, brightness or intensity.
Consequently, the present invention can be practiced directly with color M-mode images to trace borders defined in the images and using the velocity information of the images to assist the tracking of image tis~ue characteristics.
The tracing techn;que of the present invention can also be used to follow the motion of the heart wall from one 2D image to another as shown in FIGURE
10. This drawing depicts a short axis view of the heart muscle 130, with the heart wall 120 repetitively moving toward the "+" in the center of the chamber and outward again with each contraction and relaxation of the heart. The velocities of each segment of the heart wall, depicted by the arrows in the drawing around the heart muscle, are computed from Doppler velocity scanning and processing. Since the arrows are arrayed in a radial pattern, several different scanning directions may be needed to overcome the angle sensitivity of the Doppler beam.
The velocities around the heart wall are then used to CA 02239330 1998-0~-29 ~_ track the heart wall as it moves inward and outward during the cyclic beating of the heart.
The principles of the present invention can also be applied to tracking moving tissue in three dimensional ultrasonic images. In the case of three dimensional images, it will be appreciated that the tracking of moving tissue can occur in any of the three x, y, or z directions, enabling tracing from one voxel to the next in all three dimensions. The three dimensional images can be either greyscale images or M-mode images for which velocity information can be obtained.
DOPPLER ASSISTED TRACKING OF TISSUE MOTION
This invention relates to ultrasonic diagnostic imaging systems and, in particular, to ultrasonic diagnostic imaging systems in which motion depicted in ultrasonic images is tracked through the assistance of ultrasonic velocity information.
One of the primary uses of medical diagnostic ultrasound is the diagnosis of cardiac function and disease. Echocardiography revolves around the functioning of the heart and the diagnosis of abnormal conditions of cardiac performance and pathology. Of significant interest to the cardiologist in making these diagnoses is clear delineation of the heart structure itself, and in particular the motion of the well-defined heart muscle as it cyclically pumps blood. But many patients are what are known as "difficult patients,"
meaning that the patient's physiology impedes the creation of sharply defined images from strong, clear echo returns. In those and other, more normal situations, the echocardiographer is limited to diagnosing from echocardiograms of poorly defined cardiac structures. To aid in diagnoses in these situations, ultrasound systems facilitate the tracing of the outlines of cardiac structures on the displayed echocardiogram by the echocardiographer, generally with a joystick or trackball. Quantified measures and analyses can then be performed based upon the dimensions and measurement of the outlined structures. But such tracing can be a painstaking task, since the heart is constantly contracting and relaxing and valves are continually opening and closing. For a complete and accurate diagnosis, it CA 02239330 1998-0~-29 may be necessary to define the structures of the heart in every echocardiogram throughout the cardiac cycle so that the heart muscle is continually defined while it is in motion. Realtime cardiac images are produced at a nominal rate of 30 images per second or more, presenting the echocardiographer with a daunting number of images with image structures to trace. It would be desirable, of course, for the ultrasound system to outline the structures of the heart automatically, and to maintain such outlines by tracking these structures automatically as they are engaged in the dynamics of contraction and relaxation, opening and closing.
In accordance with the principles of the present invention, an ultrasonic diagnostic imaging system is provided in which moving structures in a realtime ultrasonic image sequence are automatically tracked by use of velocity information. In a constructed embodiment the velocity information is derived from ultrasonic Doppler signals. Such automatic tracking enables the heart wall to be outlined or traced continually as it moves with the motion of the beating heart cycle. A preferred method for automatically tracing a moving structure in a realtime image sequence comprises the steps of detecting the edges of structures in the images, identifying an edge which is to be tracked or traced, and tracking or tracing the identified edge during its motional sequence with the use of velocity information concerning the identified edge.
In an illustrated embodiment, the borders of the heart wall are identified and traced in an M-mode sequence by use of the velocity characteristics of the beating heart. This embodiment leads to a number of useful diagnostic tools, such as analysis of the CA 02239330 1998-0~-29 mean velocity of the heart wall or measurement of the velocity gradient across the heart muscle.
In the drawings: .
FIGURE 1 is a block diagram of an ultrasonic diagnostic imaging system constructed in accordance with the principles of the present invention;
FIGURE 2 is a two dimensional image of the left ventricle of the heart with an M-line cursor across the heart chamber;
FIGURE 3 is a greyscale M-mode display acquired from the M-line location in FIGURE 2;
FIGURE 3a represents one M-line of the M-mode display of FIGURE 3;
FIGURE 4 is a block diagram illustrating a preferred edge detection technique for the ultrasound system of FIGURE l;
FIGURES 5a-5c are plots illustrating the operation of the preferred edge detection technique of FIGURE 4;
FIGURE 6 is a color M-mode display correlated to the greyscale M-mode display of FIGURE 3;
FIGURE 7 illustrates the tracking of moving tissue structure in accordance with the principles of the present invention;
FIGURE 8 illustrates the tracing of the heart wall in a greyscale M-mode display in accordance with the present invention;
FIGURE 9 is a graph of the mean heart wall velocity produced from the color M-mode data of FIGURE 6 in accordance with the principles of the present invention; and FIGURE 10 illustrates the application of the principles of the present invention to two dimensional greyscale images.
Referring first to FIGURE 1, an ultrasonic CA 02239330 1998-0~-29 diagnostic imaging system constructed in accordance with the principles of the present invention is shown in block diagram form. An ultrasonic array probe 10 is operated under control of a beamformer 12 to transmit ultrasonic waves into the body of a patient and receives echoes returning from the body in response to those waves. The beamformer 12 is controlled by a beamformer controller 14 which controls the time of activation of the probe transducer array elements to steer and focus a transmitted beam. The beamformer is also controlled to appropriately delay and sum received echo signals to form received beams of coherent echoes along each desired scanline.
The echoes from the received scanlines may then be processed in a variety of ways. One way is by means of a B mode processor 18, which detects the echo amplitudes. The echo amplitudes may then be displayed spatially as a greyscale or B mode image of the body. Another processing technique is by means of a Doppler processor, which processes ensembles of scanlines from a common location to detect the Doppler shift of moving blood and other tissues along the scanlines. The Doppler signals can be continually received from a sample volume in the body and displayed in a spectral Doppler display. A
second means of display is a colorflow display, in which the Doppler signals are spatially overlaid for display with a greyscale image, thereby depicting flow parameters within the structures of the greyscale image.
In accordance with one aspect of the present invention, the detected and processed echo signals are displayed in a time-motion or M-mode display format. In an M-mode display a one dimensional CA 02239330 1998-0~-29 "_ image, a single scanline termed an A-line, is repetitively scanned and displayed over time. The resultant display reveals changes occurring along the scanline location with time. M-mode displays can be formed using either Doppler or greyscale scanlines.
In the system of FIGURE 1, Doppler scanlines produced by the Doppler processor 16 are stored as color M-lines in storage device 20, and greyscale scanlines produced by the B mode processor 18 are stored as greyscale M-lines in the storage device 20. The M-mode display is most useful for displaying structures in the body which are moving, such as the beating heart, as this movement can be seen in the time-displaced sequence of scanlines of the display.
The color M-lines stored in the storage device 20 may be read from the device and converted into video signals by a video processor 22, then displayed as a color M-mode display on a display 26. The greyscale M-lines may also be applied to the video processor 22 for display as a greyscale M-mode display.
In accordance with the principles of the present invention, the edges or borders of structures in the greyscale scanlines are traced by a motion tracker 30. The motion tracker 30 uses the Doppler signals to track the movement of moving structures such as blood cells and other tissue in the greyscale scanlines. The motion tracker first removes extraneous artifacts and noise from the greyscale scanlines by means of a noise eliminator 32 and, optionally, a noise eliminator 34 eliminates noise from the Doppler scanlines. The greyscale scanlines then undergo edge detection by an edge detector 36, which identifies the locations of edges or boundaries of structures along the scanlines. Finally, the edges or boundaries depicted in a sequence of scanlines are traced with an edge tracer 38. The edge tracing is initiated by user initialization, whereby a user points to one point along the edge which is to be traced. The edge tracer then traces the edge of which the identified point is a part. It does this through use of the velocity data of the Doppler information. The Doppler information provides data as to the velocity at locations along the greyscale scanline, including the boundary point identified by the user. By knowing the time interval which lapsed between the acquisition of each greyscale scanline, and the velocity and direction of motion of an identified boundary point, the system can compute the location to which the boundary point has moved by the time of the subsequently acquired scanline. The system then looks for the boundary point in the vicinity of the computed position and produces a trace connecting the points of the two scanlines. In this way the motion tracker 30 can continuously track moving structures or structure boundaries from one scanline to the next, and can trace or outline them in the M-mode display.
The trace defined by the edge tracer 38 is applied to the video processor in concert with the greyscale M-lines, where the trace is overlaid over the M-mode display to define the traced edge or border or other characteristic in the display. The trace can also be overlaid over the temporally corresponding color M-mode display.
Once the border of a structure such as the heart wall has been defined by edge tracing, a number of quantified measurements can be made of the now precisely delineated structure. Such quantified measures are performed by a motion analysis module 24 CA 02239330 1998-0~-29 , ''-as described below, in which the computation of measures such as mean velocities or velocity gradients are discussed.
An understanding of the operation of the system of FIGURE 1 may be had with reference to the following drawings. FIGURE 2 illustrates a two dimensional ultrasound image of the heart. Typically a sector image of the heart is produced with a phased array probe, although curved and linear array probes may also be used. At the top of the image is a portion of the right ventricle, marked RV in the drawing. A septum 40 separates the right ventricle from the left ventricle LV. The chamber of the left ventricle is bounded by the posterior wall 42. At the right side of the left ventricle is the mitral valve 44.
Extending from the apex of the image and across the right and left ventricles is an M-mode cursor 46.
The M-mode cursor marks the scanline of the two dimensional (2D) image which is repetitively scanned and reproduced to form an M-mode image as shown in FIGURE 3. An M-mode image is formed by repetitively scanning the same spatial scanline and reproducing each new scan of the scanline alongside the previous scanlines. Thus, the M-mode image is a time sequence of one dimensional images of a single A-line. The benefit provided by such a display is that the image reveals the motion of tissues and tissue boundaries intersected by the single scanline. In the present example, as the heart beats the left ventricle will cyclically contract, then relax. As the left ventricle contract, the septum 40 and the posterior wall will move closer to each other. And as the left ventricle relaxes these two boundaries of the left ventricle will move apart again. This cyclical CA 02239330 1998-0~-29 , ~._ motion is revealed in the M-mode display of FIGURE 3, where the septum 40 and the posterior wall are seen to periodically approach each other and move apart again. Since tissue boundaries tend to return relatively strong echoes, the boundaries of the septum 40 and the posterior wall 42 will be defined by strong echo returns 50, 52 in each scanline.
FIGURE 3a, for instance, illustrates the strong echo returns by dots on a single scanline 48, one of many such scanlines that are returned from along the location of the M-mode cursor 46. In FIGURE 3a the smaller dots between the strong echo returns 50, 52 represent relatively weaker echo returns from the blood pool within the left ventricle (between the dots 50, 52), and from the right ventricle above the septum 40 (50) and the abdominal tissue below the posterior wall 42 (52).
In accordance with the principles of the present invention, the ultrasound system of FIGURE 1 automatically traces one or more of the borders in the M-mode display of FIGURE 3. First, the preferred embodiment performs the optional step of noise elimination on the greyscale data of the M-mode display. Any of a number of noise reduction techniques may be employed, such as filtering the greyscale signals or threshold elimination of low level noise or averaging a number of consecutively received M-lines to form each M-line of the M-mode display. However, the present inventor have found a preferable noise reduction technique to be the mathematical morphologic operations of dilation and erosion. This technique is described in the article Haralick et al., "Image Analysis Using Mathematical Morphology" in IEEE Trans. Pattern Anal. and Machine Intell., vol. PAMI-9, (July 1987), which provides a CA 02239330 1998-0~-29 , ~
_g _ number of examples of mathematical operands which may be use in erosion and dilation. Dilation and erosion are performed along each M-line from the minimum depth of the scanline to the maximum depth of the scanline using the dilation equation A(i,j) = maXi-~kSi+~(A(k'i)) followed by the erosion equation of A(i,j) = mini_~k~i+~(A(k,j)) for each point i along each of j M-lines. Dilation and erosion will remove much of the unwanted noise from the greyscale data such as speckle artifacts and noise in the blood pool of the heart. Each M-line is thus smoother with decreased noise, but without loss of the sharpness of the edges of tissue defined by i5 the greyscale data.
The second step is to develop an "edge map" of the edges of tissue boundaries represented by the greyscale M-line data. A variety of techniques may be used to detect tissue edges, such as detecting transitions between successive data points above a certain level, or differentiating the M-line data. A
preferred technique is to use an algorithm which detects edges in the greyscale data by computing the variance of the mean value in a sliding window applied to each M-line from the minimum depth to the maximum depth of each scanline. The preferred algorithm is illustrated in FIGURE 4, which is a succession of summation functions ~ (60) and squaring functions ( )2 (62), ending in a comparison (subtraction) computation 64. In the preferred algorithm each summation function is of the form k=i+~
~ (A(k,j)) (2*~+1) CA 02239330 1998-0~-29 and each squaring function is of the form ( ) = (A(i,j)) The resolution of the edge map is a function of the size of the window 2~ of successive data points upon S which the algorithm is performed, with a smaller window producing a relatively large number of sharply defined edges and a larger window producing a relatively smaller number of more broadly defined edges. For example, suppose-that the data values of the data points of one M-line are as shown in FIGURE
5a, which represents the amplitudes of the data points proc~e~;ng from the minimum depth to the maximum depth. In this plot the first peak 72 represents data points returned from the septum 40 of FIGURE 2 or 3, the second peak 76 represents data points returned from the posterior wall 42 of these drawings, and the region 74 between these two peaks represents data points returned from the blood pool of the left ventricle. When a window of three successive data points is used in the algorithm of FIGURE 4 to produce the edge map, an edge map for the M-line such as that shown in FIGURE 5b is produced.
In this plot the peaks 71 and 73 are the leading and trailing edges of the septum 70 of FIGURE 5a, and the peaks 75 and 77 are the leading and trailing edges of the posterior wall 76 of FIGURE 5a. The smaller peaks 79 represent smaller transition in the data values from the septum and blood pool.
When a wider window such as five or six data point is used for the window, a more broadly defined group of edge peaks are produced as illustrated by FIGURE 5c, in which peaks 71' and 73' are the edges of the septum 72 and peaks 75' and 77' are the edges of the posterior wall. Experimentation with a number of different window sizes and M-mode data sets has CA 02239330 1998-0~-29 , ~
led the present inventors to prefer a smaller window of three data points for the edge detection algorithm of FIGURE 4.
The next step is to identify a point on an edge which is to be traced by the edge tracker. This is shown in FIGURE 8, where an arrow cursor 80 is manipulated over the M-mode display until the arrow is pointing at one point on an edge. In this example the arrow cursor 80 is pointing at the inner border of the septum 40. When the arrow cursor is positioned as desired, the user presses a control on the system to actuate the edge tracer 38, which traces the edge selected by the user.
In accordance with the principles of the present invention, the edge indicated by the user in the M-mode display is traced using a corresponding map of ultrasonic velocity data. An example of a velocity data map is shown in FIGURE 6, which represents a color M-mode display. The color M-mode display of FIGURE 6 is produced by the ultrasound system in an interleaved fashion with the M-lines of the M-mode display of FIGURE 8. For instance, a pulse is transmitted along the M-mode cursor 46 and echoes are returned from the scanline and their amplitudes detected and displayed as one M-line in the greyscale M-mode display of FIGURE 8. A succession of further scanlines are repetitively acquired from along the M-mode cursor over time as a Doppler ensemble of scanlines. The ensemble of scanlines is Doppler processed to detect the velocity of motion at each point along the M-mode cursor at substantially the same time as the time of acquisition of the greyscale M-line. The echoes used to produce the greyscale M-line can of course also be Doppler processed as one of the ensemble of Doppler lines. The ensemble of CA 02239330 1998-0~-29 Doppler lines can also be acquired both before and after the time of acquisition of the greyscale M-line. The succession of lines of velocity data lines can be mapped in the same manner as the M-mode display as indicated in FIGURE 6, in which each new line of velocity data is displayed next to the previously acquired line of velocity data, and with color coding of the direction of the motion, hence the term color M-mode display. When the M-mode cursor intersects a chamber of the heart as does the cursor 46, the color M-mode display will vary in color with the direction of the motion of blood and the heart walls. When blue is used for the color of motion away from the transducer and red for the color of motion toward the transducer, the color M-mode display will appear as alternating bands of blue and red as indicated by the alternating B and R bands indicated in the drawing. Since opposite walls of the heart are constantly moving in opposite directions relative to the transducer, i.e., both simultaneously toward the center of the chamber or both simultaneously away from the center of the chamber, the opposite walls will exhibit opposite colors on the same vertical scanline in the color M-mode display. The color M-mode display of FIGURE 6 comprises a velocity map corresponding to the greyscale M-mode display of FIGURE 8.
The velocity information of the color M-mode display is used to track the location of edges in the greyscale M-mode display as the edges appear at different depths along the M-line from one M-line to another as illustrated by FIGURE 7. In this drawing scanlines 88 and 98 represent successively acquired M-lines of a greyscale M-mode display. The large dots 81 and 83 represent detected outer and inner ATL-lS5 CA 02239330 l998-0~-29 edges of the septum 40 of FIGURE 8, and the large dots 85 and 87 represent detected inner and outer edges of the posterior wall 42 of FIGURE 8. The second greyscale M-line 98 is acquired after a time interval t~ from the time of acquisition of M-line 88, and its edge map identifies outer and inner edges 91,93 of the septum and inner and outer edges 95,97 of the posterior wall. Between the times of acquisition of these greyscale M-lines, a color Doppler M-line 100 is acquired and processed. Since the heart is contracting during these times the opposite walls of the heart are moving toward each other as depicted by the arrows on the color M-line 100. The arrows are larger in the vicinity of the heart walls at the top and bottom of the color M-line, depicting greater velocities at these locations. In the center of the heart there is little or no motion along the color M-line, as indicated by the dots without arrows in the center of the color M-line 100.
Suppose that the user has pointed the arrow cursor 80 at the inner edge of the septum as shown in FIGURE 8. When the greyscale M-line (or edge map line) 88 is so designated, the arrow cursor 80 is pointing at septum inner edge dot 83. The edge tracker 38 now produces a trace to the location of the septum inner edge on the adjacent greyscale line.
To do this the edge tracker must identify the inner edge of the septum on the adjacent greyscale M-line.
It does this by observing the velocity at the location of inner edge dot 83, which is the velocity value of velocity arrow 102 on the color M-line 100.
As the velocity arrow 102 indicates, the velocity of the inner edge of the septum is moving toward the center of the heart chamber (down in the drawing) at CA 02239330 l998-0~-29 , this time. The quantified value of the velocity arrow will provide both speed and direction. This velocity is multiplied by the time interval t~
between greyscale lines, yielding a displacement value. The system now has computed that the inner septum edge on line 98 should be displaced downward by the calculated displacement from the corresponding location of the inner edge 83 on line 88. Moving down M-line 98 by this displacement, the system finds inner edge dot 93 on the next M-line 98.
In a constructed embodiment the edge tracker 30 follows an edge from one M-line to the next by use of this technique, but with the added benefit of a window function. This window function is indicated by bracket W1 for M-line 88, which is centered around the edge dot 83 which the system is tracking. By calculating the velocity times time displacement value shown above, the system relocates the window function to the displaced bracket location W2 as shown for M-line 98. The system then identifies an edge of the edge map which is within the window W2 as the edge which is being tracked.
The system of the present invention proceeds in the above manner to track the indicated edge from one greyscale M-line to the next, using the velocity data to compute the new location of the edge in each successive M-line. The technique of the present invention may also be used to track the interior region of tissue such as inside the myocardial wall, enabling the motion of the heart muscle to be defined by a tracing. In a constructed embodiment the edge is tracked both forward and backward in time through the M-mode display. For example, if the arrow cursor 80 is pointing at edge 93, the location of edge 83 in the previous line can be computed using the interval CA 02239330 l998-0~-29 t~, and the velocity of velocity arrow 102 but with its sign (direction) reversed. The finished edge tracing will appear as a highlighted line 110 drawn over the M-mode display in FIGURE 8.
The ability to track an edge either forward or backward allows erroneous tracings to be corrected.
The greater the distance over which an edge is traced, the greater the possibility that the tracing will move to an adjacent image border. This raises the possibility of the endpoint of a long tracing being on a border other than the desired border. To overcome this condition, the user can point to the desired border at one end of the image and execute a tracing which proceeds forward in time through the image. Then the user can point to the desired border at the other end of the image and execute a tracing which proceeds backward in time through the image.
Since both traces are likely to be accurate near their starting points and less so at their end points, the two tracings can be combined on a weighted basis, with the weighting being a function of the proximity of each trace to its starting point.
The weighted average of the two traces will generally be an accurate trace of the desired border.
The edge tracing can be displayed over the color M-mode display as shown in FIGURE 6, since the color M-mode display corresponds to the greyscale M-mode display. In FIGURE 6 the user has used the ultrasound system of the present invention to trace four edges, the inner and outer borders of both the septum 40 and the posterior heart wall 42. With the borders of the heart walls thus delineated, several quantified measures of cardiac performance can be readily calculated. For instance, the color pixels of a heart wall on each vertical color Doppler line CA 02239330 l998-0~-29 , ~-between the tracings can be summed or averaged to compute the instantaneous or mean heart wall velocity. An exemplary plot 120 of the mean velocity of the heart wall which is computed in this manner is shown in FIGURE 9 in correspondence with the color M-mode display of FIGURE 6. Alternatively, the velocity gradient along each vertical color Doppler line between the tracings can be computed to produce a quantified measure of the velocity gradient across the heart wall. Other quantified measures of heart performance will be readily apparent to those skilled in the art.
It will be appreciated that tissue boundaries can also be identified in color M-mode images by transitions in color, hue, brightness or intensity.
Consequently, the present invention can be practiced directly with color M-mode images to trace borders defined in the images and using the velocity information of the images to assist the tracking of image tis~ue characteristics.
The tracing techn;que of the present invention can also be used to follow the motion of the heart wall from one 2D image to another as shown in FIGURE
10. This drawing depicts a short axis view of the heart muscle 130, with the heart wall 120 repetitively moving toward the "+" in the center of the chamber and outward again with each contraction and relaxation of the heart. The velocities of each segment of the heart wall, depicted by the arrows in the drawing around the heart muscle, are computed from Doppler velocity scanning and processing. Since the arrows are arrayed in a radial pattern, several different scanning directions may be needed to overcome the angle sensitivity of the Doppler beam.
The velocities around the heart wall are then used to CA 02239330 1998-0~-29 ~_ track the heart wall as it moves inward and outward during the cyclic beating of the heart.
The principles of the present invention can also be applied to tracking moving tissue in three dimensional ultrasonic images. In the case of three dimensional images, it will be appreciated that the tracking of moving tissue can occur in any of the three x, y, or z directions, enabling tracing from one voxel to the next in all three dimensions. The three dimensional images can be either greyscale images or M-mode images for which velocity information can be obtained.
Claims (39)
1. In an ultrasonic diagnostic imaging system, a method for tracking a characteristic of moving tissue in time sequential ultrasonic image lines comprising the steps of:
identifying said characteristic in one of said ultrasonic image lines;
determining the velocity of said characteristic;
and using said velocity to identify said characteristic in the other of said ultrasonic image lines.
identifying said characteristic in one of said ultrasonic image lines;
determining the velocity of said characteristic;
and using said velocity to identify said characteristic in the other of said ultrasonic image lines.
2. In the ultrasonic diagnostic imaging system of Claim 1, wherein said time sequential ultrasonic image lines comprise lines of an M-mode display.
3. In the ultrasonic diagnostic imaging system of Claim 2, wherein said time sequential ultrasonic image lines comprise greyscale lines of an M-mode display.
4. The method of Claim 1 or 2, wherein said characteristic comprises a tissue boundary.
5. The method of Claim 1, 2, or 4, wherein said velocity is determined by an ultrasonic Doppler technique.
6. The method of Claim 5, wherein said velocity is determined in close time proximity to the acquisition of one of said ultrasonic image lines.
7. The method of Claim 3, wherein said velocity is determined from Doppler M-line data.
8. The method of Claim 7, wherein said greyscale lines and said Doppler M-line data are acquired in a time interleaved sequence.
9. In an ultrasonic diagnostic imaging system, apparatus for tracking a characteristic of moving objects in the body comprising:
means for acquiring a plurality of temporally different echoes of an object;
means for identifying a characteristic of said object in one of said echoes;
means for determining the velocity of said characteristic; and means for using said velocity to identify said characteristic in a temporally different echo.
means for acquiring a plurality of temporally different echoes of an object;
means for identifying a characteristic of said object in one of said echoes;
means for determining the velocity of said characteristic; and means for using said velocity to identify said characteristic in a temporally different echo.
10. In the ultrasonic diagnostic imaging system of Claim 9, wherein said means for acquiring comprises means for acquiring lines of an M-mode display.
11. In the ultrasonic diagnostic imaging system of Claim 9, wherein said means for identifying includes a user controlled cursor.
12. In the ultrasonic diagnostic imaging system of Claim 9, wherein said means for determining the velocity comprises a Doppler signal processor.
13. A method for tracking a tissue boundary in an ultrasonic M-mode display of greyscale M-lines comprising the steps of:
producing a map of edge characteristics from said greyscale M-lines;
producing a map of velocity values corresponding to said map of edge characteristics;
denoting a tissue boundary to be traced; and using velocity values corresponding to the edge characteristics of said denoted tissue boundary to track said denoted tissue boundary.
producing a map of edge characteristics from said greyscale M-lines;
producing a map of velocity values corresponding to said map of edge characteristics;
denoting a tissue boundary to be traced; and using velocity values corresponding to the edge characteristics of said denoted tissue boundary to track said denoted tissue boundary.
14. The method of Claim 13, wherein said step of producing an edge map comprises detecting the variance of the mean level along each M-line.
15. The method of Claim 13, wherein said step of producing a map of velocity values comprises acquiring a sequence of Doppler M-lines.
16. The method of Claim 15, wherein said Doppler M-lines are acquired in a time interleaved sequence with said greyscale M-lines.
17. The method of Claim 13, further comprising the step of reducing noise in said greyscale M-lines.
18. The method of Claim 17, wherein said step of reducing noise comprises performing the mathematical morphologic operations of dilation and erosion on each M-line.
19. In a diagnostic ultrasound system, including means for acquiring M-mode lines of tissue of the body, apparatus for tracing a tissue boundary in an M-mode image comprising:
means for processing said M-mode lines to define the edges of tissue in said lines; and means, responsive to said defined tissue edges, for tracing a tissue boundary across a plurality of said lines.
means for processing said M-mode lines to define the edges of tissue in said lines; and means, responsive to said defined tissue edges, for tracing a tissue boundary across a plurality of said lines.
20. The diagnostic ultrasound system of Claim 19, wherein said processing means comprises means for computing the variance of the mean value of data points along said lines.
21. The diagnostic ultrasound system of Claim 19, wherein said M-mode lines comprise lines of greyscale data, and further comprising means for producing a velocity map corresponding to said M-mode lines.
22. The diagnostic ultrasound system of Claim 21, wherein said tracing means is further responsive to said velocity map for tracing a tissue boundary.
23. The diagnostic ultrasound system of Claim 22, wherein said velocity map comprises Doppler M-mode lines corresponding to said lines of greyscale data.
23. The diagnostic ultrasound system of Claim 22, wherein said velocity map comprises Doppler M-mode lines corresponding to said lines of greyscale data.
23. The diagnostic ultrasound system of Claim 19 or 22, further comprising means for reducing noise in said lines of greyscale data.
24. In a diagnostic ultrasound system, a method for tracking a characteristic of moving tissue through a temporally acquired sequence of two dimensional ultrasonic images comprising the steps of:
identifying said characteristic in one of said ultrasonic images;
obtaining the velocity of said identified characteristic; and using said velocity to identify said characteristic in another of said ultrasonic images.
identifying said characteristic in one of said ultrasonic images;
obtaining the velocity of said identified characteristic; and using said velocity to identify said characteristic in another of said ultrasonic images.
25. In the diagnostic ultrasound system of Claim 24, wherein the step of obtaining the velocity of said identified characteristic comprises obtaining said velocity by means of ultrasonic Doppler processing.
26. In the diagnostic ultrasound system of Claim 25, wherein said velocity is obtained at approximately the time of acquisition of said one of said ultrasonic images.
27. In a diagnostic ultrasonic imaging system, apparatus for diagnosing the human cardiovascular system comprising:
a processing system for producing an M-mode display of temporally acquired A-line data;
a processing system for producing an M-mode display of temporally acquired Doppler data temporally corresponding to said A-line data;
a tracing system for producing a trace of an image characteristic of said M-mode display of A-line data; and display means, responsive to said trace, for displaying said trace on said M-mode display of Doppler data.
a processing system for producing an M-mode display of temporally acquired A-line data;
a processing system for producing an M-mode display of temporally acquired Doppler data temporally corresponding to said A-line data;
a tracing system for producing a trace of an image characteristic of said M-mode display of A-line data; and display means, responsive to said trace, for displaying said trace on said M-mode display of Doppler data.
28. In the diagnostic ultrasonic imaging system of Claim 27, further comprising analysis software, responsive to said trace, for performing analysis of the information of said M-mode display of Doppler data.
29. In a diagnostic ultrasound system, a method for tracking a characteristic of moving tissue in a three dimensional voxel image comprising the steps of:
identifying said characteristic in one of the voxels of said three dimensional image;
obtaining the velocity of said identified characteristic; and using said velocity to identify said characteristic in another of said image voxels.
identifying said characteristic in one of the voxels of said three dimensional image;
obtaining the velocity of said identified characteristic; and using said velocity to identify said characteristic in another of said image voxels.
30. The method of Claim 29, wherein said three dimensional voxel image comprises a three dimensional M-mode image.
31. A method for tracking a tissue boundary in an ultrasonic color M-mode display of Doppler M-lines comprising the steps of:
producing a map of edge characteristics from said Doppler M-lines;
denoting a tissue boundary to be traced; and using velocity values corresponding to the edge characteristics of said denoted tissue boundary to track said denoted tissue boundary.
producing a map of edge characteristics from said Doppler M-lines;
denoting a tissue boundary to be traced; and using velocity values corresponding to the edge characteristics of said denoted tissue boundary to track said denoted tissue boundary.
32. The method of Claim 31, wherein said Doppler M-lines contain velocity information, and wherein said velocity information is used to track said denoted tissue boundary.
33. The method of Claim 31, wherein said Doppler M-lines contain Doppler signal intensity information.
34. In a diagnostic ultrasound system, including means for acquiring M-mode lines of tissue of the body, apparatus for tracing a tissue boundary in an M-mode image comprising:
means for initiating the tracing of a tissue boundary in said image in a first direction;
means for initiating the tracing of said tissue boundary in said image in a second direction which is opposite to the first direction; and means for combining said tracings to define the tissue boundary.
means for initiating the tracing of a tissue boundary in said image in a first direction;
means for initiating the tracing of said tissue boundary in said image in a second direction which is opposite to the first direction; and means for combining said tracings to define the tissue boundary.
35. In the diagnostic ultrasound system of Claim 34, wherein said combining means further comprises means for combining said tracings on a weighted basis.
36. In the diagnostic ultrasound system of Claim 35, wherein said combining means further comprises means for combining said tracings on a weighted basis as a function of the proximity of a trace from its point of initiation.
37. In the diagnostic ultrasound system of Claim 34, wherein said means for initiating the tracing of a tissue boundary in said image in a first direction further comprises means for initiating the tracing of a tissue boundary in a first direction from one side of said image; and wherein said means for initiating the tracing of said tissue boundary in said image in a second direction further comprises means for initiating the tracing of a tissue boundary in a second direction from the other side of said image.
38. In the diagnostic ultrasound system of Claim 34, wherein each of said means for initiating a tracing comprises means for utilizing ultrasonic motion signals to trace a tissue boundary.
39. In the diagnostic ultrasound system of Claim 38, wherein said ultrasonic motion signals comprise Doppler velocity signals.
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US08/865,340 | 1997-05-29 | ||
| US08/865,340 US5800356A (en) | 1997-05-29 | 1997-05-29 | Ultrasonic diagnostic imaging system with doppler assisted tracking of tissue motion |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CA2239330A1 true CA2239330A1 (en) | 1998-11-29 |
Family
ID=25345289
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CA002239330A Abandoned CA2239330A1 (en) | 1997-05-29 | 1998-05-29 | Ultrasonic diagnostic imaging system with doppler assisted tracking of tissue motion |
Country Status (9)
| Country | Link |
|---|---|
| US (1) | US5800356A (en) |
| EP (1) | EP0880937B1 (en) |
| JP (1) | JPH1142227A (en) |
| AT (1) | ATE223177T1 (en) |
| AU (1) | AU6979198A (en) |
| CA (1) | CA2239330A1 (en) |
| DE (1) | DE69807575T2 (en) |
| IL (1) | IL124491A0 (en) |
| NO (1) | NO982447L (en) |
Families Citing this family (92)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6086539A (en) | 1996-12-04 | 2000-07-11 | Acuson Corporation | Methods and apparatus for ultrasound image quantification |
| US6030345A (en) * | 1997-05-22 | 2000-02-29 | Acuson Corporation | Method and system for ultrasound enhanced-resolution spectral Doppler |
| EP0908137A1 (en) * | 1997-10-06 | 1999-04-14 | Technologiestichting STW | A method and apparatus for making an image of a lumen or other body cavity and its surrounding tissue |
| 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 |
| EP0961577A1 (en) * | 1997-11-18 | 1999-12-08 | Koninklijke Philips Electronics N.V. | Method for the processing of signals relating to an object having moving parts and echographic device for carrying out this method |
| JP4116122B2 (en) * | 1997-11-28 | 2008-07-09 | 株式会社東芝 | Ultrasonic diagnostic apparatus and ultrasonic image processing apparatus |
| JPH11221210A (en) * | 1998-02-09 | 1999-08-17 | Aloka Co Ltd | Ultrasonograph |
| JP2889568B1 (en) * | 1998-05-18 | 1999-05-10 | 正男 伊藤 | Vascular thickness measurement device and arteriosclerosis diagnosis device |
| US6196972B1 (en) * | 1998-11-11 | 2001-03-06 | Spentech, Inc. | Doppler ultrasound method and apparatus for monitoring blood flow |
| US6464641B1 (en) * | 1998-12-01 | 2002-10-15 | Ge Medical Systems Global Technology Company Llc | Method and apparatus for automatic vessel tracking in ultrasound imaging |
| EP1055133A1 (en) | 1998-12-15 | 2000-11-29 | Koninklijke Philips Electronics N.V. | Ultrasound process and apparatus for the determination of the location of a parietal surface in a tissue and of the absolute radius of an artery |
| FR2790657B1 (en) * | 1999-03-12 | 2001-08-03 | Iodp | METHOD FOR MEASURING THE THICKNESS OF A MUSCLE FROM AN IMAGE |
| US6193660B1 (en) | 1999-03-31 | 2001-02-27 | Acuson Corporation | Medical diagnostic ultrasound system and method for region of interest determination |
| US6213945B1 (en) * | 1999-08-18 | 2001-04-10 | Acuson Corporation | Ultrasound system and method for generating a graphical vascular report |
| US6352507B1 (en) * | 1999-08-23 | 2002-03-05 | G.E. Vingmed Ultrasound As | Method and apparatus for providing real-time calculation and display of tissue deformation in ultrasound imaging |
| US6415046B1 (en) | 1999-10-07 | 2002-07-02 | Edmund Kenneth Kerut, Sr. | Method and apparatus for the early detection of tissue pathology using wavelet transformation |
| JP3386025B2 (en) * | 1999-12-15 | 2003-03-10 | 株式会社ニコン | Image feature extraction apparatus, image feature extraction method, monitoring inspection system, semiconductor exposure system, and interface system |
| JP2003518404A (en) * | 1999-12-28 | 2003-06-10 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | Ultrasonic image processing method and apparatus for displaying an ultrasonic color-coded image of an object having a moving part |
| NL1018864C2 (en) | 2001-08-31 | 2003-03-03 | Technologiestichting Stw | Device and method for generating three-dimensional images with tissue hardness information. |
| US20050124898A1 (en) * | 2002-01-16 | 2005-06-09 | Ep Medsystems, Inc. | Method and apparatus for isolating a catheter interface |
| US20080146943A1 (en) * | 2006-12-14 | 2008-06-19 | Ep Medsystems, Inc. | Integrated Beam Former And Isolation For An Ultrasound Probe |
| WO2003071950A1 (en) * | 2002-02-27 | 2003-09-04 | Amid Srl | M-tracking for space-time imaging |
| JP4060615B2 (en) * | 2002-03-05 | 2008-03-12 | 株式会社東芝 | Image processing apparatus and ultrasonic diagnostic apparatus |
| 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 |
| IL149844A (en) | 2002-05-23 | 2010-11-30 | Hadasit Med Res Service | Method, system and apparatus for evaluation of esophageal function |
| US7314446B2 (en) * | 2002-07-22 | 2008-01-01 | Ep Medsystems, Inc. | Method and apparatus for time gating of medical images |
| US20070167809A1 (en) * | 2002-07-22 | 2007-07-19 | Ep Medsystems, Inc. | Method and System For Estimating Cardiac Ejection Volume And Placing Pacemaker Electrodes Using Speckle Tracking |
| US20050245822A1 (en) * | 2002-07-22 | 2005-11-03 | Ep Medsystems, Inc. | Method and apparatus for imaging distant anatomical structures in intra-cardiac ultrasound imaging |
| US20070083118A1 (en) * | 2002-07-22 | 2007-04-12 | Ep Medsystems, Inc. | Method and System For Estimating Cardiac Ejection Volume Using Ultrasound Spectral Doppler Image Data |
| JP4192545B2 (en) * | 2002-09-27 | 2008-12-10 | 株式会社日立メディコ | Ultrasonic diagnostic equipment |
| US6835177B2 (en) * | 2002-11-06 | 2004-12-28 | Sonosite, Inc. | Ultrasonic blood vessel measurement apparatus and method |
| KR20050084366A (en) * | 2002-12-18 | 2005-08-26 | 코닌클리케 필립스 일렉트로닉스 엔.브이. | Ultrasonic doppler system for determining movement of artery walls |
| US7727153B2 (en) * | 2003-04-07 | 2010-06-01 | Sonosite, Inc. | Ultrasonic blood vessel measurement apparatus and method |
| EP1616201A1 (en) * | 2003-04-15 | 2006-01-18 | Koninklijke Philips Electronics N.V. | Heart wall strain imaging |
| US7092749B2 (en) * | 2003-06-11 | 2006-08-15 | Siemens Medical Solutions Usa, Inc. | System and method for adapting the behavior of a diagnostic medical ultrasound system based on anatomic features present in ultrasound images |
| US20050131474A1 (en) * | 2003-12-11 | 2005-06-16 | Ep Medsystems, Inc. | Systems and methods for pacemaker programming |
| US8900149B2 (en) * | 2004-04-02 | 2014-12-02 | Teratech Corporation | Wall motion analyzer |
| US7507205B2 (en) * | 2004-04-07 | 2009-03-24 | St. Jude Medical, Atrial Fibrillation Division, Inc. | Steerable ultrasound catheter |
| US7654958B2 (en) * | 2004-04-20 | 2010-02-02 | St. Jude Medical, Atrial Fibrillation Division, Inc. | Method and apparatus for ultrasound imaging with autofrequency selection |
| US20060122505A1 (en) * | 2004-11-23 | 2006-06-08 | Ep Medsystems, Inc. | M-Mode presentation of an ultrasound scan |
| US7713210B2 (en) * | 2004-11-23 | 2010-05-11 | St. Jude Medical, Atrial Fibrillation Division, Inc. | Method and apparatus for localizing an ultrasound catheter |
| US9814439B2 (en) * | 2005-01-19 | 2017-11-14 | Siemens Medical Solutions Usa, Inc. | Tissue motion comparison display |
| US7771358B2 (en) * | 2005-05-20 | 2010-08-10 | Spentech, Inc. | System and method for grading microemboli monitored by a multi-gate doppler ultrasound system |
| US8162837B2 (en) * | 2005-06-13 | 2012-04-24 | Spentech, Inc. | Medical doppler ultrasound system for locating and tracking blood flow |
| US8070684B2 (en) * | 2005-12-14 | 2011-12-06 | St. Jude Medical, Atrial Fibrillation Division, Inc. | Method and system for evaluating valvular function |
| US20070167793A1 (en) * | 2005-12-14 | 2007-07-19 | Ep Medsystems, Inc. | Method and system for enhancing spectral doppler presentation |
| US20070196005A1 (en) * | 2006-02-23 | 2007-08-23 | White Christopher A | Feature Tracing Process for M-mode Images |
| US20100234730A1 (en) * | 2006-03-31 | 2010-09-16 | National University Corporation Kyoto Institute Of Technology | Image processing device, ultrasonic imaging apparatus including the same, and image processing method |
| US20070232949A1 (en) * | 2006-03-31 | 2007-10-04 | Ep Medsystems, Inc. | Method For Simultaneous Bi-Atrial Mapping Of Atrial Fibrillation |
| US20070299479A1 (en) * | 2006-06-27 | 2007-12-27 | Ep Medsystems, Inc. | Method for Reversing Ventricular Dyssynchrony |
| US20080009733A1 (en) * | 2006-06-27 | 2008-01-10 | Ep Medsystems, Inc. | Method for Evaluating Regional Ventricular Function and Incoordinate Ventricular Contraction |
| US8107694B2 (en) * | 2006-07-20 | 2012-01-31 | Ultrasound Medical Devices, Inc. | Method of tracking speckle displacement between two images |
| US20080021945A1 (en) * | 2006-07-20 | 2008-01-24 | James Hamilton | Method of processing spatial-temporal data processing |
| US20080021319A1 (en) * | 2006-07-20 | 2008-01-24 | James Hamilton | Method of modifying data acquisition parameters of an ultrasound device |
| US20100138191A1 (en) * | 2006-07-20 | 2010-06-03 | James Hamilton | Method and system for acquiring and transforming ultrasound data |
| US20080119735A1 (en) * | 2006-11-20 | 2008-05-22 | Sonowise, Inc. | Ultrasound imaging system and method with offset alternate-mode line |
| US20080146942A1 (en) * | 2006-12-13 | 2008-06-19 | Ep Medsystems, Inc. | Catheter Position Tracking Methods Using Fluoroscopy and Rotational Sensors |
| US8187190B2 (en) * | 2006-12-14 | 2012-05-29 | St. Jude Medical, Atrial Fibrillation Division, Inc. | Method and system for configuration of a pacemaker and for placement of pacemaker electrodes |
| US20080146940A1 (en) * | 2006-12-14 | 2008-06-19 | Ep Medsystems, Inc. | External and Internal Ultrasound Imaging System |
| US8317711B2 (en) | 2007-06-16 | 2012-11-27 | St. Jude Medical, Atrial Fibrillation Division, Inc. | Oscillating phased-array ultrasound imaging catheter system |
| US8057394B2 (en) | 2007-06-30 | 2011-11-15 | St. Jude Medical, Atrial Fibrillation Division, Inc. | Ultrasound image processing to render three-dimensional images from two-dimensional images |
| US20100185093A1 (en) * | 2009-01-19 | 2010-07-22 | James Hamilton | System and method for processing a real-time ultrasound signal within a time window |
| US20100185085A1 (en) * | 2009-01-19 | 2010-07-22 | James Hamilton | Dynamic ultrasound processing using object motion calculation |
| US9275471B2 (en) | 2007-07-20 | 2016-03-01 | Ultrasound Medical Devices, Inc. | Method for ultrasound motion tracking via synthetic speckle patterns |
| US8052607B2 (en) | 2008-04-22 | 2011-11-08 | St. Jude Medical, Atrial Fibrillation Division, Inc. | Ultrasound imaging catheter with pivoting head |
| WO2010039555A1 (en) * | 2008-09-23 | 2010-04-08 | Ultrasound Medical Devices, Inc. | System and method for flexible rate processing of ultrasound data |
| US20100081937A1 (en) * | 2008-09-23 | 2010-04-01 | James Hamilton | System and method for processing a real-time ultrasound signal within a time window |
| JP2012503501A (en) | 2008-09-25 | 2012-02-09 | シーエーイー ヘルスケア インク | Simulation of medical image diagnosis |
| US8199994B2 (en) * | 2009-03-13 | 2012-06-12 | International Business Machines Corporation | Automatic analysis of cardiac M-mode views |
| EP2238912B1 (en) * | 2009-04-07 | 2016-01-27 | Samsung Medison Co., Ltd. | Ultrasound system and method of providing color M mode image and brightness M mode image |
| US8343053B2 (en) * | 2009-07-21 | 2013-01-01 | Siemens Medical Solutions Usa, Inc. | Detection of structure in ultrasound M-mode imaging |
| WO2011027252A1 (en) | 2009-09-04 | 2011-03-10 | Koninklijke Philips Electronics, N.V. | Ultrasonic elastographic imaging of relative strain ratios |
| CN102483455A (en) | 2009-09-04 | 2012-05-30 | 皇家飞利浦电子股份有限公司 | Ultrasonic elastographic strain imaging with fade control |
| WO2012035472A1 (en) | 2010-09-16 | 2012-03-22 | Koninklijke Philips Electronics N.V. | Quantification of tissue strain in ultrasonic elastography images |
| JP2012071037A (en) * | 2010-09-29 | 2012-04-12 | Fujifilm Corp | Ultrasonic image diagnostic apparatus, ultrasonic image forming method, and program |
| US8939909B2 (en) | 2011-10-28 | 2015-01-27 | Decision Sciences International Corporation | Spread spectrum coded waveforms in ultrasound imaging |
| JP2016506809A (en) * | 2013-02-11 | 2016-03-07 | コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． | Ultrasound imaging system and method |
| US8995021B2 (en) * | 2013-06-12 | 2015-03-31 | Xerox Corporation | Black trapping methods, apparatus and systems for binary images |
| US9844359B2 (en) | 2013-09-13 | 2017-12-19 | Decision Sciences Medical Company, LLC | Coherent spread-spectrum coded waveforms in synthetic aperture image formation |
| CN105426927B (en) * | 2014-08-26 | 2019-05-10 | 东芝医疗系统株式会社 | Medical image processing devices, medical image processing method and medical image equipment |
| WO2016138257A1 (en) | 2015-02-25 | 2016-09-01 | Decision Sciences Medical Company, LLC | Acoustic signal transmission couplants and coupling mediums |
| JP6412458B2 (en) * | 2015-03-31 | 2018-10-24 | セコム株式会社 | Ultrasonic sensor |
| AU2016334258B2 (en) | 2015-10-08 | 2021-07-01 | Decision Sciences Medical Company, LLC | Acoustic orthopedic tracking system and methods |
| US20190029651A1 (en) | 2016-02-26 | 2019-01-31 | Koninklijke Philips N.V. | Clutter filters for strain and other ultrasonic deformation imaging |
| KR20190031567A (en) | 2016-08-03 | 2019-03-26 | 피아이-하베스트 홀딩스 아게 | System and method for non-invasive measurement of intra-body pressure including intra-arterial blood pressure |
| 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 |
| WO2018051265A1 (en) | 2016-09-15 | 2018-03-22 | Koninklijke Philips N.V. | Ultrasonic elastographic pre-load measurement and display |
| CN109414246B (en) * | 2016-11-09 | 2021-09-14 | 深圳市理邦精密仪器股份有限公司 | System and method for Doppler spectrum time duration |
| JP7081143B2 (en) * | 2017-12-27 | 2022-06-07 | セイコーエプソン株式会社 | Ultrasonic device and ultrasonic measurement method |
| US11154274B2 (en) | 2019-04-23 | 2021-10-26 | Decision Sciences Medical Company, LLC | Semi-rigid acoustic coupling articles for ultrasound diagnostic and treatment applications |
| WO2022073242A1 (en) * | 2020-10-10 | 2022-04-14 | 深圳迈瑞生物医疗电子股份有限公司 | Ultrasonic imaging method in combination with physiological signal and electronic device |
| JP2023549818A (en) | 2020-11-13 | 2023-11-29 | ディスィジョン サイエンシズ メディカル カンパニー，エルエルシー | System and method for synthetic aperture ultrasound imaging of objects |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4398540A (en) * | 1979-11-05 | 1983-08-16 | Tokyo Shibaura Denki Kabushiki Kaisha | Compound mode ultrasound diagnosis apparatus |
| JPS5883942A (en) * | 1981-11-12 | 1983-05-19 | 株式会社東芝 | Ultrasonic pulse doppler apparatus |
| JPH0793927B2 (en) * | 1990-11-02 | 1995-10-11 | 富士通株式会社 | Ultrasonic color Doppler diagnostic device |
| US5195521A (en) * | 1990-11-09 | 1993-03-23 | Hewlett-Packard Company | Tissue measurements |
| US5622174A (en) * | 1992-10-02 | 1997-04-22 | Kabushiki Kaisha Toshiba | Ultrasonic diagnosis apparatus and image displaying system |
| US5435310A (en) * | 1993-06-23 | 1995-07-25 | University Of Washington | Determining cardiac wall thickness and motion by imaging and three-dimensional modeling |
| JP3045642B2 (en) * | 1994-01-25 | 2000-05-29 | アロカ株式会社 | Ultrasound diagnostic equipment |
| JP3833282B2 (en) * | 1994-06-24 | 2006-10-11 | 株式会社東芝 | Ultrasonic diagnostic equipment |
| NO943696D0 (en) * | 1994-10-04 | 1994-10-04 | Vingmed Sound As | Method of ultrasound imaging |
-
1997
- 1997-05-29 US US08/865,340 patent/US5800356A/en not_active Expired - Fee Related
-
1998
- 1998-05-14 IL IL12449198A patent/IL124491A0/en unknown
- 1998-05-28 AT AT98304234T patent/ATE223177T1/en not_active IP Right Cessation
- 1998-05-28 AU AU69791/98A patent/AU6979198A/en not_active Abandoned
- 1998-05-28 NO NO982447A patent/NO982447L/en unknown
- 1998-05-28 EP EP98304234A patent/EP0880937B1/en not_active Expired - Lifetime
- 1998-05-28 DE DE69807575T patent/DE69807575T2/en not_active Expired - Fee Related
- 1998-05-29 CA CA002239330A patent/CA2239330A1/en not_active Abandoned
- 1998-05-29 JP JP10149355A patent/JPH1142227A/en not_active Withdrawn
Also Published As
| Publication number | Publication date |
|---|---|
| EP0880937A3 (en) | 1999-06-09 |
| JPH1142227A (en) | 1999-02-16 |
| NO982447L (en) | 1998-11-30 |
| EP0880937B1 (en) | 2002-09-04 |
| US5800356A (en) | 1998-09-01 |
| EP0880937A2 (en) | 1998-12-02 |
| DE69807575T2 (en) | 2003-08-07 |
| NO982447D0 (en) | 1998-05-28 |
| ATE223177T1 (en) | 2002-09-15 |
| DE69807575D1 (en) | 2002-10-10 |
| IL124491A0 (en) | 1998-12-06 |
| AU6979198A (en) | 1998-12-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5800356A (en) | Ultrasonic diagnostic imaging system with doppler assisted tracking of tissue motion | |
| US11510651B2 (en) | Ultrasonic diagnosis of cardiac performance using heart model chamber segmentation with user control | |
| US8144956B2 (en) | Ultrasonic diagnosis by quantification of myocardial performance | |
| EP1673013B1 (en) | Ultrasonic cardiac volume quantification | |
| US6976961B2 (en) | Tissue motion analysis medical diagnostic ultrasound system and method | |
| EP2030042B1 (en) | Quantification and display of cardiac chamber wall thickening | |
| JP5764058B2 (en) | Ultrasound assessment of cardiac synchrony and viability | |
| EP1778093B1 (en) | Ultrasonic diagnosis of ischemic cardiodisease | |
| EP1189074B1 (en) | Method and apparatus for locking sample volume onto moving vessel in pulsed doppler ultrasound imaging | |
| US20030171668A1 (en) | Image processing apparatus and ultrasonic diagnosis apparatus | |
| US20050288589A1 (en) | Surface model parametric ultrasound imaging | |
| US20060241457A1 (en) | Cyclical information determination with medical diagnostic ultrasound | |
| US20180192987A1 (en) | Ultrasound systems and methods for automatic determination of heart chamber characteristics | |
| EP3267896B1 (en) | Ultrasonic diagnosis of cardiac performance by single degree of freedom chamber segmentation | |
| EP1876567A1 (en) | A method of determining the time dependent behavior of non rigid moving objects, particularly of biological tissues, from echographic ultrasound imaging data | |
| US20230240645A1 (en) | Systems and methods for measuring cardiac stiffness |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| FZDE | Discontinued | ||
| FZDE | Discontinued |
Effective date: 20010529 |