Diffuse Optical Tomography Imaging Research Articles

Image quality improvements of diffuse optical tomography by using multiple polarization components

The inverse process of diffuse optical tomography (DOT) is ill-posed, which has a significant effect on the quality of the reconstructed results. In this paper, the polarization characteristics of light are introduced into the inverse solution model of DOT to increase the solution constraints and improve the imaging quality. The spectral element method is used to numerically solve the vector radiative transfer equation, which is used to simulate the propagation of polarized light in a scattering medium. The DOT reverse solution algorithm based on the polarization component is constructed so that the Stokes vector can be used for absorption coefficient and scattering coefficient reconstruction. The simulation and experimental results show that when polarized light is used for incidence, using the intensity component and the multi-polarization component to solve the DOT inverse model can improve the reconstruction accuracy of the medium absorption coefficient distribution and the DOT image quality.

Difüz Optik Tomografi’de Modellenmiş Meme Tümörü Benzeri Dokuların Teşhisi için Simülasyon Deneyleriyle Geri Çatım Algoritmalarının Uygulanması

In Diffuse Optical Tomography (DOT), data processing and reconstruction stages are crucial to obtain high-quality images. Thus, choosing suitable algorithms for the system is a critical choice. This study aims to determine an appropriate reconstruction algorithm for DOT imaging. There are several reconstruction algorithms used in DOT systems. Some algorithms have been improved for solving specific cases, and some still need to be improved. In this study, we used three algorithms for the reconstruction process: Singular Value Decomposition (SVD), Bi-Conjugated Gradient (Bi-CG), and Transpose Free Quasi Minimal Residual (TFQMR). In testing the algorithms, data of the simulation experiments have been used. The simulation experiments model the tumoral tissue within the breast. All three algorithms were produced correct images while the tumor close to the surface. In the case of the tumor that is not close to the breast surface, the tumor location on the images created by Bi-CG and SVD algorithms was not its actual location. However, the tumor location in the image created by the TFQMR algorithm was close to its actual location. Outcomes of the reconstruction algorithms were evaluated based on correctly defining the location of the tumors by using Mean Percentage Error (MPE), Mean Squared Error (MSE), and Mean Absolute Error (MAE) metrics. We have demonstrated the TFQMR algorithm is a more appropriate reconstruction technique for DOT systems. Thus, we have concluded that TFQMR can have the potential to be used in medical imaging systems.

Investigation of effect of modulation frequency on high-density diffuse optical tomography image quality.

.Significance: By incorporating multiple overlapping functional near-infrared spectroscopy (fNIRS) measurements, high-density diffuse optical tomography (HD-DOT) images human brain function with fidelity comparable to functional magnetic resonance imaging (fMRI). Previous work has shown that frequency domain high-density diffuse optical tomography (FD-HD-DOT) may further improve image quality over more traditional continuous wave (CW) HD-DOT.Aim: The effects of modulation frequency on image quality as obtainable with FD-HD-DOT is investigated through simulations with a realistic noise model of functional activations in human head models, arising from 11 source modulation frequencies between CW and 1000 MHz.Approach: Simulations were performed using five representative head models with an HD regular grid of 158 light sources and 166 detectors and an empirically derived noise model. Functional reconstructions were quantitatively assessed with multiple image quality metrics including the localization error (LE), success rate, full width at half maximum, and full volume at half maximum (FVHM). All metrics were evaluated against CW-based models.Results: Compared to CW, localization accuracy is improved by >40% throughout brain depths of 13 to 25 mm below the surface with 300 to 500 MHz modulation frequencies. Additionally, the reliable field of view in brain tissue is enlarged by 35% to 48% within an optimal frequency of 300 MHz after considering realistic noise, depending on the dynamic range of the system.Conclusions: These results point to the tremendous opportunities in further development of high bandwidth FD-HD-DOT system hardware for applications in human brain mapping.

Frequency shifting model for diffuse optical tomography

