Infrared multi-spectral camera and process of using infrared multi-spectral camera

[0001] This is related to and claims priority from U.S. Provisional Patent Application Ser. No. 60/616,800, which was filed on Oct. 6, 2004.

[0002] The present invention relates to cameras, particularly cameras that provide thermal images for analysis.

[0003] Multi-spectral cameras have long been used for many different applications. Multi-spectral camera systems typically include software for camera control, image acquisition, and image analysis, so that the imaged object can be used for some diagnostic purpose. For example, such images can be used for airborne reconnaissance and terrestrial observation, environmental characterization, and military applications such as target acquisition, camouflage penetration, and surveillance. These cameras are cryogenic, that is, they use liquid nitrogen or other coolants to reduce thermal noise present in the images so that the signal-to-noise ratio is adequate. Thus, these apparatus are not suitable for use as consumer devices.

[0004] Night-vision long infrared cameras are well-known, and are sold commercially by a number of companies. These conventional cameras are non-cryogenic cameras, that is, they are built without liquid nitrogen or other coolants. These cameras provide a visible image under low-light conditions, but they are blind with respect to spectrum. This limits conventional cameras such that they cannot be effectively used for certain applications. For example, conventional cameras cannot be used to quantitatively determine early tumor development. Rather, a middle infrared spectrum camera must be used for this application and, due to low signal-to-noise ratio of the resulting image, one must use liquid nitrogen to cool down the camera detector backplane.

[0005] It would be advantageous to provide a multi-spectral camera that is commercially available and preferably features a non-cryogenic design. It also would be advantageous to provide a process by which such a camera can be used to test for and diagnose medical conditions, such as cancer, for use in human and veterinary medicine. Understanding of genotype & phenotype of cancers had led to advances in diagnoses and treatments of cancer; nevertheless, the battle against cancer remained to be a major concern; for example, in the U.S. alone over 200,000 women and 1,500 men are diagnosed with malignant breast tumors every year. It would be advantageous to provide a convenient and reliable screening methodology with inexpensive telemedicine decision aids in households, to complement a public awareness in healthy living style.

[0006] Several active imaging modalities have already been adopted in hospitals and laboratories using, for example, targeted radiation spectrum of radio waves in f-NMI, X-rays for mammograms, and gamma-rays for PET. Modern satellite imaging & smart military ATR technology passively receive radiation from targets and are non-intrusive. Less expensive versions for ground operation can be utilized for household or lab purposes. Traditional thermal imaging known as the thermograms can be improved over the conventional practice of using one IR spectrum-integrated camera to take two pictures of a patient sitting in a dark cold room a few minutes apart, which practiced with a mixed reputation for decades. Because high-precision satellite-grade cameras operating cryogenically in seven spectral bands from visible to IR are capable of detecting a hidden tank hot engine in the desert over thousand miles above the earth, a less expensive and portable versions of IR two-color CCD cameras (mid IR band from 3 to 5 microns wavelengths and the long IR band from 8 to 12 microns in terms of semi-conductor bang-gap material), having a reduced size of about 512×512 pixels operated at a minimum resolvable temperature difference (MRTD) of about 0.02 to 0.1 degree Kelvin can be used for close-up imaging of patients by household caretakers. Since such non-intrusive imaging screening methodology is passive and therefore non-specific, a NIH/NIBIB benchmark would require anyone to conduct an experiment of human subjects of 100,000 healthy but risky population: (i) Institute Review Board (IRB), (ii) Consent Patient Form (CPF) & (iii) tabulated results, in terms of the Receiver Operation Characteristics (ROC), of which a control study for 0.1% disease rate over five years might result in 100 sick patients with probability of detection (PD) versus False Alarm/Negative Rate (FAR). These requirements were beyond available resources and commitment level, so instead, a reverse approach was chosen. Rather than taking a large group of healthy but high-risk volunteers, a volunteer group of sick patients was used, in collaboration with Thermal Scan Inc. U.S., and Pontifical Lateral University at Vatican City of Rome, Italy, and nearby hospitals. The uncertainty of initial search was avoided and the sick patients were tracked under chemotherapy using IR two-color cameras. Vatican involvement was beneficial because passive and non-intrusive screening for breast tumors was based on (i) a ten-fold higher risk factor for nulliperous women, according to recent epidemiology studies in Singapore and (ii) a smart ATR algorithm was demonstrated by a pair of satellite-grade spectral-grams in tracking the Wien spectral peak displacement law over time. A sequence of time-order snapshots documented that successful recovery history in a time-contract-animated video. Getting better is not necessarily equivalent to getting worse physiologically, even when identical cameras and a smart algorithm are used. Moreover, such a new generation of non-mercury-contact spectrum thermometer can be available for general purpose usage, for those who cannot frequently take expensive and intrusive diagnostic modalities. A physics model of two IR color spectral-grams is provided, and a brief overview of active imaging and second generation anti-angiogenesis drugs is provided for the convenience of interdisciplinary studies.

