The present invention relates generally to intravascular ultrasound imaging systems, and in particular to mechanically-scanned intravascular ultrasound (IVUS) imaging devices, systems, and methods directed to forming a cross-sectional image of a blood vessel and measuring the velocity of blood flow within the vessel.
Intravascular ultrasound imaging is widely used in interventional cardiology as a diagnostic tool for assessing a diseased vessel, such as an artery, within the human body to determine the need for treatment, to guide the intervention, and/or to assess its effectiveness. IVUS imaging uses ultrasound echoes to form a cross-sectional image of a vessel of interest. Typically, an ultrasound transducer on an IVUS catheter both emits ultrasound pulses and receives the reflected ultrasound echoes. The ultrasound waves pass easily through most tissues and blood, but they are partially reflected from discontinuities arising from tissue structures (such as the various layers of the vessel wall), red blood cells, and other features of interest. The IVUS imaging system, which is connected to the IVUS catheter by way of a patient interface module (PIM), processes the received ultrasound echoes to produce a cross-sectional image of the vessel where the transducer is located.
To establish the need for treatment, the IVUS system is used to measure the lumen diameter or cross-sectional area of the vessel. For this purpose, it is important to distinguish blood from vessel wall tissue so that the luminal border can be accurately identified. In an IVUS image, the blood echoes are distinguished from tissue echoes by slight differences in the strengths of the echoes (e.g., vessel wall echoes are generally stronger than blood echoes) and from subtle differences in the texture of the image (i.e., speckle) arising from structural differences between blood and vessel wall tissue. As IVUS imaging has evolved, there has been a steady migration towards higher ultrasound frequencies to improve the resolution in the display. But as ultrasound frequency is increased, there is diminished contrast between the blood echoes and vessel wall tissue echoes. At the 20 MHz center frequency used in early generations of IVUS, the blood echoes are very weak in comparison to the vessel wall echoes due to the small size of the red blood cell compared to the acoustic wavelength. However, at the 40 MHz ultrasound center frequency now commonly used for IVUS imaging, there is only a modest difference between blood and tissue echoes because the ultrasound wavelength at this higher frequency is closer to the dimensions of the red blood cells.
Another use of IVUS imaging in interventional cardiology is to help identify the most appropriate course of treatment. For example, IVUS imaging may be used to assist in recognizing the presence of a mural thrombus (i.e., coagulated blood attached to the vessel wall and stationary within the blood vessel) in an artery prior to initiating treatment. If a thrombus is identified in a region where disease has caused a localized narrowing of the arterial lumen, then the treatment plan could be modified to include aspiration (i.e., removal) of the thrombus prior to placing a stent in the artery to expand and stabilize the cross-sectional area of the vessel. In addition, the identification of a thrombus could lead the physician to order a more aggressive course of anti-coagulant drug therapy to prevent the subsequent occurrence of potentially deadly thrombosis. In a conventional IVUS image, however, there is very little difference in appearance between stationary thrombi and moving blood.
Another use of IVUS imaging in interventional cardiology is to visualize the proper deployment of a stent within an artery. A stent is an expandable cylinder that is generally deployed within the artery to enlarge and/or stabilize the lumen of the artery. The expansion of the stent typically stretches the vessel and displaces the plaque that otherwise forms a partial obstruction of the vessel lumen. The expanded stent forms a scaffold propping the vessel lumen open and preventing elastic recoil of the vessel wall after it has been stretched. In this context, it is important to recognize proper stent apposition; that is, the stent struts should be pressed firmly against the vessel wall. A poorly deployed stent may leave stent struts in the stream of the blood flow and these exposed stent struts are prone to initiate thrombus formation. Thrombus formation following stent deployment is referred to as “late stent thrombosis” and these thrombi can occlude the artery or break free from the stent strut to occlude a downstream branch of a coronary artery and trigger a heart attack.
In these examples of IVUS imaging, it is particularly useful to identify moving blood, and to distinguish the moving blood from the relatively stationary or static tissue or thrombi. Motion information can be helpful in delineating the interface between blood and vessel wall so that the luminal boundary can be more easily and accurately identified. Motion parameters such as velocity may be the most robust ultrasound-detectable parameters for distinguishing moving blood from stationary thrombi. In the case of stent malapposition, the observation of moving blood behind a stent strut is a clear indication that the stent strut is not firmly pressed against the vessel wall as it should be, possibly indicating a need to further expand the stent. In each of the aforementioned IVUS imaging examples, the addition of motion parameters to the traditional IVUS display of echo amplitude can improve the diagnosis and treatment of a patient.
Traditionally, IVUS catheters, whether rotational or solid-state catheters, are side-looking devices, wherein the ultrasound pulses are transmitted substantially perpendicular to the axis of the catheter to produce a cross-sectional image representing a slice through the blood vessel. The blood flow in the vessel is normally parallel to the axis of the catheter and perpendicular to the plane of the image. IVUS images are typically presented in a grey-scale format, with strong reflectors (vessel boundary, calcified tissue, metal stents, etc.) displayed as bright (white) pixels, with weaker echoes (blood and soft tissue) displayed as dark (grey or black) pixels. Thus, flowing blood and may appear very similar to soft tissue or static blood (i.e., thrombi) in a traditional IVUS display.
In non-invasive ultrasound imaging applications, Doppler ultrasound methods are often used to measure blood and tissue velocity, and the velocity information is used to distinguish moving blood echoes from stationary tissue echoes. Commonly, the velocity information is used to colorize the grey-scale ultrasound image in a format referred to as Doppler color flow ultrasound imaging, with fast moving blood tinted red or blue, depending on its direction of flow, and with stationary tissue displayed in grey scale.
Traditionally, IVUS imaging has not been amenable to Doppler color flow imaging since the direction of blood flow is predominantly perpendicular to the IVUS imaging plane. More specifically, Doppler color flow imaging and other Doppler techniques do not function well when the velocity of interest (i.e., blood flow velocity) is perpendicular to the imaging plane and perpendicular to the direction of ultrasound propagation, resulting in near zero Doppler shift attributable to blood flow. In the case of rotational IVUS, there is an added complication due to the continuous rotation of the transducer, which makes it problematic to collect the multiple echo signals from the same volume of tissue needed to make an accurate estimate of the velocity-induced Doppler shift. Various image correlation methods attempt to overcome the directional limitations of the Doppler method for intravascular motion detection, but are generally inferior to Doppler methods. Moreover, such image correlation techniques are not suitable for rotational IVUS because the rate of decorrelation due to the rotating ultrasound beam is comparable to the rate of decorrelation for the blood flow.
Accordingly, there is a need for apparatuses, systems, and methods that can produce intravascular ultrasound images using a mechanically-scanned ultrasound transducer, and that differentiate between moving blood and stationary tissue within a vessel. The apparatuses, systems, and methods disclosed herein overcome one or more of the deficiencies of the prior art.
Embodiments of the present disclosure describe a mechanically-scanned intravascular ultrasound (IVUS) imaging system that produces a rotational IVUS image with the addition of velocity data encoded as a color overlay on a grey-scale IVUS image to enhance the differentiation between moving blood echoes and stationary tissue echoes.
In one aspect, the present disclosure provides a rotational intravascular ultrasound system for imaging a vessel. The system comprises an ultrasound transducer rotationally disposed within an elongate member, and an actuator coupled to the transducer, the actuator moving the transducer through at least a portion of a revolution. The imaging system includes a control system controlling the emission of a sequence of ultrasound pulses and reception of the associated ultrasound echo signals. The control system processing the ultrasound echo signals to produce a cross-sectional image of the vessel based on both the echo amplitude and the Doppler frequency shift (indicative of the velocity of blood or other tissue within the vessel). In one embodiment, the actuator is coupled to the ultrasound transducer through a flexible drive cable extending substantially the entire length of the elongate member, the actuator continuously rotating the ultrasound transducer generally about a longitudinal axis of the elongate member.
In another aspect, the disclosure provides an ultrasound imaging system having a distal portion insertable into a vessel of a living body. The system comprising an elongate member having a longitudinal axis extending along a distal portion, the elongate member having an ultrasound transducer mounted adjacent to the distal portion such that an ultrasound pulse emitted by the transducer propagates away from the elongate member at a substantially non-perpendicular angle relative to the longitudinal axis. The ultrasound transducer being configured to receive ultrasound echo signals and convey these signals through a plurality of conductors extending between the ultrasound transducer disposed adjacent the distal portion to a connection assembly disposed adjacent an opposite proximal portion of the elongate member. The system further includes an actuator, which may be within the elongate member or external to the body, coupled to the ultrasound transducer. The actuator is configured to move the transducer through a range of positions extending over at least a portion of a revolution. In one embodiment the movement is continuous rotation about the longitudinal axis while in an alternative embodiment, the movement is an oscillatory action over a portion of a revolution. The system further includes a control system coupled to the connection assembly. The control system being configured to control the position of the ultrasound transducer and the timing of the ultrasound pulses emitted by the transducer. The control system receiving the ultrasound echo signals from the ultrasound transducer through the plurality of conductors and processing the echo signals to generate an image of the vessel. In one aspect, the vessel image includes amplitude data represented in grey-scale overlaid with colorized areas representative of the velocity data, derived by the control system from the Doppler frequency shifts detected in the ultrasound echo signals. In an alternative form, the grey-scale amplitude data is altered to reflect changes associated with the determination of the velocity data. In one form, pixels associated with velocities above a threshold value are suppressed or diminished in brightness to enhance the view of relatively stationary features of the blood vessel.
In another aspect, the present invention includes a method of imaging a vessel. The imaging method comprising positioning an elongate member having a distal portion with a longitudinal axis within the vessel, the catheter including an ultrasound transducer movably mounted within the distal portion. The method continues with emitting a sequence of ultrasound pulses from the transducer at a substantially non-perpendicular angle relative to the longitudinal axis while moving the transducer through at least a portion of a revolution with respect to the longitudinal axis. The method includes receiving the corresponding sequence of ultrasound echo signals from vessel features including blood within the vessel, processing the sequence of ultrasound echo signals to generate a single composite amplitude ray associated with a position along the arc, processing the sequence of ultrasound echoes to determine the velocity of vessel structures, and displaying a vessel image combining velocity and amplitude information. In one aspect, the display includes a color encoding of the velocity information combined with a grey-scale representation of the echo amplitude, while in an alternative aspect, pixels associated with velocities above a threshold value are suppressed or diminished in brightness in the displayed image. In a further aspect, the method can include automated determination of the vessel boundaries based on an algorithm utilizing the both velocity and amplitude information. In still a further aspect, the velocity information can be used to quantify blood flow within a vessel.
It is to be understood that both the foregoing general description and the following detailed description are exemplary and explanatory in nature and are intended to provide an understanding of the present disclosure without limiting the scope of the present disclosure. In that regard, additional aspects, features, and advantages of the present disclosure will be apparent to one skilled in the art from the following detailed description.
The accompanying drawings illustrate embodiments of the devices and methods disclosed herein and together with the description, serve to explain the principles of the present disclosure. Throughout this description, like elements, in whatever embodiment described, refer to common elements wherever referred to and referenced by the same reference number. The characteristics, attributes, functions, interrelations ascribed to a particular element in one location apply to those elements when referred to by the same reference number in another location unless specifically stated otherwise.
The figures referenced below are drawn for ease of explanation of the basic teachings of the present disclosure only; the extensions of the figures with respect to number, position, relationship and dimensions of the parts to form the preferred embodiment will be explained or will be within the skill of the art after the following description has been read and understood. Further, the exact dimensions and dimensional proportions to conform to specific force, weight, strength, and similar requirements will likewise be within the skill of the art after the following description has been read and understood.
The following is a brief description of each figure used to describe the present invention, and thus, is being presented for illustrative purposes only and should not be limitative of the scope of the present invention.
a is an illustration of a partially cross-sectional view of the rotational IVUS catheter shown in
b is an illustration of an IVUS grey-scale image according to one embodiment of the present disclosure.
a is an illustration of an IVUS velocity image according to one embodiment of the present disclosure.
b is an illustration of a hybrid color flow IVUS image according to one embodiment of the present disclosure.
For the purposes of promoting an understanding of the principles of the present disclosure, reference will now be made to the embodiments illustrated in the drawings, and specific language will be used to describe the same. It will nevertheless be understood that no limitation of the scope of the disclosure is intended. Any alterations and further modifications to the described devices, instruments, methods, and any further application of the principles of the present disclosure are fully contemplated as would normally occur to one skilled in the art to which the disclosure relates. In particular, it is fully contemplated that the features, components, and/or steps described with respect to one embodiment may be combined with the features, components, and/or steps described with respect to other embodiments of the present disclosure. For simplicity, in some instances the same reference numbers are used throughout the drawings to refer to the same or like parts.
The present disclosure describes apparatuses, systems, and methods for producing Doppler color flow ultrasound images from a mechanically-scanned ultrasound transducer, deployable with a rotational intravascular ultrasound (IVUS) imaging system, to facilitate interpretation of cross-sectional images of a blood vessel of interest and to facilitate the qualitative or quantitative measurement of blood flow in the vessel. In particular, the disclosure describes one embodiment of the apparatuses, systems, and methods that produce a rotational IVUS image with the addition of velocity data encoded as a color overlay on a grey-scale IVUS image to enhance the differentiation between moving blood echoes and stationary tissue echoes.
The catheter 100 is an elongate member shaped and configured for insertion into a lumen of a blood vessel (not shown) such that a longitudinal axis LA of the catheter 100 substantially aligns with a longitudinal axis of the vessel at any given location along the length of the vessel. In that regard, the curved configuration illustrated in
It will be understood that although the imaging system of
In a typical rotational IVUS catheter, a single ultrasound transducer element is mounted near the tip of a flexible driveshaft, which spins inside a plastic sheath inserted into the vessel of interest. The transducer element is oriented such that the ultrasound beam propagates generally perpendicular to the axis of the catheter. As the driveshaft rotates (typically at approximately 30 revolutions per second), the transducer is periodically excited with a high voltage pulse to emit a short burst of ultrasound. The same transducer then receives the returning ultrasound echoes reflected from various tissue structures, and the IVUS imaging system assembles a two dimensional display of the vessel cross-section from a sequence of several hundred of these ultrasound pulse/echo acquisition sequences occurring during a single revolution of the transducer. In a conventional IVUS imaging system, the imaging surface is nominally planar and generally perpendicular to the axis LA and the longitudinal axis of the blood vessel. An alternative configuration for a rotational IVUS catheter uses a rotating acoustic mirror combined with a stationary ultrasound transducer to produce a similar effect. In a still further configuration, an ultrasound transducer, mounted on a shaft extending from a motor or other actuator mounted in the distal portion of the catheter or guidewire may be used to mechanically scan the transducer through continuous rotation or over a portion of revolution. In response to electrical power or control signals, the motor or other actuator operates at predetermined speed to cause the shaft and transducer to rotate at a predetermined speed or oscillate at a predetermined rate. Examples of such systems are shown in U.S. Pat. Nos. 5,375,602 and 7,658,715, and in U.S. Patent Application Publication 2011/0071401, each of which is hereby incorporated by reference in its entirety.
The Doppler-enabled rotational IVUS catheter 100 includes an ultrasound transducer 118 tilted at a modest angle away from a perpendicular to the longitudinal axis LA, and it operates much as a conventional rotational IVUS catheter. As the driveshaft rotates (typically at approximately 30 revolutions per second), the transducer is periodically excited with a high voltage pulse to emit a short burst of ultrasound. During the brief time period following each transmitted pulse from the transducer 118, the echo signals from the surrounding tissue and blood are received by the transducer 118 and detected by the control system 300 by way of the PIM 200. The control system 300 then assembles a two-dimensional ultrasound image of a blood vessel cross-section from these hundreds of pulse/acquisition cycles occurring during a single revolution of the device. One cross-sectional image is produced for every rotation of the transducer 118, so the display is updated at approximately 30 frames per second to create the appearance of continuous real-time intravascular imaging. By virtue of the tilted transducer mount, the Doppler-enabled rotational IVUS catheter produces a shallow conical imaging surface 500 instead of the traditional imaging plane, and it may also acquire the Doppler frequency shift information needed to detect the blood velocity component parallel to the longitudinal axis of the catheter.
During operation of the imaging system 10, the control system 300 cooperates with the PIM 200 to generate the appropriate sequence of ultrasound transmit/receive cycles at each image angle to facilitate the extraction of velocity information from the sequence of echo signals. Specifically, as shown in
The proximal bearing 122 supports the rotational interface 111 of the rotating imaging core 110. In the illustrated embodiment, the proximal bearing 122 includes a port 124 for injecting saline into a lumen 131 of the sheath 120 and a fluid seal (not shown) to prevent the fluid from leaking out of a proximal end 132 of the sheath 120. The telescoping section 123 permits the sheath 120 to be reduced or extended in length, causing the rotating imaging core 110 to be correspondingly advanced or retracted with respect to the window segment 128 of the sheath. This telescoping configuration facilitates longitudinal pullback of the transducer 118 through a segment of vessel that is to be examined by the rotational IVUS imaging system 10.
The proximal 126 and distal 127 portions of the sheath partially or fully encase the rotating imaging core 110. The proximal portion 126 comprises a robust, flexible, cylindrical tube that extends from the telescoping section 123 to the window segment 128. The window segment 128 is structurally contiguous with the proximal shaft 126, but it is formed of different material than the proximal portion 126 of the sheath. The window segment 128 may be composed of a material (or combination of materials) that has an acoustic impedance and a sound speed particularly well-suited for conducting the ultrasound beam from the transducer 118 out into the blood vessel with minimal reflection, attenuation, or beam distortion. The tip assembly 130 extends distally from the window segment 128 and is shaped and configured to engage with a conventional coronary guidewire to enable the IVUS catheter 100 to be easily directed into a vessel of interest and to be easily removed from the guidewire.
The ultrasound transducer in a conventional rotational IVUS catheter is mounted substantially in line with the catheter axis LA so that the ultrasound beam emerges substantially perpendicular to the axis of the catheter. In practice, the transducer is frequently mounted at a slight angle in order to reduce the strength of the echo from the catheter sheath. The echo received by the transducer from the catheter sheath) will be strongest when the reflector is parallel to the transducer face and the echoes from different portions of the reflector arrive back at the transducer in phase with one another. If the transducer surface is tilted at an angle such that there is at least one wavelength of path length difference across the axial length of the transducer, then the echoes from the different portions of the sheath will tend to cancel and the echo will be reduced.
As an example of the degree of transducer tilt preferred for a conventional rotational IVUS catheter, the window or aperture width for a typical rotational IVUS catheter is approximately twelve wavelengths (for example, a 500 μm transducer length and approximately 40 μm wavelength at a 40 MHz transducer center frequency). To introduce one wavelength of round trip path length difference across the aperture would require one-half wavelength of tilt over the same width, or an angle of approximately 1/24 radians (approximately 2.5°). With optimum sheath design, the sheath reflection can be small enough that no transducer tilt is needed. Transducer tilt angles in the range of 0° to 8° are common for conventional rotational IVUS catheters.
As described above, in a conventional rotational IVUS catheter, the ultrasound beam emerges from the transducer substantially perpendicular to the longitudinal axis of the catheter, and as the imaging core rotates, the ultrasound beam sweeps out a substantially planar imaging surface to produce an ultrasound image of a cross-sectional slice through the vessel. In this traditional configuration, there is minimal Doppler frequency shift introduced by the moving blood, since the blood motion is substantially parallel to the axis of the vessel, parallel to the axis of the catheter, and accordingly perpendicular to the direction of propagation of the ultrasound beam. Since there is minimal Doppler frequency shift attributable to blood motion, it is difficult to derive a blood flow velocity estimate by Doppler methods using a conventional rotational IVUS catheter.
However, tilting the transducer mounting angle away from the traditional orientation causes a significant component of the blood velocity to be aligned with the ultrasound propagation direction such that the blood velocity (generally parallel to the catheter longitudinal axis LA) can be detected by measuring a Doppler frequency shift. For example, in
The tilted transducer orientation required to produce a significant Doppler frequency shift in the ultrasound echoes from the moving blood transforms the traditional imaging plane into a shallow conical imaging surface. But with a modest transducer tilt, only slight geometric distortion is introduced in projecting this conical imaging surface onto a planar display, and this distortion does not significantly impair the physician's ability to interpret the image. There are two competing considerations for choosing the transducer tilt angle θ for the catheter 100: (1) the larger the tilt angle, the greater (and more readily detectable) will be the Doppler component in the ultrasound echo, and (2) the larger the tilt angle, the greater will be the geometric distortion when a conical imaging surface is projected onto a planar display.
The Doppler shift measured by an ultrasound system is proportional to the cosine of the angle between the direction of the motion and the direction of propagation of the ultrasound beam. In the idealized circumstance, where the longitudinal axis LA of the catheter is aligned with the longitudinal axis of the vessel, and where the velocity of blood flow is parallel to the longitudinal axis of the vessel as well, the angle between the direction of the blood flow and the direction of the ultrasound beam is the complement of the transducer tilt angle θ. Accordingly, the Doppler shift is proportional to the sine of the transducer tilt angle θ. For a zero tilt angle, there is no significant Doppler shift, and the velocity information cannot be obtained from traditional Doppler signal processing. The theoretical maximum Doppler shift would be obtained with a transducer tilt angle of 90°, but that would preclude the possibility for IVUS imaging since the ultrasound beam would then be aligned with the axis of rotation (axis LA). At a modest tilt angle of 30°, the Doppler shift is 50% of the theoretical maximum, and a reasonable IVUS image from the shallow conical imaging surface 500 (shown in
For the intracoronary IVUS application, the Doppler velocity data is important for its role in helping to differentiate blood from tissue, hence the importance of distinguishing the Doppler shift of fast moving blood from the Doppler shift of slow moving tissue. In color flow imaging applications throughout most of the body (e.g., liver, carotid, or peripheral artery), the tissue motion is negligible, so the velocity threshold for classification of an echo as a moving blood echo can be very low. However in the case of coronary imaging, the tissue motion can be quite prominent, and it is more difficult to reliably distinguish tissue motion from blood flow. In general, it is helpful for distinguishing the Doppler shift of fast moving blood from the Doppler shift of slow moving tissue to use a larger transducer tilt angle to increase the Doppler shift for the blood with its predominantly axial velocity, while having little effect on the Doppler shift of the tissue with its predominantly radial velocity.
Although the motion of the heart muscle is quite rapid during early systole when the ventricles contract, the IVUS catheter tends to move with the heart by virtue of its capture within the coronary artery. Thus, the relative motion between the catheter and the surrounding tissue is usually significantly less than the absolute motion of the heart. An example of a fast movement of the IVUS catheter with respect to the heart would be for the catheter to shift one vessel diameter (˜3 mm) during the approximately 100 msec that constitutes the early portion of systole. The corresponding relative tissue velocity in this case would be ˜3 cm/sec. Throughout most of the cardiac cycle, and in the majority of locations throughout the epicardial arterial tree, the actual tissue velocity will be much less than this estimate. In particular, in the coronary arteries, blood flow is most significant (typically in the range of 10 cm/sec to 100 cm/sec) during diastole, the portion of the cardiac cycle when the heart motion is at its minimum (as the heart muscle gradually relaxes). Accordingly, in some embodiments, it is desirable to gate the Doppler color flow imaging with the ECG to capture blood flow measurements only during diastole, when the blood flow is maximum, and the heart motion (and relative tissue velocities) are at a minimum.
a illustrates the distal portion 127 of the Doppler-enabled rotational IVUS catheter 100 positioned within a vessel 600, which includes a lesion 601 attached to the vessel wall 601 within a lumen 602. The catheter 100 includes the transducer 118 mounted at a significant tilt angle within the housing 116. In
The choice of the transducer tilt angle for Doppler color flow rotational IVUS imaging should consider the robustness of the Doppler velocity measurement in the face of misalignment between the catheter axis and the axis of the blood vessel, as well as the ability to distinguish the Doppler shift of fast moving blood from the Doppler shift of slow moving tissue. In the course of normal clinical use, there may be a misalignment between the axis LA of the catheter 100 and the axis VA of the vessel (and the direction of blood flow). If that misalignment is comparable to the transducer tilt angle θ, then the Doppler shift in a portion of the vessel might be reduced to zero where the catheter misalignment cancels the transducer tilt angle θ. However, if the transducer tilt angle θ is significantly greater than the typical range for catheter misalignment, then the system 10 will retain a robust capability for estimating blood velocity across the entire vessel lumen. Because the human anatomy may include significant tortuosity in regions where IVUS imaging is commonly used (e.g., but not by way of limitation, the coronary arteries), it is difficult to predict the largest misalignment that can exist between the vessel axis VA and the catheter axis LA. In one example, a large misalignment that might be found in clinical practice would be the equivalent of a 1 millimeter (mm) diameter catheter traversing a 3 mm in diameter vessel lumen over a 10 mm length of vessel, corresponding to a likely maximum misalignment angle of approximately 12°. Over much of the epicardial arterial tree, however, the actual misalignment angle would be less than this maximum likely value. Nevertheless, it would be helpful for maintaining a robust Doppler signal if the transducer tilt angle θ was greater than 12°. Based on this consideration, the transducer tilt angle θ should preferably be greater than 15° to allow a small margin above the 12° maximum likely misalignment angle predicted above. More preferably, the transducer tilt angle θ should be approximately 20° to provide a greater margin of tolerance for catheter to vessel misalignment.
As illustrated in
In most respects, the grey-scale image 700 is substantially the same as that produced by the traditional IVUS method using a non-tilted or marginally tilted transducer, except for a slight geometric distortion arising from projecting the conical imaging surface onto a planar display. Since the ultrasound image produced from the conical surface 500 is typically displayed on a planar video monitor, there is a geometric distortion introduced in the conical to planar transformation. The degree of distortion can be quantified by a figure of merit which represents the discrepancy between radial and tangential distance measurements on the distorted planar display. The distortion figure of merit can be calculated as one minus the cosine of the tilt angle. A zero tilt angle produces a planar imaging surface with no distortion, while a tilt angle of 20° produces 6% distortion. A modest degree of distortion will not interfere with the qualitative interpretation of the image which requires the identification of the inner and outer borders of the vessel wall structures and assessment of the general character of the echoes from lesions within the vessel wall. Any quantitative measurements, such as lumen diameter or plaque cross-sectional area to be made from the distorted planar display can be easily corrected by applying the appropriate geometric formula to remove the conical distortion from the calculation. For the preferred range of tilt angles θ from 10° to 30°, the visual distortion ranges from 1.5% to 13%, while for the more preferred range of tilt angles θ from 15° to 25°, the visual distortion ranges from 3% to 9%.
Therefore, the choice of the transducer tilt angle θ for the Doppler Color Flow Rotational IVUS imaging system 10 may involve consideration of the following factors: (1) the robustness of the Doppler velocity measurement in the face of misalignment between the catheter axis LA and the longitudinal axis of the blood vessel, (2) the ability to differentiate Doppler shift of fast moving blood from the Doppler shift of slow moving tissue, and (3) the degree of distortion of the IVUS image when the conical image surface is projected onto a planar display (with a view to minimize such distortion). A compromise tilt angle can be chosen wherein the image distortion due to projection of the conical imaging surface onto the planar display is acceptably small, while the Doppler shift is large enough that it provides a robust blood velocity measurement, tolerant of small misalignments between the catheter axis LA and vessel axis VA, and sufficient to differentiate fast moving blood from stationary or slow moving tissue.
In the grey-scale image 700 depicted in
It is important to note that noninvasive color flow imaging systems cannot take advantage of high ultrasound frequencies, such as 40 MHz, due to the frequency-dependent attenuation of the ultrasound in tissue which severely limits the penetration depth. Furthermore, noninvasive color flow imaging systems cannot take advantage of high ultrasound frequencies, such as 40 MHz, since the high ultrasound frequency results in large Doppler frequency shifts that necessitate a high pulse repetition frequency and short period between successive ultrasound pulses, once again limiting the usable penetration depth. However, for rotational IVUS imaging, the shallow penetration depth (approximately 5 mm) permits the use of a high pulse repetition frequency adequate to capture the maximum velocity likely to be encountered in a physiological environment. For rotational IVUS imaging, the required penetration depth is only about 5 mm, and the attenuation in blood, even at a 40 MHz ultrasound frequency, is low enough to allow adequate signal to noise ratio for such a shallow penetration depth.
To improve the diagnostic value of the IVUS image 700, the Doppler enabled IVUS catheter 100 and the Doppler-capable IVUS control system 300 utilize a separate signal processing path operating in parallel to the standard imaging path to provide velocity information for the various components within the vessel 600. While the standard image processing algorithm translates the amplitude of the echo signal into grey-scale brightness on the display image 700, a parallel signal processing path extracts a velocity estimate for every pixel of the display image 700 from the information contained in the Doppler frequency shifts of the echo signals.
a depicts the image that would be obtained if the imaging system 10 was programmed to display an image of the velocity estimates extracted from the ultrasound echoes received by the transducer 118 instead an image of the amplitudes of the ultrasound echoes received by the transducer 118. The lesion representation 711 and the vessel wall representation 712 of the velocity image 710 would indicate low velocities for the relatively static lesion 601 and vessel wall tissue 602, respectively, while the relatively fast-moving blood 603 within the vessel lumen 602 would be prominently highlighted by the blood velocity representation 713.
In practice, the separate grey-scale IVUS image 700 and velocity image 710 may be difficult to interpret, but a synergistic image may be produced by combining the echo amplitude and velocity information together in a hybrid Doppler color flow image 720, as shown in
The encoder 210, which is coupled to the motor 260 that drives the rotating imaging core 110 (not shown), generates pulses at regular intervals throughout the rotation of the imaging core 110 (i.e., typically 512 pulses per revolution). Instead of each encoder pulse generating a single trigger pulse as in a traditional IVUS system, each encoder pulse triggers a sequence of multiple transmit triggers (such as, by way of not-limiting example, 2 to 16 triggers) via the sequencer 310. Each transmit trigger initiates a pulse from the transmitter 220 which passes through the ultrasound T/R switch 230 to reach the ultrasound transducer 118 by way of the rotary transformer 240. The T/R switch 230 protects the delicate circuitry of the amplifier 250 from the high-voltage transmit pulses while permitting the low amplitude echo signals to pass freely to the amplifier input. The rotary transformer 240 allows the transmit pulses and the echo signals to pass between the stationary elements of the PIM 200 and the rotating imaging core 110 that carries the transducer 118.
For rotational IVUS imaging at 30 frames per second, the orientation of the transducer 118 is constantly changing, making it difficult to collect the repeated measurements from a single direction that are preferred for creating a Doppler velocity estimate. However, given the high speed of sound propagation through tissue compared to the scanning rate, together with the short penetration depth for IVUS imaging of approximately 5 mm, there is sufficient time to include a sequence of several ultrasound transmit/receive cycles for each imaging angle within the IVUS display. The duration of this rapid sequence of pulses can be sufficiently short that several successive transmit/receive cycles will capture echoes from substantially the same tissue/blood volume such that the Doppler frequency shift can be extracted from the echoes received during this sequence of transmit/receive cycles.
For a period of time after each transmit pulse (e.g., typically approximately 10 μsec), the amplifier 250 receives a low level echo from the transducer 118 and applies the appropriate time dependent gain to produce an amplified echo signal. The amplitude of the echo signal versus time (relative to the transmit pulse) is representative of the reflectivity of the (reflecting) tissue as a function of distance from the transducer 118. In addition, information regarding the motion of the tissue is encoded in the small changes, particularly the phase shift, between one echo signal and the next within a sequence. The amplifier 250 transmits the amplified echo signals to signal processing hardware in the control system 300 for processing.
For Doppler processing, it is convenient to transform the radio frequency (RF) echo waveform into a baseband representation according to methods well-known to those skilled in the art, wherein the transducer center frequency is shifted down to DC, and the echo signal is represented as digitized samples of a complex modulation waveform comprised of in-phase (I) and quadrature (Q) components. The demodulator/digitizer 330 transforms the amplified echo signal from the amplifier 250 into a baseband representation of the signal comprising digitized samples of the I and Q components of the complex modulation waveform. This function can be performed in the analog domain using a pair of mixers, followed by a pair of analog-to-digital converters to provide digital samples of the I and Q components. Alternatively, the demodulation step can be performed in the digital domain by direct sampling of the RF echo waveform with a high speed analog-to-digital converter, followed by digital filtering to produce digital samples of the I and Q components of the complex modulation waveform.
The grey-scale analyzer 350 and the velocity computer 360 process the multiple echo signals from a single sequence as a group, and use the information contained in the sequence of echo signals to (1) detect the echo amplitude as a function of depth to generate a single ray or radial line of the grey-scale image (commonly referred to as an A-line), and (2) to calculate the Doppler-derived velocity for each position along that ray, respectively. Specifically, the grey-scale analyzer 350 uses the information contained in the sequence of echo signals to detect the echo amplitude as a function of depth to generate a single A-line of the grey-scale image with a low noise level and wide dynamic range, while the velocity computer 360 calculates the Doppler-derived velocity estimate for each position along that A-line from the small phase changes from one echo signal to the next within a single sequence. In theory, the velocity data could be used to produce a velocity image 710 of a cross-section through the vessel 600, but in practice, it is convenient to combine the echo amplitude data with the velocity data to produce a hybrid color flow image 720 combining the grey-scale IVUS image with velocity information encoded as shades of red and blue (for antegrade and retrograde flow), and with stationary and slow moving tissues displayed in shades of grey. As the rotating imaging core 110 and transducer 118 rotate within the sheath 120, the IVUS imaging system 10 builds a complete cross-sectional image 700 (commonly referred to as a B-scan image) of the artery 600 from the succession of A-lines from the grey-scale analyzer 250. The amplitude and velocity data are also combined into color-coded A-lines and scan converted in the display processor 370 for display as the hybrid color flow image 720 on the color monitor 400.
Following each of the high voltage transmit pulses 221 in the sequence, the amplifier 250 (not shown in
Each sequence of transmit trigger pulses 221 initiates the acquisition of a sequence of echo signals 251, and although the echo amplitude can be easily derived from a single echo signal from within that sequence, it is more difficult to extract the velocity information from just one echo signal. However, the velocity can be estimated by analyzing how the echo signal changes from one echo signal to the next within a sequence. For stationary or slow-moving tissue, the echo signal 251 changes very little from one echo signal 251 to the next within a sequence, since the pulses within a sequence are so closely spaced in time that there is little tissue motion over that short interval. Furthermore, the rotation of the transducer 118 over this short interval is small enough compared to the dimensions of the ultrasound beam 121 (shown in
It is advantageous under circumstances described below, to arrange the sequence of acquisitions such that the first amplified echo signal 251a in the sequence is acquired with a lower analog gain setting in the amplifier 250 compared to the gain used for the second 251b through the Nth 251n acquisitions. In this case, the first acquisition, acquired with a low analog gain setting, accurately capture the echoes from strong reflectors such as calcified plaque or metal stent struts without the distortion that arises when the amplifier is driven into saturation by a strong echo signal acquired with a high amplifier gain setting. The subsequent acquisitions within the sequence are collected using the higher gain setting to faithfully capture the low amplitude echo signals from blood, soft plaques, and other low reflectivity tissues. The first acquisition, captured using a low analog gain setting, is not particularly useful for Doppler velocity estimation since the weak echoes from blood and soft tissue are likely to be partially obscured by the quantization noise from the analog-to-digital converter.
However, this low-gain acquisition is useful for generating wide dynamic range grey-scale image data, and it also serves to initialize the reverberations within the transducer, the catheter, and the medium, thereby reducing the Doppler artifact that can arise from such reverberations or the absence thereof. Reverberations arise from acoustic signals originating from the transmit pulse(s) prior to the most recent one. For grey-scale IVUS imaging, these reverberations are generally of little consequence since they are greatly diminished over the time between acquisitions and only generate a small perturbation in the amplitude of the echo signal. However, these small perturbations can give rise to phase artifacts that may be significant for low level echoes that are of great interest with respect to blood flow, causing artifacts in the blood velocity estimates. Excluding the first echo acquisition from the Doppler velocity computation ensures that each of the subsequent acquisitions used for the velocity algorithm includes a similar pattern of reverberation that will tend to introduce zero Doppler shift, at least to the extent that the reverberation is consistent from pulse to pulse. In this case, only the transmit pulse from the first acquisition affects the reverberations in subsequent acquisitions, and the amplifier gain for the first acquisition is irrelevant to the reverberations. Acquiring one low-gain acquisition as the first in a sequence facilitates wide dynamic range grey-scale IVUS imaging without compromising the velocity measurement potential of the system.
Although the amplitude signal can be derived from just a single echo signal 251 chosen from the sequence of echo signals acquired, it is advantageous to construct a composite echo signal 258 by signal averaging or a similar technique applied to the entire sequence of echo waveforms to provide an improved signal-to-noise ratio. The envelope derived from this composite echo signal 258 exhibits improved signal-to-noise ratio compared to that derived from just a single echo waveform 251. Signal processing across the ensemble of echo waveforms within a sequence is facilitated by processing the data in the digital domain where the multiple waveforms can be readily stored, retrieved, and processed.
The amplitude and velocity information can be independently presented as separate images of echo amplitude and Doppler velocity over the field of view, but it is preferred to combine these two sets of information into the hybrid color flow image 720 (shown in
In particular, the combined amplitude and velocity data may be utilized, for example but without limitation, by the imaging system 10 to enhance suppression of blood echoes, luminal border and cross-sectional area detection, quantitative blood flow measurements, and thrombus detection. The imaging system may enhance the contrast between the blood echoes and the vessel wall by using color to highlight the moving blood or by simply deemphasizing the moving blood by reducing the brightness of the fast-moving blood elements of the image. For example, to suppress blood echoes from the final image 720, the imaging system 10 diminishes the brightness of the echoes that contain a significant velocity component so that the vessel lumen 602 (shown in
To enhance luminal border detection, the imaging system 10 uses the velocity data to improve the algorithm for automatic (computer-based) detection of the lumen border and the lumen cross-sectional area. These can be determined by manually tracing the lumen borders or by placing markers at intervals around the vessel border, but it is highly advantageous if those measurements are automatically provided by a computer algorithm that identifies the border on its own. Some IVUS imaging systems include such automated measurement algorithms, but these frequently require human intervention to improve the quality of the border detection. Such a system is described in U.S. Pat. No. 6,945,938 hereby incorporated by referenced herein in its entirety. Providing velocity information to the automatic border detection algorithm can improve the quality of the result and reduce the need for tedious manual tracing of the borders.
As described above, in the traditional IVUS imaging system, the differentiation between moving blood and stationary thrombus may be very subtle. There may a slight difference in the temporal appearance of the speckle pattern, but there is often very little difference in the echogenicity of blood versus thrombus (particularly fresh thrombus). However, velocity provides a very strong signature to differentiate moving blood from stationary thrombus, and Doppler color flow imaging by the imaging system 10 utilizing blood velocity data may greatly improve the detection of thrombus.
To estimate quantitative blood flow, the imaging system 10 numerically integrates the blood velocity data over the cross-sectional area of the vessel lumen 602 to provide a quantitative measurement of volumetric blood flow within the artery 600 (shown in
In addition, the imaging system 10 may extend the dynamic range of the grey-scale IVUS image by using the same pulse sequence used to provide the information needed for measuring Doppler frequency shift. In a healthy artery, a clear IVUS image may be obtained with relatively modest dynamic range. However, in the pathological conditions of greatest interest to the physician, a wider dynamic range is frequently needed. In diseased arteries, there are often deposits of calcium within a plaque, and these calcium deposits produce strong echoes that may tend to dominate the image and obscure the nearby low level echoes. Similarly, IVUS is frequently used for imaging arteries where metal stents have been placed, and the metal stent struts produce strong echoes which tend to obscure the vessel wall image behind the struts. The wide dynamic range feature offers a significant benefit under these pathological conditions by reducing the noise level in the image, enhancing the visibility of weak reflectors, and reducing the image artifacts arising from strong reflectors such as calcium deposits or stent struts. Extending the dynamic range of the IVUS signals can make it easier to detect the weak echoes from soft tissue while simultaneously detecting the strong echoes from metal stent struts or calcified plaques embedded in the vessel wall. Extending the dynamic range of the IVUS signals is discussed in more detail below with respect to
The amplitude detection circuitry could be as simple as calculating the envelope of a single A-line at a time and forwarding that envelope data on to the display processor. But because the phase detection circuitry used for the Doppler velocity computation preferably uses a sequence of echo signal acquisitions covering substantially the same volume of tissue, it is advantageous to use the same multiple acquisitions to improve the signal-to-noise ratio and dynamic range available for the grey-scale display. IVUS imaging produces a wide dynamic range of echoes, ranging from the weak echoes from blood or soft tissue to the strong echoes from calcified plaque or metal stent struts. One method to expand the dynamic range available for display is to increase the signal-to-noise ratio by averaging multiple echo signals together to reinforce the coherent echoes from the tissue compared to the incoherent noise present in the individual echo signals. Another method for expanding the dynamic range is to acquire echo signals with different analog gain settings, and to then splice together the low gain samples from strong reflectors with high gain samples from weak reflectors, including digital compensation for the different analog gains. Additional noise reduction may be achieved by applying a nonlinear algorithm across the sequence of echo signal acquisitions to reject isolated impulsive noise spikes.
In the embodiment shown in
A separate low gain line buffer 352 stores an echo signal, which is also in baseband I/Q format, acquired using a lower amplifier gain compared to the gain used for acquiring the other echo signals in the acquisition sequence. The echo signal acquired with lower amplifier gain captures a low distortion representation of strong echoes from calcified plaques or stent struts that might saturate the amplifiers and acquisition circuitry when a higher amplifier gain is used for the other echo signal acquisitions in the sequence. In a typical example, the low gain amplifier setting would be −12 dB relative to the high gain setting (a factor of one-quarter), with correspondingly reduced distortion from amplifier or signal acquisition stage saturation due to strong echoes. If the signal amplitude is less than one-quarter of the maximum low gain value, then the low-noise composite echo signal derived from the multiple high gain acquisitions provides the best representation of the echo signal, but for any echo amplitude above that threshold, the composite echo signal is likely to be saturated, and the low gain acquisition samples should be used instead.
The amplitude of each of these buffered signals (the composite echo signal and the low gain echo signal) can be calculated using a variety of methods known to those skilled in the art, but the method described herein is well-suited for implementation in an FPGA, since it requires relatively simple logic and small lookup tables stored in the FPGA memory, and it operates directly on echo signal waveforms captured in baseband I/Q format. Essentially identical circuitry is shown for envelope detection of the accumulated composite echo signal and the buffered low gain echo signal, but this circuitry could be shared by time-multiplexing the signals from those two separate signal paths through one set of envelope detection circuitry in order to reduce the required FPGA resources.
The first step of the envelope detection process is to convert the linear representation of the baseband I and Q values into a more compact representation (requiring fewer bits) in order to simplify the subsequent calculations and reduce the size of the required lookup tables. A block priority encoder 353 or 354 converts an I/Q sample pair into a floating point format, using a shared exponent for both samples. The block priority encoder determines which of the I and Q samples is the largest, preserving the most significant bits of that value as the mantissa and using simple logic to calculate the associated exponent (power of 2) required for the floating point representation of the original sample. The smaller of the two samples is shifted by the number of bit positions specified by the shared exponent, and the high order bits become the mantissa for the smaller of the two samples. In this illustrative example, the 12- or 14-bit I and Q samples (11 or 13 bits plus sign) are converted to floating point representations, each with a sign (not needed for amplitude calculation) plus 7-bit mantissa, and with a shared 4-bit exponent.
At this point, the benefit of the block priority encoder becomes apparent in the small size of the magnitude lookup table (LUT) 355 or 356 required for calculating the magnitude of the I/Q sample pair as the square-root of the sum of the squares of the two values. In this illustrative example, two 7-bit mantissas necessitate a modest-sized 8-kbyte (213×1 byte) LUT to provide the square-root of the sum of the squares of the two mantissas. Only 13 bits are required to address the magnitude LUT because the most significant bit of the larger mantissa is omitted since it is always a 1 and the sign bits are ignored since they don't affect the magnitude calculation. Without this block floating point approach or another efficient method, an impractically large 128-megabyte (226×2-byte) LUT would be needed to calculate the square-root of the sum of squares of a pair of 14-bit composite I and Q samples (sign bits are ignored for the magnitude calculation). Once the mantissa of the magnitude is provided by the magnitude LUT 355 or 356, the exponent (representing a common factor shared by both I and Q values) is applied through the shifter 357 or 358 to reverse the floating point conversion and restore a linear representation of the magnitude of the envelope for the corresponding signal (either the composite echo signal or the low gain echo signal).
Finally, to provide wide dynamic range grey-scale image information to the display processor 370, the dynamic range combiner 359 splices together the low-amplitude, low-noise envelope data from the composite echo signal with the high-amplitude, low-distortion envelope data from the low-gain echo signal. The result is exceptionally wide dynamic range for the grey-scale image data, facilitating the display of weak tissue and blood echoes visible above a very low noise floor, while strong echoes from stent struts or calcified plaques appear without saturation. The dynamic range combiner may be as simple as a comparator that switches between either of the two signal sources based on the strength of the echo at that particular point. If the signal is weak, then the low-noise composite echo signal is used as the source for the grey-scale information, and if the signal is strong enough to approach the full-scale maximum for the composite echo signal, the dynamic range combiner switches seamlessly to the low-distortion echo signal acquired using reduced amplifier gain.
More advanced schemes for combining these two echo signals may include a transition zone between the two signal sources, wherein low amplitude echoes are derived solely from the low-noise composite echo signal, strong echoes are derived from the low-distortion, low-gain echo signal, and over a broad intermediate range, the grey-scale information is provided by interpolation between the two signal sources using amplitude-dependent coefficients. For example, if the low-gain sample is less than one-sixteenth of the full-scale amplitude, then the low-noise composite sample should be well below its full scale limit of one-quarter of the low-gain maximum, and the composite sample with its low noise level is used alone. For samples greater than one-quarter of the full-scale amplitude of the low-gain acquisition, the composite sample is likely to be beyond its full scale maximum value and the low-gain sample may be used alone to provide a low-distortion sample of the high amplitude echo amplitude. For the samples in the transition range between one-sixteenth and one-quarter of the full-scale amplitude of the low-gain acquisition, a weighted average of the composite and low-gain samples may be used, with amplitude-dependent coefficients gradually phasing in one or the other source based on the echo amplitude at a particular point. At the output of the dynamic range combiner, the amplitude signal covers a very wide dynamic range with reduced noise compared to the noise level expected from a single acquisition. The wide dynamic range A-line is encoded as a 16-bit integer capable of encoding a dynamic range on the order of 96 dB.
In the embodiment shown in
A variety of methods are known to those skilled in the art for extracting a Doppler-derived velocity estimate from a sequence of echo signals. These methods have been applied extensively to noninvasive Doppler color flow ultrasound imaging systems, but heretofore this technology had not been thought to be applicable to rotational IVUS imaging systems for the reasons discussed above. The phase detection circuitry is designed to extract a Doppler velocity estimate for each radial position along an A-line of the image, from the sequence of echo signal acquisitions that is obtained from that angular position. The movement of blood or tissue at a particular depth along a ray of the image is encoded in the sequence of echo signals as a rate of change in phase of the echo signal at that radial position. The phase detection circuitry is designed to determine the phase change at every sample depth for each echo signal acquisition (ignoring the low-gain acquisition if that feature is implemented), to calculate the change in phase from one echo signal to the next within a sequence, to accumulate the change in phase over the series of echo signal acquisitions within the sequence, and to estimate the tissue or blood flow velocity from the rate of change in phase according to the Doppler equation.
The first step of the phase detection process is to convert the linear representation of the baseband I and Q values into a more compact representation (requiring fewer bits) in order to simplify the subsequent calculations and reduce the size of the required lookup tables. A block priority encoder 361 converts an I/Q sample pair into a floating point format, using a shared exponent for both samples. The block priority encoder determines which of the I and Q samples is the largest, and it preserves the most significant bits of that value as the larger mantissa while generating an I>Q comparison flag according to which of the two values is larger. The priority encoder uses simple logic to calculate the associated exponent (power of 2) required for the floating point representation of the larger sample, and the smaller of the two samples is shifted the number of bit positions specified by the shared exponent, with the high order bits then becoming the mantissa for that smaller of the two samples. In this illustrative example, the 12-bit I and Q samples (11-bits plus sign) are converted to floating point representations, each with a sign plus 7-bit mantissa, and with a shared 4-bit exponent (which is not needed for phase detection).
At this point, the benefit of the block priority encoder becomes apparent in the small size of the phase LUT 362 required for calculating the phase of the I/Q sample pair as the arctangent of Q/I. In this illustrative example, two 7-bit mantissas necessitate a modest-sized 8-kbyte (2{circumflex over (0)}×1 byte) LUT to provide the arctangent computation over one octant. Only 13 bits are required to address the magnitude LUT because the most significant bit of the larger mantissa is omitted as it is always 1, and the exponent is ignored since it is common to both I and Q and does not affect the ratio of Q/I. The one-octant phase angle from the phase LUT is expanded to a full four-quadrant phase angle by octant logic 363 by utilizing the two sign bits to identify one of the four quadrants and utilizing the I>Q comparison bit from the block priority encoder to identify which of the I/Q sample pair is larger. Without this block floating point approach or other efficient method, an impractically large 16-megabyte (224×1-byte) LUT would be needed to compute the arctangent of Q/I for a pair of 12-bit I and Q samples.
The Doppler velocity is estimated for each radial position along an A-line by measuring the rate of change of phase at that point. This may be accomplished by buffering the phase from one line of acquisition in the one-line phase buffer 364, and then subtracting this buffered phase data, sample by sample, from the next line of phase data as it is loaded into the one-line phase buffer replacing the prior line of buffered phase data. The subtraction operation used to calculate the phase change is performed (modulo 2π) in the delta phase block 365. The phase change as calculated covers a range of 2π, but this phase change must be properly interpreted to distinguish positive velocity (antegrade flow) from negative velocity (retrograde flow). If there is no a priori knowledge regarding the direction of flow, then an unbiased approach would be to interpret the phase change values to represent a range from −π to +π. If there is some a priori knowledge that flow is strictly one direction, then a biased approach may be to assign all phase changes to be either positive or negative, according to the assumptions about the directional nature of the blood flow. There can be an intermediate interpretation as well, for example if the velocity is predominantly in one direction, the phase change could be interpreted to represent a range from −π/2 to +3π/2.
The range of velocities that can be unambiguously interpreted is limited by the requirement that the Doppler frequency shift must not exceed one-half of the pulse repetition frequency between successive transmit pulses within a sequence (for the unbiased case described previously). In a typical example, the transmit pulses within a sequence are spaced 10 μsec apart, corresponding to a 100 kHz pulse repetition frequency, and the corresponding maximum Doppler frequency is 50 kHz. The Doppler equation can be used to translate this maximum Doppler frequency into a maximum detectable blood velocity. The Doppler equation states:
In this equation, c represents the speed of sound in blood, 1540 m/sec, α is the angle between the direction of blood flow and the direction of ultrasound propagation, that is, the complement of the transducer tilt angle (for example, α would be 70° for a typical transducer tilt angle of 20° anticipated for the Doppler-enabled rotational IVUS catheter). For rotational IVUS catheters, the typical ultrasound center frequency is 40 MHz, and for the maximum Doppler frequency of 50 kHz, the calculated maximum measurable blood velocity is ±2.80 m/sec. This range of velocities covers most clinical conditions where the device is likely to be used, with a blood flow velocity generally less than 1.0 m/sec in a relatively unobstructed coronary artery under resting conditions. This range of velocities may be extended by implementing a biased delta phase detector as described previously.
The delta phase block 365 calculates the difference in phase between corresponding samples from two successive acquisitions of phase data, and it applies the selected interpretation of phase over a 2π range. The cumulative phase change is then calculated over a sequence of acquisitions to provide a robust estimate of the rate of change of phase at each point along the corresponding A-line. This rate of phase change can be interpreted as a velocity estimate by applying a constant factor derived from the Doppler equation according to methods known to those skilled in the art.
After a complete acquisition sequence is processed, the output of the velocity computer 360 is a single line of velocity data corresponding to the single A-line of grey-scale amplitude data provided by the grey-scale analyzer 350. As subsequent acquisition sequences are processed, additional lines of grey-scale and velocity data are produced and these lines of data are used to paint a complete tomographic image of the vessel, including color encoded velocity information to aid in the interpretation of the image.
Referring back to
In some embodiments, a negative threshold and a positive threshold may be used, wherein any velocity below the negative threshold is assumed to retrograde flow, any velocity above the positive threshold is assume to be antegrade flow, and any velocity between the positive and negative thresholds is assumed to be stationary or slow-moving tissue. This velocity threshold scheme can be used to generate a simple, three level color mask, with blue tint applied to the grey-scale value for any velocity below the negative threshold, red tint for any velocity above the positive threshold, and no tint (grey) for any velocity values between these threshold values representing stationary or slow-moving tissue. In other embodiments, the color flow imaging may use the mask to define the vessel boundaries and support border detection, virtual histology, or a de-speckling algorithm to more clearly distinguish the blood from the stationary tissue.
In some embodiments, a more elaborate scheme may be used with shades of red through yellow encoding positive velocities, with shades of blue through green encoding a range of negative velocities, and with stationary and slow moving tissue receiving a neutral (grey) tint.
Another option might be to forego the color display entirely, and simply use the velocity information to identify moving blood, and then to suppress the grey-scale brightness of the blood speckle to more clearly differentiate the moving blood from the stationary or slow-moving vessel wall. Advanced algorithms might even integrate the velocity map over the cross-section of the artery to provide a quantitative measurement of volumetric flow in the artery.
In some embodiments, the velocity threshold may be chosen to separate moving blood with axial velocities in the 10 to 200 centimeters per second (cm/sec) range from moving tissue with typical velocities on the order of 3 cm/sec or less. The Doppler component will be only 30% of the axial velocity due to the transducer tilt angle, while vessel wall motion is likely to be in the direction of the ultrasound beam where it will cause a maximum clutter signal.
It is important to note that while the apparent tissue motion due to rotation of the transducer in a rotational IVUS catheter creates an obstacle to use of image correlation methods for motion detection, this apparent velocity does not produce a significant Doppler frequency shift because the apparent motion is in a tangential direction, perpendicular to the direction of propagation of the ultrasound beam. Therefore, the transducer rotation does not present a fundamental impediment to Doppler-based velocity measurement.
For application in coronary IVUS imaging, the Doppler velocity data is important for its role in helping to differentiate blood from tissue. The anatomy and physiology of the coronary arteries creates unique blood flow characteristics. In the coronary arteries, blood flow occurs predominantly during the diastolic phase of the cardiac cycle, during which tissue motion is at a minimum because the heart muscle is relaxed. In some instances, early diastole is a preferred phase of the cardiac cycle in which to use the blood velocity to assist with border detection. During early diastole, the velocity information provides maximum differentiation between the stationary tissue and the moving blood, since blood velocity is at its maximum and the heart motion is minimal.
Diastole is also a preferred time for measuring the artery cross-sectional area and diameter, while the distortion of the lumen due to compression of the heart muscle is minimized. In general, it is advantageous to use electrocardiogram (ECG) gating to identify the diastolic phase, and to select diastolic frames of the Doppler color flow IVUS image for detailed quantitative analysis.
In the peripheral arteries, however, where flow occurs predominantly during systole and where tissue motion is less significant, systolic gated frames may be more appropriate for detailed quantitative analysis.
In color flow imaging applications throughout most parts of the body (e.g., but not by way of limitation, hepatic, carotid, or peripheral artery), the tissue motion is negligible, so the velocity threshold for classification of an echo as a moving blood echo can be very low. In the case of coronary imaging, however, the tissue motion can be quite prominent because the coronary vessels overlie the heart muscle, thereby making it more difficult to distinguish tissue motion from blood flow. Although the motion of the heart muscle is quite rapid during early systole when the ventricles rapidly contract, the IVUS catheter 100 tends to move with the heart by virtue of its position within the coronary artery. Thus, the relative motion between the catheter and the surrounding tissue is usually significantly less than the absolute motion of the heart.
An example of a fast movement of the IVUS catheter 100 with respect to the heart would be for the catheter to shift one vessel diameter (approximately 3 millimeters) during the approximately 100 milliseconds that constitutes the early portion of systole. The corresponding relative tissue velocity in this case would be approximately 3 centimeters per second. Throughout most of the cardiac cycle, and in the majority of locations throughout the epicardial arterial tree, the actual tissue velocity will be much less than this maximum likely tissue velocity estimate. Thus, for this additional reason, it may be beneficial to obtain a diastolic gated flow measurement to more accurately obtain blood flow velocities.
In addition, the principles of the above described imaging system and methods can be applied to imaging systems based on other types of waves, such as electromagnetic waves (light, or radio waves), wherein the waves might be emitted/received by an angled emitter/receiver or deflected by an angled mirror, such that the waves propagate at an angle substantially tilted away from a perpendicular to the axis of the catheter. More specifically, while the above disclosure discusses application of the concepts to emitters/receivers that continuously rotate 360° in a single direction, it will be understood that the same methods can be applied to oscillating emitters/receivers. In oscillating systems, a transducer or mirror may be controlled to oscillate between 90° to 400° with preferred ranges being approximately 120° to approximately 360°. Still further, while the description above is set forth in relation to use of an ultrasound transducer, other forms of wave emitters/receivers such lasers or light sources may be controlled to take advantage of the systems and methods described above.
The above described imaging system is disclosed in a non-limiting example of at least one application for use as an intravascular ultrasound system for imaging blood vessels. It will be understood that blood vessels are only one type of structure within a living body that may be imaged by the described system in accordance with the methods set forth above. More specifically, fluid filled or surrounded structures, both natural and man-made, within a living body that may be imaged can include for example, but without limitation, structures such as: organs including the liver, heart, kidneys, gall bladder, pancreas, lungs; ducts; intestines; nervous system structures including the brain, dural sac, spinal cord and peripheral nerves; the urinary tract; as well as valves within the blood or other systems of the body. In addition to imaging natural structures, the images may also include imaging man-made structures such as, but without limitation, heart valves, stents, shunts, filters and other devices positioned within the body.
Persons of ordinary skill in the art will appreciate that the embodiments encompassed by the present disclosure are not limited to the particular exemplary embodiments described above. In that regard, although illustrative embodiments have been shown and described, a wide range of modification, change, and substitution is contemplated in the foregoing disclosure. It is understood that such variations may be made to the foregoing without departing from the scope of the present disclosure. Accordingly, it is appropriate that the appended claims be construed broadly and in a manner consistent with the present disclosure.
This application claims priority to and the benefit of U.S. Provisional Application No. 61/646,080, filed May 11, 2012, which is hereby incorporated by reference in its entirety.
Number | Date | Country | |
---|---|---|---|
61646080 | May 2012 | US |