The Diffuse Optical Tomography (DOT) imaging modality is an interesting research field for researchers since it has uncertainties in the solution space. DOT imaging modality is an unsolved scientific problem. Inverse problem solution and image reconstruction have never been of their best quality. Reconstructed images have low spatial resolution. The scattering nature of low-energized diffusive light is the obscuring effect for DOT modality. DOT has 3 functional sub-branches which are Continuous Wave, Time-Resolved, and Frequency-Domain. In this work, one new approach to Frequency Domain Diffuse Optical Tomography (FDDOT) biomedical optic imaging modality is presented to the readers. Frequency shifting data were added to the forward model problem which has source-detector couplings and several imaging voxels. 100 MHz center core light modulation frequency was selected as modulation frequency. 169 source-detector matches were used on back-reflected imaging geometry. Absorption coefficient µa was selected 0.1 cm−1. Scattering coefficient µs was selected 100 cm−1. 1-µm x, y, z cartesian grid coordinates were used in each direction for imaging tissue-like simulation media. A total of 100 frequency shifts were added to the forward model problem which has 5 Hz frequency step. 2 inclusion objects were embedded inside the imaging simulation phantom to be reconstructed. Their optical absorption coefficient µa = 0.7 cm−1. Two inclusion images were successfully reconstructed with the low contrast to noise ratio (CNR) and position error. The frequency shifting technique is first applied for FDDOT here in this work. This technique increased the total number of equations in the forward model problem; hence it is helping to solve the inverse problem more accurately. In this work, the positive effect of using frequency shifting methodology was observed. Differentiation of 2 embedded inclusions was completed and illustrated in this work. Two different scenarios were compared with each other. In the first scenario, only 100 MHz center core frequency was used for one source-detector match. In the second scenario, 100 additional shifted frequencies over 100 MHz center core frequency were added. It was seen that the second scenario which has more frequency shifts is superior to the general concept FDDOT methodology. This work shows that the frequency shift technique might be used for future FDDOT devices.

Ultrasound segmentation-guided edge artifact reduction in diffuse optical tomography using connected component analysis.

Ultrasound (US)-guided diffuse optical tomography (DOT) has demonstrated potential value for breast cancer diagnosis and treatment response assessment. However, in clinical use, the chest wall, poor probe-tissue contact, and tissue heterogeneity can all cause image artifacts. These image artifacts, appearing commonly as hot spots in the non-lesion regions (edge artifacts), can decrease the reconstruction accuracy and cause misinterpretation of lesion images. Here we introduce an iterative, connected component analysis-based image artifact reduction algorithm. A convolutional neural network (CNN) is used to segment co-registered US images to extract the lesion location and size to guide the artifact reduction. We demonstrate its performance using Monte Carlo simulations on VICTRE digital breast phantoms and breast patient images. In simulated tissue mismatch models, this algorithm successfully reduces edge artifacts without significantly changing the reconstructed target absorption coefficients. With clinical data it improves the optical contrast between malignant and benign groups, from 1.55 without artifact reduction to 1.91 with artifact reduction. The proposed algorithm has a broad range of applications in other modality-guided DOT imaging.

Three-compartment-breast (3CB) prior-guided diffuse optical tomography based on dual-energy digital breast tomosynthesis (DBT).

Diffuse optical tomography (DOT) is a non-invasive functional imaging modality that uses near-infrared (NIR) light to measure the oxygenation state and the concentration of hemoglobin. By complementarily using DOT with other anatomical imaging modalities, physicians can diagnose more accurately through additional functional image information. In breast imaging, diagnosis of dense breasts is often challenging because the bulky fibrous tissues may hinder the correct tumor characterization. In this work, we proposed a three-compartment-breast (3CB) decomposition-based prior-guided optical tomography for enhancing DOT image quality. We conjectured that the 3CB prior would lead to improvement of the spatial resolution and also of the contrast of the reconstructed tumor image, particularly for the dense breasts. We conducted a Monte-Carlo simulation to acquire dual-energy X-ray projections of a realistic 3D numerical breast phantom and performed digital breast tomosynthesis (DBT) for setting up a 3CB model. The 3CB prior was then used as a structural guide in DOT image reconstruction. The proposed method resulted in the higher spatial resolution of the recovered tumor even when the tumor is surrounded by the fibroglandular tissues compared with the typical two-composition-prior method or the standard Tikhonov regularization method.

Age-related changes in diffuse optical tomography sensitivity profiles in infancy.