[0007] To facilitate interdisciplinary collaboration, a simple but representative overview of various imaging modalities and cancer drugs usages are provided. According to oncological practice, an excessive new blood vessel can be generated, via Vascular Endothelial Growth Factor (VEGF) & other molecular signaling receptors, to supply the metabolic needs of oxygen, glucose, and other nutrition to a rapid-growth malignant tumor site, known as angiogenic blood vessel generation, an effect common to the metastasis of most malignant tumors. For example, ionizing UV & IR radiation of sunburn-damaged ATR & ATM of cells causes gene defects at CHK1 & CHK2, respectively. On the other hand, BRCA2 gene-deficient cells caused instability of chromosomes due to spindle abnormal cytokines. A common mechanism of cancer was suspected to be DNA reproduction without apoptosis programming of death. Advances with molecule-tagged cellular imaging had improved doctors' ability to diagnose and treat patients: (i) functional-Nuclear Magnetic Imaging (f-NMI) imaged the blood “hemodynamics” following the metabolic need of oxygen, of which a new algorithm improved higher-order spectral correlation when spectral lines increased the resolution under an increased magnetic field & cost of f-NMI device; (ii) in addition to oxygen, glucose is also required conjugated with an unstable isotope capable of decaying into a Positron, which is annihilated locally inside a patient with an electron within sub-mm mean free path resolution, Emitting gamma rays, 0.5 MeV, in two opposite directions which provide a direct internal radiation projection imaging Tomogram in the so-called PET for the imaging sugar “glucodynamics”; (iii) red-light tagged florescence molecular imaged in contrast with the opposite green-light-tagged florescence; (iv) an improved safety margin of X-ray dosage of mammograms. These high-energy radiation image modalities are too sophisticated costly for typical household ownership and use. Current passive IR spectral-grams study augmented these modalities, as double-blind tests when determined to be successful, could be potential supplements at labs and in homes.

[0008] Because current reverse or backward study to track the progress of chemotherapy treatments of already-sick patients involved some new anti-angiogenesis cancer drugs, they were briefly reviewed under “starving tumors of blood” with drug and chemotherapy for the convenience of cross-disciplinary imaging experts. For instance, FDA approved Genetech's anti-body “Avastin”, injected to augment the chemotherapy, to extend a patient's life from 6 to 11 months. Some of new anti-angiogenesis drugs are already in Phase III trial and could be made available for treatments, which, rather than being injected, are small enough to swallow as pills and made more than just “starving tumors of blood”, but to aim at multiple molecular targets for inhibition or competition. Bayer's Sorafenib, discovered 4 years ago for kidney cancer, can inhibit tyrosine kinases: Raf and other receptors of VEGF. Another second generation drug, SUGEN's Sutant, was discovered a year earlier and was recently acquired by Pfizer, proved its inhibition of the protein produced by KIT oncogene of stomach cancer GIST patients, and demonstrated furthermore with some positive affects on breast cancer and other cancers. Moreover, passive IR spectral-gram tracking might help monitor nano-scale targeted-drug delivery systems which had made significant promises, including in-situ laser burning of cancer cells having rich folate receptors using the vitamin-folate-guided carbon nanotube for IR-absorption generating rapidly heat killing cancer cells without harming healthy tissue. All these new drugs and nano-technology were timely for a new epidemiology study revealing an elevated increase of kidney cancer patients on the east coast of the U.S. Another large category of cancer is breast cancer, which has struck over 200,000 in 2005 so far, with a mortality rate of about 18.8%. Breast cancer deaths are less than those attributed to lung cancer, which number over 80,000 per year at a 85% mortality rate. Nevertheless, a 0.1% of risk of Ductile Carcinoma In Situ (DCIS) increased to 1% for nulliparous mothers without children according to Singapore public health statistics & biannual reports of epidemiologist Frank Speizer et al.

[0009] According to an aspect of the invention, a process of performing a medical test includes taking multi-spectral images of an area of interest of a patient. The patient can be a human being or an animal, and can be known to be healthy or known to have health issues or problems.

[0010] Taking multi-spectral images can include taking substantially simultaneous images of the area of interest using a plurality of cameras, wherein each of the cameras provides an image in respective different spectra. The cameras can be, for example, infrared cameras, and the spectra can be infrared spectra. The plurality of cameras can be, for example, two cameras.

[0011] Alternatively, taking multispectral images can include taking images using a multiple-spectrum camera. For example, the multiple-spectrum camera can be a dual-spectrum camera, such as a dual-spectrum infrared camera. The camera can include any two of a long-infrared wavelength detector, a mid-infrared wavelength detector, and a short-infrared wavelength detector. Preferably, the dual-spectrum infrared camera includes a long-infrared wavelength detector and a mid-infrared wavelength detector. The mid-infrared wavelength detector can include a detector array having carbon nanotubes. The carbon nanotubes can be detector elements. For example, the carbon nanotubes can be tuned-bandgap carbon nanotubes. Each pixel of resolution of the detector array can include a balanced Wheatstone bridge circuit including one of the tuned-bandgap carbon nanotubes.

[0012] The camera can be cooled by electrical dissipation. Alternatively, the camera can be cooled by refrigeration. For example, the camera can include a backplane, which can be refrigerated. As another alternative, the multiple-spectrum camera can be cryogenically cooled, such as by using liquid nitrogen as a coolant.

[0013] A process of performing a medical diagnosis can include performing a medical test as described above, comparing the images to spectrograms of subjects having a known health issue, and diagnosing a health status of the patient based on a correlation of the images to the spectrograms.

[0014] A process of performing a medical prognosis can include performing a medical test as described above, wherein the patient has a known health issue of a particular type, comparing the images to spectrograms of subjects having the known health issue of the particular type, and providing a prognosis for the patient based on a correlation of the images to the spectrograms.

