Diffuse optical tomography (DOT) is an attractive noninvasive imaging technique that allows three-dimensional (3-D) volumetric reconstructions of optical properties within highly scattering media (e.g., biological tissues) down to a depth of few centimeters. Typically, visible and/or near-infrared light is shined/detected due to proper arrangements of many (e.g., few tens) injection and collection points.1,2 DOT is useful for plenty of biomedical applications from neurology to oncology since optical absorption and scattering are, respectively, linked to chemical composition and microstructure of the investigated tissue.18.104.22.168.–6
Depending on the thickness and geometry of the tissue under investigation, either transmittance (sources and detectors placed on opposite sides of the medium) or reflectance (sources and detectors placed on the same side at a given distance) approaches can be employed. However, the adoption of a reflectance configuration is, in most cases, mandatory owing that many regions of the human body are not accessible using a transmittance approach.
Three different measurement strategies can be employed since sources and detectors can be operated in three different domains: (1) continuous wave (CW), (2) frequency-domain (FD), and (3) time-domain (TD). The CW regime is the simplest one in terms of instrumentation, but the information content brought by one single measurement is poor. Indeed, it does not allow disentangling absorption from scattering contributions while this is, in principle, feasible in both FD and TD techniques.2 However, only the TD approach has the added advantage of encoding different average penetration depths in the photon arrival time.7,8 Therefore, a time-gated analysis9 and/or detection10,11 of the reflectance curve allows one to separate information coming from different layers of the investigated medium. The use of the photon arrival time as an additional information to differently weight spatial regions explored by the photons also leads to ideally better image quality with respect to other techniques,12 as well as to lower dependence on moving artifacts.
An effective analysis of the reflectance TD curve can be achieved using the Mellin–Laplace transform (MLT).13 Indeed, this approach permits, on one hand, one to achieve a variable time-windowing of the diffusion curve to fully take advantage of the information content of TD measurements. On the other hand, due to the windowing of the curves, it reduces the amount of information to be processed, thus reducing calculation steps and processing time.
Unfortunately, the widespread adoption of TD technique for DOT is mainly hampered by both overall cost and dimension of the system, being based on components much more complex than those required for CW instruments.14,15 During the last two decades, picosecond pulsed lasers, time-resolved single-photon detectors, and time-tagging electronics have experienced considerable advancements, but the state-of-the-art TD-DOT systems like MONSTIR16,17 and other instruments18,19 are still large rack-based systems, limiting somewhat a widespread clinical adoption. In particular, since DOT requires a large number of detection points (e.g., few tens), the main bottleneck is still represented by the use of expensive, bulky, and delicate photocathode-based vacuum tubes like photomultiplier tubes (PMTs) or, more recently, hybrid PMTs.20
The unmet need for cost-effective and rugged detectors forced the interest of the scientific community toward microelectronic devices like single-photon avalanche diodes (SPADs).21,22 Yet, their performance in diffuse optics systems is often unsatisfactory due to the small active area (few tens of micrometers in diameter) that affects the collection efficiency of diffused light. To cope with this limitation, SPADs require the use of lasers with high average power and, also, proper solutions for reducing the light power density so as to match the maximum permissible exposure for skin.23
Very recently, silicon photomultipliers (SiPMs) were introduced in TD diffuse optics.24 These promising detectors are inexpensive and rugged, they do not require complex front-end circuits, and they feature wide active areas of few square millimeters, thus permitting to maximize the photon collection. Essentially, they bring together the advantages of vacuum tubes and microelectronics detectors.
After initial positive results in proof-of-concept studies,25,26 SiPM started to be integrated into compact and easy-to-operate detection modules27 and to replace traditional PMTs in systems designed for clinical use.28,29 Still, the state-of-the-art for SiPM technology is not yet fully miniaturized, but due to the possibility to integrate into a single silicon chip in both detectors and ancillary electronic circuits,30,31 SiPMs represent the most attractive technology for the next-generation of DOT systems. Therefore, the demonstration of SiPM suitability also for TD-DOT could open the way to novel immediate perspectives due to the possibility to conceive low-cost and reliable multichannel systems.
Even if SiPMs have already been validated for diffuse optics in general, their suitability also for DOT is not straightforward. Indeed, as highlighted in Ref. 32, the present generation of SiPMs employed as time-resolved single-photon detectors is affected by some limitations. The most critical one is the presence of a slow exponential decaying tail (known as “diffusion tail” being due to the diffusion of photogenerated carriers inside the detector) in the single-photon response shape of the device. This tail is characterized by a few nanoseconds decaying time-constant starting 1 to 2 decades below the response peak.32 As demonstrated in Ref. 33, the tails in the response function can seriously affect the performance of the device in TD diffuse optics in general, possibly preventing the capability to detect and localize absorption perturbations within the diffusive medium. Additionally, this long tail limits the dynamic range of the detector response function to about 2 orders of magnitude,32 thus possibly affecting DOT reconstructions due to the reduced range suitable for data analysis.34 As reported in Ref. 32, SiPMs are also characterized by high noise background (up to few hundred thousand counts per second, depending on the operating conditions), which limits the maximum number of photons per second that can be detected to avoid pile-up effects in TD acquisitions using time-correlated single-photon counting (TCSPC) systems,20 and also, the minimum photon bunch, which can be extracted out of the noise. As a final concern, the temporal stability of the detector response function has to be carefully considered in TD-DOT systems. Indeed, DOT can involve long acquisition times, in particular, when tomographic reconstructions are obtained by scanning the tissue surface instead of using a multichannel system.10
The aim of this paper is to address these issues and demonstrate the suitability of SiPMs for optical tomography in the TD. In particular, we have studied on tissue-mimicking heterogeneous phantoms (using absorbing inclusions at different depths), the sensitivity of SiPMs to deep perturbations into the medium and their capability to properly localize and reconstruct buried inclusions. The paper is organized as follows: in Sec. 2, we describe the experimental setup and the signal analysis technique; Sec. 3 deals with simulations and phantom measurements; finally, in Sec. 4, we summarize the main results of the paper and discuss future perspectives.
Material and Methods
Figure 1 reports the setup used for the experiments. The light source was a mode-locked laser (Fianium Ltd., United Kingdom) providing optical pulses [duration: 26-ps full-width at half maximum (FWHM)] at 820 nm at the repetition rate of 40 MHz. The light passed through a variable optical attenuator and was then injected into the sample through an optical fiber [core diameter: ; numerical aperture (NA): 0.22; Lightech srl, Italy]. The photons re-emitted from the sample were then collected at two different points by means of two optical fibers (core diameter: 1 mm, NA: 0.37; Thorlabs Gmbh, Germany). We decided to use two detection points so as to increase the quality of the reconstructed image by acquiring data from multiple points of view. The distance between injection and collection fibers was set to 30 mm. The light exiting from each collection fiber was then imaged onto a SiPM detector by a set of two lenses (achromatic doublets lens, focal length: 30 mm) providing a magnification.
Each module embeds a active area SiPM (C30742-11-050-T1, Excelitas Technologies, Canada) and all the required electronic circuits for sensing and amplifying the avalanche signal and yielding the pulse synchronous with the avalanche. Additionally, the module hosts a thermoelectric cooler controller board, which is set to keep the detector at about 15°C. This value represents the best compromise to limit both primary dark count rate and afterpulsing (i.e., noise contribution due to charge carriers trapped during the avalanche ignition). For further details about the SiPM module design and characterization, see Ref. 27. Each SiPM module was then coupled to a TCSPC board (SPC-130, Becker&Hickl GmbH, Germany) to provide the start signal for time-stamp generation. The stop pulse was given to both TCSPC boards by the “synchronization” output of the laser.
To acquire the instrument response function (IRF) of the system, we directly faced the source fiber and one detector fiber at a time, putting a thin Teflon layer in-between to ensure the isotropic illumination of the latter. The IRFs obtained for a 10-s measurement time and after the subtraction of the mean value of the background noise for both SiPM modules are reported in Fig. 2. The peak timing jitter (FWHM) is around 260 ps for both detectors (262 ps for module 1 and 255 ps for module 2) while the tail with a long time constant () starts 2.5 decades below the peak, in agreement with Ref. 32.
The homogeneous liquid phantom was a water dilution of black Indian ink (which acts as an absorber) and Intralipid (providing scattering property). The proportion of the components were chosen so as to have an absorption coefficient () of and a reduced scattering one () of at 820 nm.
As an absorbing inhomogeneity within the homogeneous medium, we used five totally absorbing inclusions [black cylinders made of black polyvinyl chloride (PVC)] of different volumes. As reported in previous studies,35,36 small totally absorbing objects properly reproduce realistic absorption inhomogeneity in a reliable and reproducible way and for a wide range of applicability. In particular, for a given background, it is possible to identify a whole equivalence class of absorbing perturbations, with different size and shape, which yield the same temporal perturbation for any choice of geometry (e.g., reflectance, transmittance, source–detector distance), position of the object within the medium, and background absorption. Each equivalence class also contains a totally absorbing object whose volume [equivalent black volume (EBV)] can be taken as representative of the whole class. Some breaking of the equivalence are observed only when black object has a very small volume () and the object is quite close (depth ) either to the source or to the detector. The black inclusions we used in this work were cylinders with height equal to the diameter with values of 3.2, 4, 5, 6.8, 8.6 mm, a volume , 50, 100, 250, , and corresponding to an equivalent , 0.087, 0.15, 0.37 and for a volume, respectively.36 In this paper, the inclusions will be named as “inclusion 1, 2, 3, 4, or 5” depending on their volume (e.g., we will refer to the as “inclusion 3”). For what concerns data presented in the following section, in most cases, there is negligible discrepancy in the equivalence relation: only in two cases (, , with ), some 10% to 20% discrepancy in the amount of perturbation as compared to the equivalent realistic inhomogeneity is observed.36
Since multiple measurements are needed to obtain a tomographic reconstruction during measurements, the injection and collection fibers [represented, respectively, as yellow star and open red circles in Fig. 3(a)] were kept fixed while the inclusion was moved. This strategy was selected in order to avoid possible artifacts due to the movement of fibers, which were instead bound to the phantom tank using a black PVC plate in which holes for fiber tips were drilled. Using a motorized three-axes translation stage, the inclusion was moved over an area of at steps of 7.5 mm in both and directions [scan area: from to 18.8 mm and from to 15 mm for and direction, respectively]. For each inclusion position [small black dots in Fig. 3(a)], we acquired the signal for 15 s (15 repetitions of 1 s each). As reference measurement, for each line, we moved the inclusion far enough () to consider its effect negligible.
For each perturbation, we first aligned the inclusion at 7.5 mm from the injection fiber in both and directions, and this point was set as origin for the scanning coordinates system. This choice was done to obtain a scanning area not centered to the inclusion position, thus being able to disentangle the reconstruction of the inclusion from possible artifacts arising in the center of the image.
Scans of all the inclusions in the plane were done at different depths (where “depth” refers to the distance from the surface of the phantom to the center of the inclusion) ranging from 10 mm up to 20 mm at step of 5 mm. Also, only for “inclusion 3,” 25 and 30 mm depth scans were acquired in order to test the depth sensitivity of the setup. The depth was automatically changed using the motorized translation stage when a scan on the plane was finished.
In the following, the diffusion approximation was used to model light propagation.37 Data analysis was performed considering two sets of measurements and , acquired on a reference medium (homogeneous medium of known absorption ) and on the unknown medium (with unknown absorption , but estimated at iteration by ). As explained in Ref. 13, the reference measurement is needed to take into account the effect of the IRF. In the previous expression of measurements, is the time, is an index on the source position, and is an index on a detector position. As in Ref. 13, we used the combination obtained with these measurements and known Green’s functions of medium and Green’s functions of estimated medium at iteration , which is linearly related to the absorption update . Green’s functions are defined as solutions of the diffusion equation for a Dirac source. In the reported experiment, was measured () and was not considered as an unknown parameter in the data analysis. After space discretization with finite volume method and time discretization with the MLTs,13,38 we ended up with this system:13), is the MLT order (an integer number, here growing from 0 to 20). The space discretization was performed on a regular mesh grid of the medium, partitioned in tetrahedra with 31119 nodes indexed by . The unknown was determined using a conjugate gradient method. Since Eq. (1) (bottom) is a development at the order 1 of , an iterative approach must be undertaken if accurate absorption determination is sought; therefore, we processed our reconstruction by step and updated the absorption  and the predicted Green’s functions () reported in Eq. (1) at each iteration up to 10.
To evaluate the performance of the system, we first carried out a simulation study in the same conditions (same geometry and same reconstruction method) as the described experiments. We generated simulated measurements by time-convoluting computed Green’s functions (for each simulated phantom) and the experimental IRF of the SiPM (within the region of interest highlighted in Fig. 2) and by adding Poisson noise. We decided to consider only Poisson noise since it is the only unavoidable noise contribution.39 Indeed, other sources of noise (e.g., laser power instability, thermal drifts, and so on) are dependent on the particular system being used. It is worth noting that a detailed simulation of the performance of this particular system is out of the scope of this paper. The scope of the presented simulations is to investigate whether or not SiPMs, with their distinctive IRF, can allow tomographic reconstructions. Therefore, simulations should represent the best results allowed by this technology. The simulated phantoms were composed of the same homogenous background medium as the experiment and the same perturbations: spheres with equivalent to those provided by the total absorbing inclusions we used for experiments. Figure 3(b) shows an example of a 3-D reconstruction result for an inclusion set at 20 mm depth obtained in simulation.
To ascertain whether the SiPM performances can alter the quantification capabilities of the system, we have also performed reconstructions on simulations using the IRF of a high-performance SPAD. Indeed, it can be considered as a state-of-the-art, quasi-ideal detector due to its fast response (, diffusion tail time constant 40). Due to the narrow jitter of the SPAD detector, its contribution to the system jitter is expected to be negligible. Indeed, the main contribution to the system jitter is the dispersion introduced by 1-mm core optical fiber, whose value can be derived from the SiPM IRF since all other contributions (e.g., laser jitter, TCSPC card jitter, and so on) are known and can be considered negligible. We therefore considered a FWHM value of the SPAD-based IRF of 258 ps, which is about the average value of the IRF jitter reported in Sec. 2.1.
Parameters of interest
Data analysis of the 3-D reconstructions (coming both from simulations and measurements) was computed to extract spatial and absorption quantification information, as reported in Table 1. To get the values of spatial localization (, , and depth), the weighted centroids of each reconstructed inclusions were calculated on a delimited region defined by all values above a threshold of 50% of the absorption variation , where is calculated by taking the average over a sphere at the expected depth of the inclusion. The expected values of the recovered and position are and for all depths.
Recovered value of x, y and depth localization and value of quantified Δμa computed for both simulations and experiments for all inclusions. “Sim ideal” corresponds to the simulations using the IRF of SPAD.
|x localization (expected: x=0)||y localization (expected: y=0)||Depth localization||Recovered Δμa|
|Inclusion||EBV (mm3)||Δμa for 1 cm3 (cm−1)||Depth (mm)||Sim (mm)||Exp (mm)||Sim (mm)||Exp (mm)||Sim (mm)||Exp (mm)||Sim ideal (cm−1)||Sim (cm−1)||Exp (cm−1)|
Results and Discussion
Figure 4 shows the reconstructed maps (in the and planes, first and second row, respectively) for the “inclusion 3” at depths ranging from 10 mm up to 30 mm (columns). The black circle reported in each map represents the volume of in which the absorption perturbation given by the inclusion is supposed to be. It is clearly noteworthy that in the plane, the position of the absorption perturbation is properly localized in depth till 25 mm, while for larger depths, it is reconstructed in a position shallower than the real one. On the other hand, in the plane, there is a proper localization of the inclusion at all depths. For both and slices at almost all depths, the reconstructed volume is larger than the expected and, more precisely, it progressively increases at higher depths. This is due to the fact that the information of deep inclusions is carried by photons with a long time-of-flight, thus providing a lower spatial resolution.41
With the aim to study the sensitivity to the amount of absorption perturbation, Fig. 5 reports the reconstructed maps (at a given depth: 15 mm) obtained using the whole set of totally absorbing objects. Regardless the volume, the position of the inclusion in both planes ( and ) is precisely reconstructed (for values, see Table 1). In addition, we clearly see that as the real of the perturbation increases, the retrieved increases as well (even if the absolute value is lower than expected) while the reconstructed volume does not, as expected. Indeed, since the inclusions depth is the same, the spatial resolution of the system is expected to be the same.
From simulations, we can conclude that SiPMs are suitable for the tomographic reconstruction of even small absorption perturbations () buried down to 20 mm (see Table 1), with good accuracy in terms of depth and lateral resolution. Indeed, for both and localization, the recovered position always differs of from the expected one, whereas for the axis, the relative error in depth localization is always (for inclusion 2 at depth 15 mm, which corresponds to an error of 16%). By contrast, the absolute quantification of the perturbation is not accurate, and its accuracy decreases upon the increase of the inclusion’s volume. However, the simulations done with the quasi-ideal detector reported in Table 1 clearly demonstrate that this lack in accuracy is not dependent on the peculiar features of the SiPM (e.g., slow diffusion tail), since similar results in quantification are recovered also using a state-of-the-art detector.
Due to the encouraging results obtained in the simulations, we have studied the suitability of the SiPM for TD-DOT with experimental measurements.
The effect of the depth (ranging from 10 to 30 mm, columns) for “inclusion 3” in phantom measurements can be seen in Fig. 6, where the (first row) and (second row) maps are reported. Up to 25 mm depth, the spatial localization (depth and ) is comparable with simulations. For the 30-mm depth measurement, the perturbation is probably too low to properly detect the inclusion, thus resulting in an unsatisfactory localization of the perturbation. For all depths, there is a larger number of artifacts than in simulations, but this is most likely due to instabilities of the experimental system (e.g., laser instability, ageing of the phantom) that were not taken into account in simulations. In any case, it is worth noting that the current setup is able to reconstruct an equivalent absorbing perturbation of (for volume), which is buried down to 25 mm in depth. Those values are typical for biomedical applications (e.g., optical mammography or brain imaging) thus making the use of SiPMs a good solution in medical imaging.
To ascertain the sensitivity on the amount of absorption perturbation, Fig. 7 shows the and maps for the five inclusions at the same depth (15 mm). It is worth noting that all inclusions are properly reconstructed in terms of depth and position.
Concerning the quantification of the absorption perturbation for experimental measurements, it is similar to that obtained from simulations for both the absolute value and for the trend, showing a lower value of upon an increase of the depth.
As can be seen in Table 1, the depth localization is quite accurate for all depths with an error always lower than 2.5 mm (up to 20-mm depth) that corresponds to a maximum error of 17%. We can also notice that experiments and simulations show a good agreement with a discrepancy in the position of the inclusion, which is always lower than 2 mm for depths (except for 30 mm).
The retrieved coordinates of the position of the inclusion recovered from the experiments seem to systematically differ by about 1 to 2 mm from the expected position. Maybe this is due to some systematic errors in the alignment of the inclusion in the experimental measurements. On the other hand, the localization presents a slighter error () for both the simulation and the experimental data. To conclude, the error in the localization is small enough not to affect significantly the capability of the system to retrieve the position of an absorbing perturbation in a real application.
The “inclusion 3”, set at a depth of 30 mm, produced a too low perturbation to be able to properly localize it not only in depth but also in the and directions. For this reason, we can state that the maximum depth of the inclusion allowing proper tomographic reconstruction is 25 mm.
Concerning the quantification capability on the retrieved , measurements and simulations are usually in good agreement (see Table 1). The only exception is for inclusion 1 and 2 at a shallower depth, which could be partially ascribed to the deviation of the equivalence relations between EBV and the realistic , which, for these particular cases, implies higher contrast for the black object in the experiments. On the other hand, there is a clear underestimation in the reconstructed upon increasing the depth. Since this trend is present also in simulations, we suppose that it is not possible to ascribe it to an experimental problem, but more to an issue connected with the model used for reconstruction. Indeed, as discussed in Sec. 3.1, simulations using an almost ideal detector (without long tail) are substantially in agreement with SiPM simulations. Thus, we can conclude that SiPM performances do not alter quantification capabilities as compared to a state-of-the-art detector.
The absolute quantification of the perturbation is a well-known problem in tomography42 in particular for deep inhomogeneities. More comprehensive studies are needed to quantify this effect properly, identify the determinants, and try to improve quantitation. Surely, the decrease in retrieved for increasing depth is a concern when imaging complex structures since fluctuations in reconstruction due to shallower inhomogeneities can hinder deeper and attenuated structures. From the data in Table 1, this effect—for the background properties presented in this study (, )—is around a factor of 3 to 4 between a depth of 10 and 20 mm. Yet, this effect is compensated by the broadening of the reconstructed inhomogeneity, which implies that the overall impact of the deeper perturbation on the image is somehow augmented. In a different study, with conventional detectors, we have shown that fixing the size of the reconstructed inhomogeneity to the true value, as could be obtained through a coregistration with another imaging modality, completely solves the depth reduction in reconstructed .43
In this work, SiPMs are proposed as potentially revolutionary detectors for TD-DOT due to their unique features, bringing together the advantages of PMTs and SPADs. In particular, our aim was to validate their performance on heterogeneous phantoms to evaluate the impact of their inherent disadvantages: a slow tail and low dynamic range in the temporal response to single photons, a high background noise, and unknown stability for long measurement times in TCSPC applications. We performed a study based both on simulations (taking into account the true detector temporal response) and experiments in order to disentangle nonidealities of the employed system from the actual detector behavior.
We demonstrated from both points of view that SiPMs clearly allow one to detect perturbations of in a volume of down to 3 cm depth into a scattering medium featuring and using a 3-cm interfiber distance. An excellent precision in lateral localization (better than 2 mm) is obtained down to a 2.5-cm depth, whereas depth localization is good up to 2 cm (better than 3 mm). Additionally, it is worth pointing out that with a perturbation placed at a depth of 2 cm, a lateral localization better than 1.5 mm and a depth localization better than 1.5 mm can be achieved with a large range of perturbation values, ranging from to . By contrast, the accuracy in quantification of the absolute value of the perturbation is much worst, with a strong underestimation of the retrieved upon increasing depth. Still, this effect is not ascribed to the SiPMs since it is replicated in simulations also for a state-of-the-art, almost ideal detector. The problem of the absolute quantification in diffuse optics measurements is well known in the literature42 and a dedicated effort is needed to properly address this issue, which is out of the scope of this work. Preliminary results show that this issue is due to the physics of diffuse optics and the retrieval algorithm and not to the specific experimental implementation.43
The demonstration of SiPM suitability for TD-DOT opens now the way to interesting perspectives. In this work, we made use of two custom-made detection modules based on SiPMs. As reported in Ref. 27, the overall dimensions of each module is while the cost in terms of components is around 1000 €. This solution allowed us to get rid of thermal drift and electromagnetic interferences and to maintain high timing resolution due to the high-performance radiofrequency electronics that can be hosted inside such detection modules.27 However, with the aim to conceive multichannel systems, most of the electronics presently hosted into the module can be shared between different channels (e.g., power supply generators and high-voltage circuitry for SiPM biasing), thus further lowering sizes and costs even with the state-of-the-art technology. More farseeing, fully integrated devices embedding the detector, its ancillary electronics, and also the TCSPC circuitry represent an incredible possibility to fabricate detection chains with millimetric size, with the strong potential to revolutionize the entire field.
The research leading to these results has partially received funding from the European Union’s Horizon 2020 research and innovation program under Grant Agreement No. 731877 SOLUS: Smart Optical and Ultrasound diagnostics of breast cancer and under Grant Agreement No. 654148 Laserlab-Europe.
Laura Di Sieno received her master’s degree in electronics engineering at the Politecnico di Milano, Italy, in 2011, and her PhD in physics at the same university in 2015. Since January 2015, she has been a postdoctoral fellow, and her activity is mainly focused on the study and application of instrumentation for time-resolved optical spectroscopy of highly scattering media using SPAD enabled in fast-gated mode.
Judy Zouaoui received her Engineer’s degree in biomedical engineering in 2013 from the Engineering School of Polytech, Lyon, France. The same year, she obtained her master’s degree in medical imaging from the University of Lyon 1 in France. Since November 2013, she has been a PhD student at CEA LETI and at the University Grenoble Alpes on the topic of multispectral time-resolved diffuse optical tomography.
Lionel Hervé has mainly worked on signal processing for medical imaging. He studied the relationship between mammographic densities and cancer occurrence on large population trials in a two-year postdoc fellowship in UCSF. Since 2005, he worked on the diffuse optical tomography modality which consists of reconstructing 3-D optical property maps from IR optical measurements. He particularly focused on the processing of time-domain optical measurements to probe deeper inside biological media.
Antonio Pifferi graduated in nuclear engineering in 1991 and received his PhD in physics in 1995 at the Politecnico di Milano. Since 2003, he has been an associate professor in the Physics Department of the Polytechnic of Milan. His research focuses on the development of new laser and instrumentation exploiting light propagation in highly scattering media with applications to optical characterization of biological tissues, optical mammography.
Andrea Farina received his master’s degree (cum laude) in electronics engineering from Politecnico di Milano, Italy, in 2005, and his PhD degree in physics at the same institution in 2009. Since December 2011, he has held a permanent position as a researcher at the Italian National Research Council (CNR) in Milan. His current research includes the development and the application of time-resolved techniques for the imaging and spectral characterization of diffusive media.
Edoardo Martinenghi received his master’s degree in electronics engineering from the Politecnico di Milano, Italy, in 2013, where he is currently working toward his PhD. Currently, his research activity concerns the design of TD diffuse spectroscopy instrumentations and their applications to biomedical diagnosis.
Jacques Derouard is a professor emeritus of physics at the University Grenoble Alpes. His current research interests are related to optics and spectroscopy methodologies applied to life sciences, more specifically fluorescence correlation spectroscopy and optics in turbid media.
Jean-Marc Dinten is a senior scientist at the Biology and Health Division in CEA-LETI. For more than 20 years, he has been widely developing medical image processing and reconstruction algorithms associated to the development of innovative x-rays and optical imaging systems. He now heads the Imaging Readout Systems Laboratory, which develops new optical imaging systems for health and biology applications.
Alberto Dalla Mora has been an assistant professor in the Department of Physics, Politecnico di Milano, Italy, since 2011. He graduated summa cum laude in electronics engineering from Politecnico di Milano in 2006, and he received his PhD degree summa cum laude in information and communication technology at the same university in 2010. He authored more than 35 papers in international peer-reviewed journals. Currently, his research interests include time-resolved diffuse spectroscopy instrumentation and applications for biomedical diagnosis.