Diffuse optical tomography uses near-infrared light spectroscopy to measure changes in cerebral hemoglobin concentration. Anatomical interpretations of the location that generates the hemodynamic signal requires accurate descriptions of diffuse optical tomography sensitivity to the underlying cortical structures. Such information is limited for pediatric populations because they undergo rapid head and brain development. The present study used photon propagation simulation methods to examine diffuse optical tomography sensitivity profiles in realistic head models among infants ranging from 2 weeks to 24 months with narrow age bins, children (4 and 12 years) and adults (20 to 24 years). The sensitivity profiles changed systematically with the source-detector separation distance. The peak of the sensitivity function in the head was largest at the smallest separation distance and decreased as separation distance increased. The fluence value dissipated more quickly with sampling depth at the shorter source-detector separations than the longer separation distances. There were age-related differences in the shape and variance of sensitivity profiles across a wide range of source-detector separation distances. Our findings have important implications in the design of sensor placement and diffuse optical tomography image reconstruction in (functional) near-infrared light spectroscopy research. Age-appropriate realistic head models should be used to provide anatomical guidance for standalone near-infrared light spectroscopy data in infants.

Neuroimaging of depression with diffuse optical tomography during repetitive transcranial magnetic stimulation

Repetitive transcranial magnetic stimulation (rTMS) is an effective and safe treatment for depression; however, its potential has likely been hindered due to non-optimized targeting, unclear ideal stimulation parameters, and lack of information regarding how the brain is physiologically responding during and after stimulation. While neuroimaging is ideal for obtaining such critical information, existing modalities have been limited due to poor resolutions, along with significant noise interference from the electromagnetic spectrum. In this study, we used a novel diffuse optical tomography (DOT) device in order to advance our understanding of the neurophysiological effects of rTMS in depression. Healthy and depressed subjects aged 18–70 were recruited. Treatment parameters were standardized with targeting of the left dorsolateral prefrontal cortex with a magnetic field intensity of 100% of motor threshold, pulse frequency of 10 per second, a 4 s stimulation time and a 26 s rest time. DOT imaging was simultaneously acquired from the contralateral dorsolateral prefrontal cortex. Six healthy and seven depressed subjects were included for final analysis. Hemoglobin changes and volumetric three-dimensional activation patterns were successfully captured. Depressed subjects were observed to have a delayed and less robust response to rTMS with a decreased volume of activation compared to healthy subjects. In this first-in-human study, we demonstrated the ability of DOT to safely and reliably capture and compare cortical response patterns to rTMS in depressed and healthy subjects. We introduced this emerging optical functional imaging modality as a novel approach to investigating targeting, new treatment parameters, and physiological effects of rTMS in depression.

Changes in Diffuse Optical Tomography Images During Early Stages of Neoadjuvant Chemotherapy Correlate with Tumor Response in Different Breast Cancer Subtypes.

This study's primary objective was to evaluate the changes in optically derived parameters acquired with a diffuse optical tomography breast imaging system (DOTBIS) in the tumor volume of patients with breast carcinoma receiving neoadjuvant chemotherapy (NAC). In this analysis of 105 patients with stage II-III breast cancer, normalized mean values of total hemoglobin ([Formula: see text]), oxyhemoglobin ([Formula: see text]), deoxy-hemoglobin concentration ([Formula: see text]), water, and oxygen saturation ([Formula: see text]) percentages were collected at different timepoints during NAC and compared with baseline measurements. This report compared changes in these optical biomarkers measured in patients who did not achieve a pathologic complete response (non-pCR) and those with a pCR. Differences regarding molecular subtypes were included for hormone receptor-positive and HER2-negative, HER2-positive, and triple-negative breast cancer. At baseline, [Formula: see text] was higher for pCR tumors (3.97 ± 2.29) compared with non-pCR tumors (3.00 ± 1.72; P = 0.031). At the earliest imaging point after starting therapy, the mean change of [Formula: see text] compared with baseline ([Formula: see text]) was statistically significantly higher in non-pCR (1.23 ± 0.67) than in those with a pCR (0.87 ± 0.61; P < 0.0005), and significantly correlated to residual cancer burden classification (r = 0.448; P < 0.0005). [Formula: see text] combined with HER2 status was proposed as a two-predictor logistic model, with AUC = 0.891; P < 0.0005; and 95% confidence interval, 0.812-0.969. This study demonstrates that DOTBIS measured features change over time according to tumor pCR status and may predict early in the NAC treatment course whether a patient is responding to NAC.

Effect and correction of optode coupling errors in breast imaging using diffuse optical tomography.