[0015] According to another aspect of the invention, a multi-spectral camera includes a long-infrared charge-coupled device, a mid-infrared detector array, and a control device that synchronizes operation of the charge-coupled device and the detector array. The mid-infrared detector array can include carbon nanotubes. The carbon nanotubes can be detector elements. For example, the carbon nanotubes can be tuned-bandgap carbon nanotubes. Each pixel of resolution of the detector array can include a balanced Wheatstone bridge circuit including one of the tuned-bandgap carbon nanotubes. Adjacent pixels of the detector array can be arranged for orthogonal polarization.

[0016] The multi-spectral camera can also include conductive members that cool the detector array by dissipation. Alternatively, the camera can include a refrigeration element that cools the detector array. As another alternative, the camera can include a cryogenic cooling element that cools the detector array. For example, the cryogenic cooling element can use liquid nitrogen as a coolant.

[0017] Preferably, the charge-coupled device and the detector array are co-axially aligned.

[0031] High-precision spy cameras on board satellites located thousands of miles away can precisely detect and image the hot engine of a car or tank, or a missile plume, with the help of aided target recognition (ATR) techniques using multiple spectrogram features of a spontaneous thermal emission. Likewise, by transferring such military technology “from tank to tumor”, we can discover a hidden ductile carcinoma in-situ (DCIS) of a patient in a close-up setting according to the angiogenesis heating effect of new blood vessels working to feed a fast-growing malignant tumor.

[0032] As shown in FIG. 1, Wien's displacement rule of Planck blackbody cavity radiations is a linear law in terms of the log-log plot of the peak radiation intensity versus the wavelengths, of which Einstein photon dispersion is a special case in a vacuum at the slope m=1. Spectrograms features abnormal {right arrow over (a)}=(a₁, a₂); body {right arrow over (b)}=(b₁, b₂) based on the Planck Radiation Spectrum Distribution whose mean values of mid IR band (3-5 micron) and long IR band (8-12 micron) wavelength.

[0033] Almost all portions of the electromagnetic wave spectrum have previously been explored for medical applications, for example, blackbody radiation spectrum (FIG. 1); a short wavelength at gamma ray: PET; and at X-rays: mammogram; and radio waves: f-NMI. All utilize some man-made radiation sources to actively probe patients with sophisticated equipment and imaging processing algorithms. A notable exception is passive thermographs, which utilize a single thermal camera to image a patient's self-emitted heat radiation in a dark cold room, similar to colorless night vision, and subsequently exam the patient again for any remnant hot spot after being cool down, as shown, for example, in FIG. 2. In this conventional practice, two recordings per session are necessary: the first recording is made soon after the patient undresses, and the second is made after some duration has passed. This procedure requires an embarrassing wait in a chilly room (typically cooled lower than 21° C.). The first image (left) was taken within 1 minute after the patient undressed while the second one was taken 10 minutes later, during which time a normal pair of breasts became blue cold (right).

[0034] FIG. 3 shows that in this example the right breast continues to emanate heat radiation in the cold dark room even after the average body temperature decreased, for example, by immersion of the patient's hands in ice water. The extra heat could be an indication of active blood vessels generated to feed a malignant tumor, known as the angiogenesis effect. However, conventional non-intrusive and affordable thermographs can only lead a medical professional to speculate as to the cause behind those remnant hot spots. According to the present invention, the single camera is replaced by two cameras in recording the temporal increase, with the dual IR spectrograms ratio to be the salient invariant feature of malignance.

[0035] The approach of the present invention is to improve the precision of two cameras and to derive from the equilibrium thermodynamics principle the invariant feature to modify its pattern recognition by unsupervised ATR, or unlabelled data classifier by Dude & Hart definition, as follows. A starting point is the feasibility of a pair of satellite-grade cameras taking dual IR spectrograms of patients under anti-angiogenesis treatments. A volunteer patient is exemplified in FIG. 4, who is tested according to the process of the invention and diagnosed subsequently to be DCIS zero stage in her right breast nipple. According to this two-camera multispectral infrared breast image (left mid-IR and right long-IR), only one instantaneous and simultaneous recording per session was sufficient for the unsupervised classification of sub-pixel super-resolution algorithm. This non-intrusive, dark-room private, and passive imaging permits one with conveniently tracking of the hemodynamics phase transition to follow the angiogenesis of the breast pre-cancer tumor development for a potential telemedicine decision of early pre-cancer intervention. Several passive and non-intrusive spectrogram snapshots are animated into a time-reversible video, and playing backward, give virtual cancer-development dynamics that might not be accurate but the demonstrated tracking capability justifies a direct screen study for pre-cancer detection.

[0036] According to Fourier's conduction law, heat always flows from a hotter inside region to a colder ambient outside region to reach an equilibrium condition no matter how deep the heat sources are hidden. In addition to different degrees of angiogenesis activity, the differences in the depths can also result in different total image intensity and image spot size. Thus, instead of dealing with the unpredictable total intensity value, it is advantageous to explore an invariant feature of the dual IR spectrograms ratio. It is assumed that an unknown mixture of benign or malignant DCIS tumors exist in an arbitrary metastasis activity state. They might be located at an unknown depth within normal breast body tissue. According to Planck's blackbody radiation law, a healthy human body emits invisible IR radiation less efficiently than does an ideal blackbody oven kept at a constant temperature T_o=37° C.+273° K with a radiation leakage pinhole. We are brown bodies, so to speak, producing dual IR spectrograms: {circumflex over (b)}=(b₁, b₂) normalized at an equivalent heat source s′₂, as shown in FIG. 1. A malignant tumor of strength s′₁ can radiate as an equivalent blackbody at a slightly elevated temperature T_o=37° C.+Δ+273° K, whose mean values of dual spectrograms results in an abnormal feature vector, â=(a₁, a₂), with reference again to FIG. 1. Assuming two IR spectral cameras, perfectly registered & calibrated, we can take two spectral images at time t, resulting in two spectrograms. The image gray scale values at the corresponding pixel p are denoted as the spectral image vector: {right arrow over (X)}′_p(t)≡(x′₁(t), x′₂(t))^T, where the prime indicates the physical unit of spectrograms in watts per cm² per pixel. According to the thermal physics model, we assume a weakly linear mixture of two isotropic sources that generates an unknown mixture of IR dual spectrograms:

