Patent Application
20100113921
1. Field of the Invention
This disclosure relates to scan patterns, and more particularly to scan patterns for estimating unknown surfaces with reflection data.
2. Description of Related Art
Imaging, through various means, can generally allow visualization of objects that are obscured from the human eye. In medical applications, for example, imaging can allow a physician to probe internal parts of a body without resorting to invasive procedures such as open surgery. The choice of which imaging modality to use can depend upon the disease or trauma the physician wishes to visualize. For example, a broken arm may, in certain situations, warrant a simple x-ray to determine the extent of bone fracture. In other circumstances, however, magnetic resonance imaging may be justified to determine if a particular disease is present within bone.
In some cases, imaging can include transforming object signal data that corresponds to object structure into a form that is generally recognizable by some attribute, such as object shape. For example, MRI signal data may consist largely of tables of binary information that, when viewed, looks nothing like the object structure being interrogated. Computer algorithms and programs can transform this information into recognizable images that a physician can use to interpret disease.
From the foregoing discussion, it should be apparent that a need exists for a system and method for object surface estimation.
An embodiment of a method for object surface estimation is presented. In one embodiment, the method includes receiving, at a sensor, one or more reflections from a surface of an object while the sensor is scanned in a pattern that localizes the surface of the object within a medium. The method may also include converting the reflections into one or more signals. Additionally, the method may include estimating a position of the surface of the object within the medium in response to the one or more signals.
In a further embodiment, the method may include locating a physiological feature on the object using a grid scan. The grid scan may include receiving signals corresponding to reflections from the object surface while the sensor is moved in a pre-defined path, and determining a distance from the sensor to the physiological feature. Additionally, the method may include performing a scan in a selected area of the surface of the object, wherein an antenna scan pattern is determined in response to information from a coordinate representation of the surface of the object. In another embodiment, the method may include generating a model of the surface of the object in response to a plurality of estimates of the position of a plurality of portions of the surface of the object within the medium.
In a further embodiment, the method may include conducting a first scan of the surface of the object with a first sensor, and conducting a second scan, substantially concurrently with the first scan, of an interior portion of the object with a second sensor. Additionally, the method may include positioning an antenna within a predetermined proximity of the surface of the object in response to the estimate of the position of the surface of the object within the medium.
In a certain embodiment, the method includes emitting a beam of light from a laser directed at the surface of the object, receiving, at a photo-detector, reflections of the light from the surface of the object, and converting the received reflections into one or more electrical signals. In an alternative embodiment, the method may include moving an antenna that receives microwave energy reflections from the surface of the object in a pre-selected pattern, the pre-selected pattern being defined to allow overlap of a first geometrically-estimated area calculated from said energy reflection on the object surface at a first antenna location with a second geometrically-estimated area calculated from the energy reflection on the object surface at a second antenna location, and generating a coordinate representation of the surface of the object. In still another embodiment, the method may include emitting an ultrasonic pulse from a transducer, the ultrasonic pulse directed at the surface of the object, receiving a reflection of the ultrasonic pulse from the surface of the object, and converting the received reflection into one or more electrical signals. Alternatively, the method may include capturing a digital image of the surface of the object with a digital camera, and estimating the location of the surface of the object in response to one or more properties of pixels comprising the digital image.
A system is also presented for object surface estimation. In one embodiment, the system includes a sensor. The sensor may receive one or more reflections from a surface of an object while being scanned in a pattern that localizes the surface of object within a medium, and convert the reflections into one or more signals. Additionally, the system may include a signal processing device coupled to the sensor, the signal processing device configured to estimate a position of the surface of the object within the medium in response to the one or more signals.
In a further embodiment, the signal processing device may locate a physiological feature on the object in response to data received from a grid scan. The grid scan may include receiving signals corresponding to reflections from the object surface while the sensor is moved in a pre-defined path, and determining a distance from the sensor to the physiological feature. In another embodiment, the sensor is further configured to scan in a selected area of the object, wherein a sensor scan pattern is determined in response to information from the coordinate representation or any other means of determining an outline of the surface of the object.
In another embodiment, the system includes a first sensor configured to conduct a first scan of the surface of the object, and a second sensor configured to conduct a second scan, substantially concurrently with the first scan, of an interior portion of the object. In a further embodiment, the signal processing device is further configured to generate a model of the surface of the object in response to a plurality of estimates of the position of a plurality of portions of the surface of the object within the medium. In a further embodiment, the system includes a positioning arm coupled to an antenna, the positioning arm configured to position the antenna within a predetermined proximity of the surface of the object in response to the estimate of the position of the surface of the object within the medium.
In one embodiment, the system includes a laser configured to emit a beam of light directed at the object, and a photo-detector configured to receive reflections of the light from the surface of the object, and convert the received reflections into one or more electrical signals. Alternatively, the system may include an antenna configured to receive microwave energy reflections from an object surface, and a positioning arm coupled to the antenna, the positioning arm configured to move the antenna in a pre-selected pattern, the pre-selected pattern being defined to allow overlap of a first geometrically-estimated area calculated from the energy reflection on the object surface at a first antenna location with a second geometrically-estimated area calculated from the energy reflection on the object surface at a second antenna location. In an alternative embodiment, the system may include an ultrasound transducer configured to emit an ultrasonic pulse directed at the surface of the object, receive a reflection of the ultrasonic pulse from the surface of the object, and convert the received reflection into one or more electrical signals. In still another embodiment, the system may include a digital camera configured to capture a digital image of the surface of the object, and the signal processing device coupled to the digital camera, the signal processing device configured to estimate the location of the surface of the object in response to one or more properties of pixels comprising the digital image.
The term “coupled” is defined as connected, although not necessarily directly, and not necessarily mechanically.
The terms “a” and “an” are defined as one or more unless this disclosure explicitly requires otherwise.
The term “substantially” and its variations are defined as being largely but not necessarily wholly what is specified as understood by one of ordinary skill in the art, and in one non-limiting embodiment “substantially” refers to ranges within 10%, preferably within 5%, more preferably within 1%, and most preferably within 0.5% of what is specified.
The terms “comprise” (and any form of comprise, such as “comprises” and “comprising”), “have” (and any form of have, such as “has” and “having”), “include” (and any form of include, such as “includes” and “including”) and “contain” (and any form of contain, such as “contains” and “containing”) are open-ended linking verbs. As a result, a method or device that “comprises,” “has,” “includes” or “contains” one or more steps or elements possesses those one or more steps or elements, but is not limited to possessing only those one or more elements. Likewise, a step of a method or an element of a device that “comprises,” “has,” “includes” or “contains” one or more features possesses those one or more features, but is not limited to possessing only those one or more features. Furthermore, a device or structure that is configured in a certain way is configured in at least that way, but may also be configured in ways that are not listed.
Other features and associated advantages will become apparent with reference to the following detailed description of specific embodiments in connection with the accompanying drawings.
Like reference symbols in the various drawings indicate like elements.
In general, analysis of biological structures using medical imaging methods can be enhanced using an accurate estimation of one or more of the surface topography, outline, and volume of the structure. In some cancer screening imaging modalities, for example, creating an estimate of the surface outline of the target organ or tissue can be an important factor in developing an overall imaging approach. Some radar-based microwave imaging techniques require an antenna to be placed a known distance from, or on the surface of, the object to be imaged. In some cases this may require a priori knowledge of the object location, size, and shape, which may be difficult to determine precisely in some circumstances. An imaging method providing the ability to scan an object using, e.g., lasers, mobile antennas, ultrasound transducers, or photographic equipment, without any prior knowledge of the object properties would be advantageous in a number of ways.
In general, methods for scanning a surface of a biological structure and estimating the location of the surface relative to the antenna locations are disclosed. One embodiment of the method is exemplified in the Tissue Sensing Adaptive Radar (TSAR) method (described below). Another embodiment of the method is described in the laser imaging method (described below). It will be understood that the methods are equally practicable in other scanning methods and analytical approaches whose objective it is to determine, for example, one or more of structure, volume, and surface features (including surface area) of an object. Additional examples may include ultrasound systems, photographic imaging systems, and the like.
Some imaging techniques use an antenna that directs energy toward a region of interest (such as a breast) and also receives scattered or reflected signals from the region, or from objects within the region. An example of such a technique is TSAR, described in U.S. patent application Ser. No. 10/942,945, which is incorporated herein by reference in its entirety.
Generally, TSAR is a low-power, ultra-wideband microwave imaging method that has been proposed for early stage breast cancer detection. Radar-based breast imaging can include illuminating the breast with short pulses of microwave energy and recording reflections. An image of the breast interior can be constructed in one embodiment by creating a synthetic array of antennas. In one embodiment, a synthetic array can be a single antenna placed at a number of locations, where a signal is collected at each point.
One operational theory of TSAR is that differences in the electromagnetic properties of healthy and diseased tissues can be expected to result in differences in the reflected signals. The reflected signals can be typically composed of a skin reflection and reflections from the breast interior, e.g., interfaces between fatty tissue, glandular tissue and tumors. The reflected signals can be analyzed and focused to create images that indicate the location of strongly scattering objects.
In some implementations of TSAR, an active ultra-wideband (UWB) monostatic radar approach is used to create an electromagnetic scattering map of the interior of a target imaging region, such as a breast. In some implementations of TSAR, a person may lie prone on a table, where a single breast can be suspended in a scan chamber through an opening in the table. The scan chamber can contain one or more antennas and an immersion medium into which the breast is immersed, either fully or partially. In another embodiment, the scan chamber may include a laser emitting device, a camera, an ultrasound transducer, or the like. Data relating to the breast and its surface, including topology, area, and other characteristics can be determined using antenna scan patterns and surface estimation algorithms, such as those disclosed below.
Preliminary scans, i.e., those that are performed to provide an imaging outline of the entire breast, may yield real coordinates of the skin surface with respect to the antenna; these data can provide information for subsequent scan patterns required for tumor detection. In one embodiment, a preliminary scan can be performed with a microwave antenna. In other embodiments, any modality that provides information for subsequent scans may be used, for example, photographic, optical (both laser or stereoscopic), acoustic, and mechanical methods of determining a surface outline may be used. The preliminary scan pattern can also aid in skin subtraction and focusing algorithms which can be used for tumor detection.
In general, a focusing algorithm can be enhanced using an estimated outline of the breast. For example, voxels included in the focusing algorithm may be limited to the voxels determined to exist inside the breast. This may lead to an increase in computational efficiency. In another example, signal propagation through different regions can be utilized in the focusing procedure. In one embodiment, time delays corresponding to the distances traveled in immersion medium and breast tissue can depend on both the properties of the materials and the distance traveled in each material. The surface estimate can be used to improve the estimates of distances traveled in each material relative to estimating the skin location at one point. Examples of this are discussed in U.S. patent application Ser. No. 12/407,671 which is incorporated herein by reference in its entirety. In some cases, for example using the skin estimation algorithm combined with a surface estimation algorithm, the outline of the inside of the skin can also be estimated. Signal propagation estimates may therefore be expanded to include propagation through the skin as well as the immersion medium and breast interior.
In some cases, a surface estimation algorithm can be used in an iterative manner. For example, a first estimated outline of a breast may indicate areas of low and high curvature where different concentrations of antenna or sensor locations may be appropriate. After additional data acquisition, a second and potentially more accurate outline may be created.
In some cases, parallel computing software and specialized hardware can be used to aid in signal processing speed.
In one embodiment, a priori knowledge of the breast location, shape, and size may be required for optimal antenna placement relative to the breast surface, and also to provide an accurate reconstruction of the processed image. In some cases this information may be necessary in order to prevent the antenna from contacting, or being placed at too great a distance from the breast during the imaging (i.e., scanning) process.
In general, generating an acceptable surface estimate of an object (e.g., a breast), can assist in the determination of certain variables for estimating the volume and other characteristics of the object. In some cases, a carefully chosen antenna scan pattern can have an effect on the generation of the surface outline.
The scan patterns described herein to irradiate the breast can be an important factor that may affect the outcome of the tumor detection efficacy, both in living samples and increasingly complex breast models used for research and development purposes. In some cases, the surface estimation process in TSAR uses no a priori information as to the shape or placement (within the scan region) of the breast, yet still operates successfully.
A second contribution generally includes accurately estimating the path lengths in the immersion medium, skin and breast interior. This estimation can benefit the skin subtraction, clutter reduction, focusing, and imaging steps in TSAR image reconstruction algorithms such as those disclosed in U.S. Patent Ser. No. 61/038,022 the contents of which are fully incorporated by reference herein.
In one embodiment of a scanning method, a first scan pattern may include a scan in a grid formation along the bottom of a scanning chamber, such as the scanning chamber described with respect to the TSAR method. The data in this grid scan may be used to approximate, among other things, the height of the breast nipple. A generic, multi-level cylindrical or coronal scan around the perimeter of the chamber can also be performed. The combined grid and perimeter scans can be used to estimate features relating to the surface of the breast, as described below. In one embodiment, the first scan may be carried out by a microwave antenna. Alternatively, the first scan may be carried out by a laser, an ultrasound transducer, a digital camera, or the like.
In one embodiment, the estimate of the breast surface obtained with the first scan may be used to identify appropriate antenna locations for the second scan. In one embodiment, the data from the second scan can be processed through a TSAR algorithm for tumor detection. Exemplary algorithms for processing TSAR data are disclosed in U.S. Patent Ser. No. 61/038,022, previously incorporated by reference.
In some cases, a method for breast surface estimation involves a detection scan in which an antenna or a laser may be scanned around the base and/or perimeter of the imaging region, and reflections of a short pulse can be recorded. The skin location can be estimated and this information can be used to create a simple signal that describes the location of the breast. The collection of simple signals obtained at the antenna and/or laser can be focused, and the resulting image can be processed to create a surface estimate.
In general, the methods provided herein for determining the surface of the breast may not require an a priori assumption regarding the breast position or shape in the array. In many respects, localization of the nipple is straightforward.
In general, a breast surface estimate may be used to define the region of interest for tumor detection, as well as in the design of imaging algorithms. In some implementations, the antennas, lasers or other sensors are scanned around the perimeter of the imaging region, resulting in separations between the antenna and breast of up to 11 cm. In this case, scanning the antenna closer to the breast may assist in increasing the tumor-to-clutter ratio in signals, as the illumination of the breast by the antenna can become more selective and the signal may experience less attenuation due to shorter propagation distances through the immersion medium. Therefore, in general, data from a surface estimate, conducted by various embodiments of a surface estimating technique or method, may be used to position the antennas at a desired distance from the breast for a second tumor detection scan. A system capable of antenna placement can be used to collect data for both the preliminary and subsequent scans. An antenna positioning system is generally described below that can be used for this purpose.
As described above, a breast surface estimation may include a first or preliminary scan of the breast. The first scan may collect information about the breast surface which may be used by a breast surface estimation algorithm to determine a model of the breast surface. The first scan may be conducted by a sensor. Various embodiments of a sensor are described below with reference to
In one embodiment, the a microwave antenna sensor may include a Balanced Antipodal Vivaldi Antenna (“BAVA”) antenna. Alternatively, a BAVA-D configuration antenna may be used. In another embodiment, a Transverse ElectroMagnetic (“TEM”) horn antenna may be used. One embodiment of an antenna sensor is described below with reference to
In an alternative embodiment, the sensor may be a laser. One embodiment of a laser is described below with reference to
One benefit of using a laser sensor is the small beam-width or footprint of the emitted light beam. This characteristic of a laser sensor may allow more precise tracking of the curvature of a breast. Another benefit may include in increased rate of sensing as compared with a microwave antenna. For example, the response time of some embodiments of a laser sensor may be between 50 μs and 50 ms, depending on the system configuration and the characteristics of the laser. This response time may greatly reduce the time required to perform the first scan when a laser sensor is implemented.
In a particular embodiment, the laser may be calibrated for operation in an oil medium. The calibration may be accomplished by scanning an a custom calibration piece, such as a cone or semi-sphere of known dimensions, from a known distance. In one embodiment the calibration may further including calculating medium dependent light propagation parameters based on a distance-voltage relationship. In one embodiment, the distance-voltage relationship may be expressed as d=aν+b, where d is the distance between the laser and the calibration piece, v is the voltage generated by the sensor, and a and b correspond to propagation parameters of the material through which the light propagates.
In an alternative embodiment, the first scan may be accomplished with an ultrasound transducer. The ultrasound transducer may include a piezoelectric transducer configured to generate a sonic pulse of a predetermined frequency in response to an electrical impulse. The sonic pulse may propagate through, e.g., the oil medium, reflect off of the breast surface, and receive a response pulse from the reflection of the breast surface. The response pulse may be converted into an electrical potential.
In another alternative embodiment, the first scan may be accomplished with a digital camera. One embodiment of a digital camera is described below with reference to
In one embodiment, the breast surface estimate can use a preliminary scan of the breast. This preliminary scan can be accomplished by moving the antenna (or multiple antennas) around the perimeter of the imaging region. A preliminary scan can include scanning the antenna(s) around the base and perimeter of the container, as illustrated in
A scan at the base of the container can include positioning the antenna parallel to the chest wall (the x direction in
In general, the data collected during this portion of the scan can be used to determine the location of the nipple and create the surface outline of the lower portion of the breast, as described below. To collect further data with the perimeter scan, the antenna can be oriented in the z direction (
The diameter of the circular scan paths may take into account certain variables such as the physical size of the tank, the antenna dimension. With an estimate of the location of the nipple and knowledge of the chest wall position (i.e. top of the tank), the extent of the perimeter scan in the z direction (
In some cases, a breast surface estimate can be created by focusing a modified version of the recorded reflections, incorporating information and methods adapted from the TSAR algorithm. The skin location at each antenna can be estimated from the reflected signal using the impulse response technique. This method can include deconvolving the reflection from a reference object placed in the immersion liquid from the breast reflection. The skin location information can be used to create a simple signal from each reflection:
where Si is the signal corresponding to antenna i and nl is the time sample at which the skin location is estimated. The signals can be passed to the TSAR focusing algorithm, which essentially can time-shift and sum the data. At focal point location r, the pixel intensity can be calculated as:
where ri is the location of antenna i, νext is the estimate of the velocity in the immersion medium, dt is the time step and M is the number of antennas being used in the current surface estimation scan. When estimating the nipple location, M corresponds to the number of antennas in the base scan. When estimating the entire surface or volume of the breast, M can correspond to all antenna locations in the base and perimeter scans. To identify a surface from the resulting focused data, thresholding can be applied to E(r).
The threshold can be set such that pixels corresponding to the breast must contain contributions from all antennas; therefore, at each location r in the focused image, E(r) should equal one (1) for the point to be considered part of the breast. A physical interpretation of this criterion in two dimensions is shown in
In some cases, two surface estimates can be created. Base scan data can be processed first, and the nipple location can be estimated from the result. This information can be used to design the perimeter scan, as described above. The algorithm can be applied to data recorded at all base and perimeter antennas, and the surface of the breast can be estimated from the nipple to the coronal plane of antennas closest to the chest wall. In most embodiments, extrapolation is not necessary to extend the surface artificially past the scan region. With a skin thickness estimate, the method described above may also be applied when determining the surface of the breast inside the skin. A second signal similar to equation (1) can be formed, where nl is simply modified to the time sample corresponding to the location and thickness of the skin layer. By taking the difference between the surface estimates created with the two signals, a representation of the skin layer can be obtained.
In some cases, multiple reflections may exist from the breast when the antenna is located such that, e.g., the chest wall and the breast are observed, i.e., detected. In some cases, determining the impulse response of the breast, i.e., using a deconvolution procedure to estimate skin location and thickness can assist in distinguishing the first and subsequent reflections for optimal use in TSAR imaging algorithms.
In one implementation, a deconvolution procedure can generally include the reflection from the object of interest, as well as a reflection from a known, calibration object, such as a metal plate. The impulse response, which describes the reflections from all material interfaces in the object of interest, can be extracted by deconvolving the reference signal from the reflection from the object of interest.
The algorithm described above is for estimation of the breast surface and calculation of a surface model in response to data received from a microwave antenna. In an alternative embodiment, a laser may be used to scan the breast surface. In such an embodiment, the algorithm described above may be adapted for use with laser data. Alternatively, the laser may be used in conjunction with commercially available surface estimation software to determine a breast model. For example, Powercrust™ software available from the University of Texas Department of Computer Science may be incorporated into, e.g., a dedicated signal processing device, for estimation of the breast model in response to data from the laser sensor.
In this example, data were generated with the Finite Difference Time Domain (FDTD) software similar to previous work with TSAR. The hemispherical model is shown in
Three simulation scenarios were performed. The first consisted of scanning the dipole in a circular path around the breast at z=5.5 cm (for orientation, refer to
The second simulation scenario consisted of scanning the antenna around the hemisphere in 9 circular paths with 10 antenna locations per row. A “path,” as used herein, implies a trajectory of an antenna on a two-dimensional plane of a coordinate system and a “row,” as used herein, implies the set of locations along a given path where measurements are performed. Alternate rows began at references corresponding to 0 and 18 degrees, and the antenna was moved in 36-degree increments until 10 antenna positions per row were acquired. The scans were performed from z=2.5 to z=6.5 cm in increments of 0.5 cm with the antennas aligned with the z direction.
The third scenario included creating an outline in the sagital plane. The antenna was positioned at x=10.0 cm (refer to
Information from every antenna location was used to create the outline of the hemisphere using the TSAR algorithm. First, skin location and thickness estimates are found using the impulse response method. The skin location and thickness is estimated for each antenna in the 9×10 array. The average antenna skin location and thickness estimates for scenario one are 1.16 cm and 1.75 mm respectively. This location estimate corresponds to a 2.57% error when compared to the center of the dipole. Next, a focusing program examined every voxel inside the imaging region to determine whether it was at a sufficient distance to be classified as outside or inside the model. Antennas from multiple rows and locations can have an influence in determining whether a specific voxel is classified as immersion medium, skin, or model interior in this method, and does not require any a priori knowledge about the shape of the model.
The sagital outline was generated using the antennas described as in scenario three, and is shown in
More realistic models may be used to test the breast surface detection and antenna placement algorithms. These models may be derived from Magnetic Resonance (MR) images. The MR images may be converted into models suitable for use with an electromagnetics simulator. In this example, several complex models are created, however similar procedures are employed. Images are segmented by applying a threshold in order to identify the interior and exterior regions of the breast. The threshold is adjusted manually for each breast model to identify the interior of the breast, including the skin.
Breast model 1 is a simplified model consisting of a 2-mm thick layer of skin bounding a non-dispersive, homogeneous fatty tissue with εr=9 and σ=0.4 S/m. The layer of skin with εr=36 and σ=4 S/m is added at the boundary between the interior and exterior of the model that is identified in the thresholding step.
In breast model 2, variations in the MR pixel intensity are used to create variations in the electrical property distribution in the breast. The pixels of the breast are linearly mapped to the electrical properties by selecting the maximum required permittivity and conductivity values (εr=36 and σ=4 S/m). A modified version of breast model 2 (2B) has constant properties assigned to the skin (εr=36 and σ=4 S/m), which are identified in the MR image using pixel intensities. Breast model 2, therefore, has thickness and property variations in the skin, while the skin on breast model 2B exhibits thickness changes but maintains a constant electrical property of εr=36 and σ=4 S/m. To increase the realistic nature of the model, a layer of muscle (εr=50 and σ=4 S/m), representing the chest wall, is attached to the breast models. The breast models are placed in a lossless immersion liquid with similar electrical properties to canola oil (εr=3.0). To illustrate the three models,
Reflections from the breast model are obtained with simulations performed using the finite difference time domain (FDTD) method. A single, resistively loaded dipole antenna of length 2 cm transmits an ultra-wideband pulse and reflections are recorded at the same antenna. The transmitted pulse is a differentiated Gaussian signal with frequency content from 1-10 GHz. The excitation of the antenna with a pulse is simulated as the antenna is scanned through locations forming grid and perimeter scans, as described below.
For each antenna location, skin location is estimated first. For breast model 1, the actual distance between the antenna feed point and the closest point on the realistic model is calculated and found to track very closely with the estimated antenna/breast separations using the IR method, with a mean error of 1.50 mm over 224 (lower 5 of 7 scanned rows×32 positions) cylindrically scanned antennas from z=30 to z=70 mm.
Next, surface estimates are constructed for the breast models using the algorithm described above. The performance of the surface estimation algorithm can be assessed by first considering a two-dimensional result. To understand the influence of the number of antennas per row, a 2D surface estimate can be created using breast model 1 and one row of the perimeter scan collected at z=103 mm and with path diameter of 18.0 cm. Results obtained using 8, 16, and 32 antennas are compared in
Considering the scan diameter of 18 cm and 8 antennas, the beamwidths of neighboring antennas do not overlap significantly, giving rise to the errors observed in
Next, a full scan of breast model 1 was performed. The base scan was performed, and the resulting location of the nipple was estimated to be z=3.4 cm, while the actual location was 3.8 cm. Although there was a 4 mm error, the location was an estimate used to determine the number of rows for the perimeter scan.
The distance between the nipple estimate and chest wall was 9.7 cm. A row separation of 2 cm was selected, which corresponded to the length of the resistively loaded dipole. This resulted in 4 rows of antennas in the perimeter scan with antenna feed points spanning from z=40 mm to z=100 mm. Based on the results of the 2D estimate, 16 antennas per row were included. With multiple rows, increased information was provided, especially with judicious selection of scan pattern. Each row was rotationally offset by 7.5°, so every third row had the same antenna positioning. This provided offset antenna locations in rows above and below a selected row, while maintaining overlap between the antenna patterns.
Another approach to assessing error is calculating the quantity of voxels in error. Taking into account the inherent slight overestimation of the skin location (relative to the antenna location) associated with the impulse response method, the larger number of underestimation error sites compared to overestimation was predictable. Although the number of error sites associated with underestimation was large, the average error value was low. Therefore, the error associated with most sites was expected to be 1 mm.
With breast model 1, the effect of an increased number of antennas was investigated. Breast surface estimates were created with varying numbers of rows and antennas per row, and errors are summarized in Table 1. First, results obtained with 4 and 9 rows are compared. The additional 5 rows were included in such a way that the array spanned from 30 mm to 110 mm in the z-direction with 1 cm separation between rows. For comparison purposes, the surface was first estimated to z=100 mm. The breast surface estimate was improved with the additional rows, especially near the chest wall. This was demonstrated by the decrease in maximum overestimation error, average overestimation error, and error sites. Next, an additional 16 antennas were added to each row of the 9×16 array. In this 9×32 array, scan locations were simply placed at the midpoint between the existing locations. This maintains a rotational offset that repeats every third row, however now the offset is 3.75°. The overestimation error did not change significantly with the added antennas, consistent with the observations in the 2D case. A slight improvement in underestimation error was noted. This again underscores the importance of a sufficient antenna scan pattern that minimizes data collection time while maximizing breast coverage by the antennas.
In all cases reported in Table 1 for breast model 1, the constant maximum underestimation error was attributed to the same base scan antenna. Enabling one antenna to have such a dominant affect on the error is undesirable. In some embodiments, understanding the effect of the threshold applied to E(r) on overall error may be an important consideration. Referring to
Surface estimates were created for breast models 2 and 2B. The base scan gave an estimate of the nipple location at z=3.9 cm, compared with the actual location of z=3.4 cm. The rows of antennas, therefore, started at z=4.0 cm and continued to 2 cm below the chest wall. With 1 cm separation between rows, this corresponded to a perimeter scan consisting of 7 rows with 16 antennas per row. Each path has diameter of 19.0 cm. The surface estimate for model 2 is shown in
The small average errors obtained with models 2 and 2B demonstrate the success of the method and show that the algorithm is robust to variations in skin thickness and properties. Further, the size of this model and its placement inside the perimeter scan can ensure a more difficult surface to estimate. The model was slightly smaller than breast model 1, the perimeter scan was 1 cm greater in diameter, and the model was located in an offset position relative to the center of the scan pattern. Several of the perimeter antennas furthest away from the model are located significantly closer to the chest wall, so the skin location estimates at these antennas correspond to the chest wall interface. However, several antennas in the outlying regions of the base scan overestimated the skin location relative to the antenna. In combination, these errors in skin location resulted in surface estimates with no overestimation errors and increased underestimation errors when compared to breast model 1. The maximum errors obtained with both models corresponded to the region in which the perimeter scan antennas are located furthest from the model.
The results so far have demonstrated the effectiveness of the surface estimation algorithm. Next, the impact of antenna scan pattern, specifically the grid scan, on results is further explored. With breast model 1, two rows in the perimeter scan are selected. These rows contain 32 equally spaced positions on the circumference of the circle and are located at z=30 mm and z=40 mm. The grid antennas, placed perpendicular to the coronally scanned antennas, were configured in 9×9 and 5×7 arrays as shown in
The above observations suggest that the grid scans can be an important factor in estimating surface outline regions not illuminated by the perimeter scan (e.g., the nipple region). In one implementation, grid layer scans can be used when the lower perimeter rows do not illuminate the nipple sufficiently (e.g., as illustrated in
To test the differences in results obtained with various grid scans, the grid size can be systematically reduced, beginning with a 9×9 array scan, which encompasses slightly more than the area of the cage scan. The grid scan patterns can include 7×7, 5×7, 5×5, and 3×3 configurations, for example. One example grid scan pattern has been included for surface outline creation comparison and is a 5×7 grid configuration that is roughly centered in the cage scan (
Surface outlines for 5×7 and 3×3 grid scans are shown in
For the example shown, the 3×3 grid scan can be adequate for localizing the z-location of the nipple if the pattern is underneath the nipple; the lowest z-location of the overall breast volume, however, can be underestimated due to large areas below z=60 mm not being illuminated. The differences in surface estimation error in this model may not substantially differ between the 5×5 and the 5×7 scan patterns. This may, however, not be the case for another realistic model. In one implementation, the 7×7 grid scan, with very few error sites and 1.4 mm maximum error, may be a superior scan pattern when considering a practical application. The difference of 32 scan positions between a 9×9 grid and a 7×7 grid may decrease the overall scan time and lead to reduced error associated with movement artifacts.
In some cases, the orientation of an antenna can be defined according to certain geometric considerations. The preliminary scan is not limited to a cylindrical scan and also a grid scan in two dimensions. In some cases, the grid scan can move an antenna in three dimensions, for example, an antenna can be scanned in a spherical or hemispherical pattern around a breast. Because the shape and size of human breasts differ among women, a grid scan pattern can be chosen that best accommodates these and other breast characteristics, allowing an optimal initial breast scan to be performed.
The error associated with the surface estimate can be analyzed and an attempt can be made to quantify the discrepancies between real and estimated breast surfaces placed in a relevant context. In some implementations, the error associated with the surface estimate may lead to errors in antenna placement during the tumor sensing scan, as well as in erroneous wave velocity estimates when focusing to obtain an image. When considering a subsequent antenna scan, an underestimation can be worse than an overestimation, as the antenna may be placed too close to the breast. When considering image focusing, a single erroneous region may mean very little when focusing at one location using a particular antenna, but may have greater impact when focusing at another location using a different antenna, as illustrated in
In some cases, the following method can be used to describe the error. Using OE1 in
Error figures can be directly influenced by the resolution of the imported model and the resolution used in focusing. In one implementation, both resolutions are set to 1 mm to enable direct comparison between simulated and estimated surfaces. Therefore, the minimum error in a region under consideration can be 1 mm. In some implementations, when calculating the average error associated with a specific surface estimate, only erroneous pixels are included. This can create a much larger average error than would be calculated using all points on the estimated breast surface, however, a more relevant average error can be computed. With the current resolution, this can imply that the minimum average error attainable is 1 mm. The number of pixels on the surface can be included and considered erroneous, or Error Sites, used to calculate the average error. This value can be an important gauge as a maximum error value may be identical for different antenna configurations however the quantity of Error Sites may vary. This may offer insight into the relative success of the surface detection algorithm.
During data collection of a surface estimate, the antenna may be placed at a selected distance from the object being imaged. In some embodiments, the distance may be significant, such as a distance of up to 14 cm. In some embodiments, the surface estimate may be acquired from an alternate imaging modality, like those mentioned above (e.g. audio, RF, optical, infrared). For imaging the interior of the breast, for example, it may be beneficial to locate the antenna close to the breast surface. The results of the surface estimate may be used to place the antennas at a desired distance from the object for a second scan, for example, a high-resolution scan, or a tumor-detection scan if the object is within tissue. An automated algorithm has been created to place the antennas on a row-by-row basis at a desired distance from the surface.
The algorithm is illustrated as a flow diagram in
The location on the breast surface closest to the antenna can be identified, in this case the location closest to the end of the antenna (z=z1+Aheight/2). Placing the end of the antenna at the desired distance can ensure that the antenna does not touch or couple with the breast; if the feed position is used to determine the separation, the antenna may be placed too close to the breast in areas of curvature such as the chest wall. A coronal slice through the surface estimate can be obtained at the selected z location. The approximate center of the breast surface can be calculated by averaging the maximum and minimum extent of the breast in the x and y directions. A line that radiates out from the estimated center of the breast can be drawn for each desired antenna location.
In one embodiment, the point where the line and edge of the breast surface intersect can be determined using pixel intensities in the surface estimate images. The point can then be translated the desired distance, dint, from the breast surface along the line. The x-y coordinates can be calculated for this point and translated down Aheight/2 to the row selected for antenna placement. The process can be repeated for each row required to span the distance between the nipple and chest. The breast surface estimate may also be used to determine the number of antennas per row, as fewer antennas may be required near the nipple. In this case, an estimate of the breast circumference may be obtained from the 2D section through the surface estimate, and divided by the desired separation between antennas in order to determine the number of antennas required in a specific row.
The antenna placement algorithm is applied to the surface estimates created with 9 rows of 32 antennas and breast model 1. For the tumor detection scan, the resistively loaded dipole described above was used and dint=2 cm, dheight=1 cm and Aheight=2 cm based on antenna dimensions. The resulting antenna positions on breast model 1 are shown in
In one general aspect, an antenna positioning system (APS) is provided. In some cases, the positioning system can be used to carry out the scan patterns described herein. In some cases, the APS can be used to scan a breast with an antenna or with multiple antennas. The system can positioned, for example, beneath a patient bed while the patient lies with her breast extending through a hole in the table and extended into a tank that is part of the APS system. The system can scan the antenna around the breast, including from the chest to the nipple and encircling the breast. The pitch of the antenna may change to provide optimal scanning of the breast.
One embodiment of an antenna positioning system 1900 is shown in
In general, TSAR reflections may be optimal for imaging purposes when the antenna is located at an optimal distance from the breast in order to record reflections. The selection of the optimal position may be accomplished by translation, rotation, vertical positioning and pitch adjustment of the antenna.
Translation may be used to position the antenna at the correct horizontal distance from the breast. The entire tank 1910 and antenna 1930 can be moved into the desired position using the translation plate 1955. A motor 1968 can be connected to a horizontal threaded rod 1970. When threaded rod 1970 turns, the tank 1910 and stand 1945 can slide along the groove in the translation plate 1955. Vertical positioning of the antenna 1930 can be provided by a motor 1975 connected to a second threaded rod 1976 that is rotatably integrated into the positioning arm 1965 that passes through the stand 1945 (
Referring now to
In some cases, to perform a scan of the breast, the antenna can be moved to multiple positions encircling the breast. In one embodiment, the entire tank can rotate on the rotation plate in order to provide this function. Before or after rotation, a translation may be required to appropriately position the antenna.
In general, the antenna positioning system can provide flexible scan sequences. In some cases a “row” of data from antenna positions can be collected that encircle the breast. In these cases, the antenna can be positioned near the breast and the system can be rotated to capture data points around the breast circumference. In some cases minor adjustments to the position of the antenna may be warranted.
In some cases the system can provide a scan sequence to collect data in “column” format. In this case, the antenna can scan over a range of vertical positions. Adjustments to positioning (e.g. translation, pitch) may be required during these types of scan.
In some cases, the system can provide scanning of specific points on the surface of the breast or within the volume of the breast. In these cases, the antenna can be moved to specific locations where data are then collected. This type of scanning may be useful, for example, when inspecting a number of suspicious areas rather than scanning the entire breast. Generally, the system can be designed in order to provide flexible scanning capabilities and is extensible to multiple antennas. In some cases, one rotation plate and one translation plate can be incorporated, as in the design shown in
In general, an antenna can be attached to the positioning system by any means. In one embodiment, a cable can connect the antenna to the measurement system by attaching to the antenna and passing through a hole on the antenna positioning arm. This arrangement may minimize cable flex during motion.
In general, the antenna positioning system can use a cylindrical coordinate system such as that shown in
In general, an antenna fixture system 1980 is provided. In some cases, the fixture system can consist of a plastic plate with a hole to permit access for the antenna cable, two pins to ensure alignment, and two threads to attach the antenna. In some cases, a minimum area of the plate can be dictated by the antenna dimension; its thickness can be variable and chosen to provide optimal functionality for a given arrangement or antenna size. In some cases the connector hole can be located at the center of the plate, and its size can be commensurate with the cable connector diameter.
In general, the antenna positioning system can be used to accurately move and place an object (e.g., an antenna) around an object, e.g., a human breast. In general, the antenna positioning system can move in up to as many degrees of freedom as is necessary to suit the particular purpose at hand. In some cases, four degrees of freedom are used: circular movement around the object (Phi), vertical movement, radial movement (Rho), and pitch angle movement (Theta). Exemplary degrees of freedom are illustrated in
It is to be understood that while the disclosure has been described in conjunction with the detailed description thereof, the foregoing description is intended to illustrate and not limit the scope of the disclosure, which is defined by the scope of the appended claims.
For example, several alternatives to the error analysis described are possible. One approach to error calculation involves determining the error per antenna scan position relative to all focus locations. Calculating the maximum distances that are influenced by over- and underestimation areas gives a number relevant to each antenna and whether its signals are hampered by unacceptably large error. A bar graph may then be generated depicting each antenna and trends may be determined to give insight into the locations of unacceptable over- and underestimation areas. A further improvement would be to determine the error relative to each focus location. A total over- and underestimation distance may be calculated using all included antennas and the distances each signal traveled through erroneous estimations. A surface plot may then be generated of each focused voxel containing this total distance.
Although the methods and variants thereof have been discussed in connection with the Tissue Sensing Adaptive Radar (TSAR) system, they are applicable to other radar-based breast cancer detection systems. The methods and variants thereof are also easily adapted to different antennas.
As described above, alternative embodiments of the first or preliminary scan system 2200 may include a laser 2202, an ultrasound transducer (not shown), and/or a photographic device such as a digital camera 2204.
Additionally, as described in
In a further embodiment, the at least one of the antenna 1930, the laser 2202, and/or the digital cameral 2204 may be coupled to a signal processing device 2206. The signal processing device 2206 may comprise a special purpose signal processing device configured with computer readable code. When executed by the signal processing device 2206, the computer readable code may cause the signal processing device to perform operations associated with the various steps and functions of the algorithms described above.
In a further embodiment, the computer readable code may be stored on a tangible computer readable medium 2208. Examples of a tangible computer readable medium 2208 may include a hard disk, a floppy disk, an optical storage disk, a flash memory device, a Random Access Memory (RAM) device, a Read Only Memory (ROM) device such as an EPROM or an EEPROM, and the like. One of ordinary skill in the art will recognize other suitable signal processing devices 2206.
In one embodiment, the signal processing device 2206 may include a Programmable Logic Chip (PLC), a Digital Signal Processor (DSP) device, a Field Programmable Gate Array (FPGA), a microprocessor, a computer processor, or the like. Alternatively, the signal processing device 2206 may comprise a computer suitably programmed to execute the operations of the algorithm described above.
By way of comparison,
As illustrated in
Similarly, the laser surface reconstruction may be more precise than the microwave surface reconstruction, as shown in
Other aspects, advantages, and modifications are within the scope of the claims. In light of the described embodiments, one of ordinary skill in the art will recognize methods for adapting, e.g., an ultrasound transducer and/or a digital camera for surface estimation in accordance with the described system and methods.
This application claims priority to U.S. Provisional Application No. 61/058,135 filed Jun. 2, 2008, the entire text of which disclosure is specifically incorporated by reference herein without disclaimer.
| Number | Date | Country | |
|---|---|---|---|
| 61058135 | Jun 2008 | US |