In diffuse optical tomography (DOT) and spectroscopy (DOS) using handheld probes, tissue curvature can cause bad fiber-to-tissue contact. Understanding and minimizing image artifacts caused by these coupling errors would significantly improve DOT and DOS image quality. In this work, we utilized Monte Carlo simulations and experiments with gelatin-Intralipid phantoms to systematically study the influence of source or detector (optode) coupling errors. Optode coupling errors can increase the amplitude and decrease the phase of the measured diffuse reflectance, creating artifacts in the reconstructed absorption maps, such as hot spots on the edges. We propose an outlier removal algorithm that can correct these image artifacts, and we demonstrate its performance using simulations, phantom experiments, and breast patient data acquired with bad probe contact due to a dense or small breast. Further, we designed and implemented a new resistance-type thin-film force sensor array that provides real-time optode coupling feedback and guides the outlier removal to minimize optode coupling errors. Our approaches and study results have significant implications for reducing image artifacts arising from handheld probes, which are commonly used with mobile and wearable DOT and DOS devices.

Cap-Based Transcranial Optical Tomography in an Awake Infant.

Although Blood Oxygenation Level Dependent (BOLD) functional MRI (fMRI) is widely used to examine brain function in adults, the need for general anesthesia limits its practical utility in infants and small children. Functional Near-Infrared Spectroscopy - Diffuse Optical Tomography (fNIRS-DOT) imaging promises to be an alternative brain network imaging technique. Yet current versions of continuous-wave fNIRS-DOT systems are restricted to the cortical surface measurements and do not probe deep structures that are frequently injured especially in premature infants. Herein we report a transcranial near infrared optical imaging system, called Cap-based Transcranial Optical Tomography (CTOT) able to image whole brain hemodynamic activity with 3 seconds of data acquisition time. We show the system is capable of whole brain oxygenation mapping in an awake child, and that tomographically reconstructed static CTOT-derived oxy- and deoxygenated blood volumes are spatially correlated with the time-averaged BOLD fMRI volumes. By removing time bottlenecks in the current system, dynamic CTOT mapping should be possible, which would then enable evaluation of functional connectivity in awake infants.

Development of digital breast tomosynthesis and diffuse optical tomography fusion imaging for breast cancer detection

Diffuse optical tomography (DOT) non-invasively measures the functional characteristics of breast lesions using near infrared light to probe tissue optical properties. This study aimed to evaluate a new digital breast tomosynthesis (DBT)/DOT fusion imaging technique and obtain preliminary data for breast cancer detection. Twenty-eight women were prospectively enrolled and underwent both DBT and DOT examinations. DBT/DOT fusion imaging was created after acquisition of both examinations. Two breast radiologists analyzed DBT and DOT images independently, and then finally evaluated the fusion images. The diagnostic performance of each reading session was compared and interobserver agreement was assessed. The technical success rate was 96.4%, with one failure due to an error during DOT data storage. Among the 27 women finally included in the analysis, 13 had breast cancer. The areas under the receiver operating characteristic curve (AUCs) for DBT were 0.783 and 0.854 for readers 1 and 2, respectively. DOT showed comparable diagnostic performance to DBT for both readers. The AUCs were significantly improved (P = 0.004) when the DBT/DOT fusion images were used. Interobserver agreements were highest for the DBT/DOT fusion images. In conclusion, this study suggests that DBT/DOT fusion imaging technique appears to be a promising tool for breast cancer diagnosis.

Improving vascular imaging with co-planar mutually guided photoacoustic and diffuse optical tomography: a simulation study

Diffuse optical tomography (DOT) and photoacoustic tomography (PAT) are functional imaging modalities that provide absorption coefficient maps of the tissue. Spatial resolution of DOT is relatively low due to light scattering characteristics of the tissue. On the other hand, although PAT can resolve regions of different absorptions with a high spatial resolution, measuring the absolute value of optical absorptions using PAT is challenging due to unknown light fluence distribution in the tissue. Development of image guidance techniques using a priori information of imaging target structure has been shown to increase the accuracy of DOT. PAT is one such method that can be used as a complementary modality to serve as a guide for DOT image reconstruction. On the other hand, estimated fluence map provided by DOT can be used to quantitatively correct PAT images. In this study we introduce a mutually-guided imaging system for fast and simultaneous optical and photoacoustic measurements of tissue absorption map, where DOT is guided by the PAT image and vice versa. Using the obtained absorption map of the tissue, we then estimate the tissue scattering map. We conducted this study using a series of simulations on digital phantoms and demonstrated the effectiveness of the proposed method.