[0037] $\begin{array}{cc}{\overrightarrow{X}}_p^{\prime }\left(t\right)=\left[A_p^{\prime }\left(t\right)\right]{\overrightarrow{S}}_p^{\prime }\left(t\right);\left[A_p^{\prime }\left(t\right)\right]\equiv C_p^{\prime }\left(t\right)\left[\begin{array}{cc}{a}_1^{\prime }& b_1^{\prime }\\ a_2^{\prime }& b_2^{\prime }\end{array}\right];C_p^{\prime }\left(t\right)={\int }_{{\lambda }_1}^{{\lambda }_2}\underset{\mathrm{tumor}}{\int \int \int }\exp \left(\frac{-\overrightarrow{z}}{c_p^{\prime }\left(\lambda ,\overrightarrow{z}\right)}\right)R\left(\lambda ,T\left(t\right)\right)d^{3}zd\lambda ,& \left(1\right)\end{array}$
where the Planck emission distribution R is integrated over the long IR regime λ₂ to the mid IR band λ₁. These spectral intensities suffer diffusion and conduction loss with an exponential decay function c′_p(λ, {right arrow over (z)}) of the wavelength and the medium property sampled through by the radiation source at the depth |{right arrow over (z)}|. Equating the conservation law of energy, |{right arrow over (X)}′_p|≡√{square root over (x′₁²+x₂^′2)}=C′_p(s′₁+s′₂), dimensionless spectrograms are introduced without the prime {right arrow over (X)}_p≡{right arrow over (X)}′_p(x′₁²+x₂^′2)^−1/2 and normalized sources {right arrow over (S)}_p≡{right arrow over (S)}′_p(s′₁+s′₂)⁻¹; Eq. (1) can be rewritten with the unit magnitude mixing matrix as {right arrow over (X)}_p(t)=C_p(t)[{right arrow over (A)}_p(t)] S_p(t); or explicitly using matrix vector multiplication:

[0038] $\begin{array}{cc}\left(\begin{array}{c}{x}_1\left(t\right)\\ x_2\left(t\right)\end{array}\right)=C_p\left(t\right)\left[\begin{array}{cc}{a}_1& b_1\\ a_2& b_2\end{array}\right]\left(\begin{array}{c}{s}_1\left(t\right)\\ s_2\left(t\right)\end{array}\right)& \left(2\right)\end{array}$
where a₁=cos θ_a; a₂=sin θ_a and b₁=cos θ_b; b₂=sin θ_b. Since the spectral decay factor C_p(t) of an arbitrary pixel p cannot be computed from the first principle due to patients' variable physiques, the unknown tumor depth and thus the intensity should be eliminated by taking the ratio of spectrograms intensities r_x(t)≡x′₁(t)/x′₂(t). Moreover, to be further invariant to the imaging environment, the intensity ratio is inverted in terms of sources ratio ρ_s(t)≡s′₁(t)/s′₂(t):

[0039] $\begin{array}{cc}{r}_x\left(t\right)=\frac{a_1{\rho }_s\left(t\right)+b_1}{a_2{\rho }_s\left(t\right)+b_2};{\rho }_s\left(t\right)=\frac{b_2{r}_x\left(t\right)-b_1}{-a_2{r}_x\left(t\right)+a_1}{\Delta }^{-1}& \left(3\right)\end{array}$
If the mixing matrix were known, the inversion would be straightforward for a nontrivial determinant of two different unit feature vectors ∥Δ∥≡a₁b₂−a₂b₁≠0. However, in general, this unknown matrix inversion belongs to an ill-posed class of single-pixel blind source separation which consequently has many possible inverse solutions. Among all of these solutions, one should choose the dynamic equilibrium solution, which, by definition, would be realized most often experimentally.

[0040] Two unknown mixing angles of matrix [A_p] of Eq(3) remain to be determined by imposing two physics equilibrium laws:

