Abstract:
A hybrid system is developed using thermal and electrical impedance imaging methods together. The innovation of the approach relies on the frequency dependence of the tissue&#39;s electrical impedance which facilitates the acquisition of multiple thermal images with currents at different frequencies injected to the region of the body under inspection. Proposed method without current application (in passive mode of operation) provides images which are obtained by standard thermal imagers. On the other hand, the application of the electrical current (in active mode of operation) increases the temperature contrast on the body surface depending on the electrical property of tissue. Therefore, the technique while increasing the thermal contrast provides frequency dependent conductivity distribution data which can be used as a basis for the detection of the breast carcinoma. The sensitivity of the technique is increased by an infrared camera with dual band (MWIR/LWIR) imaging capability.

Description:
BACKGROUND OF THE INVENTION 
       [0001]    The present invention is mainly concentrated on the diagnosis of the breast carcinoma using a hybrid technique incorporating both thermal and electrical impedance imaging methods. The innovation of the approach relies on the frequency dependence of the tissue&#39;s electrical impedance which facilitates the acquisition of multiple thermal images with currents at different frequencies injected to the region of the body under inspection. This new technique improves the sensitivity of medical imaging over the standard thermography as follows. 
         [0002]    Current injection at a single frequency increases the temperature gradients in the unhealthy regions enhancing the visibility (detection probability) of the cancerous tissues with a thermal camera, 
         [0003]    Current injection at multiple frequencies along with thermal imaging with a thermal camera results in generation of multiple thermal images according to the frequency dependencies of the tissue conductivities since heat generation (measured by the thermal camera) due to current injection is related with the tissue conductivity. 
         [0004]    As a result, the technique allows the acquisition of the frequency dependent tissue conductivity distribution along with enhanced temperature gradients in the body regions with unhealthy tissues which, when interpreted simultaneously, are expected to considerably increase the sensitivity of medical thermal imaging. 
         [0005]    The invention can be used in all clinical applications of the thermal imagers: 
         [0006]    oncology (breast, skin, etc.), 
         [0007]    pain management control, 
         [0008]    vascular disorders, 
         [0009]    arthritis/rheumatism. 
         [0010]    neurology, 
         [0011]    surgery (open heart, transplant, etc.), 
         [0012]    ophthalmic. (cataract removal), 
         [0013]    tissue viability, 
         [0014]    thyroid (diagnosis and treatment), 
         [0015]    dentistry, 
         [0016]    monitoring efficacy of drugs and therapies, 
         [0017]    respiratory (allergies, SARS), 
         [0018]    sports and rehabilitation medicine. 
       PRIOR ART ABOUT THE INVENTION (PREVIOUS TECHNIQUE) 
       [0019]    Various imaging modalities are used for diagnostic purposes and monitoring the treatment processes. These modalities provide images of different physical properties of the body tissues. Table 1 shows a list of these modalities and corresponding physical properties: 
         [0000]    
       
         
               
             
               
               
               
             
           
               
                 TABLE 1 
               
             
             
               
                   
               
               
                 Medical Imaging Methods and Corresponding Physical Properties 
               
             
          
           
               
                   
                 Imaging Method 
                 Physical Property 
               
               
                   
                   
               
               
                   
                 X-ray Computerized 
                 Attenuation coefficient 
               
               
                   
                 Tomography 
               
               
                   
                 Ultrasound 
                 Reflection coefficient 
               
               
                   
                 Magnetic Resonance Imaging 
                 Proton density, transverse and 
               
               
                   
                   
                 longitudinal decay constant 
               
               
                   
                 Nuclear Techniques (PET) 
                 Capacity of holding radioactive 
               
               
                   
                   
                 materials 
               
               
                   
                 Electrical Impedance 
                 Electrical Impedance 
               
               
                   
                 Imaging 
               
               
                   
                 Infrared Thermal Imaging 
                 Emissivity and temperature 
               
               
                   
                   
                 distribution 
               
               
                   
                   
               
             
          
         
       
     
         [0020]    For the diagnosis of breast carcinoma, mammography is the standard test for breast screening. However, this method has accuracy problems especially for smaller tumors and it is not comfortable for patients since the breast needs to be squeezed to minimize the path length between source and detector. If the dimension of a tumor is lower than 1.6 cm, it cannot be reliably sensed by this technique, since X-rays passing through the tumor are not detectable. Statistical studies show that there 10-25% false negative (tumor is not detected following a mammogram) ratio in mammography. Another concern is the exposure of the patient to ionizing radiation (X-ray) which is potentially harmful for the human tissue. 
         [0021]    Ultrasound is another important imaging method that is used to evaluate the breast. However, this technique is not successful in detecting micro formations and monitoring deep regions of the breast. 
         [0022]    Magnetic resonance imaging (MRI) is another imaging modality used in breast imaging. However, it should be used as a complementary tool instead of utilizing as a self-sufficient screening tool. Although it is more sensitive than mammography, it has a higher false-positive rate which may result in unnecessary biopsies. At the same time, due to the relatively long data acquisition time, this technique is not the first choice by the patients. Currently, there is no gold standard in detecting cancerous tissue in the breast. To minimize diagnostic errors, multi-modal techniques are proposed that employ more than one technique to complement each other. 
         [0023]    Now imaging techniques have also been investigated. Electrical Impedance Tomography (EIT) and Infrared Thermal Imaging (origins of this invention) were developed for this purpose. However, these techniques are not also sufficient as a sole screening tool. They have been used as an adjunct to the conventional medical imaging techniques. 
         [0024]    In EIT, changes in the electrical properties of the tissues are used in the diagnosis of breast cancer. However, in EIT applications, the current applied on the surface of the tissue has the highest intensity beneath the current drive electrodes. When the electric current is applied from the surface, the current cannot be focused homogeneously in the desired body location. Therefore, desired amount of signal cannot be provided from the inner regions of the body. In addition, the number of measurements are limited by the number of current drive electrodes. This results in low-resolution images. To enhance the signal strength, the current intensity can be increased: however, the applied current should be below the medical safety limits. One way to improve the resolution is to increase the number of electrodes. However, this requires a more demanding hardware, since the signal level decreases and the interaction between the electrodes increases. 
         [0025]    Thermal imaging is another technique that may be used to screen breast cancer, it also gives complementary information about developing cancerous tissues. Abnormal breast thermogram shows a significant risk indicator for the breast cancer. Several studies report infrared imaging as an important risk marker. Thermal imaging gives information about metabolic activity and tumor induced angiogenesis instead of structural knowledge such as size of the tumor and micro calcifications. Infrared imaging may also show indications almost 10 years before the development of breast cancer, unlike other imaging modalities. When this important ability combines with its proven high specificity and sensitivity rates, thermal infrared imaging becomes an invaluable tool in breast cancer diagnosis. 
         [0026]    In thermal imaging, the radiation emitted from infrared region of the electromagnetic spectrum is measured by an infrared thermal imaging camera and temperature distributions corresponding to the measured thermal emission are displayed. Infrared imaging is the only imaging modality which can image the metabolic aspect of the breast. Thermal image is an indicator of abnormal pathophysiology and physiology. 
         [0027]    Thermal Imaging provides images related to two different physical properties of the tissue: temperature and emissivity. The method has the following properties: 
         [0028]    generates surface images instead of tomographic image, 
         [0029]    can be used as real-time imaging modality, 
         [0030]    portable, 
         [0031]    relatively inexpensive, 
         [0032]    relies on continuous investment for its development, based on other important uses (for industrial, military, and space research) 
         [0033]    it has no harmful effect on the human body due to the passive mode of operation, 
         [0034]    offers contactless and non-invasive imaging (no contact with the patient). 
         [0035]    Diagnostic capabilities of the method vary depending on the rumor size and the distance to the surface. In the early stages of the disease, diagnosis of smaller tumors is limited, especially for distant masses from the surface. 
         [0036]    Nyirjesy and associates implemented a study which compares the accuracy of clinical examination, thermography and mammography in the breast cancer detection. 16778 patients were participated in that study. There were 8757 asymptomatic patient, 4715 confirmed carcinoma patients and 3305 histologically diagnosed benign breast cancerous patients. In that study, the diagnostic performances of a mammography, thermography and clinical examination results are summarized (Table 2). 
         [0000]    
       
         
               
             
               
               
               
             
               
               
               
               
             
           
               
                 TABLE 2 
               
             
             
               
                   
               
               
                 Accuracy rates of different methods in the 
               
               
                 diagnosis of breast cancer (Nyirjesy 1986). 
               
             
          
           
               
                   
                 Specificity (%) 
                 Sensitivity (%) 
               
               
                   
                   
               
             
          
           
               
                   
                 Mammography 
                 73 
                 80 
               
               
                   
                 Thermography 
                 85 
                 88 
               
               
                   
                 Clinical examination 
                   
                 75 
               
               
                   
                   
               
             
          
         
       
     
         [0037]    In Thomassin et al. study, 130 malignant tumors ranging in diameter 3-5 mm were detected out of 4000 confirmed breast cancer cases. 10% of them were diagnosed by mammography, 50% of the patients were diagnosed by thermography and rest of them which could not be detected either mammography or thermography were diagnosed (40%) by multi modal imaging (mammography and thermography). As a result, many studies were reported that abnormal breast thermogram shows a significant risk indicator for the breast cancer. 
         [0038]    Infrared thermal imaging system can obtain high spatial resolution thermal images of the tissue without contact. Although it has higher sensitivity and specificity, smaller cancerous tissues cannot be detected from deeper locations since the temperature contrast on the breast surface due to metabolic heat generation of the cancerous tissue may not be sufficient to be sensed. Tumor tissue cannot produce adequate amount of temperature difference on the skin surface. While, the noise equivalent temperature difference (NETD) of the state-of-the-art thermal infrared cameras is as low as several mK, detection may not be possible especially for smaller tumors located at deeper regions. The thermal sensitivity can be increased by developments in the detector technology, but this requires a high cost and long-lasting R&amp;D work. 
         [0039]    A. W. Guy implemented an active thermographic study with tissue mimicking phantoms. In that study, electromagnetic fields at different frequencies are applied by various microwave sources as plane wave, aperture, slot and dipole sources, and the heating pattern is observed due to the electric field. 
         [0040]    M. Mital and E. P. Scott have studied on the diagnosis of embedded tumors with an infrared camera. They placed a heat source in a cylindrical biological phantom and imaged the temperature distribution on the phantom surface. The heater (resistance 32Ω) was embedded in the agar to simulate the tumor tissue&#39;s temperature increase. Power dissipation of the heat source was adjusted by changing the voltage of the heater. Infrared camera had 80 mK thermal resolution and the emissivity of the agar phantom (1% solution of agar) was 0.97. Genetic algorithm was implemented to be able to determine the heat generation rate and location of the source from the temperature distributions. 
         [0041]    F. J. Gonzalez presented a thermal simulation study of breast tumors. Woman breast was modeled with as 9 cm radius hemisphere and the spherical model of the tumor tissue was embedded in it. In that study, the smallest size of cancerous tissue is explored for each depth using state-of-the-art thermal infrared cameras. The results showed that 3 cm tumors can be detected from 7 cm depths and smaller tumors can be diagnosed if they are located closer to the skin surface. 
         [0042]    L. Jiang et. al. conducted simulation study on the thermography of the woman breast under elastic deformation. They modeled the woman breast by taking the gravity-induced elastic deformation into account. This was the first study adding this parameter beside the metabolic, heat generation, blood perfusion, thermal and electrical boundary conditions, and thermal and electrical conductivities of breast tissue. Numerical breast models were developed using the 3D finite element method by including the thermal and elastic properties of the breast tissue. Temperature differences due to tumor tissue were obtained numerically for both static and dynamic thermal states. The developed modeling technique establishes a connection between breast thermography and pathological and physiological factors due to cancerous tissue. 
         [0043]    Most of the clinicians and radiologists are well experienced on the structural anatomy and medical education but they are not experts on thermodynamics, thermokinetics and physical and biological basis of infrared imaging. This may appear as a problem for the effective utilization of thermography for medical imaging. However, technological developments in digital processing, image acquisition and user friendly software facilitate the use of infrared imaging. Multispectral images could be obtained with the improvements in narrow-band filtering. With dual band thermal cameras, two different images for the same body can be acquired at different wavelengths. Software advances in smart processing (artificial intelligence) make the post image processing of raw data possible. Artificial neural networks, automated target recognition and threshold algorithms are the most commonly used ones. Data summation is used in artificial neural network algorithm to obtain pattern recognition. On the other hand, threshold and automated target recognition algorithms use a normative data base. These improvements have raised the accuracy of the system and have reduced the need for skillful operators. 
         [0044]    A period of heating and then cooling was proposed (e.g. with the aid of the fan) in order to improve the diagnostic performance (Keyserlingk, 2000; Ohashi, 2000). Kachzmarek and Nowakowski (2004) proposed optical heating of tissues in the IR band. Arabi and Muttan (2008) have applied both optical and microwave energy for this purpose. Such non-destructive testing methods which are used in industry is called the active dynamic thermography, while different energy types are applied for this purpose (electromagnetic, vibration with ultrasound, laser, hot air, or microwave). The terminology of Active Dynamic Thermography is relatively new in the medical imaging literature. 
         [0045]    Cyril Smith&#39;s patent, named as Thermal Imaging using Electro-Magnetic Induction, is on active mode operation of thermal infrared cameras. In this study, the aim is determining the hidden and embedded objects in the environment (suitcases, soil, etc.) by thermal infrared camera with higher performance. 
         [0046]    System is based on, 
         [0047]    high-frequency (1-2 GHz) electromagnetic waves applied to the body, electromagnetic waves radiated from the body. 
         [0048]    interference pattern formed on an electrically conductive layer. 
         [0049]    In this method, a resistive metal layer is placed in front of a thermal camera. Thermal camera continuously displays the temperature distribution on this layer. Time-varying current (1-2 GHz) is applied in short time periods into the coil placed behind a layer. Magnetic fields created by the coil induce current in conductive objects, and magnetic fields created by the induced currents change the current distribution on the metal layer and therefore the temperature distribution. Hidden/buried object can be determined in a wide area by changing the position of the camera/metal sheet/coil assembly. This system has no clinical application. 
         [0050]    Zombo et. al. developed a system, named as “Hand held magnetic induction thermography system”. By moving a coil/IR camera system manually, the object to be displayed is scanned. Eddy-currents are induced in the object with the 100-1000V, 10-100 A current (10-1000 kHz) injected to the coil. Current and temperature distribution change in the defective areas of the material. This difference can be distinguished in the thermal images, and error location can be determined. Different current distributions in the object can be obtained by changing the coil axis, and the image contrast can be changed. No medical application was predicted. 
         [0051]    Herbert L. Berman and colleagues proposed an active thermographic system which is described in the patent US 20070213617 A1. The idea of this patent is based on the study of Kachzmarek and Nowakowski&#39;s (2004) optical heating of tissues. Breast tissue is heated by a halogen lamp (optic stimulation) for 30 seconds from 50 cm distance. In that study, this topic is embedded into various innovative ideas (dual-band imaging, the contact of the breast with the glass which is transparent to infrared wavelengths, heating with heater or it is the only active thermographic system which has a medical application. The additional temperature in the tissue is generated by the light applied from the surface. Another system, namely, “Infrared Defect Detection via Broad-Band Acoustics” was patented by Raulerson and etc. In this study, heating the object with wide-band (1 kHz-1 MHz or 10 kHz to 50 kHz) acoustic energy was proposed to improve the performance for industrial applications of thermal cameras (e.g. fault finding equipment). 
         [0052]    It is possible to access a large number of patents about active thermography, however, they refer to methods for different industrial applications and signal processing methods. The only patent that refers to medical applications is the Berman and colleagues&#39; patent which was taken in 2007 where active thermography was implemented using halogen light sources. 
         [0053]    In the present invention, external energy is applied to the body in a different manner. Electrical current is applied from the surface with electrodes yielding an additional heat source in the tissue. Since electrical conductivity and permittivity of tissue alter with frequency; multi-frequency current application provides a means to image electrical properties for different frequencies by measuring different temperature distributions on the body surface. Therefore, the technique is expected to be able to detect anomalies in terms of both temperature and tissue conductivity differences as may be displayed by the breast carcinoma. 
         [0054]    In 2009, Carlak et. al. showed that low frequency electrical currents applied from the surface provide significant performance improvements in the breast cancer diagnosis. In that study, simulation of a 2-dimensional electrical conductivity model of the breast slice was implemented, and it was shown that a 1 cm tumor located 4.5 cm away from the electrodes can be distinguished with the state of the art thermal imagers. 
         [0055]    In 2010, the same group was conducted a similar study with a 3-dimensional (cube-shaped) electrical conductivity breast model. It was shown that operating at a frequency of 10 kHz, a 40 cm 3  tumor in the breast tissue (heated 300 seconds) can be distinguished from a 1.5 cm depth from the surface at the axis of the camera. 
       SUMMARY OF THE INVENTION 
       [0056]    Invention is about a hybrid technique that incorporates thermal and electrical imaging methods for the diagnosis of cancerous tissues. In this system there are three main ideas:
       Using the electrical properties of the tissue for diagnostic aims,   obtaining different information/different images depending on the operation frequency of the applied current,   Imaging the surface temperature and emissivity parameters by measuring the infrared emission.       
 
         [0060]    Proposed method without current application (i.e., in passive mode of operation) provides images which are obtained by standard thermal imaging systems. However, when the electrical current is applied (in active mode of operation) the temperature contrast on the surface changes depending on the electrical properties of the tissue. Electrical current density below the medical safety limits is applied to the patient for this purpose. 
         [0061]    Since tissues exhibit different electrical properties, it is possible to obtain useful information on the anatomy of the human body and tissues&#39; health status by imaging, for example, the electrical conductivity distribution of the tissues. These images can be used to identify organs and evaluate their health status. The method, for example, can be employed to detect breast carcinoma. Table 3 shows the conductivity properties of several types of breast tissues for different operation frequencies. Note the relative change and the contrast between the tumor tissue and healthy breast tissue conductivities at different frequencies. 
         [0000]    
       
         
               
             
               
               
               
               
             
           
               
                 TABLE 3 
               
             
             
               
                   
               
               
                 Electrical conductivity (mS/cm) values of 
               
               
                 breast tissues at different frequencies 
               
             
          
           
               
                 Breast Tissue 
                 10 kHz 
                 100 kHz 
                 10 MHz 
               
               
                   
               
               
                 Tumor (center) 
                 3.98 
                 4.28 
                 7.24 
               
               
                 Tumor (surrounding tissue) 
                 1.62 
                 1.70 
                 2.70 
               
               
                 Fatty with tumor cells 
                 0.65 
                 0.67 
                 1.08 
               
               
                 Tissue far from tumor center 
                 0.28 
                 0.31 
                 0.37 
               
               
                 Healthy breast tissue 
                 0.22 
                 0.22 
                 0.25 
               
               
                   
               
             
          
         
       
     
         [0062]    Since the electrical conductivity and permittivity of the tissues change with frequency different images can be reconstructed for each operation frequency yielding more than one image for the same body. The time course of the temperature distribution is recorded using thermal infrared cameras. The resultant images are a function of electrical and thermal properties of the tissues and may provide valuable diagnostic information. 
         [0063]    This method utilizes 
         [0064]    the enhancement in the temperature difference generated by the cancerous tissue formation by current injection through surface electrodes placed on the inspected region of the body, 
         [0065]    tissue conductivity imaging by thermography using the dependence of the heat generation on the tissue conductivity, 
         [0066]    the different frequency dependencies of the cancerous (strong dependence) and healthy (weak dependence) tissues in a properly selected frequency interval. 
         [0067]    The technique enhances the temperature difference generated by cancerous tissue by injecting currents within safety limits to the inspected region of the body as well as measuring the conductivity distribution through thermal imaging. Comparison of the conductivity distribution data obtained through current injection at different frequencies yields information on the presence/absence and location of possible cancerous tissues in the body. 
         [0068]    In order to increase the sensitivity of the technique, a thermal camera with dual band (MWIR/LWIR) imaging capability is utilized. For this purpose, the dual band camera includes a single sensor capable of imaging in both bands (MWIR and LWIR). Due to the low cost and the narrow band spectral response nature of the Quantum Well Infrared Photodetector (QWIP) technology, a dual band MWIR/LWIR QWIP seems to be feasible for this approach although detectors based on the HgCdTe, Type II superlattice and quantum dot infrared photodetector (QDIP) technology can also be utilized. The thermal image is first recorded in the MWIR (3-5 (m) band and then in the LWIR (8-12 (m) window in two different wavelength intervals centered around (1 and (2. The ratio of the spectral exitances (M O ) of two blackbodies at these wavelengths can be expressed by the following equation 
         [0000]    
       
         
           
             
               
                 
                   M 
                   λ 
                 
                  
                 
                   ( 
                   
                     λ 
                     1 
                   
                   ) 
                 
               
               
                 
                   M 
                   λ 
                 
                  
                 
                   ( 
                   
                     λ 
                     2 
                   
                   ) 
                 
               
             
             = 
             
               
                 
                   λ 
                   2 
                   5 
                 
                 ( 
                 
                   
                      
                     
                       hc 
                       
                         
                           λ 
                           2 
                         
                          
                         
                           k 
                           B 
                         
                          
                         
                           T 
                           obj 
                         
                       
                     
                   
                   - 
                   1 
                 
                 ) 
               
               
                 
                   λ 
                   1 
                   5 
                 
                 ( 
                 
                   
                      
                     
                       hc 
                       
                         
                           λ 
                           1 
                         
                          
                         
                           k 
                           B 
                         
                          
                         
                           T 
                           obj 
                         
                       
                     
                   
                   - 
                   1 
                 
                 ) 
               
             
           
         
       
     
         [0069]    where T obj  is the temperature to be measured, h is the Planck&#39;s constant, k B  is the Boltzmann constant, and c is the speed of the light. Therefore, the temperature can be measured accurately by a dual band sensor if the emissivity of the object does not significantly change between these wavelength intervals. The other advantage of utilizing dual band imaging is benefiting from the individual advantages of both MWIR and LWIR bands for medical imaging as well as the possibility of thermal imaging performance enhancement by processing the image data obtained in two different thermal imaging bands. In order to facilitate all of the above possibilities, the dual band camera should be able to operate in the following modes selectable by the operator; 
         [0070]    i. single band imaging mode in the MWIR window 
         [0071]    ii. single band imaging mode in the LWIR window 
         [0072]    iii. dual band imaging mode in the MWIR and LWIR windows 
         [0073]    In the third (iii) mode two alternatives are possible. The camera can yield image data in both windows (MWIR and LWIR) simultaneously or sequentially on a frame by frame basis in which case one frame containing the image data in one window will follow the other. 
         [0074]    The infrared sensor can be of mid (such as 320×256 and similar formats) or large (such as 640×512 and similar formats) format while lower formats with degraded resolution or higher formats with increased cost are also possible. 
       AIMS FOR THE DEVELOPMENT OF THE INVENTION 
       [0075]    The aims for the development of the medical electro-thermal imaging method and the system thereof, are 
         [0076]    Procuring an imaging system without using ionizing radiation which is potentially harmful for the human health, 
         [0077]    Developing a method that can detect and define the tumor type in the same examination process, 
         [0078]    Developing a system that can detect the tumors in deeper areas of the tissue, and the micro calcifications in the tissue, 
         [0079]    Obtaining a comfortable examining procedure, which has relatively short diagnosing time, 
         [0080]    The advantages of the invention over the current technology and other thermal imaging methods can be summarized as follows: 
         [0081]    Using three physical properties of the tissue, thermal properties, emissivity property and electrical properties (electrical conductivity and permittivity and possibly magnetic permeability), for more reliable imaging of the examining area over the passive thermography which uses two physical property of the tissue, thermal and emissivity properties, 
         [0082]    Comfortable electric current insertion using surface electrodes instead of embedded electrodes, 
         [0083]    Applying electric current normal to the surface using electrodes to obtain better resolution in the normal direction, 
         [0084]    Attainability of different thermal images of the area regarding to the form and the number and configuration of the electrodes, 
         [0085]    Achieving high resolution images using dual band (MWIR and LWIR) active thermography, 
         [0086]    Easy reconfiguration of the location of the electrodes during examination, 
         [0087]    Obtaining as many thermal images as the inserted electric current frequency variety for obtaining larger amount of data in order to increase the reliability of the assessment, 
         [0088]    The advantages of the invention over other active thermography methods: 
         [0089]    Since the electrical current is applied from the surface of the material (tissue), it can be penetrated vertically in the tissue. By applying electrical currents, higher temperature contrast can be provided at deeper regions. 
         [0090]    In the invention, it is possible to change the number, position and shape of the electrodes easily. When the electrodes are located closer to the tumor, the current density increases in the tumor region, and the temperature contrast on the surface gets higher. 
         [0091]    In passive thermography methods, electrical properties of tissues are not expected to considerably affect the thermal images. On the other hand, active thermography method (with the electrical current application) creates thermal images (temperature distributions) depending on the electrical properties of tissues (electrical conductivity, electrical permittivity and possibly magnetic permeability). Since these features vary with the injected current frequency in the tissue, different temperature distributions can be obtained correspondingly. 
         [0092]    Applying current increases the temperature contrast. 
         [0093]    By using dual-band infrared camera, two separate images in two different wavelength intervals can be obtained for a frequency value of the applied current. These images may have different details according to each other and the individual advantages of single band MWIR and LWIR imaging can be combined in this way. At the same time, skin temperature can be measured much more reliably when compared with single band imaging as explained above. Harold Szu has a patent on dual-band infrared camera usage in medical imaging. However., dual-band infrared thermal imaging with the method proposed in our invention has not been implemented yet. 
       TECHNICAL PROBLEMS SOLVED BY THE INVENTION 
       [0094]    As mentioned above, decrease in the thermal contrast depending on the tumor depth is the fundamental problem of the passive thermography systems. The diagnostic performance of the passive thermography systems is not sufficient especially for smaller tumors located at deeper regions. This problem is solved by the present invention. 
         [0095]    Our numerical analysis show that 1.5 mm tumors, smaller than the critical tumor size (2-3 mm), can only be detected by passive thermography from 3 mm depth. On the other hand, by using the active thermography method, mentioned in the invention, the same tumor can be detected from 1 cm depth from the breast surface. If the tumor size enlarges to 3.5 mm, cancerous tissue at 1.5 cm depth generates more than 20 mK temperature contrast on the breast surface and can be sensed with the state-of-the-art thermal imagers by using active thermography. 
       The Differences of the Invention (from the other known Imaging Modalities) 
       [0096]    Invention, a hybrid modality, is the combination of electrical impedance imaging with thermal infrared imaging which we name Medical Electro-Thermal Imaging. This new functional imaging modality utilizes both techniques and minimizes individual drawbacks of these two imaging systems. Based on the innovations mentioned above, imaging performance is expected to be improved by increasing the sensitivity and specificity of EIT and Infrared Imaging Systems. This attribute has not been existed in any imaging modality. 
         [0097]    The invention is based on the idea that tissue electrical properties, which contain important information for the diagnosis, affect the temperature distribution and this effect can be acquired via contactless measurements. For this purpose, the current is applied to the patient within medical safety limits, and the temperature distribution depending on the electrical properties is measured using an infrared thermal camera. 
     
    
     
       DESCRIPTIONS OF FIGURES 
         [0098]    The names of the drawings presented for better understanding of medical electro-thermal imaging method and the system thereof are listed as follows: 
           [0099]      FIG. 1 . Schematic view of the medical electro-thermal imaging system 
           [0100]      FIG. 2 . Application of the medical electro-thermal imaging method using the medical electro-thermal imaging system. 
           [0101]      FIG. 3 . Flow chart of the medical electro-thermal imaging method. 
           [0102]      FIG. 4 . Electromagnetic model of the tissue exposed to the electric current. 
           [0103]      FIG. 5 . Thermal model of the tissue. 
           [0104]      FIG. 6 . Temperature distribution of the healthy and the tumor tissue (realistic breast phantom) taken in passive mode. 
           [0105]      FIG. 7 . Ratio of thermal images taken in active and passive mode. 
           [0106]      FIG. 8 . Image showing the difference of healthy tissue and tumor in active mode of operation. 
       
    
    
       [0107]    The numerals referred to in the following description correspond to the following, 
         [0108]      1 . Infrared camera, 
         [0109]      2 . Electrodes, 
         [0110]      3 . Cancerous tissue (tumor), 
         [0111]      4 . Breast tissue, 
         [0112]      5 . Control unit, 
         [0113]      6 . Current source, 
         [0114]      7 . Surface of the breast tissue, 
         [0115]      8 . Thermal image. 
       DETAILED DESCRIPTION OF THE INVENTION 
       [0116]    Theory 
         [0117]    We present below the theoretical analysis of thermal contrast enhancement with current injection. Electromagnetic problem of the method is modeled and the schematic of the electromagnetic problem is shown in  FIG. 4 . The electrical model of the body is represented using permeability μ=μ 0 , electrical conductivity κ and permittivity ε. Sinusoidal currents are applied using two electrodes attached on the body surface at points A and B. Applied currents generate an electric field in the conductive body. The steady-state electric field {right arrow over (E)}=−jw{right arrow over (A)}−∇Ø can be calculated using the following coupled partial differential equations, 
         [0000]      ∇ 2   {right arrow over (A)}−jw μ(κ+ jw ε) {right arrow over (A)} −μ(κ+ jw ε)∇Ø=0
 
         [0000]      ∇.[(κ+ jw ε)∇Ø]+∇(κ+ jw ε). jw{right arrow over (A)}= 0
 
         [0118]    and boundary conditions 
         [0000]    
       
         
           
             
               κ 
                
               
                   
               
                
               
                 
                   ∂ 
                   ∅ 
                 
                 
                   ∂ 
                   n 
                 
               
             
             = 
             
               { 
               
                 
                   
                     
                       I 
                        
                       
                           
                       
                        
                       on 
                        
                       
                           
                       
                        
                       A 
                     
                   
                 
                 
                   
                     
                       
                         - 
                         I 
                       
                        
                       
                           
                       
                        
                       on 
                        
                       
                           
                       
                        
                       B 
                     
                   
                 
                 
                   
                     
                       0 
                        
                       
                           
                       
                        
                       otherwise 
                     
                   
                 
               
             
           
         
       
     
         [0119]    where {right arrow over (A)} is the magnetic vector potential, Ø is the scalar potential, and I is the current applied from the surface. 
         [0120]    Thermal problem is also modeled to obtain the temperature distribution inside the tissue. 
         [0121]    Schematic of the Bio-Heat problem (including an external heat source due to current application) is shown in  FIG. 5 . Pennes Bio Heat Equation is used to describe the effects of metabolic generation and blood perfusion over the energy balance. It explains the thermal interaction between tissues and perfused blood in detail: 
         [0000]    
       
         
           
             
               
                 ρ 
                  
                 
                     
                 
                  
                 
                   C 
                   h 
                 
                  
                 
                   
                     ∂ 
                     T 
                   
                   
                     ∂ 
                     t 
                   
                 
               
               + 
               
                 ∇ 
                 
                   · 
                   
                     ( 
                     
                       
                         - 
                         k 
                       
                        
                       
                           
                       
                        
                       
                         ∇ 
                         T 
                       
                     
                     ) 
                   
                 
               
             
             = 
             
               
                 Q 
                 b 
               
               + 
               
                 Q 
                 met 
               
             
           
         
       
     
         [0122]    where, ρ is the density (kg/m 3 ), C h  is the specific heat (J/kgK), T is the absolute temperature (K), k is thermal conductivity (W/mK), Q b  is the heat source due to blood perfusion and Q met  is the metabolic heat generation (W/m 3 ). 
         [0123]    The first term on the right hand side is the source due to blood perfusion which can be expressed as 
         [0000]        Q   b =ρ b   C   hb   W   b ( T   b   −T )
 
         [0124]    where, ρ b  is blood mass density (kg/m 3 ), C hb  is the blood specific heat rate (J/kgK), W b  is the blood perfusion rate (1/s), and T b  is the blood temperature (K) which is approximated to the core temperature of the body, and the temperature of the venous blood is approximated to T which is the unknown temperature value (local tissue temperature (K)). 
         [0125]    Law of conservation of energy states that the heat lost from the skin surface is in a constant equilibrium with the heat supplied by the vascular flow to the skin in the steady state. Thus, heat transfer from the front skin surface (by both convection and radiation to the surrounding air and surfaces at specified temperatures) should be considered as the boundary conditions: 
         [0000]        Q   conv   =h   h   A   s ( T   S   −T   ∞ ) 
         [0000]        Q   rad   =eσA   s ( T   s   4   −T   sur   4 ) 
         [0126]    where, h h  is the convection heat transfer coefficient (W/m 2 .K), A s  is the surface area through which the convection heat transfer takes place, T s  is the surface temperature, T ∞  is the temperature of the air, e is the emissivity of a skin (0.95), σ is the Stefan-Boltzmann constant (W/m 2 . K 4 ), and T sur  is the temperature of the walls, ceiling and floor. In this study, T sur  is assumed to be equal to the air temperature (T sur =T ∞ ). 
         [0127]    Note that, the boundary condition at the front skin surface can also be written as: 
         [0000]      − k∇T=h   h ( T−T   ∞ )
 
         [0128]    Here, h h  (W/m 2 K) represents the overall heat transfer coefficient due to the combined effect of radiation and convection. 
         [0129]    To set the boundary condition at the rear surface of the breast, the temperature of the thoracic wall can be assumed to be the core temperature of the body (i.e., 310 K). 
         [0130]    Due to the applied external current sources, a new term should be added to the right-hand side of the Pennes Bio Heat equation: 
         [0000]    
       
         
           
             
               
                 ρ 
                  
                 
                     
                 
                  
                 C 
                  
                 
                   
                     ∂ 
                     T 
                   
                   
                     ∂ 
                     t 
                   
                 
               
               + 
               
                 ∇ 
                 
                   · 
                   
                     ( 
                     
                       
                         - 
                         k 
                       
                        
                       
                         ∇ 
                         T 
                       
                     
                     ) 
                   
                 
               
             
             = 
             
               
                 Q 
                 b 
               
               + 
               
                 Q 
                 met 
               
               + 
               
                 Q 
                 ext 
               
             
           
         
       
     
         [0131]    The external heat term Q est  is calculated using the following Joule Heat Equation: 
         [0000]    
       
         
           
             
               Q 
               ext 
             
             = 
             
               
                 1 
                 κ 
               
                
               
                 
                    
                   J 
                    
                 
                 2 
               
             
           
         
       
     
         [0132]    where J is the electrical current density and κ is the electrical conductivity of the tissue. 
         [0133]    Medical electro-thermal imaging method comprises the steps of, 
         [0134]    a. Preparation of the patient and the imaging room according to the thermal imaging standards, 
         [0135]    b. Recording the passive mode thermal image(s) of the body in two different windows (MWIR and LWIR) using dual band infrared camera ( 1 ), 
         [0136]    c. Selecting the frequency of electric current that is attached to the breast tissue ( 4 ) using control unit ( 5 ), 
         [0137]    d. Placing the electrodes to the previously specified locations on the surface of the breast tissue ( 4 ), 
         [0138]    e. Applying the electric current for predetermined time period to the targeted breast tissue ( 4 ), 
         [0139]    f. Recording the active thermal image ( 8 ) of the stimulated area (in both bands (MWIR and LWIR)) in the heating period, 
         [0140]    g. Recording the active thermal image ( 8 ) of the stimulated area (in both bands (MWIR and LWIR)) in the cooling period until the temperature of the breast tissue ( 4 ) becomes stable, 
         [0141]    h. Going to the item c and repeating the same procedure for different frequency values of the applied electric current, 
         [0142]    i. Eliminating position differences of the active ( 8 ) and passive thermal images and making necessary matching processes (image registration), 
         [0143]    j. Obtaining the absolute temperature value of the breast and cancerous tissue ( 3  &amp;  4 ) using the dual band property of the thermal camera ( 1 ), 
         [0144]    k. Comparing and determining the ratio/difference of the active ( 8 ) and passive thermal images obtained under different frequencies using the image processing algorithm which is embedded to the control unit to detect the cancerous tissue ( 3 ), 
         [0145]    l. Carrying out steps c through k until a healthy diagnosis or detection is made. 
         [0146]    The proposed hybrid imaging method is denoted in  FIG. 1 .  FIG. 2  shows the application of the system in practice. In this imaging method, electrical current within medical safety limits is applied to the breast tissue from the surface, by means of electrodes ( 2 ) placed on the fingers. The temperature distribution on the breast surface ( 4 ) is recorded in real time with an infrared camera ( 1 ). Due to the current application, the presence of tumor in the breast causes higher temperature contrast on the breast surface. This temperature difference is determined in the thermal image ( 8 ) obtained by the infrared camera ( 1 ).  FIGS. 6, 7, and 8  are thermal images of the realistic breast phantoms which mimic the healthy and cancerous breast tissue. These thermal images were obtained using a QWIP single band (LWIR band) thermal infrared camera. 
         [0147]    The preparation of the room and patient according to thermal imaging standards is an essential part of the medical electro-thermal imaging method for making correct diagnosis about the breast tissue. Imaging should be implemented in a controlled manner. Patients disrobe to the waist and let the surface of the breasts to cool to room temperature (18-22° C.) for 10-15 minutes. 
         [0148]    After the preparation procedure is completed, passive thermal image of the tissue is recorded. In this passive mode, since there are no external sources, the camera displays the surface radiance distribution due to the internal sources (metabolic heat generation and blood fusion). The passive thermal image is used for the comparison with an active thermal image ( 8 ) which is obtained after the medical electro-thermal imaging method is applied. The passive thermal image of the target area is given in the  FIG. 6 . 
         [0149]    After placement of the electrodes, the current is injected to the breast tissue ( 4 ). Since electrodes ( 2 ) are non-invasive type, there is no need for the insertion of the electrodes into the tissue. Electric current is injected from the surface of the breast tissue ( 7 ). The electric current is applied to the tissue ( 4 ) until the temperature of the tissue ( 4 ) become stable or along predetermined period and then active thermal image ( 8 ) is recorded using an infrared camera ( 1 ). 
         [0150]    Dual band sensor usage yields to measure the absolute temperature. Thermal imaging performance is improved by processing the images which are obtained in two different thermal imaging bands (MWIR and LWIR bands). 
         [0151]    Different thermal images can be recorded for different electrode locations. If the places of electrodes ( 2 ) are changed with respect to the location of the tumor tissue ( 3 ) higher temperature contrasts may be achieved by moving the electrodes ( 2 ). 
         [0152]    The parameters that affect the temperature contrast are the tumor dimension, amplitude, duration and frequency of the stimulation current, and location (depth) of the tumor tissue. Increasing the strength of the applied current causes higher contrast in the thermal images. Tumor dimension is also another important parameter affecting the temperature contrast. As the size of the tumor increases, higher temperature contrasts are obtained, and consequently, they can be diagnosed from deeper regions of the breast. An example to the ratio of passive and active ( 8 ) thermal images in other words contrast of the active ( 8 ) and passive thermal images given m the  FIG. 7 , and ratio of the healthy ( 8 ) and the tumor tissue ( 8 ) which are taken in active mode of operation is given in the  FIG. 8 . 
         [0153]    Medical electro-thermal imaging system comprises; 
         [0154]    Infrared camera ( 1 ) which can display temperature distribution in the focused area, 
         [0155]    Electrodes ( 2 ) that drive electric current to the breast tissue ( 4 ), 
         [0156]    Control unit ( 5 ) that regulates the electric current provided by current source ( 6 ), that records image by the data provided by infrared camera ( 1 ) and makes comparison between the passive and active thermal images ( 8 ). 
         [0157]    Current source ( 6 ) which provides the electric current that is inserted to the breast tissue ( 4 ).