|Publication number||US6061589 A|
|Application number||US 08/926,621|
|Publication date||9 May 2000|
|Filing date||9 Sep 1997|
|Priority date||1 Jul 1994|
|Also published as||DE69532367D1, DE69532367T2, EP0694282A2, EP0694282A3, EP0694282B1|
|Publication number||08926621, 926621, US 6061589 A, US 6061589A, US-A-6061589, US6061589 A, US6061589A|
|Inventors||Jack E. Bridges, Allen Taflov, Susan C. Hagness, Alan Sahakian|
|Original Assignee||Interstitial, Inc.|
|Export Citation||BiBTeX, EndNote, RefMan|
|Non-Patent Citations (10), Referenced by (172), Classifications (11), Legal Events (6)|
|External Links: USPTO, USPTO Assignment, Espacenet|
This application is a continuation-in-part of Ser. No. 08/641,834 filed May 2, 1996, now U.S. Pat. No. 5,829,439 which was a continuation-in-part of Ser. No. 08/492,998 filed Jun. 21, 1995, U.S. Pat. No. 5,704,355 which was a continuation-in-part of application Ser. No. 08/269,691 filed Jul. 1, 1994, now abandoned.
Breast cancer is one of the leading causes of death for women. About one out of eight or nine women are expected to develop tumors of the breast, and about one out of sixteen to twenty are expected to die prematurely from breast cancer. Mammography or other X-ray methods are currently most used for detection of breast cancers. However, every time a mammogram is taken, the patient incurs a small risk of having a breast tumor induced by the ionizing radiation properties of the X-rays used during the mammogram. Also, the mammography X-ray process is costly and sometimes imprecise. Accordingly, the National Cancer Institute has not recommended mammograms for women under fifty years of age, who are not as likely to develop breast cancers as are older women. However, while only about twenty two percent of breast cancers occur in women under fifty, data suggests that breast cancer is more aggressive in pre-menopausal women. Furthermore, women under forty are getting the disease in increasing numbers--about eleven thousand annually now-and no one knows why.
Mammograms require interpretation by radiologists. One radiologist has said "I generally can spot cancers between five and ten millimeters in diameter. The prognosis is excellent then." However, about ten to fifteen percent of tumors of this size are not detected. One study showed major clinical disagreements for about one-third of the same mammograms that were interpreted by a group of radiologists. Further, many women find that undergoing a mammogram is a decidedly painful experience.
Thus, alternative methods to detect breast cancers are needed, especially non-invasive methods that do not entail added risks, that can detect tumors as small as two millimeters in diameter, that are not unduly unpleasant to the patient, and that can be used for early detection. Such a screening system is needed because extensive studies have demonstrated that early detection of small breast tumors leads to the most effective treatment. While X-ray mammography can detect lesions of approximately five mm or larger, the accuracy may range between 30% and 75%, depending on the skill of the diagnostic radiologist. Repeated X-ray examinations, however, are not encouraged because these may become carcinogenic. These considerations, in addition to cost considerations, have led physicians to recommend that women wait until the age of fifty before having routine mammograms.
One solution would be a non-ionizing, non-invasive, and low cost detection or screening method to detect very small malignant breast tumors. It could greatly increase without hazard the number of patients examined and would identify those patients who need diagnostic X-ray examinations, where the added hazards and costs could be justified.
About one in eight women develop breast cancers and about one in sixteen die prematurely from this disease. Despite strong encouragement, less than half of the millions of women who should be are routinely screened. Some of the reasons are cost and discomfort experienced during mammography. Other concerns are the additional risks associated with ionizing radiation, especially for routine exams for women under fifty. However, while only twenty two percent of breast cancers occur in women under fit, data suggest that breast cancer is more aggressive in pre-menopausal women. A screening procedure need only identify breasts with abnormalities. The precision and imaging requirements associated with diagnostic purposes and treatment monitoring, while desirable, need not apply.
Further, mammography fails to detect between five and twenty-five percent of malignant breast cancers. While many tumors are detected by mammography, the method is not capable of determining whether the tumor is benign or malignant. When a tumor is detected, a biopsy must be made to determine the nature of the tumor. Less than twenty percent of the tumors detected in mammograms are found to be malignant, but the biopsy is both painful and costly. Thus, a complementary method is needed that uses a different modality, such as microwaves, that may be beneficial in reducing the number of false negatives and false positives.
There are several generic cancer detection methods: sonic, chemical, nuclear and non-ionizing electromagnetic. The sonic, chemical and nuclear (such as MRI) techniques have been under study for some time and, while some interesting approaches are being followed, none have been publicized as being available in the near future for low cost screening.
The method described herein is not like the known non-ionizing electromagnetic methods. Studies have considered the use of electromagnetic, non-ionizing methods to detect or image portions of the human body. An excellent summary of such activity is presented in a publication entitled "Medical Applications of Microwave Imaging", edited by L. E. Larsen and J. H. Jacobi, IEEE Press 1986.* These activities include microwave thermography, radar techniques to image biological tissues, microwave holography and tomography, video pulse radar, frequency modulation pulse compression techniques for biological imaging, microwave imaging with diffraction tomography, inverse scattering approaches, and medical imaging using an electrical impedance. The publications in this book contain about five hundred citations, some of which are duplicates. The technology cited not only includes electromagnetic disciplines, but also notes related studies in sonic imaging and seismic imaging. To update these data, the IEEE transactions on Medical Imaging, Biomedical Engineering, Microwave Theory and Techniques and Antennas and Propagation have been reviewed. Also surveyed was the publication Microwave Power and Engineering. This update has indicated little significant progress in the aforementioned electromagnetic techniques that would be important to detect breast and prostate cancers. In these publications breast and prostate cancer detection systems based on the concepts described in this specification were not presented.
Many important reasons exist for this lack of progress. In the case of microwave thermography, adequate depth of penetration, along with the required resolution, may not be realized, except for large cancers. In the case of holography, reflections at the skin-air interface tend to mask the desired returns from breast tumors beneath the skin. Further, illuminating the entire volume of a breast either requires excessive power (with possible biological hazards) or acceptance of poor signal-to-noise ratios. In the case of through-the-body electromagnetic techniques, such as tomography, the attenuation characteristics of the body are such that long wavelengths are usually used, with an attendant loss of resolution. Imaging by determining perturbations in body impedance caused by the presence of tumors as sensed by multi-electrode arrays have been either inadequate in sensitivity or subject to false alarms. A major difficulty with some of the multi-electrode or multi-antenna systems is that matrix methods are used to process the measured data into an image. With such methods, small errors in the measurements or assumptions are often enlarged during processing. Other problems, such an ambiguous or inconclusive results, may occur with the computational matrix method itself.
More recently, two microwave methods have been proposed to detect breast cancers. Both employ no more than a few narrow-band, fixed frequencies. One such method applies unfocused 900 megahertz energy directly to the breast via a resonant, open faced cavity applied directly to the breast. Some promising preliminary clinical results have been claimed if the data from adjacent breasts is compared. However, inconclusive results may occur in the vicinity of the nipple, where substantial variations in skin thicknesses occur. The other method proposes scanning a breast with a microwave beam via a dielectric slab pressed against the breast. The so-called beam of this method was relatively broad because it was developed by an open-ended wave guide pressed against the slab. The waveguide did not embody focusing features. As a result, reflections from many incidental scatterers not near a possible tumor could be expected. Such reflections could mask the desired returns. Further, it is nearly impossible to press such a plate uniformly against the breast without developing some air gaps that can cause massive reflections. Studies demonstrate that the use of such a plate is disadvantageous and results in excessive reflections and reverberations from the skin that may mask any desired return.
The present method is not like thermography, which uses the passive microwave or infrared emissions generated by malignant tumors which exhibit elevated temperatures with respect to normal breast tissue. Such radiated emissions must first pass through the normal breast tissue, then through the skin to a sensor placed on or external to the breast. Systems that have used such passive emissions have been clinically evaluated. A substantial number of analyses have been conducted, including a few that have considered but have not resolved the perturbing effects of the skin. The results are viewed as less efficacious than other detection modalities, such as X-ray mammography.
The present method is also unlike hyperthermia methods that are designed to heat malignant tumors in situ, preferentially over normal breast tissues. Typically, an antenna or an array of antennas are placed near or over the breast. Fixed frequency, microwave energy generally below three GHz is then directed into the breast through the skin, through normal breast tissue and thence to the tumor. In some cases, electronically controllable phased arrays have been considered. The heating of the tumor has been optimized by invasively emplacing a sensor in the tumor; the sensor provides feedback signals to the electronic controls of the phased array to adjust the phase of each antenna for best results.
In summary, many, if not all of the past microwave methods, such as those used in thermography, hyperthermia or microwave imaging or detection, have experienced difficulties because breast tissues are not homogeneous; they are heterogeneous. A principal heterogeneity is associated with the skin, especially near a breast nipple. In addition, the dielectric parameters of the breast of one human may be significantly different than the dielectric parameters of other humans. Another problem is the necessity to focus the microwave energy into as small a spot or voxel as possible so as to be able to resolve very small tumors. This is best done by resorting to the smaller wavelengths or higher frequencies, frequencies well above those typically used for most of the reported therapeutic or imaging electromagnetic systems. However, such use of higher frequencies increases the difficulty of coping with skin-introduced heterogeneity effects.
The methods and apparatus described herein envision an illuminator or antenna that propagates focused, short duration pulses of low power microwave beams into the breast or other tissues. The focusing may be effectively achieved either physically by lenses, reflectors, or phased arrays, or artificially by synthetic aperture methods. The pulses may be generated by pulse sources or synthetically by swept-frequency, Fourier inversion methods. When these effectively focused beams encounter a tumor, more energy is returned than from normal breast tissue. This occurs because a malignant tumor has significantly greater dielectric parameters than are exhibited by normal tissues. The backscatter returns from a possible tumor are captured by a collector that is also a part of the antenna and may be effectively focused in the region of the expected tumor. By carefully focusing the power into a small volume or voxel within the breast and scanning the focal point from the skin to chest wall and from side to side, tumors can be detected and imaged.
Unwanted returns from heterogeneity in the breast are suppressed by several methods. One method is the use of a wide-aperture, confocal illuminator and collector. Typically, many of the features of the illuminator and collector may share the same position or function. Such a wide-aperture, con-focal design tends to average out minor variations in the dielectric parameters as well and to suppress returns from sources not near the focal point. Another method is to illuminate the breast or other tissues with short duration pulses, whether generated synthetically or in real time, to isolate the returns from scatters adjacent to a possible tumor and to compensate for propagation losses.
Specifically, the preferred apparatus of the disclosed system employs short duration pulses (either in real time or synthetically generated) in combination with a multi-antenna array (realized either physically or synthetically). The perturbations introduced by the skin-related interfaces and any other heterogeneity in the breast or other tissue are detected and are used to suppress the unwanted effects of such perturbations. Such suppression may include a determination of the skin's thickness and its dielectric parameters. Another method uses electronically controllable phased or synthetic arrays. Such arrays, in combination with signal processing, can develop the approximate dielectric parameter of each breast or other tissue segment.
Such methods also can be used to determine whether or not a tumor is malignant or benign by noting the amplitude of the returns; malignant tumors usually exhibit much larger returns than benign tumors. All of the aforementioned techniques that are useful to suppress heterogeneity are also useful to help resolve whether tumors are benign or malignant.
This invention and improvements described in this application were first proposed to two government Federal agencies. Four proposals were rejected because the group of experts that reviewed them did not believe the concept would work. A fifth proposal was funded if only extensive computer simulations were conducted to demonstrate feasibility. Although the proposals incorporated technology that was proven in other areas, such as conventional radar, video-pulse radar or confocal microscopy, the review panels did not believe such technology would be viable when applied to the human breast. Some of the reasons given were that the human breast is opaque to microwaves, that the dielectric difference between normal and malignant tissues is too small for reliable detection, that resonant scattering enhancement could not occur for a tumor immersed in normal breast tissue, that the backscatter effects from heterogeneity would mask returns from malignant tumors, and that the microwave method lacked resolution needed to detect small tumors or microcalcifications. This issue of viability was put to rest only after extensive computer-aided studies that required some 250 hours of Cray computer time.
The principal object of this invention is to provide an improved microwave breast cancer radiation and reception antenna that is not only practical, but also exhibits improved performance under safe operating conditions. This is done by the use of co-sited resistance-loaded cross-polarized dipoles that are placed directly upon the breast. One of the dipoles is used to illuminate the breast with a very broad-band microwave pulse; the other dipole is used to receive backscattered returns. The cross-dipole arrangement functions so that spurious reverberating waveforms on the illuminating dipole that could mask the desired backscatter returns are substantially suppressed. Such reverberation arises from any mismatch between the wave impedances of the illuminating dipole and of the breast tissue itself. Additional reverberation from reflections from the distal ends of the illuminating dipole are suppressed by resistively loading the dipoles. This arrangement permits the use of low-noise, high-gain amplifiers in the backscatter reception channel, thereby increasing the dynamic range without resorting to increased power levels for the power source for the illuminating dipole. The combination of the foregoing features permits detection of the backscattered waveforms in a format that allows software processing to detect and image the tumors.
To realize the above, this improved sensor antenna:
a) is practical and uses printed circuit antennas and simple coaxial cables instead of hard-to-form complex antennas or higher mode waveguides;
b) partially rejects the applied pulse that reduces the reliance passive means such as directional couplers or active means, such as the use of time gating, to separate the desired returns from the applied pulses, thereby permitting the use of a practical preamplifier to enhance the dynamic range;
c) partially rejects backscatter from planar surfaces, such as skin or the chest wall, thereby permitting simpler processing of the sensed backscatter waves;
d) suppresses reverberation effects by a combination of optimized spatially distributed resistances for the antenna elements;
e) enhances the response from asymmetrical tumors that are near planar surfaces; and
f) combines functions a) through e) into a simple sensor element.
The following figures are used to explain the concepts and design of the invention as applied to microwave detection of cancers and other abnormalities in heterogeneous tissues:
FIG. 1A, is a conceptual view of an active millimeter wave (microwave) breast cancer detection system, with a patient;
FIG. 1B illustrates displays for plural focal lengths, with the generalized system of FIG. 1A;
FIG. 2 is a simplified block diagram, with two supplemental charts, that illustrates the principal functions of a mmw breast cancer detection system;
FIG. 3 is a graph of relative dielectric constants of muscle, fat, breast tissue and breast cancer as reported by various investigators;
FIG. 4 is a graph of conductivity (mho/m) of muscle, fat, breast tissue and breast cancers as reported by various investigators;
FIG. 5 is a graph of attenuation, wavelength, and depth of penetration in normal breast tissue as a function of frequency, based on the data presented in FIGS. 3 and 4;
FIG. 6 is a block diagram of a mmw breast cancer detection system that employs a "passive" signal separation technique in combination with a conventional heterodyne receiver to detect tumor-scattered returns;
FIG. 7 is a block diagram of a breast cancer detection system that is a modification of the system shown in FIG. 6; in FIG. 7 the amplitude and phase of the backscattered signal is developed by means of a phase coherent detection method. Such characterization of the returns permits use of a frequency swept waveform that generates a synthetic pulse;
FIG. 8 presents a view of the face of an array of nine double-ridged waveguides that can be used as part of the illuminator in an electronically controllable phased array that is applied directly to a breast or other tissue;
FIG. 9 illustrates the electronically controllable portion of an illuminator that employs phased array techniques;
FIG. 10 illustrates how a source that generates very short duration pulses can be used to create a synthetic aperture that progressively excites each aperture afforded by an antenna array;
FIG. 11 is a time-of-flight diagram for five aperture antennas of an array that are focused on a voxel near the center of the array;
FIG. 12 is a time-of-flight diagram for five aperture antennas of an array that are focused on a voxel near one side of the array;
FIG. 13 is a diagram that illustrates how the times-of-flight of FIG. 11 can be processed by adding a predetermined time delay to the return from each of the antennas to enhance returns from a center voxel;
FIG. 14 is a diagram that illustrates how the times-of-flight of FIG. 12 can be used to suppress the returns from an off-center voxel by adding a predetermined time delay to the return from each of the antennas to suppress the returns from a voxel near one side of the array;
FIG. 15 illustrates how composite returns from a near center voxel and from an off-center voxel can be combined in real time to show how the returns from a voxel at the focus is enhanced over a voxel not near the focus;
FIG. 16 illustrates the perturbations introduced by a layer of skin between a region near an antenna aperture and normal breast tissues for a propagating wave that impinges at right angles upon the surface of the skin;
FIG. 17 illustrates the perturbation introduced by a layer of the skin between a region near an antenna aperture and normal breast tissues for a propagating wave that impinges at an acute angle with respect to the surface of the skin;
FIG. 18 is a simplified perspective view that illustrates the conceptual design of a broad band aperture;
FIG. 19 illustrates simple apparatus that can be used to complement currently availabale tumor detection and location methods;
FIG. 20 is a perspective illustration of a Maltese Cross antenna used to explain how it suppresses interaction between two crossed dipole antennas;
FIG. 21A is a simplified circuit diagram used to explain operation of the Maltese Cross antenna of FIG. 20;
FIG. 21B is a further simplified circuit diagram for the Maltese Cross antenna;
FIG. 22 is a conceptual design, partly in perspective, of a complete microwave breast cancer detection system that employs the Maltese Cross antenna;
FIG. 23 is a perspective illustration, similar to FIG. 8, of an electronically controlled array of Maltese Cross antennas;
FIGS. 24A and 24B illustrate the geometry of a bow-tie antenna that embodies the use of continuous resistive loading;
FIG. 25 is a chart that shows the time history of the spurious scatter from a typical bow-tie antenna; and
FIG. 26 illustrates another antenna configuration that helps suppress spurious responses.
The use of non-ionizing electromagnetic microwave or millimeter waves offers several advantages over x-ray mammography in detecting incipient breast cancers. (To simplify this discussion both microwave and millimeter wavelength regimes will be referred to as mm waves, or mmw.) Non-ionizing electromagnetic systems can be operated at sufficiently low levels so as to preclude biological hazards. A contrast ratio of the order of 20:1 is potentially usable for mm waves in tissue, whereas there is less than a few per cent range of densities for X-rays for soft tissue. The tissue-mm wave interaction also exhibits additional phenomena that can be drawn upon to enhance the performance. For example, when the diameter of a highly conducting spherical tumor (e.g., an incipient cancer) is of the order of a wavelength in breast tissue, a resonance effect occurs that increases the effective scattering cross-section of the tumor. If the tumor is non-spherical, then the polarization of the scattered waves may be different than that of the impinging waveform. In some cases, side-scattered or forward-scattered energy can also be utilized. For purposes of this specification, tissue-mm waves are defined in terms of wavelength in a medium having a dielectric constant like that of breast tissue, not air. Thus, the operating frequency for an electromagnetic wave source used in the described system is preferably in the range of one to thirty GHz.
Other than the use of millimeter wave and microwave thermography to detect breast cancers, there has been little activity toward use of such mm wave approaches to detect breast cancers. As noted earlier, some of the problems that have to be overcome are formidable. First, simply flooding the torso of a female with mmw energy introduces numerous problems. How does one single out the scattered returns from a three millimeter circumference tumor from the immensely larger scattered returns from the torso? How is the defocusing effect of an air-skin interface overcome? Is the breast tissue sufficiently transparent, at mmw frequencies, to propagate energy into and out of the breast? Are the dielectric properties of tumors sufficiently different from normal breast tissue for effective detection of small (e.g., three mm circumference) incipient cancers?
To understand the invention and its novel features, the basic concept will first be briefly described. The ability of the millimeter wave electromagnetic energy to penetrate normal breast tissues will be demonstrated. The special equipment and operating conditions of prior related developments will be described, particularly as used to realize high resolution simultaneously with good penetration. Finally, the antennae of the present invention will be described, particularly in connection with FIGS. 20-26.
FIGS. 1A and 1B illustrate the basic concepts. FIG. 1A illustrates, on a conceptual basis, possible prototype equipment. The patient 21 arranges one of her breasts 22 to contact an illuminator 23 as shown. Mm waves are generated within the equipment housing 24. These mm waves are then propagated into the selected breast 22 as a refracted or reflected electromagnetic mm wave that is focused at a predetermined point or volume (voxel) within the breast. This is done by means of an interface and focusing apparatus, as described hereinafter. Further apparatus is used to cause the focal point of the beam to scan different small volumes or voxels within the breast. When this happens, the scattered mmw energy from any tumor present in the breast becomes much larger than other scattering sources, since the dielectric properties of a tumor are radically different those of normal breast tissue. The scattered returns may be collected as backscattered power via the same interface and focusing apparatus that is used to propagate the mmw power into the breast. The collected scatter waves can then be processed by either analog or digital methods to form an image of the tumor.
A stepped FM sweep, similar to pulse compression in "Chirp" radar, can be used to synthesize a time domain response to isolate shallow from in-depth scattering and to mitigate the effects of heterogeneity in the dielectric characteristic within the breast. A functional goal of the combined confocal and time-domain features of the present system is to isolate the returns from tumors from spurious returns generated by heterogeneity in adjacent normal tissues.
The amount of collected backscattered energy and its time-of-flight or round trip time delay can be presented in a three dimensional format, as shown generally in FIG. 1B. For illustrative purposes, it is assumed that the impinging energy can be selectively and sharply focused into three vertical planes that are parallel to the patient's chest, wherein the X-Y planes at maximum depth 25, medium depth 26, and shallow depth 27 are shown. The three coordinates show the backscatter returns 28, the "X" coordinate 30 and the "Y" coordinate 29. Small vertical lines 31 are shown for numerous combinations of X and Y coordinates.
The amplitudes or heights of most of these lines 31 are proportional to the non-target returns that can arise from, for example, the tissues that surround the rib cage. Note that a large return 32 exists in the center of the medium depth display 26. This large return 32 is assumed to arise from a tumor that is at the focus of the impinging energy in the medium depth plane 26. In the center of the shallow plane 27 there is a somewhat smaller response 33 caused by a tumor intercepting and only scattering a small portion of the impinging beam. Note that in the center of the deeper plane 25, the return 34 is smaller because, it is assumed, the focused energy has largely been scattered by the tumor in the medium depth plane 26 before it arrives at the center of the deeper plane 25.
FIG. 2 illustrates the different functions that are employed to detect and image tumors in a system 40. A microwave or millimeter wave source 41 generates a frequency sweep as a linear function of time. Such a sweep can be processed to generate a synthetic pulse. This output power from the microwave source 41 is applied to an input port 45 of a directional coupler or circulator 43 a via waveguide 44. The output power from a port 46 of the circulator 43, arrow 49, is directed to an illumination/collector device 47 via a waveguide 48.
The illumination/collector device 47, hereinafter called an illuminator, consists of an array of waveguides wherein the time delay of the power injected or collected by each guide can be electronically controlled or processed. The distal end serves as an antenna that both radiates and collects power. By electronically controlling or processing the time delay for each of the waveguides in illuminator 47, the power (arrow 49) that is applied to the illuminator 47 can be electronically focused into a preselected spot within the patient's breast 50. The illuminator 47, via electronic means, also collects backscattered power from the same focal point. Each antenna/waveguide can be connected to port 46 of the directional coupler 43 via the waveguide 48 and a waveguide switching arrangement, described hereafter, in the illuminator 47.
The requisite directing action of the directional coupler or circulator 43 can be realized by several "passive" means, such as a balanced bridge circuit or magic tee, a directional coupler, or a circulator; see Ramo et al. Fields and Waves in Communication Electronics, John Wiley and Sons, New York, (1965) sections 11.17, 11.8 and 9.16. "Active" means of separating the applied power of source 41 from the tumor-scattered power from illuminator/collector 47 are possible in the time domain. For example, very short duration pulses of mmw energy can be applied and the returns separated by time gating methods. Other "active" methods currently employed in some modern radar systems can be used, such as pulse compression, Chirp or frequency modulation radar; see Skolnik, Introduction to Modem Radar Systems, McGraw-Hill (1980).
More detailed explanations are presented in subsequent portions of this application. When the radiation from the excited antenna/waveguide encounters a tumor, significantly more power is backscattered or returned because the malignant tumor has a larger dielectric constant than the normal breast tissue adjacent the tumor. The return power (arrow 52) from the tumor and other scatterers is collected by the excited antenna/waveguide in the illuminator 47. Via waveguide 48, the backscatter power from the excited antenna/waveguide is injected into the port 46 of the circulator 43. At the port 54 of the circulator 43, the collected return power is separated from the input power that was applied to the port 45 by the unique properties of the circulator or directional coupler 43.
The output of port 54 is injected into a detector 56 via a waveguide 55. A portion of the swept-frequency mmw signal from the microwave source 41, is also applied, via a waveguide 57, to a detector 56. The detector 56 senses the amplitude of the backscatter or return power. The scatter detector 56 also ascertains the time of arrival of each of the amplitude components of the return power by comparing the difference in the frequency of the power from the source 41 via cable 57 with a given amplitude component of the return power from the tumor via waveguides 48 and 55. The greater the frequency difference, the farther away the tumor is from the illuminator 47. The backscatter detector may also include means to compensate the returns for the attenuation and divergence experienced by the incident and scattered energy. By these or other well known radar techniques, such as FM Chirp or pulse compression radar, the "A Scope" display 60 shown in the upper right-hand corner of FIG. 2 is developed.
The foregoing discussion assumes monostatic operation, where the illuminating antenna is also the collecting antenna, a configuration common to many radar and sonar systems. In this case, the same waveguide 48 that carries the illuminating waveform 49 also serves as the return pathway for the backscattered waveforms 52. However, performance can be greatly improved and design requirements greatly simplified if the single antenna is replaced by a pair of perpendicular bisecting dipole antennas [often called cross polarized or Maltese Cross antennas]. One of the cross-polarized antenna pair is used to illuminate the breast and collect the co-polarized [co-pol] backscatter returns 52 via cable 48. The co-pol returns would be processed as described earlier. The other dipole pair, which is oriented perpendicular to the co-polarized pair, collects the cross polarized [cross-pol] returns. These cross-pol returns 77 are applied directly to the backscatter detector 56. The time history of the cross-pol returns would be similar to the time history of the co-pol returns as illustrated by waveform 60, except that the backscatter returns from planar surfaces such as the skin (64) and the chest wall (66) would be suppressed, as discussed later Additionally, any reverberation, such as 63, in the circuits that supply power to the illuminating antenna will also be suppressed in the cross polarized returns. Rather than controlling the time-delay in the microwave band by changing the time delay in each of a plurality of waveguides, the backscatter detector 56 can convert the time histories of the returns collected by each antenna into a format suitable for digital processing.
The time and the amplitude, in the "A scope" display 60, are developed from the returns from one of the waveguides in the illuminator. Equipment generated reflections 63, reflections 64 from the skin interface, reflections 65 from the tumor, and returns 66 from the chest wall are stored as real time responses in a signal processor 70, which is fed by a cable 59. Signals 63,64 and 66 are known and can be discarded by a time gating process in the signal processor 70, using data supplied from detector 56 via the cable 59. The tumor returns 65 are further isolated by a time-gating function in the signal processor 70; they are then stored in an imaging/display subsystem 72, fed via a cable 71 from the processor 70. A cable 73 supplies an output from display 72 to a scan control subsystem 74. Other antenna waveguides may be separately and sequentially connected to cable 48 by a switching arrangement in the illuminator 47. These are actuated by signals from the scan control subsystem 74, which is connected by a cable 75 to the illuminator subsystem 47. Additional time-amplitude displays 60 are developed and the related time-amplitude data is stored in the signal processor 70. In FIG. 2 the signal processor 70, the display 72 and the scan control 74 are all part of a subsystem 76.
The signal processor 70 of FIG. 2 acts as an electronically controllable lens that is capable of positioning an apparent focal point anywhere within the breast 50 of the patient 51. Consider a round, spherical lens that is focusing sunlight into a spot. The focusing occurs, in part, because the on-axis rays that originate from the center of the lens are time-delayed more than the rays which are deflected from the circumference of the lens. The computer in processor 70 views each of the detected backscattered returns such as illustrated at 60 as one "ray". To obtain the potential backscatters for a scatterer that is, say, 30 mm below the center face of the illuminator 47, a similar process is employed. To each of the time-amplitude returns of display 60, a delay time is added. The most delay time is added to the returns from the waveguides (antennas) near the center of the face of the illuminator 47. Progressively less time delay is added, the further away the guides are from the center of the face of the illuminator. After the appropriate delays have been added to determine if a tumor is in a voxel 10 mm below the center of the illuminator 47, such delayed returns are summed together. If a scatterer or tumor is in the selected voxel, the scattered returns from the tumor add coherently. Other returns, such as from the skin or chest wall, will not add coherently. This process is further described hereinafter.
Large returns, such as returns 63 from interfaces in the equipment or returns 64 from the skin, that might overload the dynamic range of the processor 70, can be gated out prior to summing. Further refinement is possible, such as compensating for the attenuation and dispersion of the observed responses as a function of time or distance from the antenna. If the propagation time is not well known, the delay times can be varied electronically in a manner that converges onto a best estimate. Another method is to develop the approximate depth and dielectric constant of the skin, so as to further enhance the accuracy and resolution.
In the case of the cross-pol returns, a similar treatment is used as noted in the preceding discussion. Ideally, the undesired cross-pol returns from the planar surfaces of the skin 64 and chest wall 66 are completely suppressed. However, the suppression of these undesired returns cannot always be complete. The cross-pol suppression not only relies on he complete symmetry of the geometry of the crossed dipole arrangement, but also on uniformity of the dielectric material surrounding the crossed dipoles. Such cannot be expected under practical conditions, therefor time gating and other means to suppress the unwanted returns may be needed.
The system as discussed for FIG. 2 requires further explanation, especially regarding the applicator and electronic focusing features. Further discussion and means to compensate for heterogeneity of the breast will be made subsequently. To show that the breast tissue is reasonably translucent, published data on the dielectric properties of normal and malignant breast tissues will be discussed in connection with FIGS. 3, 4, and 5. FIGS. 6 and 7 provide further functional details of the microwave source 41, the directional coupler 43, and the backscatter detector 56. FIG. 8 illustrates a construction in which the faces of waveguides are emplaced in direct contact with the breast. FIG. 9 illustrates an electronically controllable illuminator that simultaneously excites the waveguides in the illuminator and that is best used with the synthetic pulse approach. FIG. 10 illustrates a system that generates and applies a pulse directly, rather than synthetically.
FIGS. 3, 4 and 5 provide data that demonstrate that non-lactating breast tissue has different dielectric properties than either tumors or muscle tissues. Moreover, the attenuation of mm waves in such breast tissue is not excessive in the 1 to 15 GHz region and hence permits reasonable operating conditions for "passive" power and signal directors. Additional attenuation can be tolerated by the use of "active" power and signal directors such that operation up to thirty GHz is possible.
FIG. 3 summarizes data on relative permittivity, scale 101, as a function of frequency, scale 102. These data demonstrate that the relative dielectric properties of low-water-content tissues and normal breast tissues are significantly lower than for high-water-content tissues and tumors, either human or non-human. The low-water content data for curve 103 were developed by Chaudhary (1984) for human breast tumors. Johnson (1972) developed the data for curve 105 for fat, bone and low-water content tissue. Edrich (1976) generated the data for cattle fat, shown at 107. Burdette (1986) generated in vivo data for canine fat, illustrated in a curve 109. The high-water-content data for curve 104 was developed by Chaudhary (1984) for human breast tumors. Johnson (1972) developed data for muscle and high water content tissues, shown in a curve 106. Rogers (1983) generated the data, shown in a curve 108, for mouse tumors. Edrich (1986) collected data for canine muscle, illustrated at 110. Burdette (1986) provided in vivo data, shown in a curve 112, for canine muscle tissue. Note that in the case of muscle or tumor tissues, the relative dielectric constant is of the order of forty or more, depending on the frequency. In the case of low-water-content tissues, such as breast or fat, the dielectric constant is in the order of five to ten, as measured for in vitro studies. The in vivo measurement of Burdette (1986), shown in curve 109, shows an approximate increase by a factor of two in the relative permeability over the data developed by Johnson (1972), curve 105. The in vitro breast tissue measurements by Chaudhary (1984), curve 103, fall somewhat in between the in vitro values developed by Johnson (1972) and the in vivo measurements of Burdette (1986).
FIG. 4 presents similar data on the conductivity of both low and high-water-content tissues. The conductivity in mhos/meter, scale 121, is the ordinate and the frequency, curve 122, is the abscissa. The low-water-content tissues are human breast tissues, shown in a curve 123 derived from Chaudhary (1984). The low-water-content fat and bone of the curves 125 is from Johnson (1972). Cattle fat, shown at 127, is from Edrich (1986). The high-water-content tissues of the curve 124 are human breast tumors, data by Chaudhary (1984). High-water-content muscle tissue is in a curve 126, data by Johnson (1972). Mouse tumors, shown in a curve 128, are from Rogers (1983). Rat muscle data 130 is derived from Edrich (1986). Canine fat data are presented in a curve 131 from Burdette (1980), and canine muscle data in a curve 132 taken from Burdette (1980). Curves 126,128 and 132 are essentially coincident. Note that conductivity, as a function of frequency, tends to increase substantially above 6 GHz and that the 40 to 90 GHz measurements 127 and 130 of Edrich tend to fall in line with the trends established by measurement make up to 10 GHz.
Based on the data presented in FIGS. 2 and 3, FIG. 5 shows the depth of penetration 140, wavelength 142, and attenuation 144 as a function of frequency 146 for the propagation of millimeter waves in non-lactating breast tissue. Above ten GHz, some uncertainty associated with the trend extrapolation is suggested by the range of possible values of the penetration depth 140 or attenuation 144. A value of nine was used for the relative dielectric constant and the extrapolated values of Chaudhary (relative to the data developed by Johnson) from FIG. 4 were used for the conductivity. From these data, it is seen that the breast tissue behaves as a lossy dielectric for frequencies substantially exceeding five GHz, wherein ω=2πF and ε=εo εr (permittivity of free space)×(relative dielectric constant), σ is the conductivity, μ is the permeability, f is the frequency, λ is the wavelength, and δ is the depth of penetration (see Ramo (1965) page 334 Sec. 6.05).
Since ωε>>σ, the approximate lossy dielectric equations are as follows:
δ=2[σ(μ/ε).sup.1/2 ].sup.-1 (2)
This defines the generic feasibility of the system to be described hereinafter. There are two requirements that must be met. First, the total path loss attenuation (in and out) should be substantially less than the dynamic range, typically in the order of 100 dB, wherein the dynamic range is defined in dB as equal to: 10 log[(largest signal power)/(smallest detectable signal power)]. Second, the wavelength in the irradiation apparatus (illuminator 47) and in the breast of the patient should be sufficiently small so that small tumors can be resolved. This, for the system discussed here, requires that, preferably, the wavelength in illuminator 47 and in the breast tissue should not exceed two or three times the circumference of the smallest tumor. If an operating frequency of 15 GHz is chosen for a passive power and signal detector, it is seen that the path loss is about 5 dB/cm, or 50 dB total path loss, in and out, for a 5 cm path length. The wavelength at 15 GHz is about 0.6 cm, which is about equal to the diameter of the smaller tumors.
As an introduction to some of the basic circuit functions; FIG. 6 illustrates a functional block diagram of a microwave breast cancer detection and imaging system 200 that employs a conventional heterodyne receiver. System 200 comprises the following subsystems: a millimeter power generator subsystem 241, a passive power and signal director 243, a focusing illuminator subsystem 228, a heterodyne receiver 253 employed for signal detection, a scanner control 248, and a signal processing and display subsystem 255. There is a connection 254 from the unit 253 to the display unit 255 and a connection 251 from the scan control 248 to the display unit 255.
Electromagnetic wave energy flows, via the power and signal director 243, from the power generation system 241 to the illuminator subsystem 228. The illuminator subsystem 228 includes features that physically or effectively focus, by synthetic methods, mmw energy into a voxel and electronic or mechanical features that allow the focal point to scan over a focal or voxel point 262 within the breast of a patient. The scanner control 248, through a mechanical connection or electrical cable, controls the location of the focal point 262 in three dimensions such that the focal point 262 is progressively positioned into each voxel (smallest volume element) of the breast under consideration. When the focal point encounters a tumor, the scattered returns are collected by the illuminator 247 of the subsystem 228 and then the scattered power (arrow 245) is supplied to the power and signal director 243 and then to the detection subsystem 253.
The scattered power 245 is separated from the impinging power 244 (supplied to the illuminator) by means of a directional coupler or circulator 224 within the power and signal director 243. A discussion of each of the aforementioned subsystems follows.
The power generation and control subsystem 241 is comprised of two functional blocks: an electromagnetic wave power source 220 connected via a cable 221 to an isolator and power splitter 222. Power splitter 222, in turn, is connected, via a cable 242, to the power input port 225 of the circulator 224 in the power and signal director 243. The power output and backscattered input port 227 of the circulator 224 is connected, via a cable 246, to the matching network 226 at the input of the illuminator subsystem 228. The output port 234 of the circulator 224 is connected to the signal processor subsystem 253 via a cable 252.
Many of the functions of these components are obvious. The isolator/power splitter 222 electrically isolates the power source 220 from any load variations that might be introduced by the circulator 224 or other components of the illuminator subsystem 228. The function of the circulator or directional coupler 224 is to extract the backscattered returns from the applied power. Otherwise, the high level of the power applied to the illuminator subsystem 228 would tend to overload the signal detector 253. Thus, the electromagnetic input signal injected into port 225 is directed out of port 227 and thence to illuminator 228. The backscattered returns that are applied to port 227 appear at port 234, wherein the amplitude of the applied power is greatly suppressed.
The performance requirements for the signal detection system 253 are not too stringent. The simplest version may use a simple heterodyne receiver as an RF voltmeter to measure the output of the circulator 224 at port 234. A reference signal from the power generation subsystem 241 can be supplied to the hetrodyne receiver 253 via a conductor 258 to stabilize the local oscillators in the receiver.
Other versions of the invention, such as the system 300 shown in FIG. 7, offer additional signal processing options. The system 300 of FIG. 7 illustrates the use of two synchronous receivers or detectors in a modification of the system of FIG. 6 in which only the signal detection subsystem is changed, with subsystem 253 of FIG. 6 replaced by a dual subsystem 353 that includes two hybrid tee synchronous detectors 361 and 362. The power splitter 222 (FIG. 6) provides a reference signal, on conductor 258 (FIG. 7) to system 353 as well as to the power and signal director 243. Each of the hybrid tee devices 361 and 362 forms a product between the applied input signal and the composite backscatter returns. However, one of the reference waveforms is shifted ninety degrees with respect to the other reference waveform; the following relationships result, where:
ω is the angular frequency of the millimeter waves;
θ is an arbitrary reference fixed phase angle;
λ is the wavelength;
χ is the electrical path length from the scatterer to the hybrid tee;
β is the propagation phase constant and β=2π/λ;
χβ is the accumulated phase shift.
The output from each of the hybrid tees 361 and 362 is the product of the returned, scattered waveform and the reference waveform. Considering just the low frequency components of such products, the output of each of the hybrid tees is as follows:
A cos [ωt+θ-2χβ] cos [ωt+Θ]=A/2[cos(-2χβ)]+high freq. component (3)
A cos [ωt+θ-2χβ] cos [ωt+θ-π/2]=A/2[sin(-2χβ)]+high freq. component (4)
The dual receiver system depicted in FIG. 7 draws the reference waveforms from the isolator-power splitter 222 via cable 258 of FIG. 6. An attenuator-power splitter 339 is used to reduce the amplitude of the waveform presented to the two phase shifters 340 and 341 via appropriate cables or other conductors 352 and 351. The output waveform of phase shifter 341 is advanced or retarded ninety degrees relative to phase shifter 340 to provide the desired quadrature relationship. The quadrature reference waveforms from circuits 340 and 341 are applied, via cables 342 and 343, to the hybrid tees 362 and 361, respectively. The output of port 234 of the circulator 224 supplies the power from the backscattered returns, via cable 252, to a power splitter-isolator 335. This circuit 335 diverts the backscatter return signal equally into cables 336 and 337, thus supplying the backscatter signals to the hybrid tees 361 and 362, respectively. These tees 361 and 362 each form a product between the reference waveform (from conductors 343 and 342, respectively) and the backscatter signals (on lines 336 and 337, respectively). The low frequency output from these two devices 361 and 362, on cables 345 and 346, provides critical inputs to the signal processor and display subsystem 255 (FIG. 6). Other variations of the above technique may be used to improve the signal-to-noise ratio, such as modulating the reference waveforms with another frequency well above the highest frequency of interest in the detected backscattered return. This removes the output signal well away from the troublesome shot noise that occurs at very low frequencies.
The scanner control subsystems (74 in FIG. 2, 248 in FIG. 6, control how the breast of the patient is scanned. In the case of electronic scanning, FIG. 2, scanner control 74 controls the X and the Y positions by methods described in conjunction with FIGS. 9 and 10. The scanner control (74 or 248) is connected to the illuminator (47 or 247) and, via the display unit 72, to the signal processing unit (70 or 255) in each of the described systems of FIGS. 2 and 6. Other methods, particularly techniques that synthesize large aperture antennas, could also be used; see FIGS. 10-13 and 15, as described hereinafter.
In any of the described systems the signal processing and display subsystem (e.g., subsystem 255 in FIG. 6) can employ any number of processing or display methods so as to suitably display the scattered returns.
The parameters needed to determine the spot size (voxel size) are the diameter D of the antenna aperture, the focal distance R, the spot diameter d, and the wavelength X of the millimeter wave in the media. See Kay (1966) and Smith (1966) for more complete development of relationships. The spot size becomes:
d=2R/80 /D (5)
As was noted earlier in regard to FIG. 5, reasonable penetration losses of about five dB/cm occur for wavelengths of the order of six mm. Thus, if tumors in the order of three mm in circumference are to be resolved, the beam width or effective spot size should not exceed the tumor circumference by much more than a factor of three. Preferably, for improved spatial resolution, the wavelength should not exceed the tumor circumference by a factor of three. To achieve a spot diameter of six mm, the ratio of the focal distance R to the aperture diameter D should be about 0.5.
Another design consideration is the depth of field Δ, as related to the aforementioned parameters and the apparent angle of resolution Φ. Thus the depth of field becomes:
Δ=[R.sup.2 Φ]/[DħRΦ], (6)
Again, to obtain good spatial discrimination, the focal distance R should be small compared to the aperture diameter D.
However, short focal lengths cannot be easily developed if the dielectric constant between the media that form an interface are greatly different. This would be the case if an attempt is made to propagate millimeter wave power in air and thence into the breast. As seen in FIG. 3, the dielectric constant of breast tissue is of the order of nine, and such a large value (relative to a value of one for air) causes substantial reflection and refraction of the incident power at any air-breast interface. More importantly, the apparent R/D ratio is reduced; that can lead to a radical increase in the spot size. This difficulty is mitigated by making the properties of dielectric materials within the illuminator similar to those for the human breast or skin and placing the illuminator directly on and in contact with the breast.
The foregoing can be better understood by referring to FIG. 8, which shows a nine aperture waveguide module 450. Nine double-ridged waveguides 451 are used. Each of these waveguides is filled with a dielectric material 452 that approximates the relative dielectric constant of either the normal breast tissue or of the ski. In the case of a screening system, only four waveguide apertures might be used. The combination of waveguide apertures is pressed directly against the breast. By proper timing or phasing of the signals to each of the waveguides 451, the focal point can be positioned within the breast without the need for mechanical movement.
FIG. 9 illustrates one way that the phase or timing of the signal applied to the waveguides may be controlled to position the focal point in a medium 418 which contains the aperture antennas 410A, 410B, 410C and 410D and focal points 412 and 413. For initial explanation, the thickness of the skin is neglected. A source 400 of microwave power applies equi-phased power via wave guides 401 and 402 to two power splitters 403 and 404. The outputs of the splitters are applied, via wave guides 405A, 405B, 405C and 405D, to the circulators or directional couplers 406A, 406B, 406C and 406D. The forward power through these devices is transferred via the guides 407A, 407B, 407C and 407D to the variable time delay or phase control devices 408A, 408B, 408C and 408D. The return power is transferred via wave guides 416A, 416B, 416C and 416D to a subsystem 417 that collects the returns in a format suitable for additional processing by subsystem 76 of FIG. 2. The time delay in each device 408A-408D may be controlled by changing the magnetic field bias applied to a ferrite element within each of the devices. Such bias may be supplied via the cables 414A, 414B, 414C and 414D from the time delay control subsystem 415. The outputs from wave guides 405 are controlled by the signal processing and display subsystem 72 of FIG. 2. Via wave guides 409A-409D, the time delayed or phased controlled power is supplied to the aperture antennas 410A-410D. A portion of the outputs from these antennas reaches the desired focal point 412 via pathways 411A, 411B, 411C and 411D. At point 412, the phases of the rays shown are nearly identical.
Assuming a time delay of t1, t2, etc. for each of the delay control elements 408A-408D and path lengths (411) d1, d2, etc., then t1 +d1 /v=d4 /v for constructive addition where v is the velocity of propagation in the medium 418. To meet this requirement, t1 =(d4 -d1)/v. Other time delays can be calculated in the same way.
Other methods of control are possible by controlling the phase of the signals applied to each aperture instead of by the timing devices. In this case, the relative phase between the signals applied to apertures 411C and 411D can be redefined by noting the following, where ω is the radian frequency [2πf and θ12 is the phase difference between 411C and 411D, such that
θ.sub.12 =ωt.sub.2 -t.sub.1 ]. (7)
The described confocal arrangement permits the reflected or back-scattered signals from voxel 412 in medium (breast 418) to return by the same pathways as the applied wave form. These backscatter signals are collected by the aperture antennas 410 and progress back through the time delay devices 408 to the circulators 406. These, in turn, supply data on the backscatter returns to subsystem 417.
The discussion of FIG. 2 also covers USC of the illuminator and focusing arrangements described above in combination with an "active" or time domain method of separating the applied power from the scattered power. Other such active or time domain methods utilize a "Chirp radar" to produce added resolution in depth and additional clutter suppression. The amount of resolution possible is dependent of the bandwidth of the source 400, the waveguides, and the antenna. For example, if a ridged guide antenna is used as illustrated in FIG. 8, a swept frequency band might extend from 3 to 8 GHz. The resulting 6 GHz bandwidth would generate range cell resolutions in the other of 10 mm. If broader bandwidth may be increased to about 15 GHz. The swept frequency band might extend from 1 GHz to 15 GHz. This produces a bandwidth of 14 GHz and a range solution in the order of 4 mm. If a time-gating, pulse system is used, pulse duration in the order of 130 picoseconds will produce a resolution in the order of 10 mm; and pulse duration of about 46 picosecond will result a resolution in the order of 4 mm. This added resolution would not only discriminate the incident power, but also may suppress backscatter returns from the different dielectric interfaces, such as the skin or the chest wall.
Active methods are of particular interest because these methods may be functional with total path losses up to 100 dB. Such path losses might be difficult to overcome with a passive system, since it may be difficult to reduce clutter levels below 50 to 70 dB in the applied power. Since some of the clutter can be reduced by considering only the returns in just one voxel, active systems may be more viable over a wider dynamic range. Also, shorter wavelengths with greater resolution can be used, since active systems can accept greater path losses, possibly as much as might be experienced by a system with an operating frequency as high as 30 GHz.
Swept frequency methods can be considered. For example, an FM Chirp radar method that has been used in weapons detection systems effectively separates desired returns from those generated by system discontinuities. A version of this would be attractive in conjunction with the confocal illumination method to separate the effects of near surface discontinuities or heterogeneities from the returns at greater depth. Linear FM pulse compression radar (PCR) techniques might also be considered. These have been described by Jacobi (7) for biological imaging applications. The theoretical resolution of a PCR is given by R=C/2B, where C is the is the velocity of propagation in the media, and B is the bandwidth of the transmitted wave form. Assuming a mid-band frequency of 8 GHz, a 5 GHz sweep and a medium with a dielectric constant of nine, a range resolution of one cm is indicated. However, a 2.5 GHz sweep may be more readily realized and could produce an in-tissue resolution of two cm. To realize this performance, the FM sweep must be highly linear, a pulse compression filter developed for this application and the dispersion effects of the dielectric compensated.
A stepped or swept frequency input impedance Fourier inversion alternative exists. This option transforms data developed from the frequency domain measurements to the time domain via digital processing, thereby eliminating the need for a pulse compression filter. This can be implemented by using either the phased array of FIG. 9 or the synthetic aperture technique described hereinafter in connection with FIGS. 10-15. The output signals from the circuit shown in FIG. 7 on lines 346 and 345 can be viewed as a complex input impedance, S(jω), at a radian frequency of ω(ω=2πf) to the illuminator. As the frequency is stepped from a low frequency to a higher frequency, the complex input impedance for each frequency is stored in a digital computer. If the frequency is swept or stepped over a band similar to that noted for the PCR system, similar spatial resolutions can be realized. Via digital processing, the complex input impedance data can be used to develop the complex Fourier component for each stepped frequency. Using inverse Fourier transformation, this series of spectral components can be processed to develop the real time response to an applied pulse or step function. The transformed data is then in the form of an amplitude vs. time response, similar to a radar A scope display, as if an impulse or stepped function had been applied at port 234 of the circulator or directional coupler 224 in FIG. 7. Initially, as illustrated in FIG. 2, the returns from system discontinuities, such as from connectors and the interface with the antenna in the illuminator, will be displayed. Then, the reflections from the near-surface (e.g. skin) anomalies in the breast will be displayed, along with the reflections from the deeper anomalies occurring at the longer times. The stepped frequency option offers the opportunity to include a standard correction at each frequency increment for a typical dispersion characteristic for normal breast tissues and could also include compensation for other factors such as path loss or system dispersion in the ferrite phase shifters. Some of the more modern network analyzers include a built-in stepped or swept frequency to time domain processing option.
The underlying mathematical basis is as follows. The genera Fourier transformations are: ##EQU1## Where
F(f) is the impulse or step response,
S(jω) is the Fourier transformation of F(f),
ω=2πf, and t is time.
FIG. 10 shows an alternative way of generating very short duration rise-time pulses. As opposed to the use of synthetic pulse techniques (FM pulse compression radar) to artificially generate the pulse, pulse sources and digital oscilloscopes are commercially available that can generate and resolve pulses with rise or fall times in the order of 14 picoseconds. Such suggests effective operating bandwidths of up to 30 GHz. The main limitation is that the dynamic range is more limited, in the order of 40 to 60 dB, unless mitigation measures are employed. In the case of the synthetic pulse method, commercially available network analyzers can replace and supplement the circuit functions illustrated in FIG. 7. These instruments can generate a synthetic pulse and yet have dynamic ranges better than 100 dB.
The functional block diagram of system 500 in FIG. 10 shows a combination of a 14 picosecond pulse source 502 that supplies pulses to a power/signal director or circulator 243 via a waveguide 504. In many digital oscilloscope test configurations, the signal director function 243 is just a passive three-way pad or power splitter. The r)ad equally divides the power applied from one port 225 to the other two ports 227 and 234, typically with a 6 dB loss from the injection port 225 to either of the two output ports. Because the power that is applied to the digital scope 528 through port 234, via a waveguide 509, cannot exceed a specified value to avoid overload or burnout, the power that can be applied to the illuminator 547 of system 500 will be similarly limited.
One solution is to continue to use the power/signal director function 243, such as illustrated in FIG. 6. Other devices, such as a directional coupler, could be used. These will suppress the power applied to the scope 528 by about 15 to 30 dB relative to the power applied to the illuminator 547, thereby increasing the level of backscatter returns that can be applied to the digital scope by about 15-30 dB. This improves the signal-to-noise ratio of the system. The signal-to-noise ratio can be further increased by time averaging successive backscatter returns from the same voxel. The improvement is proportional to the square root of the number of returns that are summed together.
A more optimum method would be to use a pulse source and digital sampling scope in combination with a Maltese Cross antenna and other components described hereinafter. This arrangement suppresses the amplitude of the applied pulse in the cross-pol dipole circuits. This allows peak pulse powers in the source 502 that are substantially greater than the allowable burnout threshold for the input to the digital sampling scope 528. In addition, this arrangement permits the use of low-noise, high-gain preamplifiers, which can substantially increase the dynamic range.
The power/signal director 243 supplies power to the illuminator 547 and to the digital sampling oscilloscope 528 and to the illuminator 547 via the waveguides 509 and 508 respectively. The illuminator 547 is comprised of two rotary waveguide switches 510 and 516 that provide a selectable, single path interconnection to any one of an array of twenty-five aperture antennas 520; only five antennas 520A through 520E are shown. Switch 510 directs power to any one of five switches like switch 516 via the waveguide terminals 512A through 512E. Switch 510 may be connected to any one of its five terminals in response to signals from the scan control 74 via cable 75; see FIG. 2. As illustrated, terminal 512B of switch 510 supplies power via a waveguide 514 to switch 516B. Other switches and related waveguides that connect to switch 510 are not shown. Switch 516, in response from signals from the scan control, may direct power into any one of five antenna apertures, 520A, 520B, 520C, 520D and 520E via the waveguides 518A through 518E. As shown, switch 516 is connected to the terminal of switch 516 that supplies power to antenna aperture 520B via waveguide 518B.
Each of the twenty-five aperture antennas 520 (only five are shown) is excited sequentially, via the illustrated switches, from the power/signal director 243. Backscatter returns go to the digital scope 528 via waveguide 509. Scope 528 also receives a reference signal from source 502 via a connection 522, a delay circuit 524, and a further connector 526. That is, backscatter of tumors and other sources in the breast are collected by one antenna aperture at a time and, via the pad or circulator 243, are directed to the scope 528. The scope 528 displays and stores the backscattered returns as is suggested in the A scope display 60, FIG. 2. The scope 528 also supplies the data stored in the scope, via a cable 530, to the signal processing subsystem 70. This subsystem, if needed, can time average the returns as needed to improve the signal-to-noise ratio. It further electronically manipulates the time-history from each of the aperture antennas to form an image of the scatterers in a breast by means of the display subsystem 72.
FIGS. 11, 12, 13, 14 and 15 illustrate how the signal process)r 70 of FIG. 10 can electronically manipulate the time-amplitude backscatter returns that may be simultaneously developed by the synthetic pulse system of FIGS. 2 and 7 in combination with the electronically controllable phased array of FIG. 9. In the case of the real-time pulse system of FIG. 10, FIGS. 11-15 also illustrate how the sequentially developed returns can be similarly manipulated.
FIG. 11 is a round-trip, time-of-flight-distance diagram 600. Shown are the round-trip times-of-flight 602A through 602E for selected propagation pathways from the aperture antennas 520A through 520E. Assuming that the velocity of propagation is known, the length of the time of each flight line or pathway 602A-602E also represents twice the distance between the aperture antennas 520A-520E and a scatterer at a point (voxel) 604. FIG. 12 shows a similar time-of-flight distance diagram 610 for a different point (voxel) 614. Assuming that the velocity of propagation is known, each of the round-trip time-of-flight pathways 610A-610E also represents twice the distance from each of the aperture antennas 520A-520E to a different scatterer positioned at point 614, FIG. 12.
FIGS. 13 and 14 show how the observed round-trip time-of-flight data can be processed such that the returns from point 604 (FIG. 11 and FIG. 13) are enhanced and returns elsewhere, whereas returns from point 614 (FIGS. 12 and 14) are suppressed. On the assumption that the propagation velocity is known, the processor 70 (FIG. 2 or FIG. 10) selects a location or voxel, such as 604, and then adds a predetermined additional time 620A-620E (FIGS. 13 and 14) to the round-trip time histories of each of the backscatter returns. The time added, 620A through 620E, is equal to the estimated round-trip time-of-flight from the most distant aperture less the estimated round-trip time-of-flight from the nearest aperture plus a fixed buffer time 622. These estimated times are based on twice the distance between the aperture and the selected voxel divided by the assumed propagation velocity.
FIG. 13 presents data based on assuming the correct propagation velocity. Chart 619 of FIG. 13, beginning at time 624, presents lines whose lengths each include one of the predetermined additional times 620A-620E. These predetermined times are added linearly to the observed round-trip times 602A-602E as well as a fixed buffer time 622. If the propagation velocity is correctly chosen, the sum of the times (or line lengths in the chart 619, FIG. 13) will be equal to the interval 628 between the starting point 624 and the voxel point 604. When this occurs, the returns from the selected voxel 604 for each of the apertures add constructively, thereby suppressing returns from voxels not at the selected synthetic focal point.
FIG. 14 illustrates a chart 640 that is similar to chart 619, except the returns from voxel 614 (FIG. 12) are added to the estimated round-trip time-of-flights from voxel 604 (FIG. 11). The sum of these times (line lengths) are not equal. Some have small values, whereas others are greater. As a consequence, the returns from voxel 614 do not add constructively if estimated time-of-flights are based on the position of voxel 604. Note that the dash line 630 in FIGS. 13 and 14 has the cross-sectional configuration of a convex lens. In the case of a lens that focuses rays from the sun into a spot, there is just one focal point. To focus the lens at some other point, the lens must be moved, or its curvature altered. Likewise, the various estimated time delays 620A-620E must be changed to focus the synthetic lens other than at the position of voxel 604.
FIG. 15 illustrates how the time histories of return (backscatter) signals from each of the antennas 520A through 520E can be combined to enhance the returns from the presumed tumor in the preselected voxel 604 over those from the nearby voxel 614, in the case where both voxels 604 and 614 are present; see FIGS. 11 and 12. For illustrative purposes, the return from the scatterer at each of the points 604 and 614 is assumed to be in the form of a one-cycle sine wave. It is also assumed that the path attenuation and divergence are normalized.
In FIG. 15, the time histories 642A through 642B for each of the antennas 520A through 520B, respectively, are shown at the top of the chart; the combined output of the antennas is illustrated by the curve 642 at the bottom of FIG. 15. Each time history is plotted as a function of time. The individual waveforms from the scatterer at point 604 (FIG. 11) are shown as single-cycle sine waves 648A through 648E, all of which are coincident in time. In a similar fashion, the waveforms from the scatterer at voxel 614 (FIG. 12) are shown as the single-cycle sine waves 650A through 650E, which are not coincident in time. See also FIGS. 13 and 14.
The curve 642 at the bottom of FIG. 15 summarizes the five individual time history curves 642A-642E. As before, curve 642 is plotted as a function of time, beginning at a time 624 taken as time zero. The returns 648A-648E (from voxel 604) all occur at about the same time, and coherently combine to afford a pulse 668 having an amplitude that is about five times greater than the amplitude of the return (648A-E) from each antenna. Because the returns 650A through 650E (from voxel 614) do not occur at the same time, and hence do not add coherently, they form a series of returns of smaller amplitude. These and the backscatter returns from other voxels are used to develop a three-dimensional image as suggested in FIG. 1B.
It cannot be assumed that the propagation velocity in the breasts of different women will be the same. One procedure would be to assume a likely overall average value and then observe the resulting images, especially for a known scatterer, such as the chest wall. Next, additional values for the overall average propagation can be selected and the resulting effect on image clarity noted, especially for known scatterers. Incremental changes in the propagation velocity should be changed in the direction that produces the greatest improvement. Such a procedure may lead to dithering the changes around a point of improved clarity or amplitude of returns. To fine tune the results, the propagation time from each aperture to a specific voxel, an interesting abnormality or regions of voxels can be changed or dithered, again initially choosing incremental changes in a direction that leads to the greatest improvement.
Recent studies demonstrate that the applicator-skin-breast tissue interface can introduce unexpected difficulties. The relative dielectric constant of the skin is in the order of 36 with a loss tangent of about 0.4 (as noted by Stuchly, 1980). The thickness of the skin of the breast can be expected to vary, with the most variation in the vicinity of the nipple. In addition, the skin of breast from different patients can also be expected to vary. The wave impedance of the skin is approximately one-half of that for the normal breast tissue. The velocity of propagation through the skin is also about one-half of that for normal breast tissue. The two-to-one difference in wave impedance results in substantial reflections. The slower velocity results in in-skin wavelengths of the order of 5 mm. Skin thicknesses as small as 1.2 mm can result in quarter wavelength transmission or refection enhancements, depending on the nature of the interface. These effects tend to reduce the signal transmitted into and out of the normal breast tissues. In addition, reflections may tend to obscure the desired returns.
Methods and means are included here to overcome such difficulties. The reason is that the combined effect of both the thickness and the dielectric properties of the skin, especially the skin of the human breast around the nipple, has not been fully appreciated. This lack of appreciation is compounded because very little has been published on the dielectric properties of human skin.
It has been thought that the skin would not introduce a significant perturbation. However, simply matching the dielectric constant or wave impedance in the antennas within an illuminator to the dielectric constant or wave impedance of the normal breast tissue, or the skin of the breast, may not be adequate for high resolution and accurate imaging. On the other hand, if detection of just abnormalities is the goal for mass screening applications, such a consideration becomes less crucial and may possibly be ignored.
The methods described in connection with the systems depicted in FIGS. 2, 7, 9 and 10 indicate that the backscatter introduced by various interfaces, including that of the skin, could be made evident in the time-history of the backscattered returns. By so doing, the primary effect of the backscatter from the skin interfaces may be eliminated by gating or other processing. However, if more precision and resolution is needed, additional steps can be taken. One method would be to estimate the dielectric properties and thickness of each area of the skin that is in immediate contact with each aperture-antenna. The smaller wavelengths are used to help resolve the aperture-skin interface returns from the skin-breast/tissue interface returns. The reflection coefficients introduced at each of the aperture-skin interfaces are first developed. From such reflection data, the approximate dielectric properties, wave impedance, and propagation velocity within the skin can be estimated. The round-trip time between the aperture-skin and the skin-tissue interfaces can be determined from the time-history of the returns from each aperture. Knowing this round-trip time, the distance from the aperture-skin interface to the skin-breast/tissue interface can then be estimated. A possible next step would be to estimate the dielectric properties of the breast tissues voxels that are nearest to each of the apertures. This can be clone by developing the reflection coefficient from the observed return from the skin-breast/tissue interface.
This general method can be better understood by referring to the diagram 700, FIG. 16, which shows the aperture-skin interface 703 and skin-breast/tissue interface 705 between the dielectric or wave impedance properties within the aperture region 702, the skin region 704, and the region of the normal breast tissue 706. FIG. 16 depicts the case for perpendicular incidence of an electromagnetic wave 708 that impinges upon the skin 704 of the breast from the aperture region 702. On encountering the interface 703, a reflected wave 712 and a transmitted wave 710 are produced. The transmitted wave 710 passes through the skin 704 and is partially attenuated, such that the remaining portion 714 encounters the skin-breast/tissue interface 705. As before, a reflected wave 718 and a transmitted wave 716 are produced. The reflected wave 718 returns back through the skin 704 such that the remaining attenuated portion 720 encounters the aperture-skin interface 703. As before a transmitted wave 722 and a reflected wave 724 are produced. The transmitted wave 722 is sensed by the aperture antenna and becomes part of the time-amplitude backscatter return.
To estimate the approximate dielectric properties of the skin, the reflection coefficient, ρ, is the ratio of the amplitude of the reflected wave as noted in the time-amplitude return to the impinging wave (also known from calibration measurements) in terms of the wave impedance n1 of the source region 702 (the aperture region).
The wave impedance of the adjacent skin region n2 can be estimated as follows: ##EQU2##
ρ=(n.sub.2 -n.sub.1)/(n.sub.2 +n.sub.1),
n.sub.2 =(μ.sub.2 /ε.sub.2).sup.1/2
See the prior text immediately preceding equations (1) and (2).
The unknown value of n2 may be estimated in terms of the known values of n1 and ρ.
The transmission coefficient, τ, is the ratio of the transmitted wave to the impinging wave as follows:
τ=(2n.sub.2)/(n.sub.2 +n.sub.1) (11)
from the calculated value of n2 and the known value of n1.
Similarly, the transmitted and reflected waves at the other interfaces indicated in FIG. 16 can also be estimated and used to determine the approximate value of the dielectric parameters of the normal breast tissue near different specific apertures. See Ramo (1965, Chapter 6) for more complete discussions on reflection and transmission coefficients.
Knowledge of the dielectric parameters and thickness of the skin is helpful to compensate for refraction and reverberation effects that may be introduced from non-orthogonal propagation pathways. FIG. 17 illustrates the case where the impinging wave 768, in the plane of the paper, strikes the interface 763 at an angle 771. Media 762 and 766, corresponding to media 702 and 706 in FIG. 16, have a much smaller dielectric constant than the media 764 (704 in FIG. 16) that represents the skin. That is, the incident wave 768 strikes the interface 763 between media 762 and 764; the skin, medium 764, is assumed to have a much higher dielectric constant than the aperture region 762. A reflected ray 772 leaves the interface 763 at the same angle 771 to the normal to the surface that is formed by the incident ray 768. The transmitted ray 770, however, is refracted so that its angle 771' to the normal is smaller than the related angle 771 for the incident ray 768. The transmitted ray proceeds along the path 774 until the skin-breast tissue interface 765 is encountered. At interface 765, a reflected ray 778 and a transmitted ray 776 are formed. Reflected ray 778 proceeds along a path 780 until it encounters interface 763, where a reflected ray 784 and a refracted ray 782 are again developed and progress respectively along the paths 792 and 794. A key point is that the expected pathway of the incident wave 768 beyond interface 763, without accounting for the refraction effects of the skin, would be different than the actual pathway. In addition, reverberation effects within the skin layer may also occur.
The lossy dielectric approximation that is used above with equations (1) and (2), also applies for the foregoing interface relationships. This is justified because the measured (as published, Stuchly 1980) loss-tangent for skin is substantially less than one. Computer aided iterative approaches can be employed to develop more precise values As the first step, the data developed by the aforementioned approximate lossy dielectric method is used. Then various values for the conductivity or dielectric constant are varied in an iterative maimer. This is done until the best match between the theoretically predicted response (based on the best iterated estimate for the dielectric parameters) and measured response is realized.
Several different applications for the above discussed apparatus and methods are possible. If detection and imaging of the smallest incipient malignant tumors is desired, the number of antenna apertures should be large, the overall dimension of the array of the apertures should be large, and the greatest bandwidth, the highest possible frequency and the smallest rise time pulses should be used. These requirements are needed to assure high resolution and precise location. However, if the goal is just to detect an abnormality for mass screening, then the number of apertures can be reduced and the area of the array can be reduced. The maximum frequency and bandwidth can be reduced and longer rise time pulses can be tolerated.
To determine whether or not a palpable tumor is benign or malignant, the number of apertures can be reduced and the array can be reduced, to perhaps just one aperture. In such cases as this, the beam may be characterized as tubular, and the necessary information developed from a very short rise time pulse, from either a synthetic or a real time system.
Several types of aperture antennas can be considered. To take best advantage of the synthetic pulse or the real time step function or short duration pulse, a broad band radiating antenna is desirable. The abrupt termination of the double-ridge guides (FIG. 8) into high dielectric constant material, such as the skin of the human breast, may in some cases not produce the desired frequency response, especially at the lower frequencies. One broader band radiation structure is the tapered, dihedral horn 801 illustrated in the isometric view of FIG. 18. The overall arrangement of horn 801 comprises two conductive sheets 802 and 804 each constituting an isosceles triangle. The bases of these triangles are arranged to form an aperture 805 and the apexes are terminated near a position 807, with one sheet's apex connected to the plus phased point of a balanced feed waveguide 806 and other sheet's apex connected to the negative phased point of the balanced feed waveguide 806. A balun 808 converts the balanced system for coupling to the more commonly used unbalanced waveguide 810, such as a 3 mm diameter coaxial cable. The aperture antenna 801 can then be encapsulated in a gel or plastic that has the required dielectric properties to match either those for the skin or for the breast tissue. The dihedral arrangement shown in FIG. 18 can be used in an array configuration similar to that shown in FIG. 8.
FIG. 19 illustrates an apparatus 900 that may use the dihedral antenna of FIG. 18 (or an equivalent broad band antenna) for the purpose of determining whether a tumor that was detected by mammography or by palpitation is benign or malignant. Precise imaging may not be needed, and therefore the wide apertures and short focal length needed to develop a compact focal point may not be required. Fewer antennas may be needed, along with reduced overall size of the aperture. Such simplified arrangements could also be used to help resolve indeterminate images that are sometimes found in mammograms. In many cases, benign tumors contain less water than malignant tumors. Further, many benign tumors do not encourage the development of small blood vessels, such as those that usually surround a malignant tumor. The presence of blood vessels that surround a malignant tumor increases the backscatter, because blood has a high water content.
A non-conductive template 906 with plural slots 905 is placed on the breast 902 of a human or animal patient (not completely shown). The slots 905 of template 906, when contacting the breast, cause the skin of the breast to pucker up within each slot 905. The system 900 includes an illuminator 922 that may contain a dihedral aperture antenna, such as that illustrated in FIG. 18. The shape of the plastic molding 920 that contains the dihedral antenna is keyed to fit into the slots 905 in template 906 and to contact the skin of the breast firmly. The system 900 also includes a unit 910 that includes operating components similar to those shown in FIG. 2, including a microwave source, a directional coupler or equivalent, a backscatter detector, a signal processor, and a display. The system 900 of FIG. 19 is designed to detect the antenna-skin and the skin-tissue interfaces, as well as backscatter. The backscatter from each of the slots 905 in template 906 is measured by placing the illuminator 920 into each slot. The returns are stored and processed as described for FIGS. 10, 11, 12, 13, 14 and 15. The switching functions shown in FIG. 10 are replaced by manually moving the illuminator 922 to each slot position while at the same time noting which slot positions are employed. The material of template 906 should be a plastic or other material that has dielectric properties similar to those of skin or of normal breast tissue when a dihedral antenna, such as shown in FIG. 18, is used. Alternatively, the template 906 could be metal, if an aperture antenna such as shown in FIG. 8 is used.
The antenna arrangement shown in FIG. 18 can be improved by replacing the single dipole dihedral horn 805 of FIG. 18 with a Maltese Cross dipole arrangement 850 as shown in FIG. 20. In the antenna arrangement 850, there are two "bow-tie" antennas. The first antenna 852A and 852B is a horizontal receiving antenna. The other "bow-tie" antenna 853A and 853A is a vertical illumination antenna. Both of these two "bow-tie" antennas are supported on a planar substrate 851 that in turn is mounted upon a housing 855. The central apexes of the two bow-tie antennas are connected, respectively, to twin line balanced cables. The cable connections are not shown in FIG. 20. The apexes of the vertical illumination antenna 853A-853B are connected to one of the two cables and the apexes of the horizontal bow-tie antenna are connected to the other cable. The two cables are converted by a balun housed within housing 855 to unbalanced coaxial cables 857 and 856 respectively.
FIGS. 21A and 21B illustrate conceptually how the cross-polarized dipoles of FIG. 20 suppress co-polar returns. The waveform in the co-polar horizontal illumination antenna elements 852A and 852B are suppressed in the cross polarized vertical receiving antenna 853A and 853B. Balanced lines 874A and 874B from a balanced receiver source 856 are connected to antenna elements 852A and 852B, respectively. Lines 873A and 873B from a balanced source 855 are connected to both elements 853A and 853B of the other antenna. Mutual coupling impedances 830A, 830B, 830C and 830D are shown; they can provide a mechanism for the power in the illuminating antenna to spill over into the cross-polar receiving antenna circuit. If the geometry and distribution of the dielectric media surrounding the axes of the two dipole antennas is symmetrical, the waveform in the co-polar antenna 852A, 852B is completely suppressed in the cross-polar antenna 853A, 853B. However, should an asymmetry occur, such as the presence of a tumor that is off axis, then the cross-polar antenna will collect the backscatter from the tumor free of any wave forms in the co-polar circuits. On the other hand, if the asymmetry is in the dielectric media, spillover will occur. This is conceptually suggested by the Wheatstone Bridge circuit 890 of FIG. 21A. So long as the impedances 892, 893, 894 and 895 are identical, no interaction between sources 891 and 896, which correspond to the two antennas, will occur.
Other antenna or signal processing arrangements are possible to suppress certain common features, such as planar surfaces. Alternatively, a different arrangement can be used to enhance the detection of a unique geometric feature, such as the edges of a hockey-puck-shaped buried anti-personnel land mine. An example of this (Ogus, U. and L. Gurel 1997) is the use of three identical bow-tie co-polar antennas, wherein the center of the bow-tie is impulsively excited. The backscatter response is collected by the left and right antennas, and these two responses are fed into a hybrid network that forms the difference between these two waveforms. These arrangements reject the co-polar, bi-static cross-talk that have even left to right symmetry while enhancing the response to edge-type asymmetries near one of the outer antennas.
A conceptual system 1000, shown in FIG. 22, utilizes a manually operated synthetic aperture arrangement that employs the improved crossed dipole sensor of FIGS. 20 and 21A. A broad-band pulse source 1001 provide a series of picosecond rise time pulses at a one MHz rate. Via a coaxial cable 1002, the pulse source 1001 supplies those pulses (arrow 1003) to a cable bundling box 1004. Cable 1002 is combined with another coaxial cable 1013 in a cable bundle 1005 that carries the pulses 1003 to a balun in an applicator 1006. The balun in applicator 1006 converts the unbalanced mode to the balanced mode to supply energy to the co-polar antennas of a Maltese Cross antenna 850 located at the bottom of the keyed section 1008.
The breast or other tissue to be examined is scanned by manually moving the applicator 1006 such that the Maltese Cross antenna 850 periodically contacts the breast 1010 in a manner similar to that described for FIG. 19. The position of the keyed section 1008 and Maltese Cross antenna 850 is confined by a series of slots 1031 in a guide plate 1011. The location of section 1008 in the guide plate 1011 may be controlled by a series of detents. Fiber optical sensors (not shown) in the guide plate may provide additional control via a cable optic bundle 1030 connected to the digital processor 1022 of system 1000.
Cross-polarized backscatter is collected by the cross-pol antenna in device 850. The balanced voltages are converted to unbalanced return voltages 1012 suitable for use with coaxial cables by the balun in unit 1006. Via the additional cable in cable bundle 1005, the backscatter voltages 1012 are applied to the cable bundling box 1004 and thence carried by the separated coaxial cable 1013 to a preamplifier 1015. The amplified signals are impressed upon a wave shaping network 1017, which supplies a modified wave form 1018 via a cable 1019 to a digital sampling scope 1020.
The digital sampling scope 1020 converts the time history of the backscatter returns, that takes place over about a 1000 picosecond interval, to a similar time history that takes place over a microsecond duration. For each position of the applicator, the time histories of about five thousand returns are integrated to improve the signal to noise ratio and to produce a integrated time history of the return every five milliseconds. The digital sampling scope 1020 is synchronized to the repetition rate from signals from the pulse source 1001.
The location of the pulse shaping network 1017, if needed, should be after the preamplifier 1015. If the location of the pulse shaper is before the preamplifier, it will reduce the amplitude by between 10 to 20 dB. In marginal situations, this may suppress the desired signal into the noise level.
The output of the digital sampling scope 1020 is applied to the digital processor 1022 via a cable 1021. The digital processor develops a three dimensional image of the interior of the breast as described for FIGS. 10 through 15. The principal difference is that the cross-polar returns 77 shown in FIG. 2 are supplied directly to the sampling scope and not through the directional coupler. This may require an additional switching circuit similar to that shown in FIG. 10, except that the cross-pol antennas would be switched in conjunction with the co-pol switches, but the returns would be applied directly to the digital sampling scope.
Signals from digital signal processor 1022, via cable bundle 223, are sent to a display unit 1024. Data from past exams is stored in a digital store 1025 and can be accessed, using a connector cable 1026 to determine if major changes have taken place.
The manually operated synthetic aperture system shown in FIG. 22 can be converted into a electronically controlled antenna array. This can be done by replacing the guide plate 1011 and the applicator assembly 1006 by the antenna array shown in FIG. 23. Pulsed signals 1003 are supplied to electronically controlled switches (not shown) in a housing 1054 via a given switch configuration, a connection is made via each of the baluns (not shown) in a housing 1052 to the co-pol or the cross-pol antenna. The Maltese (ross antennas are formed on a substrate 1052 in a manner similar to that described for FIG. 20. An opening 1055 in the substrate 1052 is made to avoid compressing he nipple of the person being tested.
Computer analyses have demonstrated the benefit of using the cross-dipole techniques just described. Under idealized conditions, the cross-talk between the co-pol and the cross-pol returns is completely rejected. Any perturbation introduced by planar surfaces, such as skin or the muscles that surround the chest wall, are also suppressed. At the same time usable returns are realized from a spherical off-axis tumor or an asymmetrical axis tumor near the chest wall.
However, this idealized situation may not occur, and additional techniques must be employed. One is the use of a resistively loaded bow-tie antenna shown in FIGS. 22A and 22B. FIG. 22A illustrates the front view 1100 of the bow-tie antenna. It consists of two triangular conducting surfaces, 1101 and 1102. These surfaces can be formed by silk-screening a mixture of metallic and graphite particles on a substrate or lossy dielectric backing that has the same dielectric properties as normal breast tissue. The surface resistance of the each section of the bow-tie 100, as expressed in ohms/square, varies from a very low value near the apexes 1104 and 1005 to a very large value near the distal ends 1106 and 1107. The apexes are connected to a balance transmission line, such as discussed for FIGS. 21A and 21B. The side view antenna 1100 in FIG. 24B shows the both sections 1111 and 1112 of the bow-tie antenna in direct contact with the skin 1108 that is over the normal breast tissue 1109. The bow-tie elements are silk-screened onto a lossy dielectric substrate 1103. This substrate is on the side of the bow-tie that is opposite to the breast.
When a pulse is applied to the balanced transmission line that feeds the bow-tie antenna 1100, two effects occur that must be mitigated, assuming that these antennas will also collect the backscatter in the monostatic mode. FIG. 25 shows the time history of the backscatter where the amplitude 1120 is plotted as a function of time 1122 when a short duration pulse is applied to the antenna shown in FIG. 24A and 24B. The first reflection or backscatter 1125 occurs because the wave impedance of the antenna, when combined with skin and breast tissue, is different from the wave impedance of the transmission line for the antenna. This can be gated out as discussed earlier, but large values can complicate the design of the gating circuits. The second effect arises from the backscatter from the distal ends 1106 and 1107 of the bow-tie antenna. This second effect is more serious because this distal end backscatter (1124 in FIG. 25) can occur in a time interval when the backscattered returns from a possible tumor might arrive. Both effects can be suppressed by the use of the graded resistance method discussed above. Additionally, the mismatch between the transmission line and antenna can be further mitigated by varying the angle at the apex of the antenna.
Studies have been made regarding the benefits of the graded resistance loading for bow-tie antenna (Slager, Smith and Maloney, 1994). For bow-tie antennas in air that use all metal antennas, the backscatter from 1123 and 1124 are -20 dB and -14 dB, respectively, below the peak value of the incident pulse. If a graded resistance design is employed, both values are reduced to about -40 dB to -50 dB below the peak value of the incident pulse. These reflection levels are high and could mask the desired returns.
To reduce these backscatter signals to a more tractable level for breast cancer detection, a graded resistance antenna was analyzed similar to that shown in FIGS. 24A and 24B. That antenna was placed upon the breast, except that the dielectric 1103 (FIG. 24B) was air. The results were disappointing. Next, the air dielectric was replaced by a material that had dielectric characteristics similar to those found for normal breast tissue. This resulted in the distal end reflections 1124 being reduced to -106 dB below the exciting pulse.
Such a large reduction was unexpected. The key feature for this redaction was the use of a dielectric material 1103 similar to normal breast tissue.
As defined in FIG. 24A, the following resistive loading taper in ohms per square was used: where R0 =3.5×10-5, R1/2 =110, z/h is the normalized distance along the bow-tie from the apex to the distal end and R1/2 is the surface resistance when z/h=0.5. ##EQU3## ohms/square. Depending on the height, h, of the dipole, values of R1/2 can range between 10 to 200 ohms per square.
Other methods are possible to reduce the distal end reflection, such as the use of material that has a fixed value of R0 but where serrations are employed as illustrated in FIG. 25, which otherwise is the same as FIG. 24A. In any case, the key design criterion is backing the bow-tie antenna with material that has dielectric characteristic similar to that found for the human breast.
Combining the features of the above resistively loaded dipole into the Maltese Cross configuration is beneficial. While the crossed-dipole arrangement can radically suppress the undersized returns, 1125 and 1124 shown in FIG. 25 under idealized conditions, such cannot be expected in practice. The incorporation of the resistively loaded dipole into the Maltese Cross arrangements illustrated in FIGS. 20-23 provides an extra margin of safety to suppress reverberation effects in general and more specifically the backscatter from the distal ends of the bow-tie antenna.
Throughout the discussion of pulse-type systems, methods are also incorporated to suppress reverberation between interfaces within the system. Examples of the use of waveguides that are somewhat lossy include the waveguides 509 and 519 in FIG. 10 and the flexible waveguide or coaxial cable 912 of FIG. 19. In addition, lumped or distributed lossy or reactive elements can be located near or within system interfaces, such as the switches shown in FIG. 10, to further suppress unwanted returns from incidental interfaces. Further electrical parameters of the material used for the waveguides or the mouths of the apertures can be progressively changed to suppress unwanted reflections or to change the radiation pattern of individual aperture antennas.
Other applications of the above described techniques are possible. For example, benign tumors can be distinguished from malignant tumors by noting the amplitude of the backscatter returns. The dielectric constant of malignant tumors is in the order of 40 to 55 whereas a lipoma (a generally benign tumor) has a dielectric constant in the order of 20 (Foster 1981).
Persons skilled in the art can also modify the above methods to measure selected volumes of human or animal bodies, such as to determine the thickness of the fat layer immediately beneath the skin. Alternatively, the methods may be used to detect tumors or other types of abnormalities that are embedded in muscle or other high water content tissue. The dielectric constant of normal muscle tissue is large, in the order of 40 to 50, and its conductivity is also large. This results in depths of penetration that are too small at frequencies suitable for penetrating breast tissues. As a consequence, the operating frequency must be reduced to the 1 GHz lower level as opposed to the approximately 3 GHz lower level noted for breast tissues. Fortunately, the wavelength at a given frequency in muscle is about one-half of that for a similar frequency in normal breast tissue.
The illuminating waveform can be different from the described impulse or swept frequency waveforms used here. For example, the amplitude of the higher frequency spectral components of the illuminating wave can be enhanced to compensate for dispersion effects. Alternatively, the spectral distribution of the illuminating waveform can be altered to improve the signal-to-noise ratio or to enhance the detection of a specific form of breast cancer. The tapered, dihedral-horn dipole of FIG. 18 can be used as part of a dual dihedral-horn to form a cross-dipole arrangement. Further, the spatial distribution of the resistivity of the conducting portions of these dipoles can also be tapered to suppress reverberation. The way in which the resistivity of the conducting portions is varied is not limited to the techniques considered for FIGS. 24A, 24B and 26, so long as the resistivity is low at the feed point and continuously increases to a large value at the distal ends of the antenna.
Throughout the foregoing discussion of pulse-type systems, methods are also incorporated to suppress reverberation between interfaces within the system. For example, waveguides that are somewhat lossy may be used, such as in the case of waveguides 508 and 519 in FIG. 10 or the flexible waveguide or coaxial cable 912 of FIG. 19. In addition, lumped or distributed elements can be located near or within system interfaces, such as the switches shown in FIG. 10, to further suppress the unwanted returns from such incidental interfaces. Further, electrical parameters of the material used for the waveguides or the mouths of the apertures can be tapered (progressively changed) to suppress unwanted reflections or to change the radiation pattern of individual aperture antennas.
Throughout the foregoing specification and in the appended claims the terms millimeter waves, or mm, waves or mmw have been used to generically represent the wavelengths of the electromagnetic waves that propagate in the human breast tissue. Since the relative dielectric constant of the human breast is in the order of 9 to 20, the free-space wave length will be reduced by a factor of three or more. Thus, the in-tissue wavelengths in the breast, over a frequency range of 2 to 60 GHz, will range from about 30 mm to 1 mm. For higher water content tissues, such as muscle, the frequency range is reduced by a factor of about three. In addition, the backscatter returns can be mathematically characterized in the time domain or the frequency domain.
As opposed to certain microwave hypothermia cancer treatment technology, none of the technology presented here is intended to heat significantly any portion of the breast. This requirement limits the power deposition density onto the surface of the breast to less than 10 milliwatts/cm2 or the volumetric heating rate in any portion of the breast to less than 0.8 milliwatts per gram of tissue as averaged over a time period of about six minutes. To further assure minimal thermal effects, the input power is turned off if the scanning system falters for any reason.
Other usages are as follows: The term impedance refers to the ratio of the voltage to the current or to the electric field to the magnetic field at a specified location. This term impedance is qualified as "electrical" or "wave" respectively, depending on whether voltages and currents or electromagnetic fields are concerned. The term wave guide is used in the generic sense and includes both cables and higher mode wave guides with just a single transverse field. The terms effective aperture and effective focal point are used in the generic sense wherein apertures and focal points can be created physically or synthetically (such as often done in synthetic aperture radar).
The effective focal point is not a point but rather is defined here as a region where the illuminating energy is most concentrated in the breast. The effective focal point is further defined as the region or volume where this energy concentration occurs as affected by the heterogeneity of dielectric characteristics of the normal breast tissues, the in-tissue wavelength, the size and distance of the illuminating globular aperture or the geometry and number of apertures used in a phased array. The focal point positioning may be either mechanical or electronic as in the case of a phased array.
The terms "detect" or "detection" are also used in the generic sense, and may mean simply indicating the presence of a tumor or more broadly providing data that permits imaging the location, size and geometry of the tumor. Detecting, identifying, imaging or locating a tumor also means noting the presence of any abnormality. The terms "power and signal director" or "input power and signal separation" are also used in a generic sense. Both passive and active techniques not only enhance detection by suppressing the direct effects of impinging power waves, but also can reduce false signals or clutter. Such are introduced by imperfect matches between impedances or by non-tumor scattering sources, such as the breast/lung interface. In the appended claims, and in this specification, the term "time delay" includes phase delay unless otherwise indicated.
The following references are of utility in understanding the foregoing specification:
Burdette, E. C., et. al. (1980): In vivo measurement techniques for determining dielectric properties at VHF through microwave frequencies, IEEE Trans., MTT, Vol MTT-28, No. 4 April, pp. 414-427.
Burdette, E. C., et. al. (1986): In situ permittivity at microwave frequencies: perspective, techniques, results, medical applications of microwave imaging, Medical Applications of Microwave Imaging, Larsen, L. E. and J. H. Jacobi, IEEE Press pp. 13-40.
Chaudhary, S. S., et. al. (1984): Dielectric properties of normal and malignant human breast tissues at radiowave and microwave frequencies, Indian Journal of Biochemistry and Biophysics, Vol. 21, February pp. 76-79.
Edrich, J., et. al. (1976): Complex permittivity and penetration depth of muscle and fat tissues between 40 and 90 GHz, (1976) IEEE Trans., MTT, vol. MTT-24, May pp 273-275.
Johnson, E. C., et. al. (1972): Nonionizing electromagnetic wave effects in biological materials and systems, Proceedings of the IEEE, Vol. 60, No. 6, June pp. 694-695.
Kay, A. F. (1966): Millimeter wave antennas, Proceedings of the IEEE, Vol. 54, No. 4, pp. 641-647.
Larson, E. E. and J. H. Jacobi, Eds. (1986): Medical Applications of Microwave Imaging, IEEE Press, Institute of Electrical and Electronic Engineers, New York.
Ramo, S., et. al. (1965): Fields and Waves in Communication Electronics, John Wiley and Sons, New York.
Rogers, J. A., et. al. (1983): The dielectric properties of normal and tumor mouse tissue between 50 MHz and 10 GHz, British Journal of Radiology, vol. 56, May, pp. 335-338.
Smith, W. J., (1966): Modern Optical Engineering, McGraw-Hill, New York, N.Y.
Stuchly, M. A. et. al. (1980): Dielectric properties of biological tissues, Journal of Microwave Power, 15(1), pp. 19-26.
Terris, M. K. et. al. (1995): Prediction of prostate cancer volume using prostate-specific antigen levels, transrectal ultrasound, and systematic sextant biopsies, Urology, January Volume 45 No. 1 pp. 75-80.
Rogers, J. A. et. al. (1983): The dielectric properties of normal and tumor mouse tissue between 50 MHz and 10 GHz. The British Journal of Radiology, 56, pp. 335-338.
Foster, K. R. et. al. (1981) Dielectric properties of tumor and normal tissues at radio through microwave frequencies, Journal of Microwave Power, 16(2), pp. 107-119.
Shlager, K. L., et. al. Optimization of bow-tie antennas for pulse radiation IEEE trans. Antennas and Propagation, 1994, 42 (7) pp. 975-982.
Oquz, U. and L. Gurel Subsurface scattering calculations via the 3D FDTD method employing PML-ABC for layered models. Paper 97.9 at 1997 IEEE AP-S International Symposium Montreal Quebec July 1997.
|1||*||ALbert K. Y. Lai, Albert L. Sinopoli, and Walter D. Burnside, A Novel Antenna for Ultra Wide Band Applications, IEEE Transactions, on Antennas and Propagation, vol. 40, No. 7, Jul. 1992.|
|2||ALbert K. Y. Lai, Albert L. Sinopoli, and Walter D. Burnside, A Novel Antenna for Ultra-Wide-Band Applications, IEEE Transactions, on Antennas and Propagation, vol. 40, No. 7, Jul. 1992.|
|3||*||James G. Maloney and Glenn S. Smith, Optimization of a Conical Antenna for Pulse Radiation: An Efficient Design Using Resistive Loading, IEEE Transactions on Antennas and Propagation, vol. 41, No. 7, Jul. 1993.|
|4||*||Kurt L. Shlager, Glenn S. Smith and James G. Maloney, Optimization of Bow Tie Antennas for Pulse Radiation, IEEE Transactions on Antennas and Propagation, vol. 42, No. 7, Jul. 1994.|
|5||Kurt L. Shlager, Glenn S. Smith and James G. Maloney, Optimization of Bow-Tie Antennas for Pulse Radiation, IEEE Transactions on Antennas and Propagation, vol. 42, No. 7, Jul. 1994.|
|6||Motohisa Kanda, The Effects of Resistive Loading of "TEM" Horns, IEEE Transactions on Electromagnetic Compatibility, vol. EMC-24, No. 2, May 1982.|
|7||*||Motohisa Kanda, The Effects of Resistive Loading of TEM Horns, IEEE Transactions on Electromagnetic Compatibility, vol. EMC 24, No. 2, May 1982.|
|8||*||Motoshisa Kanda, IEEE Transactions on Antennas and Propagation, vol. AP 28, No. 1, Jan. 1980.|
|9||Motoshisa Kanda, IEEE Transactions on Antennas and Propagation, vol. AP-28, No. 1, Jan. 1980.|
|10||*||T. T. WU and R. W. P. King, The Cylindrical Antenna with Nonreflecting Resistive Loading, IEEE Transactions on Antennas and Propagation, The imperfectly conducting cylindrical antenna, Tech Rept 440, Cruft Lab., Harvard University, Cambridge, Mass., Mar. 1964.|
|Citing Patent||Filing date||Publication date||Applicant||Title|
|US6543933 *||16 Jan 2002||8 Apr 2003||Her Majesty The Queen As Represented By The Minister Of National Defence Of Her Majesty's Canadian Government||Non-invasive 3-D intracranial thermography system|
|US6751500||3 Sep 2002||15 Jun 2004||Medrad, Inc.||Apparatuses and methods for extravasation detection|
|US6813515||4 Jan 2002||2 Nov 2004||Dune Medical Devices Ltd.||Method and system for examining tissue according to the dielectric properties thereof|
|US6927691 *||25 Nov 2002||9 Aug 2005||Spatial Dynamics, Ltd.||Dielectric personnel scanning|
|US7006676 *||17 Mar 2000||28 Feb 2006||Medical Optical Imaging, Inc.||Method and apparatus for detecting an abnormality within a host medium utilizing frequency-swept modulation diffusion tomography|
|US7047058||30 Jan 2002||16 May 2006||Medrad, Inc.||Apparatuses, systems and methods for extravasation detection|
|US7122012||26 Jul 2002||17 Oct 2006||Medrad, Inc.||Detection of fluids in tissue|
|US7184824||18 Oct 2004||27 Feb 2007||Dune Medical Devices Ltd.||Method and system for examining tissue according to the dielectric properties thereof|
|US7298536||7 Sep 2005||20 Nov 2007||Cascade Microtech, Inc.||Fiber optic wafer probe|
|US7340292 *||14 Feb 2004||4 Mar 2008||University Of Florida Research Foundation, Inc.||Breast cancer detection system|
|US7355420||19 Aug 2002||8 Apr 2008||Cascade Microtech, Inc.||Membrane probing system|
|US7420381||8 Sep 2005||2 Sep 2008||Cascade Microtech, Inc.||Double sided probing structures|
|US7454242||17 Sep 2004||18 Nov 2008||Elise Fear||Tissue sensing adaptive radar imaging for breast tumor detection|
|US7474257 *||7 Nov 2005||6 Jan 2009||The United States Of America As Represented By The Secretary Of The Navy||Multistatic adaptive pulse compression method and system|
|US7492172||21 Apr 2004||17 Feb 2009||Cascade Microtech, Inc.||Chuck for holding a device under test|
|US7492175||10 Jan 2008||17 Feb 2009||Cascade Microtech, Inc.||Membrane probing system|
|US7505811||18 Nov 2002||17 Mar 2009||Dune Medical Devices Ltd.||Method and apparatus for examining tissue for predefined target cells, particularly cancerous cells, and a probe useful in such method and apparatus|
|US7570063 *||3 Jul 2002||4 Aug 2009||Wisconsin Alumni Research Foundation||Space-time microwave imaging for cancer detection|
|US7591792||26 Jul 2002||22 Sep 2009||Medrad, Inc.||Electromagnetic sensors for biological tissue applications and methods for their use|
|US7647089||26 Jan 2006||12 Jan 2010||Wisconsin Alumni Research Foundation||Surface identification using microwave signals for microwave-based detection of cancer|
|US7656172||18 Jan 2006||2 Feb 2010||Cascade Microtech, Inc.||System for testing semiconductors|
|US7671784||10 Sep 2007||2 Mar 2010||L-3 Communications Cyterra Corporation||Computerized tomography using radar|
|US7681312||31 Jul 2007||23 Mar 2010||Cascade Microtech, Inc.||Membrane probing system|
|US7688062||18 Oct 2007||30 Mar 2010||Cascade Microtech, Inc.||Probe station|
|US7688091||10 Mar 2008||30 Mar 2010||Cascade Microtech, Inc.||Chuck with integrated wafer support|
|US7688097||26 Apr 2007||30 Mar 2010||Cascade Microtech, Inc.||Wafer probe|
|US7720532||23 Mar 2005||18 May 2010||Dune Medical Ltd.||Clean margin assessment tool|
|US7723999||22 Feb 2007||25 May 2010||Cascade Microtech, Inc.||Calibration structures for differential signal probing|
|US7725151 *||27 May 2004||25 May 2010||Van Der Weide Daniel Warren||Apparatus and method for near-field imaging of tissue|
|US7740803||22 Sep 2003||22 Jun 2010||Sonoptix (Uk) Limited||Biochemical sensor with thermoelastic probes|
|US7750652||11 Jun 2008||6 Jul 2010||Cascade Microtech, Inc.||Test structure and probe for differential signals|
|US7759953||14 Aug 2008||20 Jul 2010||Cascade Microtech, Inc.||Active wafer probe|
|US7761983||18 Oct 2007||27 Jul 2010||Cascade Microtech, Inc.||Method of assembling a wafer probe|
|US7761986||10 Nov 2003||27 Jul 2010||Cascade Microtech, Inc.||Membrane probing method using improved contact|
|US7764072||22 Feb 2007||27 Jul 2010||Cascade Microtech, Inc.||Differential signal probing system|
|US7809425||17 Jul 2006||5 Oct 2010||Dune Medical Devices Ltd.||Method and apparatus for examining a substance, particularly tissue, to characterize its type|
|US7809427||11 Feb 2005||5 Oct 2010||Wisconsin Alumni Research Foundation||Time domain inverse scattering techniques for use in microwave imaging|
|US7876114||7 Aug 2008||25 Jan 2011||Cascade Microtech, Inc.||Differential waveguide probe|
|US7876115||17 Feb 2009||25 Jan 2011||Cascade Microtech, Inc.||Chuck for holding a device under test|
|US7888957||6 Oct 2008||15 Feb 2011||Cascade Microtech, Inc.||Probing apparatus with impedance optimized interface|
|US7893704||20 Mar 2009||22 Feb 2011||Cascade Microtech, Inc.||Membrane probing structure with laterally scrubbing contacts|
|US7898273||17 Feb 2009||1 Mar 2011||Cascade Microtech, Inc.||Probe for testing a device under test|
|US7898281||12 Dec 2008||1 Mar 2011||Cascade Mircotech, Inc.||Interface for testing semiconductors|
|US7899515||29 Mar 2006||1 Mar 2011||Dune Medical Devices Ltd.||Electromagnetic sensors for tissue characterization|
|US7904145||7 May 2007||8 Mar 2011||Dune Medical Devices Ltd.||Clean margin assessment tool|
|US7940069||15 Dec 2009||10 May 2011||Cascade Microtech, Inc.||System for testing semiconductors|
|US7956752||13 Jan 2006||7 Jun 2011||Matthew Henderson||Transponder bolt seal and a housing for a transponder|
|US7969173||23 Oct 2007||28 Jun 2011||Cascade Microtech, Inc.||Chuck for holding a device under test|
|US8013623||3 Jul 2008||6 Sep 2011||Cascade Microtech, Inc.||Double sided probing structures|
|US8019411||12 Feb 2007||13 Sep 2011||Dune Medical Devices Ltd.||Probes, systems, and methods for examining tissue according to the dielectric properties thereof|
|US8032211||1 May 2007||4 Oct 2011||Dune Medical Devices Ltd.||Probes, systems, and methods for examining tissue according to the dielectric properties thereof|
|US8050740 *||15 Sep 2004||1 Nov 2011||Wisconsin Alumni Research Foundation||Microwave-based examination using hypothesis testing|
|US8069491||20 Jun 2007||29 Nov 2011||Cascade Microtech, Inc.||Probe testing structure|
|US8089417||30 May 2008||3 Jan 2012||The Royal Institution For The Advancement Of Learning/Mcgill University||Microwave scanning system and miniaturized microwave antenna|
|US8095204||8 Aug 2003||10 Jan 2012||Interstitial, Llc||Apparatus and method for diagnosing breast cancer including examination table|
|US8116845||9 Feb 2006||14 Feb 2012||Dune Medical Devices Ltd.||Tissue-characterization probe with effective sensor-to-tissue contact|
|US8147423||1 Mar 2007||3 Apr 2012||Dune Medical Devices, Ltd.||Tissue-characterization system and method|
|US8195282||1 May 2007||5 Jun 2012||Dune Medical Devices Ltd||Method and apparatus for examining tissue for predefined target cells, particularly cancerous cells, and a probe useful in such method and apparatus|
|US8295920||25 Oct 2004||23 Oct 2012||Medrad, Inc.||System for detecting fluid changes and sensoring devices therefor|
|US8319503||16 Nov 2009||27 Nov 2012||Cascade Microtech, Inc.||Test apparatus for measuring a characteristic of a device under test|
|US8352015||27 May 2008||8 Jan 2013||Kyma Medical Technologies, Ltd.||Location tracking of a metallic object in a living body using a radar detector and guiding an ultrasound probe to direct ultrasound waves at the location|
|US8376947 *||26 Mar 2008||19 Feb 2013||Bioquantetics, Inc.||Application of image-based dynamic ultrasound spectrography (IDUS) in detection and localization of breast microcalcifcation|
|US8410806||20 Nov 2009||2 Apr 2013||Cascade Microtech, Inc.||Replaceable coupon for a probing apparatus|
|US8436766||5 Nov 2010||7 May 2013||Technology Service Corporation||Systems and methods for suppressing radar sidelobes using time and spectral control|
|US8451017||18 Jun 2010||28 May 2013||Cascade Microtech, Inc.||Membrane probing method using improved contact|
|US8568404||19 Feb 2010||29 Oct 2013||Covidien Lp||Bipolar electrode probe for ablation monitoring|
|US8717430||26 Apr 2010||6 May 2014||Medtronic Navigation, Inc.||System and method for radio-frequency imaging, registration, and localization|
|US8721565||4 Aug 2005||13 May 2014||Dune Medical Devices Ltd.||Device for forming an effective sensor-to-tissue contact|
|US8724864||22 May 2007||13 May 2014||Medfield Diagnostics Ab||System and method relating to examination of an object|
|US8736486 *||23 Dec 2008||27 May 2014||Interstitial, Llc||Synthetic aperture radar system|
|US8882759||18 Dec 2009||11 Nov 2014||Covidien Lp||Microwave ablation system with dielectric temperature probe|
|US8907682||30 Jul 2010||9 Dec 2014||Sensible Medical Innovations Ltd.||System and method for calibration of measurements of interacted EM signals in real time|
|US8923951 *||6 Jan 2011||30 Dec 2014||Mmtc, Inc.||Diagnostic devices using thermal noise|
|US8928524||15 Jul 2013||6 Jan 2015||Technology Service Corporation||Method and system for enhancing data rates|
|US8977340||9 Feb 2011||10 Mar 2015||Dartmounth College||System and method for collection and use of magnetic resonance data and microwave data to identify boundaries of interest|
|US8989837||14 Apr 2010||24 Mar 2015||Kyma Medical Technologies Ltd.||Methods and systems for determining fluid content of tissue|
|US9044158||31 Mar 2011||2 Jun 2015||Indiana University Research And Technology Corp.||Wearable microstrip antennas for skin placement for biomedical applications|
|US9182481||18 Mar 2009||10 Nov 2015||Radio Physics Solutions Ltd.||Remote detection and measurement of objects|
|US9226979||31 May 2012||5 Jan 2016||Dune Medical Devices Ltd.||Method and apparatus for examining tissue for predefined target cells, particularly cancerous cells, and a probe useful in such method and apparatus|
|US9265438||1 Dec 2009||23 Feb 2016||Kyma Medical Technologies Ltd.||Locating features in the heart using radio frequency imaging|
|US9289550||20 Mar 2006||22 Mar 2016||Bayer Healthcare Llc||Apparatus and method for detecting fluid extravasation|
|US9316734 *||21 Oct 2013||19 Apr 2016||The Curators Of The University Of Missouri||Free-hand scanning and imaging|
|US9335407||17 Sep 2010||10 May 2016||Radio Physics Solutions Ltd||Detection of objects|
|US9429638||1 Apr 2013||30 Aug 2016||Cascade Microtech, Inc.||Method of replacing an existing contact of a wafer probing assembly|
|US9526460||26 Dec 2011||27 Dec 2016||Dune Medical Devices Ltd.||Tissue-characterization probe with effective sensor-to-tissue contact|
|US9572511 *||4 Mar 2010||21 Feb 2017||Sensible Medical Innovations Ltd.||Methods and systems for monitoring intrabody tissues|
|US9572512||12 Feb 2015||21 Feb 2017||Kyma Medical Technologies Ltd.||Methods and systems for determining fluid content of tissue|
|US9713437 *||27 Jan 2014||25 Jul 2017||Cianna Medical, Inc.||Microwave antenna apparatus, systems, and methods for localizing markers or tissue structures within a body|
|US9746552||2 Oct 2015||29 Aug 2017||Radio Physics Solutions Ltd.||Remote detection and measurement of objects|
|US9750425||12 Jun 2009||5 Sep 2017||Dune Medical Devices Ltd.||Graphical user interfaces (GUI), methods and apparatus for data presentation|
|US9788752||20 Nov 2015||17 Oct 2017||Zoll Medical Israel Ltd.||Implantable dielectrometer|
|US20020120189 *||23 Jan 2002||29 Aug 2002||Clarbruno Vedruccio||Electromagnetic analyzer of anisotropy in chemical organized systems|
|US20030036674 *||26 Jul 2002||20 Feb 2003||Bouton Chad Edward||Electromagnetic sensors for biological tissue applications and methods for their use|
|US20030036713 *||26 Jul 2002||20 Feb 2003||Chad Bouton||Detection of fluids in tissue|
|US20030088180 *||3 Jul 2002||8 May 2003||Van Veen Barry D.||Space-time microwave imaging for cancer detection|
|US20030104783 *||4 Dec 2001||5 Jun 2003||Esion-Tech, Llc||Adaptive electromagnetic interference rejection system and method|
|US20030138378 *||18 Nov 2002||24 Jul 2003||Dune Medical Devices Ltd.|
|US20030178034 *||25 Nov 2002||25 Sep 2003||Spatial Dynamics, Ltd.||Dielectric personnel scanning|
|US20040097811 *||8 Aug 2003||20 May 2004||Smith Randall M.||Apparatus and method for diagnosing breast-cancer including examination table|
|US20040167399 *||14 Feb 2004||26 Aug 2004||University Of Florida||Breast cancer detection system|
|US20040254457 *||27 May 2004||16 Dec 2004||Van Der Weide Daniel Warren||Apparatus and method for near-field imaging of tissue|
|US20050107693 *||17 Sep 2004||19 May 2005||Elise Fear||Tissue sensing adaptive radar imaging for breast tumor detection|
|US20050107718 *||18 Oct 2004||19 May 2005||Dan Hashimshony||Method and system for examining tissue according to the dielectric properties thereof|
|US20060008226 *||7 Sep 2005||12 Jan 2006||Cascade Microtech, Inc.||Fiber optic wafer probe|
|US20060058606 *||15 Sep 2004||16 Mar 2006||Davis Shakti K||Microwave-based examination using hypothesis testing|
|US20060097909 *||7 Nov 2005||11 May 2006||Blunt Shannon D||Multistatic adaptive pulse compression method and system|
|US20060148092 *||22 Sep 2003||6 Jul 2006||Hayat Sindi||Biochemical sensor with thermoelastic probes|
|US20060183995 *||26 Jan 2006||17 Aug 2006||Bond Essex J||Surface identification using microwave signals for microwave-based detection of cancer|
|US20060241409 *||11 Feb 2005||26 Oct 2006||Winters David W||Time domain inverse scattering techniques for use in microwave imaging|
|US20060253107 *||23 Mar 2005||9 Nov 2006||Dune Medical Devices Ltd.||Clean margin assessment tool|
|US20060264738 *||17 Jul 2006||23 Nov 2006||Dune Medical Devices Ltd.||Method and apparatus for examining a substance, particularly tissue, to characterize its type|
|US20070032739 *||4 Aug 2005||8 Feb 2007||Dune Medical Devices Ltd.||Device for forming an effective sensor-to-tissue contact|
|US20070032747 *||9 Feb 2006||8 Feb 2007||Dune Medical Devices Ltd.||Tissue-characterization probe with effective sensor-to-tissue contact|
|US20070060816 *||12 Sep 2005||15 Mar 2007||Simpkin Ray A||Imaging system|
|US20070073144 *||12 Sep 2005||29 Mar 2007||Simpkin Ray A||Imaging system|
|US20070179397 *||12 Feb 2007||2 Aug 2007||Dune Medical Devices Ltd.||Probes, systems, and methods for examining tissue according to the dielectric properties thereof|
|US20070230767 *||11 Jul 2005||4 Oct 2007||Intellectual Property Bank Corp.||Method for Displaying Bioinformation Using Millimeter-Wave Band Electromagnetic Wave, Device for Acquiring and Displaying Bioinformation|
|US20070255169 *||1 May 2007||1 Nov 2007||Dune Medical Devices Ltd.||Clean margin assessment tool|
|US20070293752 *||12 Sep 2005||20 Dec 2007||Industrial Research Limited||Synthetic Focusing Method|
|US20080021343 *||1 May 2007||24 Jan 2008||Dune Medical Devices Ltd.||Probes, systems, and methods for examining tissue according to the dielectric properties thereof|
|US20080143581 *||19 Dec 2006||19 Jun 2008||Bellsouth Intellectual Property Corporation||Scanning and projection systems and methods|
|US20080169961 *||10 Sep 2007||17 Jul 2008||L-3 Communications Cyterra Corporation||Computerized Tomography Using Radar|
|US20080214953 *||1 Mar 2007||4 Sep 2008||Dune Medical Devices Ltd.||Tissue-characterization system and method|
|US20090012391 *||12 Sep 2005||8 Jan 2009||Simpkin Ray A||Imaging system|
|US20090024026 *||12 Sep 2005||22 Jan 2009||Simpkin Ray A||Imaging system|
|US20090060366 *||29 Oct 2007||5 Mar 2009||Riverain Medical Group, Llc||Object segmentation in images|
|US20090062637 *||29 Mar 2006||5 Mar 2009||Dune Medical Devices Ltd.||Electromagnetic Sensors for Tissue Characterization|
|US20090121877 *||13 Jan 2006||14 May 2009||Matthew Henderson||Transponder bolt seal and a housing for a transponder|
|US20090187109 *||15 Dec 2008||23 Jul 2009||Dune Medical Devices Ltd.|
|US20090247869 *||26 Mar 2008||1 Oct 2009||Bioquantetics, Inc.||Application of image-based dynamic ultrasound spectrography (IDUS) in detection and localization of breast micro-calcification|
|US20090253978 *||12 Jun 2009||8 Oct 2009||Dune Medical Devices Ltd.||Graphical User Interfaces (GUI), Methods And Apparatus For Data Presentation|
|US20090299175 *||27 May 2008||3 Dec 2009||Kyma Medical Technologies||Location tracking of a metallic object in a living body|
|US20090309786 *||23 Dec 2008||17 Dec 2009||Stolpman James L||Synthetic aperture radar system|
|US20100010335 *||9 Jul 2009||14 Jan 2010||Electronics And Telecommunications Research Institute||Method and apparatus for diagnosing cancer using electromagnetic wave|
|US20100056907 *||20 Aug 2009||4 Mar 2010||Sensible Medical Innovations Ltd.||Methods and devices of cardaic tissue monitoring and analysis|
|US20100067770 *||22 May 2007||18 Mar 2010||Medfields Diagnostics Ab||System and method relating to examination of an object|
|US20100069744 *||10 Mar 2006||18 Mar 2010||Ray Andrew Simpkin||Imaging System|
|US20100256462 *||4 Sep 2008||7 Oct 2010||Sensible Medical Innovations Ltd.||Method and system for monitoring thoracic tissue fluid|
|US20100283578 *||16 Jun 2008||11 Nov 2010||Matthew Henderson||Transponder Bolt Seal and a Housing for a Transponder|
|US20110025295 *||30 Jul 2010||3 Feb 2011||Sensible Medical Innovations Ltd.||System and method for calibration of measurements of interacted em signals in real time|
|US20110130800 *||14 Apr 2010||2 Jun 2011||Kyma Medical Technologies Ltd||Microwave Monitoring of Heart Function|
|US20110152853 *||18 Dec 2009||23 Jun 2011||Prakash Manley||Microwave Ablation System With Dielectric Temperature Probe|
|US20110160549 *||4 Sep 2008||30 Jun 2011||Saroka Amir||Method, system and apparatus for using electromagnetic radiation for monitoring a tissue of a user|
|US20110319746 *||4 Mar 2010||29 Dec 2011||Sensible Medical Innovations Ltd.||Methods and systems for monitoring intrabody tissues|
|US20120179048 *||6 Jan 2011||12 Jul 2012||Mmtc, Inc.||Diagnostic devices using thermal noise|
|US20140111374 *||21 Oct 2013||24 Apr 2014||The Curators Of The University Of Missouri||Free-hand scanning and imaging|
|US20140187930 *||22 Nov 2013||3 Jul 2014||Sony Corporation||Microwave transmission device and microwave transmission system|
|US20140309522 *||27 Jan 2014||16 Oct 2014||Larry W. Fullerton||Microwave antenna apparatus, systems, and methods for localizing markers or tissue structures within a body|
|EP1647306A2||26 Jul 2002||19 Apr 2006||Medrad, Inc.||Electromagnetic sensors for biological tissue applications|
|EP1788940A1 *||12 Sep 2005||30 May 2007||Industrial Research Limited||Synthetic focusing method|
|EP1788940A4 *||12 Sep 2005||9 Sep 2009||Ind Res Ltd||Synthetic focusing method|
|EP1788946A1 *||12 Sep 2005||30 May 2007||Industrial Research Limited||Imaging system|
|EP1788946A4 *||12 Sep 2005||4 Nov 2009||Ind Res Ltd||Imaging system|
|EP1788947A1 *||12 Sep 2005||30 May 2007||Industrial Research Limited||Imaging system|
|EP1788947A4 *||12 Sep 2005||9 Sep 2009||Ind Res Ltd||Imaging system|
|EP1834667A2||26 Jul 2002||19 Sep 2007||Medrad, Inc.||Electromagnetic sensors for biological tissue applications|
|EP2418599A2||12 Dec 2007||15 Feb 2012||Dune Medical Devices Ltd.||Graphical user interfaces (gui), methods and apparatus for data presentation|
|EP2465428A1 *||12 Jan 2011||20 Jun 2012||Medielma S.r.l.||Electromagnetic detection system and method for locating tumors/calcifications in tissues|
|EP2526857A2||25 Oct 2004||28 Nov 2012||Medrad, Inc.||Systems for detecting fluid changes and sensor devices therefor|
|EP2595694A1 *||21 Jul 2011||29 May 2013||Kyma Medical Technologies Ltd||Implantable radio-frequency sensor|
|EP2595694A4 *||21 Jul 2011||9 Apr 2014||Kyma Medical Technologies Ltd||Implantable radio-frequency sensor|
|EP3181045A1||4 Mar 2010||21 Jun 2017||Sensible Medical Innovations Ltd.||Methods and systems for monitoring intrabody tissues|
|WO2002025773A1 *||17 Sep 2001||28 Mar 2002||Medipeda Ab||Microwave antenna and imaging system for imaging human tissue|
|WO2006091123A1 *||20 Feb 2006||31 Aug 2006||Otkrytoe Akzionernoe Obschestvo 'sistema Venture'||Microwave radiothermotomograph|
|WO2007136334A1 *||22 May 2007||29 Nov 2007||Medfields Diagnostics Ab||System and method relating to examination of an object|
|WO2008082722A2 *||10 Sep 2007||10 Jul 2008||L-3 Communications Cyterra Corporation||Computerized tomography using radar|
|WO2008082722A3 *||10 Sep 2007||2 Oct 2008||L3 Communications Cyterra Corp||Computerized tomography using radar|
|WO2010100649A1||4 Mar 2010||10 Sep 2010||Sensible Medical Innovations Ltd.||Methods and systems for monitoring intrabody tissues|
|WO2011067623A1||1 Dec 2009||9 Jun 2011||Kyma Medical Technologies Ltd||Locating features in the heart using radio frequency imaging|
|WO2012080300A1 *||14 Dec 2011||21 Jun 2012||Medielma Srl||Electromagnetic detection system and method for locating tumors/calcifications in tissues|
|WO2014149183A3 *||27 Jan 2014||30 Jul 2015||Cianna Medical, Inc.||Microwave antenna apparatus, systems and methods for localizing markers or tissue structures within a body|
|WO2015101921A1 *||29 Dec 2014||9 Jul 2015||University Of Manitoba||Imaging using gated elements|
|International Classification||A61B5/05, A61B6/00|
|Cooperative Classification||A61B6/502, A61B5/05, A61B6/0435, A61B5/0507|
|European Classification||A61B5/05M, A61B6/50D, A61B6/04A6, A61B5/05|
|27 Mar 1998||AS||Assignment|
Owner name: INTERSTITIAL, INC., ILLINOIS
Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:BRIDGES, JACK E.;TAFLOVE, ALLEN;HAGNESS, SUSAN C.;AND OTHERS;REEL/FRAME:009069/0649
Effective date: 19980311
|6 Nov 2000||AS||Assignment|
Owner name: INTERSTITIAL, LLC, ILLINOIS
Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:INTERSTITIAL, INC.;REEL/FRAME:011231/0143
Effective date: 20000825
|15 May 2001||CC||Certificate of correction|
|14 Oct 2003||FPAY||Fee payment|
Year of fee payment: 4
|5 Nov 2007||FPAY||Fee payment|
Year of fee payment: 8
|1 Nov 2011||FPAY||Fee payment|
Year of fee payment: 12