This application is a 35 U.S.C. §371 national stage entry of PCT/IB2009/055438, which has an international filing date of Dec. 1, 2009 and which is a continuation-in-part of U.S. patent application Ser. No. 12/127,544 filed on May 27, 2008, which is assigned to the assignee of the present patent application, the contents of which are incorporated herein by reference in their entireties.
The present invention relates generally to methods and systems for medical diagnostic imaging, and specifically to radio frequency (RF)-based imaging.
Narrowing of the coronary arteries due to atherosclerosis is commonly treated by implantation of a stent, using a catheter, to hold the artery open. In a large fraction of cases, however, the treated artery closes up again due to in-stent restenosis, necessitating further treatment. Accurate assessment of such restenosis generally requires re-catheterization. A number of non-invasive techniques have been proposed, such as in U.S. Pat. No. 6,729,336, in which an electromagnetic wave transmitter is used to excite a stent, and an acoustic sensor detects stent acoustic oscillations.
RF imaging is best known in the context of radar systems, but RF diagnostic imaging systems have also been developed for medical applications. For example, U.S. Patent Application Publication 2008/0169961, whose disclosure is incorporated herein by reference, describes computerized tomography using radar, which may be used for generating an image of living tissue. As another example, U.S. Pat. No. 7,454,242, whose disclosure is incorporated herein by reference, describes tissue-sensing adaptive radar imaging for breast tumor detection.
Various antenna designs have been proposed for RF imaging of body tissues. For example, U.S. Pat. No. 6,061,589, whose disclosure is incorporated herein by reference, describes a microwave antenna for use in a system for detecting an incipient tumor in living tissue, such as that of a human breast, in accordance with differences in relative dielectric characteristics. A composite Maltese Cross or bow-tie antenna construction is employed to irradiate the living tissue and to collect backscatter or other scatter returns.
Embodiments of the present invention that are described hereinbelow provide improved devices and methods for detecting features inside a living body using RF imaging techniques. Although some of these embodiments are directed specifically to detection of features in the heart, and specifically in the coronary arteries, the principles of these embodiments may similarly be applied in imaging, detection and tracking of features elsewhere in the body.
There is therefore provided, in accordance with an embodiment of the present invention, diagnostic apparatus, including an antenna, which is configured to direct radio frequency (RF) electromagnetic waves into a living body and to generate signals responsively to the waves that are scattered from within the body. Processing circuitry is configured to process the signals so as to locate a feature in a blood vessel in the body.
In disclosed embodiments, the apparatus includes an ultrasound transducer, and the processing circuitry is configured to guide the ultrasound transducer to direct an ultrasonic beam toward the feature. In one embodiment, the feature located by the processing circuitry includes a stent, and the ultrasound transducer is configured to generate a Doppler signal responsively to a flow of blood through the stent.
Additionally or alternatively, the apparatus includes a tracking unit, which is configured to track respective coordinates of the antenna and of the ultrasound transducer, and the processing circuitry is configured to guide the ultrasound transducer responsively to the respective coordinates. The apparatus typically includes position transducers fixed respectively to the ultrasound transducer and to the antenna, wherein the tracking unit is configured to track the respective coordinates responsively to position signals exchanged between the position transducers and the tracking system.
Further additionally or alternatively, the apparatus includes a display, wherein the processing circuitry is configured to guide the ultrasound transducer by driving the display to present to an operator of the ultrasound transducer an indication of a direction in which the ultrasound transducer should be aimed.
In a disclosed embodiment, the blood vessel is a coronary artery. The processing circuitry may be configured to track a cyclical motion of the feature over multiple cycles of the heart.
Typically, the processing circuitry is configured to locate the feature responsively to a difference in a dielectric constant of the feature relative to surrounding tissue.
There is also provided, in accordance with an embodiment of the present invention, diagnostic apparatus, including an antenna, having a front surface configured to brought into contact with an outer surface of a living body so as to direct radio frequency (RF) electromagnetic waves into the body and to generate signals responsively to the waves that are scattered from within the body. A dielectric gel is adapted to be spread between the outer surface of the body and the front surface of the antenna. Processing circuitry is configured to process the signals so as to locate a feature in the body.
Typically, the body has a first dielectric constant, and the gel has a second dielectric constant that is chosen to match the first dielectric constant. In disclosed embodiments, the gel has a dielectric constant that is between 30 and 75. The gel may be adhesive.
In one embodiment, the gel is water-based and includes an additive selected from a group of additives consisting of an alcohol, a salt, a sugar, and glycerin. Alternatively, the gel includes silicone and an additive having a dielectric constant greater than 70.
There is additionally provided, in accordance with an embodiment of the present invention, diagnostic apparatus, including an antenna, which has a front surface and is configured to direct radio frequency (RF) electromagnetic waves from the front surface into a living body and to generate signals responsively to the waves that are scattered from within the body, and which includes an array of antenna elements, each antenna element including a planar element at the front surface of the antenna and a cavity behind the planar element. Processing circuitry is configured to process the signals so as to locate a feature in the body.
In disclosed embodiments, the front surface of the antenna includes a printed circuit board, and the planar element of each antenna element includes a conductive radiator printed on the printed circuit board. The printed circuit board may include multiple conductive vias surrounding the radiator for isolating the antenna elements from one another.
There is further provided, in accordance with an embodiment of the present invention, diagnostic apparatus, including an antenna, including an array of antenna elements, which are configured to direct radio frequency (RF) electromagnetic waves into a living body and to generate signals responsively to the waves that are scattered from within the body. Excitation circuitry is coupled to apply a RF excitation waveform at multiple different frequencies to different transmitting antenna elements, selected from the array, according to a predetermined temporal pattern. Processing circuitry is coupled to receive the signals from different receiving antenna elements, selected from the array, and to process the signals at the different frequencies due to the different transmitting and receiving antenna elements so as to locate a feature in the body.
In some embodiments, the excitation circuitry includes a driver circuit, which is configured to generate the RF excitation waveform with a variable frequency, and a switching matrix, which is configured to select sets of the antenna elements in alternation, each set including at least one transmitting antenna element and one receiving antenna element, and for each selected set, to couple the driver circuit to excite the at least one transmitting antenna element at a selected frequency while coupling the processing circuitry to receive the signals from the at least one receiving antenna element. In a disclosed embodiment, the driver circuit and the switching matrix are coupled to select pairs of one transmitting antenna element and one receiving antenna element, and to excite the transmitting antenna in each pair at each of a plurality of frequencies in accordance with the predetermined temporal pattern.
In some embodiments, the apparatus includes a signal conditioning unit, which is configured to cancel a background component of the signals that arises from direct coupling between the transmitting and receiving antenna elements before the processing circuitry receives the signals. The signal conditioning unit may include an amplitude and phase modulator, which is coupled to receive the RF excitation waveform from the driver circuit, to modify a phase and amplitude of the received waveform so as to generate an anti-phased signal matching the background component, and to add the anti-phased signal to a signal received from the at least one receiving antenna element in order to cancel the background component.
In disclosed embodiments, the processing circuitry is configured to transform the signals received at the different frequencies due to the different transmitting and receiving antenna elements into a three-dimensional (3D) image, and to process the 3D image in order to find a location of the feature. In one embodiment, the processing circuitry is configured to compute a weighted sum of the signals received at the different frequencies due to the different transmitting and receiving antenna elements, using respective weights provided for a plurality of voxels in the 3D image, to determine values of the voxels in the 3D image.
There is moreover provided, in accordance with an embodiment of the present invention, a method for diagnosis, including directing radio frequency (RF) electromagnetic waves into a living body and generating signals responsively to the waves that are scattered from within the body. The signals are processed so as to locate a feature in a blood vessel in the body.
There is furthermore provided, in accordance with an embodiment of the present invention, a method for diagnosis, including spreading a dielectric gel between an outer surface of a living body and a front surface of an antenna. The front surface of the antenna is brought into contact, via the dielectric gel, with the outer surface of the living body so as to direct radio frequency (RF) electromagnetic waves into the body and to generate signals in the antenna responsively to the waves that are scattered from within the body. The signals are processed so as to locate a feature in the body.
There is also provided, in accordance with an embodiment of the present invention, a method for diagnosis, including providing an antenna, which has a front surface and which includes an array of antenna elements, each antenna element including a planar element at the front surface of the antenna and a cavity behind the planar element. Radio frequency (RF) electromagnetic waves are directed from the antenna elements via the front surface of the antenna into a living body and generating signals, using the antenna elements, responsively to the waves that are scattered from within the body. The signals are processed so as to locate a feature in the body.
There is additionally provided, in accordance with an embodiment of the present invention, a method for diagnosis, including defining a temporal pattern specifying a sequence of multiple different frequencies and spatial channels associated with an array of antenna elements. Radio frequency (RF) electromagnetic waves are directed at the multiple different frequencies into a living body from multiple different transmitting antenna elements that are selected from the array in accordance with the temporal pattern. Signals are generated responsively to the waves that are scattered from within the body and are received at multiple different receiving antenna elements that are selected from the array in accordance with the temporal pattern. The signals from the different receiving antenna elements at the different frequencies are processed so as to locate a feature in the body.
The present invention will be more fully understood from the following detailed description of the embodiments thereof, taken together with the drawings in which:
Embodiments of the present invention that are described hereinbelow use radar imaging techniques to identify and locate features in the body. Features are thus identified based on the difference in their complex dielectric constant (referring to both permittivity and conductivity) relative to the dielectric constant of the surrounding tissue. These techniques are particularly useful in detecting and tracking metal objects in the body, but they may also be used to locate features of other kinds, including both introduced objects, such as plastic objects, and naturally-occurring features, such as calcifications, and even pockets of air or other gases. The term “feature,” as used in the context of the present patent application and in the claims, should therefore be understood as referring to any item or location in the body having a distinct dielectric constant.
Some embodiments of the present invention are directed to locating features in the heart, and particularly in the coronary blood vessels. In these embodiments, an antenna directs RF electromagnetic waves toward the heart and receives the waves that are scattered from within the body. Processing circuitry processes the signals generated by the antenna due to the received waves in order to locate the feature or features of interest, and possibly to track the movement of such features over the course of the heart cycle.
The radar-based location of a feature may be used in guiding the beam of an ultrasound transducer toward the feature. In one such embodiment, the antenna and processing circuitry find the location of a stent in an artery and guide the ultrasound transducer to direct its beam toward the stent. The ultrasound transducer may operate in Doppler mode in order to measure the flow of blood through the stent and thus non-invasively assess possible restenosis in the stent.
In the embodiments that are described hereinbelow, the antenna comprises an array of antenna elements, with a front surface that is brought into contact with the outer surface (i.e., the skin) of the patient's body. A dielectric gel may be spread between the body surface and the front surface of the antenna in order to match the dielectric constants and thus improve the penetration of the RF waves into the body. Additionally or alternatively, the antenna elements may comprise a cavity and possibly other features to enhance the efficiency of coupling of electromagnetic energy from the antenna elements into the body while reducing loss and crosstalk between the elements.
In the disclosed embodiments, excitation circuitry applies a RF excitation waveform at multiple different frequencies to different transmitting antenna elements in the array. Meanwhile, the processing circuitry receives signals from different receiving antenna elements. The selection of transmitting and receiving antennas, as well as the selection of excitation frequency, follows a predetermined temporal pattern, which may be implemented by a switching matrix connected to the antenna elements.
As a result of this scheme of excitation and reception, the processing circuitry receives and processes signals from multiple spatial channels (corresponding to different pairs of antennas) at multiple different frequencies for each channel. Taken together in the time domain, these multi-frequency signals are equivalent to short pulses of RF energy. To reconstruct a three-dimensional (3D) image of the interior of the body and find the location of a feature or features, the processing circuitry applies a spatial transform to the set of received signals. The transform may, for example, comprise an inverse spherical Radon transform or an algebraic approximation of such a transform.
Despite measures that are taken to reduce coupling between antenna elements within the array, this sort of direct coupling still generates a strong background component, which tends to mask the signals due to scattered waves from the body. (The term “direct coupling,” as used in the context of the present patent application and in the claims, refers to short-range passage of RF waves between antenna elements by paths other than through the region of interest in the patient's body, including coupling that occurs within the array and near-field reflections.) In order to reduce this masking and enhance the dynamic range of the signals, in some embodiments a signal conditioning unit is used to adaptively cancel the background component out of the signals that are passed to the processing circuitry. To improve visibility of moving features, such as features in the heart, the signal conditioning unit or another element of the processing circuitry may even be configured to cancel all parts of the signals that do not vary over time.
Aiming the ultrasound probe correctly under these circumstances is difficult: The stent is small and is typically embedded in clutter in the ultrasound image due to other anatomical features; and the stent and surrounding features of the heart are in continual motion. Therefore, system 20 uses radar imaging in order to find the location of the stent and guide operator 22. For this purpose, an antenna 32 directs RF waves into the patient's chest. For good dielectric matching, in order to enhance the penetration of the RF waves into the body, a layer of a dielectric gel 33 is spread between the front surface of the antenna and the patient's skin. The gel may also have adhesive qualities, in order to aid in holding the antenna in place during the procedure.
Typically, gel 33 has a dielectric constant that is between 30 and 75. This value defines the effective dielectric constant of the antenna in its near-field. It is chosen to be close to the effective dielectric constant of the tissue layers in the path to the target region in the body. The desired dielectric constant may be achieved by increasing or decreasing the concentration of certain additives in the gel. For example, if a water-based gel is used, the additive may be an alcohol (such as ethanol), salt, sugar, or glycerin. Alternatively, a silicone gel may be used with an additive such as barium, having generally a dielectric constant greater than 70.
Antenna 32 is coupled by a cable or wireless link to a radar control unit 34. The control unit comprises processing circuitry 36, which drives the antenna to emit the RF waves into the body and processes the signals generated by the antenna due to reception of scattered waves from the body. Based on the received signals, circuitry 36 forms a 3D radar image of the interior of the body, and specifically, in the present embodiment, finds the location of the stent. These functions of circuitry 36 are described in detail hereinbelow.
In order to guide the direction of ultrasound transducer 24, the position coordinates (location and orientation) of the transducer and of antenna 32 are registered in a common coordinate frame. For this purpose, system 20 includes a tracking subsystem, comprising a tracking transmitter 38, which generates a field that is detected by sensors 40 and 42 on the antenna and on the ultrasound transducer, respectively. Transmitter 38 may, for example, generate a magnetic field, and sensors 40 and 42 may be magnetic sensors, as in the trakSTAR™ system distributed by Ascension Technology Corporation (Milton, Vt.). Alternatively, sensors 40 and 42 may be replaced by transmitting elements, which generate fields that are detected by a fixed sensor.
Further alternatively, other types of tracking devices may be used, such as optical, ultrasonic or mechanical position sensing devices, as are known in the art. For the sake of generality, the term “position transducer” is used herein to refer to the elements that are attached to ultrasound transducer 24 and antenna 32, such as sensors 40 and 42, for the purpose of finding their coordinates, regardless of the specific choice of position sensing technology. Additionally or alternatively, ultrasonic transducer 24 and antenna 32 may be mechanically fixed in a common frame of reference. A variety of alternative configurations are described in the above-mentioned U.S. patent application Ser. No. 12/127,544 and may similarly be used with the elements of system 20.
In the system configuration shown in
In an alternative embodiment (not shown in the figures), ultrasound transducer 24 may be held and manipulated by a robot arm, which is guided automatically by processing circuitry to aim the transducer in the desired direction.
Although
A driver circuit 58 generates signals, at multiple different frequencies, for exciting the transmitting antennas and demodulates the signals received by the receiving antennas. Typically, the signals are in the range of about 400 MHz to about 4 GHz, although higher and lower frequencies outside this range may also be used. A signal conditioning unit 56 between the driver circuit and switching matrix 54 amplifies the outgoing and the incoming signals and also cancels background components in the received signals. This functionality is also described below.
Front end 52 outputs the demodulated received signals (as intermediate-frequency or baseband signals) to a digital data acquisition unit 60, which samples and digitizes the signals. Unit 60 typically comprises a high-resolution analog/digital converter, such as a 14-bit converter, with suitable sampling circuits as are known in the art.
A target detection, measurement and tracking unit 62 receives and processes the digital samples. Unit 62, as described in detail hereinbelow, processes the sampled signals in order to generate a 3D radar image of the interior of the chest of patient 26. Within this image, elements having a dielectric constant that is different from that of the surrounding tissue, such as a metal stent in a coronary artery, stand out. On this basis, unit 62 identifies and measures the location coordinates of the stent relative to antenna 32.
Since the heart is in constant motion, unit 62 may also track and model the motion of the stent in order to more precisely guide ultrasound transducer 24. The direction of motion of the stent during the heart cycle also gives an indication of the direction of the stent axis (along which the ultrasound transducer should be aimed): Since the axis of the stent is oriented along the coronary artery in which the stent is implanted, and the coronary artery runs along the heart wall, the stent axis will typically be perpendicular to the direction of motion of the heart wall, and hence to the axis of motion of the stent in the radar image. As noted above, operator 22 is guided to aim ultrasound transducer 24 toward the stent in a direction along, or at least close to, the stent axis.
Although the present embodiment relates specifically to identification and tracking of a stent, the techniques and circuits that are described here may be used, by the same token, in locating and tracking other features in the coronary blood vessels, such as calcifications, as well as features elsewhere in the body.
A tracking unit 64 communicates with tracking transmitter 38 and receives position signals from position sensors 40 and 42. The tracking unit processes these signals in order to compute the coordinates of ultrasound transducer 24 and antenna 32 in the frame of reference of transmitter 38. The tracking unit may be a commercially-available device, such as in the above-mentioned FASTRAK system.
A guidance processor 66 receives the position (location and orientation) coordinates from tracking unit 64 and the position coordinates of the stent from feature detection subsystem 50. Guidance processor 66 registers the coordinates of the stent in the coordinate frame of the tracking unit or, equivalently, registers the coordinates of ultrasound transducer 24 in the coordinate frame of antenna 32, in which the stent coordinates have been found. The guidance processor is then able to compute the geometrical skew and offset between the present viewing axis of the ultrasound transducer and the desired viewing axis, which will intercept the stent along (or close to) the stent axis. Based on the computed skew and offset, the guidance processor may drive guidance display 44 to show operator 22 the required correction. For this purpose, the guidance display may show, for example, target crosshairs and directional arrows, or any other suitable sort of directional indication. Alternatively, console 28 may use the computed skew in adjusting the Doppler velocity readings to account for the angle of measurement relative to the flow.
Processing circuitry 36 typically comprises a combination of dedicated hardware circuits (such as in front end 52 and digital data acquisition unit 60) and programmable components. The front end circuits are described in detail hereinbelow. Target detection, measurement and tracking unit 62 and guidance processor 66 typically comprise programmable processors, such as a general-purpose microprocessor or a digital signal processor, which are programmed in software to carry out the functions that are described herein. Alternatively or additionally, these elements of circuitry 36 may comprise dedicated or programmable digital logic units such as an application-specific integrated circuit (ASIC) or a field-programmable gate array (FPGA). Although units 62 and 64 and processor 66 are shown, for the sake of conceptual clarity, as separate functional blocks, in practice at least some of the functions of these different blocks may be carried out by a single processor. Alternatively, the functions of a given block may be divided up among two or more separate processors.
Each antenna element 48 comprises a planar element comprising a conductive radiator 70, which is printed on a circuit board 72. This circuit board serves as the front surface of antenna, which is brought into contact with the patient's body. Circuit board 72 comprises multiple conductive vias 74 surrounding each radiator 70 for isolating antenna elements 48 from one another. The antenna elements are enclosed from behind by a case 76.
The printed (rear) side of board 72 includes an excitation transmission line 82 feeding the center point of radiator 70 through a conductive via. Transmission line 82 comprises a fifty-ohm microstrip, with a micro-miniature coaxial (MMCX) connector (not shown) for connecting to front end 50. Alternatively, other types of radiator shapes and feed lines may be used.
A conductive cavity 84 is attached to the component side of board 72 behind each radiator 70 in order to reduce antenna reverberations from back-lobe scattering and to increase the overall gain. (Antenna element 48 as shown in
wherein C0 is the speed of light, εr is the permittivity of the interior of the waveguide, and a is the largest transverse dimension of the waveguide. In the present example, with a frequency limit of 4 GHz, the depth of waveguide 86 is 15 mm. The waveguide creates an imaginary characteristic impedance, causing back-lobe radiation from radiator 70 to reflect from the cavity in phase with the back-lobe waves. This reflection enhances the external buffering of the antenna and attenuates non-TE and TM modes, and therefore reduces interference and noise.
Cross-coupling between antenna elements 48 can cause interference, which reduces the dynamic range and may saturate the receiver circuits. This cross-coupling is reduced in antenna 32 by appropriately setting the distance between the antenna elements in the array and by surrounding radiators 70 with conductive vias 74, as noted above. The vias serve as an electric wall that prevents internal waves from traveling between elements. They also create a conductive continuity between ground plane 80 on the component side of board 72 and the top conductive transverse plane of waveguide 86 located on the print side.
Antenna 32 is sealed against liquids and gels, thus preventing unwanted materials from reaching the print side and cavities of the antenna elements. Case 76, including cavities 84, can be constructed from a molded plastic with a suitable conducting coating. Additionally or alternatively, the antenna elements may be printed on the molded plastic after coating. Although switching matrix 54 is shown and described herein as a part of processing circuitry 36, it may alternatively be incorporated into antenna 32 or mounted adjacent to the antenna, thereby performing the switching alongside the patient and reducing the weight and rigidity of the cable from the antenna to control unit 34.
Although antenna 32 is shown here as a unitary assembly containing antenna elements 48, the antenna elements (of similar design to that shown in
Switching matrix 54 is designed for high isolation between channels, typically better than 40 dB over the entire frequency range of antenna 32. Switches 90, 92 and 94 are digitally controlled by a digital output control module 104 and allow fast (non-mechanical) switching. This fast switching is required in order to allow the entire waveform sequence of different antenna pairs and frequencies to be completed in a short frame time, as described hereinbelow. For this purpose, matrix 54 is typically configured to achieve a switching time of less than 1 μs.
Driver circuit 58 comprises a broadband signal generator 98, which generates the RF excitation waveform to drive the transmitting antenna elements, and a receiver 114, which receives and demodulates the signals generated by the receiving antenna elements. Signal generator 98 and receiver 114 are both synchronized and sweep their frequencies according to a predefined frequency plan, which is shown in
Receiver 114 is a tuned super-heterodyne receiver, which is able to adjust its bandwidth and gain according to the received signal. The receiver demodulates the received signals coherently, in synchronization with local oscillator 96, in order to extract both the amplitude and the phase of the signals. The complex ratio between the transmitted and received signals, as measured by detection, measurement and tracking unit 62 for each antenna pair at each selected frequency, indicates the frequency response along corresponding paths through the region of interest. This region includes the chest, thoracic cavity, beating heart and the stent itself.
Despite the measures described above for reducing coupling between different antenna elements 48 in antenna 32, the signals received from antenna 32 by signal conditioning unit 56 may still include a strong background component due to the direct coupling between the transmitting and receiving antenna elements. This background component raises the noise level due to transmitter nonlinearity and impurities in the transmitted signal and can even cause receiver 114 to saturate. It is therefore desirable to reduce the level of the background component that reaches the receiver in order to enhance the dynamic range of the radar image.
For this purpose, signal conditioning unit 56 comprises an amplitude and phase modulator, referred to here as an IQ modulator 108, which receives the sampled RF excitation waveform from coupler 100. The IQ modulator modifies the phase and amplitude of the sampled signal, under the control of an analog output control module 110, so as to generate an anti-phased signal matching the background component that is to be canceled. The amplitude and phase values of IQ modulator 108 are periodically updated and are then kept constant per frequency and per channel until coupling values change significantly and need updating. In other words, IQ modulator 108 outputs a signal that is equal in amplitude to the background component but 180° out of phase. A coupler 106 adds this anti-phased signal to the received signal from switching matrix 54 and thus cancels the background component without degrading the actual radar signal from the body. An amplifier 112 amplifies the signal following background cancellation for input to receiver 114.
Each sub-frame 122 begins with a settling time 126 (typically a few hundred microseconds) for locking the amplitude and phase of signal generator 98. Following this initial delay, switching matrix 54 selects different channels 124 in sequence. Each channel uses one transmitting antenna element and one receiving antenna element, up to a total of n channels (for example, one hundred such channels in the example shown in
In alternative embodiments (not shown in the figures), other sorts of channel configurations may be used. For example, in monostatic configurations, a selected antenna element may serve as both transmitter and receiver, as opposed to bistatic or multistatic configurations, in which each antenna either transmits or receives. As another option, antenna elements may simply transmit and receive broadband RF pulses, rather than multiple narrowband pulses as in the embodiment described above.
Based on the collected samples of the received signals, detection, measurement and tracking unit 62 detects small reflecting volumes within the region of interest (ROI) in the patient's body. As noted above, the corresponding reflections arise at the boundaries of media having different dielectric properties. The information provided by coherent detection of the signals over the broad range of frequencies covered in each frame is equivalent mathematically to the temporal information that would be provided by reflection of a single short pulse. The locations of the reflectors may be found by integrating over the propagation paths of the reflected waves, using an inverse spherical Radon transformation, for example.
In an embodiment of the present invention, detection, measurement and tracking unit 62 implements a first-order approximation of the inverse spherical Radon transform: For each voxel (x, y, z) in the ROI and for each frequency f and pair of antenna elements, a complex weight W(x,y,z,f,pair) is pre-calculated, either using an empirical calibration procedure or mathematical modeling. The weight is, in effect, the normalized complex amplitude (with conjugated phase) of the reflection that would be received at the receiving antenna in the pair from a point object at location (x,y,z) when irradiated by the transmitting antenna with a wave of frequency f. Because the body tissue through which the waves propagate is inhomogeneous, the weights may be adjusted, either empirically or by model calculation, to account for the specific tissue layers (skin, fat, muscle, lungs, etc.) through which the waves pass.
The set of weights thus derived defines a sort of matched filter. Detection, measurement and tracking unit 62 applies this filter to the matrix of complex signals Sig(f,pair) that it receives in any given frame in order to compute the reflection intensity V for each voxel, as a weighted sum over the received signals:
The inventors have found that this simplified approximation of the inverse spherical Radon transform is both robust and computationally efficient.
Front end 52 drives antenna 32 to emit and receive RF waves over multiple frequencies and spatial channels (antenna pairs), at a scanning step 130, as described above. Detection, measurement and tracking unit 62 collects samples of the received signals and applies the weights defined in equation (1) to transform the signal values to voxel intensities, at an image reconstruction step 132. To improve the clarity of the image, the processing circuitry may apply additional image processing operations, such as subtracting the mean voxel value from all voxels in the image. The mean value may be smoothed over multiple successive images using a recursive filter. Unit 62 then identifies the coordinates of the target feature, i.e., the stent or another strong reflector, in the 3D image, at a target identification step 134.
Features in the coronary arteries (or elsewhere in the heart), such as the stent, move regularly with the heart rhythm, as well as with chest movement due to respiration. In order to guide the ultrasound transducer, detection, measurement and tracking unit 62 tracks the motion of the target in the successive images, at a target tracking step 136. For example, unit 62 may apply a Kalman filter, as is known in the art, to estimate the motion trajectory of the target.
Guidance processor 66 registers the coordinates of ultrasound transducer 24 with antenna 32, at a coordinate registration step 138. The processor uses the position coordinates provided by position sensors 40 and 42 at this step, as explained above. Based on these coordinates, the processor registers the ultrasound beam in the coordinate frame of the 3D image that was reconstructed at step 132.
Guidance processor 66 drives guidance display 44 to guide operator 22 in aiming the ultrasound beam toward the target, at an aiming step 140. It could be possible but would probably be impractical for a human operator, to move ultrasound transducer 24 continually back and forth in synchronism with the motion of the heart. (Such tracking could be feasible for a robot driven by the processing circuitry.) To alleviate this difficulty, the guidance processor selects a single location within the trajectory of motion found at step 136 and guides the operator to aim at the selected location. Console 28 measures the flow through the stent at this location, at a flow measurement step 142, and thus provides an indication of the extent of any restenosis.
To select the target location at step 140, guidance processor 66 may, for example, find the center of mass of the trajectory found at step 136, and then choose a point that is displaced from the center toward the end of the trajectory that has the greater dwell time, which is the diastolic end. Blood flow through the coronary arteries occurs mainly during diastole, so that the diastolic end of the trajectory will give the strongest Doppler signal. Furthermore, aiming the ultrasound transducer toward the end of the trajectory with the greater dwell permits the ultrasonic beam to capture the stent for a longer part of each heart cycle and thus improves the signal/noise ratio.
It will be appreciated that the embodiments described above are cited by way of example, and that the present invention is not limited to what has been particularly shown and described hereinabove. Rather, the scope of the present invention includes both combinations and subcombinations of the various features described hereinabove, as well as variations and modifications thereof which would occur to persons skilled in the art upon reading the foregoing description and which are not disclosed in the prior art.
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/IB2009/055438 | 12/1/2009 | WO | 00 | 9/12/2012 |
Publishing Document | Publishing Date | Country | Kind |
---|---|---|---|
WO2011/067623 | 6/9/2011 | WO | A |
Number | Name | Date | Kind |
---|---|---|---|
4240445 | Durney et al. | Dec 1980 | A |
4344440 | Aaby et al. | Aug 1982 | A |
4557272 | Carr | Dec 1985 | A |
4632128 | Paglione et al. | Dec 1986 | A |
4640280 | Sterzer | Feb 1987 | A |
4641659 | Sepponen | Feb 1987 | A |
4774961 | Carr | Oct 1988 | A |
4825880 | Stauffer et al. | May 1989 | A |
4926868 | Larsen | May 1990 | A |
4958638 | Sharpe | Sep 1990 | A |
4986870 | Frohlich | Jan 1991 | A |
5003622 | Ma et al. | Mar 1991 | A |
5109855 | Gunter | May 1992 | A |
5394882 | Mawhinney | Mar 1995 | A |
5404877 | Nolan | Apr 1995 | A |
5540727 | Tockman et al. | Jul 1996 | A |
5549650 | Bornzin et al. | Aug 1996 | A |
5668555 | Starr | Sep 1997 | A |
5704355 | Bridges | Jan 1998 | A |
5766208 | McEwan | Jun 1998 | A |
5807257 | Bridges | Sep 1998 | A |
5829437 | Bridges | Nov 1998 | A |
5865177 | Segawa | Feb 1999 | A |
5967986 | Cimochowski et al. | Oct 1999 | A |
6019724 | Gronningsaeter et al. | Feb 2000 | A |
6061589 | Bridges et al. | May 2000 | A |
6064903 | Riechers et al. | May 2000 | A |
6093141 | Mosseri et al. | Jul 2000 | A |
6144344 | Kim | Nov 2000 | A |
6193669 | Degany et al. | Feb 2001 | B1 |
6233479 | Haddad et al. | May 2001 | B1 |
6330479 | Stauffer | Dec 2001 | B1 |
6454711 | Haddad et al. | Sep 2002 | B1 |
6480733 | Turcott | Nov 2002 | B1 |
6526318 | Ansarinia | Feb 2003 | B1 |
6592518 | Denker et al. | Jul 2003 | B2 |
6604404 | Paltieli et al. | Aug 2003 | B2 |
6729336 | Da Silva et al. | May 2004 | B2 |
6730033 | Yao et al. | May 2004 | B2 |
6755856 | Fierens et al. | Jun 2004 | B2 |
6940457 | Lee et al. | Sep 2005 | B2 |
7020508 | Stivoric et al. | Mar 2006 | B2 |
7122012 | Bouton et al. | Oct 2006 | B2 |
7130681 | Gebhardt et al. | Oct 2006 | B2 |
7191000 | Zhu et al. | Mar 2007 | B2 |
7197356 | Carr | Mar 2007 | B2 |
7266407 | Li et al. | Sep 2007 | B2 |
7267651 | Nelson | Sep 2007 | B2 |
7272431 | McGrath | Sep 2007 | B2 |
7280863 | Shachar | Oct 2007 | B2 |
7454242 | Fear et al. | Nov 2008 | B2 |
7474918 | Frants et al. | Jan 2009 | B2 |
7479790 | Choi | Jan 2009 | B2 |
7493154 | Bonner et al. | Feb 2009 | B2 |
7529398 | Zwirn et al. | May 2009 | B2 |
7570063 | Van Veen et al. | Aug 2009 | B2 |
7591792 | Bouton | Sep 2009 | B2 |
7697972 | Verard et al. | Apr 2010 | B2 |
7719280 | Lagae et al. | May 2010 | B2 |
7747302 | Milledge et al. | Jun 2010 | B2 |
7868627 | Turkovskyi | Jan 2011 | B2 |
8295920 | Bouton et al. | Oct 2012 | B2 |
8352015 | Bernstein et al. | Jan 2013 | B2 |
20020032386 | Sackner et al. | Mar 2002 | A1 |
20020045836 | Alkawwas | Apr 2002 | A1 |
20020049394 | Roy et al. | Apr 2002 | A1 |
20020050954 | Jeong-Kun et al. | May 2002 | A1 |
20020147405 | Denker et al. | Oct 2002 | A1 |
20030036713 | Bouton et al. | Feb 2003 | A1 |
20030088180 | Van Veen et al. | May 2003 | A1 |
20030100815 | Da Silva et al. | May 2003 | A1 |
20030219598 | Sakurai | Nov 2003 | A1 |
20040015087 | Boric-Lubecke et al. | Jan 2004 | A1 |
20040077952 | Rafter et al. | Apr 2004 | A1 |
20040249257 | Tupin et al. | Dec 2004 | A1 |
20040261721 | Steger | Dec 2004 | A1 |
20050038503 | Greenhalgh et al. | Feb 2005 | A1 |
20050107693 | Fear et al. | May 2005 | A1 |
20050245816 | Candidus et al. | Nov 2005 | A1 |
20060009813 | Taylor et al. | Jan 2006 | A1 |
20060265034 | Aknine et al. | Nov 2006 | A1 |
20070016032 | Aknine | Jan 2007 | A1 |
20070016050 | Moehring et al. | Jan 2007 | A1 |
20070100385 | Rawat | May 2007 | A1 |
20070123770 | Bouton et al. | May 2007 | A1 |
20070135721 | Zdeblick | Jun 2007 | A1 |
20070152812 | Wong et al. | Jul 2007 | A1 |
20070191733 | Gianchandani et al. | Aug 2007 | A1 |
20070263907 | McMakin et al. | Nov 2007 | A1 |
20080027313 | Shachar | Jan 2008 | A1 |
20080030284 | Tanaka et al. | Feb 2008 | A1 |
20080036668 | White et al. | Feb 2008 | A1 |
20080097199 | Mullen | Apr 2008 | A1 |
20080129511 | Yuen et al. | Jun 2008 | A1 |
20080167566 | Unver et al. | Jul 2008 | A1 |
20080169961 | Steinway et al. | Jul 2008 | A1 |
20080183247 | Harding | Jul 2008 | A1 |
20080200802 | Bhavaraju et al. | Aug 2008 | A1 |
20080224688 | Rubinsky et al. | Sep 2008 | A1 |
20080269589 | Thijs et al. | Oct 2008 | A1 |
20080283282 | Kawasaki et al. | Nov 2008 | A1 |
20080294036 | Hoi et al. | Nov 2008 | A1 |
20080316124 | Hook | Dec 2008 | A1 |
20090048500 | Corn | Feb 2009 | A1 |
20090076350 | Bly et al. | Mar 2009 | A1 |
20090153412 | Chiang et al. | Jun 2009 | A1 |
20090187109 | Hashimshony | Jul 2009 | A1 |
20090203972 | Heneghan et al. | Aug 2009 | A1 |
20090227882 | Foo | Sep 2009 | A1 |
20090240132 | Friedman | Sep 2009 | A1 |
20090240133 | Friedman | Sep 2009 | A1 |
20090281412 | Boyden et al. | Nov 2009 | A1 |
20090299175 | Bernstein et al. | Dec 2009 | A1 |
20090322636 | Brigham et al. | Dec 2009 | A1 |
20100056907 | Rappaport et al. | Mar 2010 | A1 |
20100081895 | Zand | Apr 2010 | A1 |
20100256462 | Rappaport et al. | Oct 2010 | A1 |
20100265159 | Ando et al. | Oct 2010 | A1 |
20100312301 | Stahmann | Dec 2010 | A1 |
20100321253 | Ayala Vazquez et al. | Dec 2010 | A1 |
20110004076 | Janna et al. | Jan 2011 | A1 |
20110022325 | Craddock et al. | Jan 2011 | A1 |
20110060215 | Tupin et al. | Mar 2011 | A1 |
20110068995 | Baliarda et al. | Mar 2011 | A1 |
20110130800 | Weinstein et al. | Jun 2011 | A1 |
20120029323 | Zhao | Feb 2012 | A1 |
20120098706 | Lin et al. | Apr 2012 | A1 |
20120104103 | Manzi | May 2012 | A1 |
Number | Date | Country |
---|---|---|
1834588 | Sep 2007 | EP |
2007-061359 | Mar 2007 | JP |
2008-515548 | May 2008 | JP |
2008-530546 | Jul 2008 | JP |
2008-542759 | Nov 2008 | JP |
2010-512190 | Apr 2010 | JP |
WO 2006130798 | Dec 2006 | WO |
WO 2007017861 | Feb 2007 | WO |
WO 2008070856 | Jun 2008 | WO |
WO2008148040 | Dec 2008 | WO |
WO2009031149 | Mar 2009 | WO |
WO2009031150 | Mar 2009 | WO |
WO2009152625 | Dec 2009 | WO |
WO2011067623 | Jun 2011 | WO |
Entry |
---|
International Search Report for International Application No. PCT/IB2009/055438, date of mailing: Jul. 20, 2010. |
Haude et al., Intracoronary Doppler—and Quantitative Coronary Angiography-Derived Predictors of Major Adverse Cardiac Events After Stent Implantation, Mar. 6, 2011, Circulation, vol. 103(9), pp. 1212-1217. |
Beyer-Enke et al., Intra-arterial Doppler flowmetry in the superficial femoral artery following angioplasty., 2000, European Radiology, vol. 10, No. 4, pp. 642-649. |
Ringer et al., Follow-up of Stented Carotid Arteries by Doppler Ultrasound, Sep. 2002, Neurosurgery, vol. 51, No. 3, pp. 639-643. |
Kantarci et al., Follow-up of Extracranial Vertebral Artery Stents with Doppler Sonography., Sep. 2006, American Journal of Roentgenology, vol. 187, pp. 779,787. |
Ghosh, et al., Immediate Evaluation of Angioplasty and Stenting Results in Supra-Aortic Arteries by Use of a Doppler-Tipped Guidewire, Aug. 2004, American Journal of Neuroradiology, vol. 25, pp. 1172-1176. |
Miura et al. “Time Domain Reflectometry: Measurement of Free Water in Normal Lung and Pulmonary Edema,” American Journal of Physiology—Lung Physiology 276:1 (1999), pp. L207-L212. |
International Search Report and Written Opinion, mailed Nov. 26, 2013 for Application No. PCT/IB2013/00663 filed Feb. 15, 2013. |
International Patent Application PCT/IB2009/055438, “Locating Features in the Heart Using Radio Frequency Imaging”, Filed on Dec. 1, 2009. |
International Application PCT/IB2009/055438 Search Report dated Jul. 20, 2010. |
International Application PCT/IB2010/054861 Search Report dated Apr. 5, 2011. |
International Search Report for International Application No. PCT/IB2011/053244, date of mailing, Dec. 2, 2011. |
Ascension Technology Corporation, “TrakSTAR Adds Versatility to Ascension's New Product Line: Desktop Model Joins driveBAY Tracker for Fast Guidance of Miniaturized Sensor”, USA, Apr. 7, 2008. |
Claron Technology Inc., “MicronTracker 3:A New Generation of Optical Trackers”, Canada, 2009. |
Immersion Corporation, “Immersion Introduces New 3D Digitizing Product-MicroScribe G2; Faster Data Transfer, USB Compatibility, New Industrial Design”, Press Release, San Jose, USA, Jul. 1, 2002. |
Polhemus, “Fastrak: The Fast and Easy Digital Tracker”, USA, 2008. |
Czum et al., “The Vascular Diagnostic Laboratory”, The Heart & Vascular Institute Newsletter, vol. 1, USA, Winter, 2001. |
Lal et al., “Duplex ultrasound velocity criteria for the stented carotid artery”, Journal of Vascular Surgery, vol. 47, No. 1, pp. 63-73, Jan. 2008. |
Larsson et al., “State Diagrams of the Heart—a New Approach to Describing Cardiac Mechanics”, Cardiovascular Utrasound 7:22 (2009). |
Bell et al., “A Low-Profile Achimedean Spiral Antenna Using an EBG Ground Plane”, IEEE Antennas and Wireless Propagation Letters 3, pp. 223-226 (2004). |
Supplementary European Search Report for Application No. EP 10834292.4 (PCT/IB2010/054861) dated Dec. 4, 2014. |
Li, J.C. et al., “Microwave Imaging of Cerebral Edema”, Proceedings of the IEEE, IEEE, NY, US, vol. 70, No. 5; May 1, 1982, pp. 523-524. |
Guido Biffi Gentili et al., “A Versatile Microwave Plethysmograph for the Monitoring of Physiological Parameters”, IEEE Transactions on Biomedical Engineering, IEEE Service Center, Pitscataway, NJ, US, vol. 49, No. 10, Oct. 1, 2002. |
Pedersen P C et al., “Microwave Reflection and Transmission Measurements for Pulmonary Diagnosis and Monitoring”, IEEE Transactions on Biomedical Engineering, IEEE Service Center, Piscataway, NJ, US, vol. BME-19, No. 1, Jan. 1, 1978; pp. 40-48. |
International Search Report for International Application No. PCT/IB2011/053246 date of mailing, Dec. 13, 2011. |
Extended Search Report for European Application No. 11809360.8, date of mailing, Mar. 11, 2014. |
Written Opinion of the International Searching Authority, mailed Jul. 20, 2010, for International Application No. PCT/IB2009/055438. |
International Preliminary Report on Patentability, mailed Jun. 14, 2012, for International Application No. PCT/IB2009/055438. |
Christine N. Paulson et al. “Ultra-wideband radar methods and techniques of medical sensing and imaging” Proceedings of Spie, vol. 6007, Nov. 9, 2005, p. 60070L. |
S. I. Alekseev et al. “Human Skin permittivity determined by millimeter wave reflection measurements”, Bioelectromagnetics, vol. 28, No. 5, Jul. 1, 2007, pp. 331-339. |
Supplementary European Search Report and European Search Opinion, dated Jun. 13, 2013, for European Application No. 09851811.1. |
Notice of Reasons for Rejection, mailed Apr. 28, 2014, for JP 2012-541588. |
Notice of Reasons for Rejection, mailed Mar. 31, 2015, for JP 2012-541588. |
Written Opinion of the International Searching Authority, mailed Dec. 2, 2011, for International Application No. PCT/IB2011/053244. |
International Preliminary Report on Patentability, mailed Jan. 31, 2013, for International Application No. PCT/IB2011/053244. |
Supplementary European Search Report and European Search Opinion, dated Feb. 26, 2014, for European Application No. 11809359. |
International Search Report and Written Opinion, mailed Feb. 26, 2015, for International Application No. PCT/IL2014/050937. |
Notice of Reasons for Rejection, mailed Apr. 17, 2015, for JP 2013-520273. |
Number | Date | Country | |
---|---|---|---|
20120330151 A1 | Dec 2012 | US |