[0041] (1) According to Einstein's theory of photons in a vacuum, light consisting of photons propagates with a constant speed of c_o, ε=ω=c_ok=c_oh/λ, and an increased photon energy ε₁ would result in a shortened wavelength ε₁/ε₂=λ₂/λ₁ inversely linearly proportional to an arbitrary reference state ε₂ and λ₂. However, in a real-world environment, such an energizing phenomenon could not happen in a vacuum and the Einstein formula must be modified according to the medium. In fact, Wien observed⁷ early that in Planck's every measurement of the radiation emitted from a blackbody cavity, all the spectral peaks at every equilibrium temperature fall on a linear negative slope, -m, on a log-log plot of the intensity versus the wavelength. This is known as Wien's displacement rule, a power law, shown in FIG. 1:
log ε₁−log ε₂=−m(log λ₁−log λ₂); ε₁/ε₂=(λ₁/λ₂)^−m,   (4)
where 1>m>0 is universal for all blackbody temperatures, which is consistent with Einstein's photon in the vacuum at Wien's power index m=1. Although Wien's index is universal for all blackbody cavity radiators at any temperature, a malignant tumor inside a human body is not as efficient as the ideal cavity radiator and further it cannot exist alone without a feeder source. For example, the infrared (IR) spectrum of a malignant tumor might be calibrated to be a brown body radiator m=½, that is, λ₁=λ₂(ε₁/ε₂)^−1/m=λ₂/(ε₁/ε₂)²; if the activity energy increases by 40%, a factor about ε₁=1.41ε₂≅√{square root over (2)}ε₂, the wavelength will be shortened by a factor 2, shifting from a long IR λ₂ (8-12 □m) at the ground state s₂ toward a mid-IR λ₁ (3-5 □m) at the excited state s₁. The local temperature raises T₁-T₂ due to the increased energies ε₁-ε₂ depending on the tumor's specific heat capacity, which can be estimated theoretically by integrating over the spectral density of tumor excited states s₁: n₁=dλ₁/ds₁ of which each degree of freedom contributes about K_BT/2, about 1/80 eV, at a warm room temperature. Nevertheless, such a change is often minutia and imperceptible to the eye; however, a pair of modern satellite cameras can detect the miniscule change by analyzing the dual infrared (IR) spectrograms images. Wien's displacement rule of the spectral peaks of Planck blackbody radiation distribution, see FIG. 1 and Eq. (4), states that hotter sources have their peaks shifted linearly and self-similarly from a long IR regime toward a middle or shorter IR regime. It has been demonstrated that Wien's spectrum shifting rule could be a salient feature of a decrease or increase of angiogenesis effect. Thus, as computed from Eq. (1), the differential slope rule of peak radiation is exactly the finite difference rule:

[0042] $\begin{array}{cc}\frac{d\left({\varepsilon }_1/{\varepsilon }_2\right)}{d\left({\lambda }_1/{\lambda }_2\right)}=-{m\left({\lambda }_1/{\lambda }_2\right)}^{-m-1}=-m\frac{\left({\varepsilon }_1/{\varepsilon }_2\right)}{\left({\lambda }_1/{\lambda }_2\right)}& \left(5\right)\end{array}$

[0043] (2) Thermodynamic equilibrium occurs at the real and non-negative Helmholtz free energy at the minimum, H_p=E_p−T_oS, illustrated in FIG. 5. In the thermal equilibrium of an open dynamic system at temperature T_o, the Helmholtz free energy H=E−T_oS should be at the non-negative minimum of which the approximation linear internal energy E (Taylor expanded near the equilibrium value) can intersect at zero, one, and two points the entropy S which has a simple convex function and maximum at equal source component s₁=s₂=0.5. A specific T_o is chosen so that only one intersection provides a unique answer. The single pixel radiation information energy E_p of an open system at a local equilibrium temperature T_o should be subtracted the Shannon-Claudius entropy valid only for a closed-equilibrium system at a maximum entropy: S/K_B=−s₁ logs₁−s₂ logs₂ normalized for two component case, s₂=1−s₁. It is assumed that the unknown internal energy E_p is analytic and expanded in a Taylor series with respect to the output feature vector {right arrow over (S)}_p near the correct inverse solution {right arrow over (S)}^(o)_p=[W]_αβX_β where the mixing matrix [A_p]⁻¹=[W].

[0044] $\begin{array}{cc}{H}_p=E_{\mathrm{po}}+{\mu }_{\alpha }\left(s_{\alpha }-S_{\mathrm{P\alpha }}^{\left(o\right)}\right)+\left(s_{\alpha }-S_{\mathrm{P\alpha }}^{\left(o\right)}\right){\varepsilon }_{\mathrm{\alpha \gamma }}\left(s_{\gamma }-S_{P\gamma }^{\left(o\right)}\right)+K_B{T}_o\left(\sum _{i=1}^m{s}_i\log s_i+\left({\mu }_o-1\right)\left(\sum _{i=1}^m{s}_i-1\right)\right)& \left(6\right)\end{array}$
where a repeated Greek index for the summation, and

[0045] ${\mu }_i=\frac{\partial E}{\partial s_i}{❘}_{s_i=s_i^{\left(o\right)}}$
was the Lagrange constraint, and

[0046] ${\varepsilon }_{\mathrm{lj}}=\frac{{\partial }^{2}E}{\partial s_j\partial s_i}{❘}_{s_i=s_i^{\left(o\right)}}\cong {\mathrm{\varepsilon \delta }}_{\mathrm{ij}}\ge 0,$
an assumed isotropic negative curvature for the convergence. To solve for s₁, one computes from the analytical assumption the Taylor series expansion (i) the linear contribution of information radiation energy from Eq. (6) to give the partial differential slope to be exactly equal to a finite difference of the Lagrange components:

[0047] $\begin{array}{cc}{E}_p-E_{\mathrm{po}}=\frac{\partial E_p}{\partial s_1}\left(s_1-s_1^{\left(o\right)}\right)=\left[\begin{array}{cc}{\mu }_1& {\mu }_2\end{array}\right]\left[\begin{array}{c}{s}_1-s_1^{\left(o\right)}\\ s_2-s_2^{\left(o\right)}\end{array}\right]=\left({\mu }_1-{\mu }_2\right)\left(s_1-s_1^{\left(o\right)}\right);\Rightarrow \left({\mu }_1-{\mu }_2\right)=\frac{\partial E}{\partial s_1}& \left(7\right)\end{array}$
where use was made of s₂^(o)≡ŵ^T₂{right arrow over (X)}=1−s₁^(o); (ii) one computes the partial differential of information energy with respect to the malignant source and obtains by the chain rule:

[0048] $\begin{array}{cc}\frac{\partial E}{\partial s_1}=\frac{\partial {\varepsilon }_1/{\varepsilon }_2}{\partial {\lambda }_1/{\lambda }_2}\frac{d{\lambda }_1/{\lambda }_2}{d{s}_1}{❘}_{s_1=s_1^{\left(o\right)}}=-\mathrm{mn}\frac{{\varepsilon }_1/{\varepsilon }_2}{{\lambda }_1/{\lambda }_2};& \left(8\right)\end{array}$
where the wavelength density of the malign states

[0049] $n\equiv \frac{d{\lambda }_1/{\lambda }_2}{d{s}_1}{❘}_{s_1=s_1^{\left(o\right)}}$
is related to the malignant tissue heat capacity; (iii) the minimum Helmholtz free energy at the isothermal equilibrium. Setting the partial differentiation of H to zero:

[0050] $\frac{\partial H}{\partial s_j}={\mu }_j+2\varepsilon \left(s_j-s_j^{\left(o\right)}\right)+K_B{T}_o\left(\log s_j+1+{\mu }_o-1\right)=0;$
and imposing the probability percentage normalization

[0051] $\sum _{i=1}^m{s}_i=1$
to eliminate μ_o. Finally, one obtains, at equilibrium, the solution s_j=s_j^(o)≡[W_jα]{right arrow over (X)}_α, the McCulloch & Pitts sigmoid logic, similar to artificial neural networks of isothermal brains,

[0052] $\begin{array}{cc}{s}_j={\left[1+\sum _{k=1,k\ne j}^m\exp \left(\left[{\mu }_j-{\mu }_k\right]/K_B{T}_o\right)\right]}^{-1}\equiv \sigma \left({\mu }_j\right),& \left(9\right)\end{array}$

[0053] In the two components case, the exact probability formula of the malignance s₁ has been derived from Eqs. (7,8,9):

[0054] $\begin{array}{cc}{s}_1={\left[1+\exp \left(-\mathrm{mn}\frac{{\varepsilon }_1/{\varepsilon }_2}{{\lambda }_1/{\lambda }_2}\right)\right]}^{-1}& \left(10\right)\end{array}$
The percentage of malignant source is mainly predicted in terms of the measured peak value of mid IR ε₁ at the wavelength λ₁ and the peak of long IR ε₂ at λ₂. Although the universal constant m for a blackbody can approximate our brown body, a realistic value m does not expect to vary appreciable from patient-to-patient. Also, a patient's tissue heat capacity is unlikely to change rapidly, in terms of the density of malignant source s₁ with respect to the wavelengths. Even without yet sufficient statistics of measurements, one can already verify the validity of the tumor formula in two limiting cases. (i) The weak source limit: mid IR ε₁<<ε₂ yielded ε₁/ε₂=0 and s₁=0.5 meaning the malignant tumor has 50% chance, of which the uncertainty can be resolved by subsequent observations further tracking the source ratio over days; (ii) The strong source limit: in the opposite limit ε₁>>ε₂ for strong mid IR and negligible long IR we have the certainty of the malignant tumor s₁=1.

[0055] The minimum H=0 occurs at E=T_oS where the approximation of information radiation energy E intersects at the convex entropy function S at only one point where the mixture temperature T_o was determined, as shown in FIG. 6. Given the percentage of sources s₁ & s₂=1−s₁, the vector parallelogram can determine those unknown abnormal spectral features and normal body spectral features to be added up to the data vector {right arrow over (X)}. It is clear that the minimum Helmholtz free energy H=0 implies in a closed system E=0, or, ε₁/ε₂=0 the maximum Shannon entropy at half chance of malignancy s₁=0.5 and half chance of benign status s₂=0.5. It is determined for

[0056] $\frac{{\lambda }_1}{{\lambda }_2}\le 0.5\mathrm{and}\frac{\partial {\lambda }_1}{\partial s_1}\ge \frac{\partial {\lambda }_2}{\partial s_1},$
one can experimentally estimate the inverse spectral density of malignant states