A wide field-of-view, modular, high-density diffuse optical tomography system for minimally constrained three-dimensional functional neuroimaging.

The ability to produce high-quality images of human brain function in any environment and during unconstrained movement of the subject has long been a goal of neuroimaging research. Diffuse optical tomography, which uses the intensity of back-scattered near-infrared light from multiple source-detector pairs to image changes in haemoglobin concentrations in the brain, is uniquely placed to achieve this goal. Here, we describe a new generation of modular, fibre-less, high-density diffuse optical tomography technology that provides exceptional sensitivity, a large dynamic range, a field-of-view sufficient to cover approximately one-third of the adult scalp, and also incorporates dedicated motion sensing into each module. Using in-vivo measures, we demonstrate a noise-equivalent power of 318 fW, and an effective dynamic range of 142 dB. We describe the application of this system to a novel somatomotor neuroimaging paradigm that involves subjects walking and texting on a smartphone. Our results demonstrate that wearable high-density diffuse optical tomography permits three-dimensional imaging of the human brain function during overt movement of the subject; images of somatomotor cortical activation can be obtained while subjects move in a relatively unconstrained manner, and these images are in good agreement with those obtained while the subjects remain stationary. The scalable nature of the technology we described here paves the way for the routine acquisition of high-quality, three-dimensional, whole-cortex diffuse optical tomography images of cerebral haemodynamics, both inside and outside of the laboratory environment, which has profound implications for neuroscience.

Optimal breast cancer diagnostic strategy using combined ultrasound and diffuse optical tomography.

Ultrasound (US)-guided near-infrared diffuse optical tomography (DOT) has demonstrated great potential as an adjunct breast cancer diagnosis tool to US imaging alone, especially in reducing unnecessary benign biopsies. However, DOT data processing and image reconstruction speeds remain slow compared to the real-time speed of US. Real-time or near real-time diagnosis with DOT is an important step toward the clinical translation of US-guided DOT. Here, to address this important need, we present a two-stage diagnostic strategy that is both computationally efficient and accurate. In the first stage, benign lesions are identified in near real-time by use of a random forest classifier acting on the DOT measurements and the radiologists' US diagnostic scores. Any lesions that cannot be reliably classified by the random forest classifier will be passed on to the second stage which begins with image reconstruction. Functional information from the reconstructed hemoglobin concentrations is employed by a Support Vector Machine (SVM) classifier for diagnosis at the end of the second stage. This two-step classification approach which combines both perturbation data and functional features, results in improved classification, as denoted by the receiver operating characteristic (ROC) curve. Using this two-step approach, the area under the ROC curve (AUC) is 0.937 ± 0.009, with a sensitivity of 91.4% and specificity of 85.7%. In comparison, using functional features and US score yields an AUC of 0.892 ± 0.027, with a sensitivity of 90.2% and specificity of 74.5%. Most notably, the specificity is increased by more than 10% due to the implementation of the random forest classifier.

Small-animal 360-deg fluorescence diffuse optical tomography using structural prior information from ultrasound imaging

.We demonstrate dual modality of free-space fluorescence diffuse optical tomography (FDOT) and handheld ultrasound (US) imaging to reveal both functional and structural information in small animals. FDOT is a noninvasive method for examining the fluorophore inside an object from the light distribution of the surface. In FDOT, a 660-nm continuous wave diode laser was used as an excitation source and an electron-multiplying charge-coupled device (EMCCD) was used for fluorescence data acquisition. Both the laser and EMCCD were mounted on a 360-deg rotation gantry for the transmission optical data collection. The structural information is obtained from a 6- to 17-MHz handheld US linear transducer by single-side access and conducts in the reconstruction as soft priors. The rotation ranges from 0 deg to 360 deg; different rotation degrees, object positions, and parameters were determined for comparison. Both phantom and tissue phantom results demonstrate that fluorophore distribution can be recovered accurately and quantitatively using this imaging system. Finally, an animal study confirms that the system can extract a dual-modality image, validating its feasibility for further in vivo experiments. In all experiments, the error and standard deviation decrease as the rotation degree is increased and the error was reduced to 10% when the rotation degree was increased over 135 deg.

Imaging System Based on Silicon Photomultipliers and Light Emitting Diodes for Functional Near-Infrared Spectroscopy

