Conventional mammography screening is the only modality that has been shown to reduce the chance of death from breast cancer in randomized control trials. Depending on the breast tissue composition, however, the chances of detecting early breast cancer are different when using conventional screening mammograms.
The breast is made up of a mixture of tissues, including fibrous connective and glandular tissues, as well as fatty tissue. Radiologists classify breast density using a four level density scale, L1 through L4, where L1 describes the lowest breast density and L4 describes the highest breast density. Having breasts with relatively dense tissue (i.e., level 3 and 4) not only may increase the risk of getting breast cancer, but may also increase the difficulty of detecting breast cancer when using two-dimensional (2D), x-ray based screening mammograms. This is because the low contrast (conspicuity) between fibrous or glandular tissue and cancerous tissue is only of about 1 percent at radiological frequencies both generally appear white at x-ray frequencies. Women with breasts consisting of dense tissue may use adjunctive (additional) imaging test to effectively increase the conspicuity, in order to increase the likelihood of detecting early cancers.
One adjunctive imaging test is the bilateral whole breast sonography (i.e., ultrasound imaging). In women with dense breast (heterogeneous dense and extremely dense), screening ultrasound may detect additional cancers to those discovered with conventional mammography, similar to that seen with other adjunctive tests. Nevertheless, the number of unnecessary biopsies resulting from ultrasound imaging increases to an unacceptable level when it is used within a general population (i.e., dense and not dense breast) as a screening test. For this reason, ultrasound is only used as an adjunctive imaging test for women with an extremely dense breast.
Another adjunctive test is the three-dimensional (3D) mammography, also known as Digital Breast Thomosynthesis (DBT). Contrary to 2D Conventional Mammography (CM), DBT collects multiple views of the compressed breast, thus enabling a 3D reconstruction. DBT enhances the likelihood of finding tumors by removing overlapping breast tissues, i.e., those located at different 2D slices. The radiation dose to the breast resulting from the use of 2D mammography and DBT is increased by a factor of 2 as compared to CM alone, although the radiation dose is still below the U.S. Food and Drug Administration (FDA) limit.
Another adjunctive imaging modality is the bilateral breast magnetic resonance imaging (MRI), although MRI requires an intravenous injection of a contrasting agent and introduces an elevated cost.
Existing modalities used for early detection of breast cancer suffer from several limitations. First, existing modalities are not optimal to analyze all four levels of breast density (L1-L4). Second, they require the use of either ionizing radiation (e.g., X-ray based), which may have potential carcinogenic effects, or potentially hazardous contrasting agents (e.g., MRI-based). Third, current CM requires the use of adjunctive imaging modalities for dense breasts (L3-L4), and adjunctive imaging may result in additional cost, unnecessary anxiety and biopsies. Fourth, current imaging modalities and their associated data analytic tools do not take into consideration the particular specifics of the patients, such as positive genetic mutations (e.g., BRCA1 and BRCA2), breast density, family history, among others. Fifth, current modalities are not capable of providing low cost and size equipment that features high contrast and resolution. Sixth, in the case of multimodal systems that provide additional information concerning the breast tissue composition (such as density, compressibility, viscosity, dielectric constant, and conductivity) to 2D/3D mammography or MRI, such multimodal systems are not capable of collecting all data in a single and quick session, and in a co-registered fashion.
X-ray based technologies such as CM and DBT are most often used to detect cancerous lesions within the breast. Unfortunately, these systems both suffer from the aforementioned low radiological contrast between healthy breast tissue and cancerous tissue. As a result, these technologies tend to produce a large number of false positives when used for early detection.
Nearfield Radar Imaging (NRI) is a less common technology for breast cancer detection that uses non-ionizing microwave radiation to assess the breast tissue. NRI is an appealing technology for breast cancer detection because, at microwave frequencies, the contrast between healthy breast tissue and cancerous tissue is on the order of 10 percent. Unfortunately, standalone NRI systems typically struggle to accurately detect cancerous lesions due to the heterogeneous distribution of tissues within the breast, having a wide dynamic range in the relative dielectric constants that ranges from 5 (fatty) up to 55 (fibrous) for healthy tissues, and about 60 for cancerous tissues at around 1 GHz.
The described embodiments relate to the use of a hybrid, multimode, breast cancer detection system. The detection system may implement a fusion of different imaging modalities, for example Digital Breast Tomosynthesis (DBT) together with Microwave Nearfield Radar Imaging (NRI). In such a hybrid system, the DBT reconstruction may be used in order to form a prior distribution of tissues in the breast that can be used as a starting point for the NRI inversion process, although other imaging modalities may be used to provide this starting point. A linear linearized matrix of the non-linear in nature sensing problem may be formulated and solved, by applying for example the Born Approximation, to provide the complex permittivity of the breast tissues, although other techniques like iterative Born or the Contrast Source method can be used in the process of inverting the non-linear problem.
The described embodiments may detect cancerous tissues, or morphologically atypical tissues in general, with respect to healthy tissue. Morphologically atypical tissue may include, but is not limited to, cancerous tissue, tumorous tissue, precancerous tissue, and invasive tissue.
In addition, the reconstruction process may be formulated as a sparse recovery problem, such that certain compressive sensing (CS) techniques can be applied. See, for example, PCT/US16/25274 (entitled “Compressive Coded Antenna/Meta-Antenna,” filed on Mar. 31, 2016, the contents of which are incorporated by reference in their entirety) for examples of such CS techniques.
Imaging results, using real DBT data and synthetic NRI data demonstrate a tremendous potential for this technique in terms of finding tumors surrounded in fibroglandular tissue.
In the described embodiments, a multimode hybrid cancer detection system may include a fusion of two or more of (i) NM, (ii) DBT, (iii) Ultrasound Imaging (USI) and (iv) Thermoacoustic Imaging (TAI). In such hybrid embodiments, the fusion of the two or more modes is accomplished by co-registering the various modes' scans, so that the scans of all modes are captured with respect to the same physical configuration of the breast under study. The fusion is further accomplished by evaluating the captured multimodal scan data jointly rather than independently. An example of a hybrid DBT and NM system may be seen in PCT/US2014/042842, which is incorporated herein by reference in its entirety.
The captured multimodal scan data is collected while the breast is under clinical compression. Then, this data is used to reconstruct the pixel-based biophysical parameters of the breast—which may be given, for example, in terms of the porosity, fluid saturation, and solid matrix composition by using forward models for each modality and solving a joint non-linear inverse problem. The inversion makes use of a biophysical model that relates the biophysical parameters and the constitutive parameters for the different technologies.
The constitutive parameters may include the following:
The biophysical parameters (porosity, fluid saturation, and solid matrix composition) uniquely identify different types of tissues inside of the breast for each pixel in which it is discretized: fatty tissue, fibrous tissue, cancerous tissues, calcifications, and so on. The multi-modal system is more capable of finding tumors at earlier stages of the disease due to at least the following reasons:
In one aspect, the invention is a cancer detection system, comprising at least two imaging systems, each of which implements an imaging modality different from others of the at least two imaging systems, and each of which provides sampled image data based on its modality. The cancer detection may further include a processor and a memory with computer code instructions stored thereon, the memory is operatively coupled to the processor such that the computer code instructions, when executed by the processor, may cause the system to implement:
In one embodiment, the at least two imaging systems may include a Digital Breast Tomosynthesis (DBT) system and a Microwave Nearfield Radar Imaging (NRI) system.
In another embodiment, the at least two imaging systems includes two or more of (i) a Digital Breast Tomosynthesis (DBT) system, (ii) a Microwave Nearfield Radar Imaging (NRI) system, (iii) a UltraSound Imaging (USI) system, and a (iv) Thermoacoustic Imaging (TAI) system.
In another embodiment, the information from each modeling unit may include biological tissue parameters. The biological tissue parameters may include one or more of (i) electrical permittivity, (ii) permeability (iii) conductivity, (iv) elastic bulk modulus, (v) density, (vi) attenuation, (vii) thermodynamic heat capacity (viii) volumetric expansion coefficient, and (ix) radiological X-ray absorption.
In one embodiment, each modeling unit may include a biophysical model, a constitutive model, a forward model, and a field simulating model.
In another embodiment, the classifier may further classify tissues corresponding to the image data based on an unmixed version of the set of joint biophysical properties.
In another embodiment, the classifier may utilize a machine learning procedure to classify the tissues corresponding to the image data. The classifier may utilize a Quadratic Discriminant Analysis procedure (or other supervised learning classifier know in the art) to classify the tissues corresponding to the image data.
In one embodiment, the at least two imaging systems may reside on a mechatronic system that is integrated with a Digital Breast Tomosynthesis (DBT) system, such that all captured image data is co-registered.
In another aspect, the invention is a method of detecting cancer, comprising performing an imaging modality that is different from others of the at least two imaging systems, using each of at least two imaging systems. The method further includes providing, from each of the at least two imaging systems, sampled image data that is based on the image system's modality. The method may further include using a processor and a memory with computer code instructions stored thereon, producing modeled image data, for each imaging modality, based on a common set of biophysical parameters, reconstructing, based on information received from each modeling unit, a set of joint biophysical properties, and revising the common set of biophysical parameters based on the set of joint biophysical properties. The method may further comprise comparing the sampled image data from each of the imaging systems to the corresponding modeled image data to determine a difference between the sampled image data and the modeled image data and determining when the difference is less than a threshold difference, and classifying tissues corresponding to the image data as healthy or cancerous, based on the set of joint biophysical properties.
One embodiment may further include sequentially activating each of the imaging systems while a test subject remains clinically advantageous position.
Another embodiment may further include mechanically rotating sensors of the two or more imaging systems, in conjunction with the activating, to accomplish co-registration of the two or more imaging systems.
Another embodiment may further include classifying tissues corresponding to the image data using a machine learning procedure.
One embodiment may further include classifying tissues corresponding to the image data using a Quadratic Discriminant Analysis procedure.
Another embodiment may further include implementing, for each imaging modality, a forward model that simulates fields corresponding to the imaging modality.
Another embodiment may further include repeatedly revising the common set of biophysical parameters until the difference between the sampled image data and the modeled image data is less than a threshold difference.
In one embodiment, the set of joint biophysical properties include one or more of (i) electrical permittivity, (ii) permeability (iii) conductivity, (iv) elastic bulk modulus, (v) density, (vi) attenuation, (vii) thermodynamic heat capacity (viii) volumetric expansion coefficient, and (ix) radiological X-ray absorption.
In another embodiment, processing modeled image data for each modality is accomplished with a biophysical model, a constitutive model, a forward model, and a field simulating model.
In another embodiment, performing the imaging modality using the each of at least two imaging system further includes the at least two imaging systems using at least two of (i) a Digital Breast Tomosynthesis (DBT) system, (ii) a Microwave Nearfield Radar Imaging (NRI) system, (iii) a UltraSound Imaging (USI) system, and a (iv) Thermoacoustic Imaging (TAI) system.
In another aspect, the invention is a cancer detection system, comprising at least two imaging systems, each of which implements an imaging modality different from others of the at least two imaging systems, and each of which provides sampled image data based on its modality. The cancer detection may further include a processor and a memory with computer code instructions stored thereon, the memory is operatively coupled to the processor such that the computer code instructions, when executed by the processor, may cause the system to implement:
The foregoing will be apparent from the following more particular description of example embodiments of the invention, as illustrated in the accompanying drawings in which like reference characters refer to the same parts throughout the different views. The drawings are not necessarily to scale, emphasis instead being placed upon illustrating embodiments of the present invention.
A description of example embodiments of the invention follows.
The teachings of all patents, published applications and references cited herein are incorporated by reference in their entirety.
The described embodiments are directed to a breast cancer detection system that uses a multimodal imaging configuration. The described embodiments may utilize a fusion of two or more imaging modes, including for example (i) Digital Breast Tomosynthesis (DBT), (ii) Microwave Nearfield Radar Imaging (NRI), (iii) Ultrasound Imaging (USI) and Thermoacoustic Imaging (TAI). The described embodiments may evaluate the captured multimodal scan data jointly rather than independently. The described embodiments may further utilize co-registration of the two or more imaging modes, which ensures that the scans of all modes are captured with respect to the same physical configuration of the breast under study, i.e., while the breast is under clinical compression. The co-registration avoids the necessity of registering the images from different independently-operating sensing modalities, which may lead to misfits in the joint inversion of the biophysical parameters.
In standalone imaging systems, the constitutive properties of the target are reconstructed by inverting the forward model operator. For instance, electromagnetic imaging, like NRI, uses the electric and magnetic fields E(r, w) and H(r, w) to estimate the electromagnetic constitutive properties permittivity, magnetic permeability, and electrical conductivity [{dot over (o)}(r,ω,σ), {circumflex over (μ)}(r,ω), {circumflex over (σ)}(r,ω)]=Oem−1{E(r,ω), H(r,ω)}; acoustic imaging, like USI, uses the acoustic pressure p(r,ω) to estimate the acoustic constitutive properties density, the attenuation factor, and compressibility [{circumflex over (ρ)}(r), {circumflex over (Q)}(r,ω), {circumflex over (κ)}(r,ω,Q)]=Oac−1{p(r,ω)}; and x-ray imaging, like 2D mammography and 3D DBT, uses the intensity I(r,ω) to estimate the absorption coefficient [{circumflex over (μ)}(r,ω)]=Odbt−1{/(r,ω)}.
In a fused multimode imaging system, a joint inversion operator is Oac/em/dbt−1{·} is used for the same purpose:
[ô(rω,σ), {circumflex over (μ)}(r,ω), {circumflex over (σ)}(r,ω),ô(r,ω,σ), {circumflex over (μ)}(r,ω), {circumflex over (σ)}(r,ω), {circumflex over (μ)}o(r,ω)]=Oac/em/dbt{p(r,ω), E(r,ω), H(r,ω), H(r,ω)}
This inversion is now jointly performed; and as a result the combined reconstruction is more informative than any one of the sensor alone (since it provides complementary information), more reliable (since it can exploit redundancy in the multiple images), more timely, more accurate and/or less expensive. Therefore multimode imaging has the potential to enable the detection of tumors with better sensitivity and/or specificity, and to allow a better classification of objects since it has access to more features.
Unfortunately, the latter approach does not consider that there is a single underlying physical property that relates all the constitutive properties. Specifically, the constitutive properties of each sensing modalities are related to the biophysical parameters pn through the nonlinear biophysical models Gem, Gae, Gdbt. In other words, each pixel in the imaging region is made of a mixture of fibrous-connective or glandular tissue, adipose tissue, and in some cases cancerous tissues that is specified by the biophysical parameters vector Pn (note that this tissue mixture is equivalent to use other equivalent parameters like porosity, fluid saturation, and solid matrix composition). This suggests that a direct inversion over the biophysical parameters should be more robust that an inversion over the constitutive parameters. Unfortunately, the biophysical models are non-linear, and ill-posed, which are difficult to invert. The described embodiments provide a solution to this problem by incorporating a robust mathematical formulation and inversion method to jointly invert the nonlinear forward operator and biophysical models in a unified step, so that an enhanced overall sensitivity and specificity are achieved.
Another useful opportunity provided by fused multimode imaging systems is that the information provided by each sensor can be combined in order to classify a pixel or an image region as healthy or tumor. Radiologists often make this decision based on the pixel intensity features of the image (e.g., a strong signal in a DBT image may be indicative of cancer in a fatty breast), morphologic features (e.g., specular morphology in a DBT image may be indicative of cancer in a dense breast), and functional features (e.g., a strong signal is achieved in an MM or PET device after injecting a contrast agent). Notwithstanding, radiologists may not use all the information contained in the reconstructed features when diagnosing a tissue as healthy or tumor. The described embodiments apply novel data analytics and machine learning to improve the sensitivity and specificity of the fused system by incorporating, in addition to intensity features, morphologic features and functional features, patient-specific factors when training the classifier and in the decision (i.e., breast density, family history and genetic testing (e.g., BRACA)) in order to outperform the state of the art.
The multi-modal measurements produced by the multimodal sensing system of the described embodiment is formulated through the nonlinear relationship y=f(x), where x is a vector that is related to the constitutive parameters of the sensing modality, this is ò(r,ω,σ), μ(r,ω), σ(r,ω), ò(r,ω,σ), μ(r,ω), σ(r,ω), μa(r,ω); and f(·) is also a nonlinear function of the constitutive parameters that describes the measurement process. The process recovers the vector x from the set of measurements y. Without any prior knowledge about the object of interest, the unknown vector x can take any value permitted by the laws of physics; and, in this case, it is difficult to accurately reconstruct x due to the ill-posed and nonlinear nature of the problem. However, if one introduces additional a priori known information to the problem (i.e., the object is constructed from a mixture of R different tissues, which is determined by the mixture values z1, z2, z3 contained in the vector Pn), one can recover the same problem in a lower dimensional space z1, z2, z3 by considering that the constitutive properties are related to the mixture through the following mapping x=h(z1, z2, z3)—note that h{·}=[Gem{·}, Gae{·}, Gdbt, {·}]. Since the problem is now resolved in the lower dimensional space, the ill-posedness of the problem is reduced and the reconstruction becomes more stable.
The use of a mechatronic system 104, that mechanically translates the NRI/USI/TAI probes (or a subset of them), enables the collection of a large number of measurements, thus reducing the ill-posedness of the collected data and enabling noise reduction by, in some embodiments, averaging consecutive measurements.
The a-priori information used to define a near-to-optimal first guess of the true biological parameters (porosity, fluids saturation, and solid matrix composition) may be obtained, in some embodiments, by inverting the biophysical model of the reconstructed x-ray absorption value at every pixel using a single-modality DBT imaging method. In other embodiments, the first guess of the true biological parameters may be obtained by an imaging modality other than the DBT system (e.g., NRI, USI or TAI).
The fusion of multiple modalities allows the co-registered classification of the biological tissues in terms of nine parameters: (i) electrical permittivity, (ii) permeability (iii) conductivity, (iv) elastic bulk modulus, (v) density, (vi) attenuation, (vii) thermodynamic heat capacity (viii) volumetric expansion coefficient, and (ix) radiological X-ray absorption. These nine features may be used to classify pixels as cancerous or healthy tissues by using basic machine learning (supervised, unsupervised and/or deep learning) classifiers. Additional spatial features may be added to the machine learning classifier. This knowledge may facilitate a breast tissue properties database, which may be used for better refining the biophysical models used during the inversion. A reconstructed vector of unmixed tissues may be used by a machine learning procedure (e.g., simple Quadratic Discriminant Analysis) in order to classify tissues under test as healthy or cancerous.
The joint inversion of the biophysical parameters reduces the dimensionality of the problem, leading to a more suitable inversion when compared with a joint inversion of the constitutive parameters.
The NRI/UST/TAT or the NRI/UST modalities avoid the use of ionizing radiation like the DBT imager, which leads to an imaging technology capable of generating high resolution images (sub-milliliter) showing high contrast between fibrous and cancerous tissues, while keeping the data collection time short (e.g., under 20 seconds).
The multimodal cancer detection system of the described embodiments may operate in a two-step fashion. In the first step, the breast is placed under clinical compression, and DBT measurements are recorded using low-dosage X-rays according to the procedures of the DBT imaging system. In the second step, the mechatronics system 104, which includes the NM, USI and TAT probes immersed in a bolus fluid, is mechanically scanned with respect to the same breast as that breast remains under the same clinical compression observed for the DBT measurements.
A set of biophysical models (gac/s, gta, gcm, gdbt) are used to relate the biophysical parameters and the traditional constitutive parameters of the different sensing modalities: 1) gac/s is a function that relates the acoustic parameters ρ(s), Q(r), {circumflex over (κ)}(r,Q) (density, attenuation, and bulk modulus) with the biophysical parameters; 2) gta is a function that relates the thermoacoustic dependent parameters, μ(r), σ(r), {circumflex over (ε)}(rσ), α((r), ρ(r), Q(r), {circumflex over (κ)}(r,Q) (electric permeability, conductivity, dielectric constant, ratio between heat capacity and compressibility, density, attenuation and bulk modulus) with the biophysical parameters; 3) em is a function that relates the electromagnetic parameters μ(r), σ(r), {circumflex over (ε)}(r,σ) (permeability, conductivity, and dielectric constant) with the biophysical parameters; and 4) as described above, dbt tht is a function that relates the X-ray parameters μa(r) and the biophysical parameters.
Once the constitutive parameters are known, a set of forward models (FWACS, FWTA, FWEM, FWDBT) (acoustics, thermoacoustics, electromagnetic and x-ray DBT) are used to synthetically predict the measured data: a) pressure ps(r) of the pp and s-waves for the acoustic model; b) electromagnetically induced pressure ps
If the synthetic data is similar to the measured data (this condition is given by a quadratic data misfit norm as well as a norm-1, norm-2 and/or norm-1,2 regularization term), then the method, which may be iterative, is stopped and Pn=[φn, Svn, Cvn] is used to compute all the constitutive parameters for the different technologies. If this condition is not satisfied, a non-linear inversion method (which may be based for instance on Born approximations, iterative born approximations, contrast source methods, Rytov methods, Eikonal inversions using norm-1, norm-2 and/or norm-1,2 regularization terms) can be used to derive the next iteration biophysical parameters Pn=[(φn, Svn, Cvn]. Since the inversion may be done at different scales for the different methods, an upscaling/downscaling technique is used to reach an unified-scale biophysical parameters. Multiple iterations of the above processing may be performed until convergence is achieved.
The final biophysical parameters and the nine constitutive parameters (electrical permittivity, permeability conductivity, elastic bulk modulus, density, attenuation, thermodynamic heat capacity volumetric expansion coefficient, and radiological X-ray absorption) may be used (e.g., by supervised, unsupervised and deep learning techniques) in order to classify pixels into the different types of tissues (classes) inside of the breast, which may include (among others) of the following tissues: fatty tissue, fibrous tissue, cancerous tissues, and calcifications.
In some embodiments, the biophysical models (ac/sta, emdbt) and the upscaling and downscaling techniques may be derived from Monte-Carlo simulations, as shown and described in
An electrical schematic of an example cancer detection system according to the described embodiments is shown in
The communications interface 610 buffers and formats the image data 606, 608 into a form suitable for transfer to a system bus 612. A processor 614 coordinates with the communications interface 610 to accept the image data and store the information 608 into a memory 616. The system may also include support electronics/logic 618, a network interface 620 for communicating with an external network 622, and a user interface 624 for communicating user information between a system user and the system bus.
The memory 616 also includes instruction code for execution by the processor 614 to perform system operations. The instruction code may include instructions for performing the processing such as image data modeling, joint non-linear inversion, scaling and machine learning, as described herein, and an operating system for coordinating and managing the compressive sensing image processor 626.
minimize ½∥Ax−y∥R
subject Re(diag(∈b)x+∈b)≥1
Im(diag(∈b)x+∈b)≥0
Although there are some artifacts in the image, the algorithm is able to locate the cancerous lesion.
In the numerical simulations, the breast geometries were excited by 17 transmitting and receiving antennas operating in a multistatic configuration. Each transmitting antenna operated at 11 frequencies linearly spaced from 500 MHz to 1500 MHz, for a total of 3179 complex measurements. Note that redundant measurements were used in the optimization routine. The healthy breast geometry simulations were used in order to generate the adjusted measurements ŷ and to compute the Jacobian matrix A required by the optimization procedure. The imaging region was constrained to 9654 positions in the breast, where the grid size of each pixel was 2 mm. In order to consider the problem in the most ideal scenario possible, random noise was not added to the measurements. As a result, the measurements were only corrupted by noise introduced into the problem when it was linearized via the Born Approximation. This noise was estimated to have 12.5% the energy of the adjusted measurement vector, i.e., η≅0.125∥ŷ∥12. In addition, the difference between the measurements of the unhealthy breast, y, and the measurements of the healthy breast, {tilde over (y)}=f(h(v1, v2, v3)), had approximately 12:69% the energy of the adjusted measurement vector, i.e. ∥y−{tilde over (y)}∥l
The mixture proportions are not exactly recovered, which is to be expected given that the true solution vector has an error of 0.125∥ŷ∥l
It will be apparent that one or more embodiments described herein may be implemented in many different forms of software and hardware. Software code and/or specialized hardware used to implement embodiments described herein is not limiting of the embodiments of the invention described herein. Thus, the operation and behavior of embodiments are described without reference to specific software code and/or specialized hardware it being understood that one would be able to design software and/or hardware to implement the embodiments based on the description herein.
Further, certain embodiments of the example embodiments described herein may be implemented as logic that performs one or more functions. This logic may be hardware-based, software-based, or a combination of hardware-based and software-based. Some or all of the logic may be stored on one or more tangible, non-transitory, computer-readable storage media and may include computer-executable instructions that may be executed by a controller or processor. The computer-executable instructions may include instructions that implement one or more embodiments of the invention. The tangible, non-transitory, computer-readable storage media may be volatile or non-volatile and may include, for example, flash memories, dynamic memories, removable disks, and non-removable disks.
While this invention has been particularly shown and described with references to example embodiments thereof, it will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the scope of the invention encompassed by the appended claims.
This application claims the benefit of and priority to U.S. Provisional Application No. 62/413,873, filed on Oct. 27, 2016, U.S. Provisional Application No. 62/412,671, filed on Oct. 25, 2016 and U.S. Provisional Application No. 62/248,041, filed on Oct. 29, 2015. The entire teachings of the above applications are incorporated herein by reference.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/US2016/059548 | 10/28/2016 | WO | 00 |
Number | Date | Country | |
---|---|---|---|
62413873 | Oct 2016 | US | |
62412671 | Oct 2016 | US | |
62248041 | Oct 2015 | US |