[0057] $n=\frac{\partial {\lambda }_1/\lambda 2}{\partial s_1}={\frac{1}{{\lambda }_2}\left[\frac{\partial {\lambda }_1}{\partial s_1}-\frac{{\lambda }_1}{{\lambda }_2}\frac{\partial {\lambda }_2}{\partial s_1}\right]}_{s_1=s_1^{\left(o\right)}}\approx \frac{1}{{\lambda }_2}\frac{\partial {\lambda }_1}{\partial s_1}{|}_{s_1=s_1^{\left(o\right)}}$
Once the minimum free energy shown in FIG. 5 is used to determine the percentage of malignancy s₁, the vector parallelogram shown in FIG. 6 determines the unknown unit feature vectors â,{circumflex over (b)} followed finally by the invariant source ratio ρ_s(t) from Eq. (3). FIG. 7 shows two independent classes discovered with DCIS zero stage near the ring around the right nipple, but not in the left nipple, when the unsupervised classification algorithm (Equations (1) through (10)) was applied to the image of the right breast. Independent classes represented good thermal classes, since most large heat classes came from inside of the breast. The marked area on the right breast indicated the existence of a DCIS of stage #0 (confined) to stage #1 (local spread). When two-color IR spectral-grams were augmented with the help of X-ray based mammography, which could detect micro-calcification—areas of cells of a few millimeters or more in diameter, which had been destroyed by cancer. That is, in the unsupervised classification images of the right breast, red means class of high probability (1) and blue means class of low probability (0). The broken ring of small red pixel dots less than millimeter size each and connected right outside quadrant, marked with the cyan circle, sharing the same texture heat supply of shallow capillary blood vessels as the rest but should not be there since the nipple did not usually have the abnormal isolated dotted characteristics unless a stage zero ductile carcinoma in situ (DCIS) is present. This dual band infrared image serves as merely a telemedicine super-resolution decision aid to doctors, which would require an intelligent data basis tracking over months or at least weeks to be ascertained by other intrusive means. Independent classes represent usual thermal diffusion Gaussian classes since most large heat classes come from the normal blood vessels of the breast. It has been derived for an open system that this generalized information theory, min H, could capture both neural network sigmoid logic as open dynamic system isothermal partition function and also the Hebbian unsupervised learning rules, Δ[Wij]=X_iμ_j. Given input data, the output is not a desired output, rather the internal Lagrange variable, whose sigmoid squashed output was the desired feature vector. By the dimensionality analysis, the synaptic weights were volts mediated by mini-volts neuro-transmittents and then for the physical power energy the internal Lagrange variable must be amperes representing the dendrite ion channel pico-amperes mediated by house-keeping glial cells. Two passive IR spectrum image data {right arrow over (X)}_p(t) the unknown feature vector {right arrow over (S)}_p(t) was extracted without external teachers, as the percentage of mid IR band versus long IR band in proportion to malignant versus benign tumors. This passive tracking of tumor-shrinkage by spectrograms might reduce the check-up frequency of X-ray mammograms.

[0058] Several remarks are in order: the study was perfromed (i) to gain the confidence in dual use of military ATR technology to public health, before many more resources were committed for comprehensive controlled studies in terms of the ROC. The study was done with volunteer patients under auxiliary passively IR spectral-grams imaging during their return visits for drug chemotherapy treatments; (ii) during treatment, the passive IR spectral-grams might provide doctors or caretakers a real-time insight for prescribing an appropriate chemotherapy dosage; and (iii) to reduce hospital return checkup frequency with a potential detection of any deadly recurrence of cancer here and elsewhere after the initial recovery. It would be beneficial to supplement that active mammogram imaging with more frequently passive IR spectrograms in-between the hospital treatment and checkups. (iv) The final sequence was animated in a time-reverse video to document the getting-better-to-complete-recovery history, which, when played backward, simulates a video of a high-risk patient getting sick as an earliest possible detection by passive spectrograms. When the predictions of these two spectrograms images were compared with the oncologist prognoses, the results consistently gave us the confidence of unsupervised ATR performance with IR dual spectrograms. Of course, the physiology change of getting better is different than that of getting sick, but the utility of IR dual spectrograms for passive screening is advantageous. (v) Modern satellite imaging is more reliable, and is passive in order to be stealthy, which qualities are suitable for screening public health status because it is a non-intrusive procedure. However, ordinarily satellite cameras are precision instruments operated cryogenically using liquid nitrogen coolant in seven or more spectrum color bands (from visible light to invisible IR) and require a supercomputer for processing. According to the present invention, the dual IR spectrograms reduced the number of satellite cameras to two cameras. Resulting spectral images are analyzed with a personal computer having an unsupervised classifier, to automatically extract the necessary features without the inconvenience of an expert-in-the-loop to adjust the threshold. In this preliminary study, commercial-off-the-self (COTS) spectral cameras were adopted. Initially, blind-controlled studies of a healthy 10,000 but risky 1% patients over several years were avoided in order to plot the results in terms of the Receiver Operation Characteristics (ROC) of the probability of detection of 100 sickness incidents versus the false negative rate. Rather than relying on the total intensity thermographs and its associated variation over time, as shown in FIGS. 2 and 3, applications in multi-spectral remote sensing on Landast seven multispectral band images were demonstrated. The unsupervised classification method described in Equations (1) through (10) were demonstrated to be capable of discovering small man-made objects located sparsely in a desert when the objects exhibit similarly-shaped spectral intensities as they would if located in a city area. According to the present invention, results of the same algorithm were shown to apply to unsupervised classification of the multi-spectral IR breast images for early breast cancer detection and tracking. Moreover, the design of the present invention enables satellite-precision cameras to be affordable and portable, not only for hospitals, but also for laboratory and household use. An electrically cooled dual-spectrum IR camera using an optically co-axial unit-frame is provided according to a biomimetic fovea. The imaging backplane houses both the long IR wavelength Charged Coupled Device (CCD) and a single quantum detector capability at the mid-IR wavelength in terms of one dimensional (1-D) Carbon Nanotubes (CNT) per pixel, which has a minimum occlusion about a nano-size in front of the CCD. The 1-D nature of CNT produces a reduced thermal noise, about ½ K_BT compared to 3/2 K_BT. Thus, electrical diffusion cooling preferably is used rather than liquid nitrogen to keep a steady backplane environment to maintain the minimum resolvable temperature difference (MRTD) similar to that of cryogenic mid-IR camera, about 0.02 degree Kelvin.

[0059] A nano-robot can be used to assemble one-dimmensional quantum mechanical band-gap material, such as carbon nanotube, at the back-plane. For example, see U.S. Pat. No. 6,862,924 to Ning Xi as an example of a device used for such nanomanipulation. The middle infrared detectors are designed to be located above the long infrared CCD, which uses the x-y plane row-sum column-sum read-out.