We built a fiber-less prototype of an optical system with 156 channels each one consisting of an optode made of a silicon photomultiplier (SiPM) and a pair of light emitting diodes (LEDs) operating at 700 nm and 830 nm. The system uses functional near-infrared spectroscopy (fNIRS) and diffuse optical tomography (DOT) imaging of the cortical activity of the human brain at frequencies above 1 Hz. In this paper, we discuss testing and system optimization performed through measurements on a multi-layered optical phantom with mechanically movable parts that simulate near-infrared light scattering inhomogeneities. The baseline optical characteristics of the phantom are carefully characterized and compared to those of human tissues. Here we discuss several technical aspects of the system development, such as LED light output drift and its possible compensation, SiPM linearity, corrections of channel signal differences, and signal-to-noise ratio (SNR). We implement an imaging algorithm that investigates large phantom regions. Thanks to the use of SiPMs, very large source-to-detector distances are acquired with a high SNR and 2 Hz time resolution. The overall results demonstrate the high potentialities of a system based on SiPMs for fNIRS/DOT human brain imaging applications.

Application of deep neural networks to improve diagnostic accuracy of rheumatoid arthritis using diffuse optical tomography

A set of deep neural network models for rheumatoid arthritis (RA) classification using a highway network, a convolutional neural network and a residual network is proposed based on the data of diffuse optical tomography (DOT) utilising near-infrared light, which ensures early diagnosis of pathophysiological changes resulting from inflammation. A numerical model of the finger is used to generate images to overcome the inherent problem of insufficient clinical DOT images available. The proposed deep neural network models are applied to automatically classify simulated DOT images of inflamed and non-inflamed joints and transfer learning is also used to improve the performance of the classification. The results demonstrate that all three deep neural network methods improve the diagnostic accuracy as compared to the widely applied support vector machine (SVM), especially for high inter-subject variability databases. In cases of distinct modelled severity of disease, residual network achieved the highest accuracy (> 99 %), and both of highway and convolutional neural networks reached 99 %, respectively. However, as the severity of the modelled disease is reduced, this accuracy is reduced to 75.2 % for residual networks. The results indicate that transfer learning can improve the performance of deep neural network methods on RA classification from DOT data and highlight their potential as a computer aided tool in DOT diagnostic systems.

Utilisation of spatial derivative measurements in circular diffuse optical tomographic imaging to improve image resolution and contrast

The resolution of images recovered using diffuse optical tomography without spatial prior information is inherently limited due to the diffusive nature of light transport in scattering dominated biological tissue. Several studies have previously reported an improvement in depth sensitivity in near infrared spectroscopy by considering the subtraction of two neighbouring measurements. A nearest neighbour, spatial derivative based algorithm is presented, in which the difference between neighbouring source-sharing source – detector pairs is considered for tomographic reconstruction. Using a circular, inward-facing optode geometry, it is shown that deeper abnormalities in oxyhaemoglobin, deoxyhaemoglobin and water can be recovered with better localisation and contrast as compared to reconstructing using absolute values of multi-wavelength data. Relative improvements of up to 25 % in quantitative accuracy of recovered physiological values are demonstrated at better resolution for cancerous tissue in a heterogeneous breast model, when using this proposed spatial derivative based parameter reconstruction.

Mathematical and numerical challenges in optical screening of female breast.

Diffuse optical tomography (DOT) is an emerging imaging technique which uses light for diagnostic purposes in a non-invasive and non-ionizing way. In this paper, we focus on DOT application to female breast screening, where the surface of the breast is illuminated by light sources and the outgoing light is collected on the surface. The comparison of measured light data with the equivalent field obtained from a relevant mathematical model yields the DOT inverse problem whose solution provides an estimate of the optical coefficients of the tissue. These latter, in turn, can be related to clinical markers for cancer detection. The goal of this work is to propose a mathematical and computational approach tailored to the concept of a DOT imaging device able to perform fast and accurate screenings at an affordable cost. Namely, we address two original points about the crucial issue of the solution of the severely ill-conditioned DOT inverse problem: (a) a computational approach based on Green's functions which do not require the exact knowledge of the tissue geometry, proposed here in the declination of the Method of Fundamental Solutions, which allows to enforce correct boundary conditions; (b) the elastic net regularization technique that shares the desirable properties of both the ℓ2 - and ℓ1 -norm penalization approaches and opens the possibility for sparsity recognition in the optical coefficients field and refinement procedures.