[0060] This architecture is similar to that of a human visual retina, which detects blue in front, green in the middle, and red behind, but read out along the z-direction. Similar to human eyes, the detector has almost single-photon detection capability using a Wheatstone bridge with 4-armed balanced circuitry to read out along z-direction pixel-by-pixel, which in turn drives an electrical current provided by a battery only when one or two of the arms, made of carbon nanotubes, receive middle infrared photons and break the balance.

[0061] Thus, the multi-spectral camera shares co-axially identical infrared optic lenses which focus naturally at multiple focal planes for different spectral wavelengths, similar to a multi-color fovea architecture. This is possible without usual color filter loss, because the 1-D quantum detector carbon nanotube (CNT) is on the order of 1 nanometer in diameter, which has a minimum occlusion of less than 1% from other radiation detectors and each can orthogonal to the other and selectively tuned for different spectral wavelengths. For example, two layers preferably consist of one mid-IR at 3 to 5 micron wavelengths and the other at a longer optical path provided by a standard un-cooled long-IR Focal Plane Array (FPA) with an intercept at 8 to 12 micron wavelengths. FIG. 8 shows an example of a two-spectrum fovea design for a 3-D FPA. The front FPA is carbon nanotube 1-D quantum detector for mid-IR wavelengths, and the back FPA is standard un-cooled CCD device for long-IR imaging. As shown in FIG. 8, in the architecture of the fovea design without filter loss, the occlusion of the CNT for long IR is less than 1%.

[0062] Orthogonal polarizations are arranged in a 2×2 pixels as shown in FIG. 9. FIG. 9 shows that unit cells of 2×2 pixels cover four orthogonal polarizations along a 1-D CNT for the electrical field direction. In addition, each pixel can detect a single photon by monitoring the balance of a Wheatstone bridge circuit, as shown in FIG. 10. As shown, a single-pixel single-photon read out by fovea cone single photon detector logic “negate the converse” implemented circuit, for example, a Wheatstone bridge balance circuit, is used per pixel. In the figure, (1) a band gap CNT is exposed on the x-y plane in a specific polarization direction; (2)-(3)-(4) a conductor CNT is in balance when no impinging photons are present; (5) a gain biased voltage is provided; and (6) a capacitor is provided for charge accumulation read-out along the z-axis. Thus, room temperature or non-cryogenic operation is possible for single-photon signal-to-noise ratio (SNR).

[0063] There are two types of CNT: conductors and semi-conductors. Nano-robotic assembly is possible at specific orientations and locations. For example, as shown in FIG. 11, two major categories of CNT include the one-dimensional quantum conductor known as the Armchair (n=m are CNT unit cell two chiral vectors); and band-gap semiconductors called Zigzag (n=0 or m=0 without being divisible by 3). Further, there are single-wall CNTs and multi-wall CNTs.

[0064] In general, a signal of 1000 to 100 photons provides statistically stable data, and therefore the usual SNR factor of 5 orders of magnitude is achieved at non-cryogenic cooling or room temperature operation. This is estimated as follows:

• [0065] CNT bandgap at Mid IR 3 to 5 micrometer
• [0066] Signal photon

[0067] $\Delta E=\mathrm{\hslash \omega }=h\frac{c}{\lambda }=0.414\mathrm{eV}\iff 0.248\mathrm{eV}$

• [0068] Between room temperature T=300° K; Liquid Nitrogen T=77° K
• [0069] Gaussian noise energy

[0070] $K_{B}T=\frac{1}{40}\mathrm{eV}=0.025\mathrm{eV};0.006\mathrm{eV}$

[0071] $\mathrm{Johnson}\mathrm{shot}\mathrm{noise}\mathrm{whose}\mathrm{mean}=\mathrm{variance}$ $1D\text{:}\frac{1}{2}{K}_{B}T<\mathrm{dark}\mathrm{current}<3D\text{:}\frac{3}{2}{K}_{B}T\mathrm{at}\mathrm{room}\mathrm{temperature}$ $0.0125\mathrm{eV}<\mathrm{dark}\mathrm{current}<0.0375\mathrm{eV}$ ${\mathrm{SNR}}_{\mathrm{room}}=\frac{0.4\iff 0.2}{0.01}=40\iff 20\mathrm{if}1\text{-}D\left(\mathrm{otherwise}\text{:}13\iff 7\mathrm{for}3\text{-}D\right)$ $\mathrm{at}77°K\Rightarrow \mathrm{noise}0.006\times 3/2\cong 0.01\left(1\%\mathrm{eV}\right)$ ${\mathrm{SNR}}_{\mathrm{cryogenic}}=\frac{0.4\iff 0.2}{0.01}=40\iff 20\mathrm{if}3\text{-}D$

[0072] FIG. 12 shows a carbon nanotube-based IR detector array. In the array array, each pixel includes a multi-walled carbon nanotube with a properly tuned bandgap for detection of a selected infrared spectrum.

[0073] A nano-robot can assemble one-dimensional quantum mechanical bandgap material such as a carbon nanotube at the backplane. Preferably, this is designed to be above the long infrared charge coupled device (CCD), which uses the x-y plane row-sum column-sum readout.

[0074] Particular exemplary embodiments of the present invention have been described in detail. These exemplary embodiments are illustrative of the inventive concept recited in the appended claims, and are not limiting of the scope or spirit of the present invention as contemplated by the inventor.