Multi-Wavelength Photo-Magnetic Imaging System for Photothermal Therapy Guidance.

BACKGROUND AND OBJECTIVES: In photothermal therapy, cancerous tissue is treated by the heat generated from absorbed light energy. For effective photothermal therapy, the parameters affecting the induced temperature should be determined before the treatment by modeling the increase in temperature via numerical simulations. However, accurate simulations can only be achieved when utilizing the accurate optical, thermal, and physiological properties of the treated tissue. Here, we propose a multi-wavelength photo-magnetic imaging (PMI) technique that provides quantitative and spatially resolved tissue optical absorption maps at any wavelength within the near-infrared (NIR) window to assist accurate photothermal therapy planning. STUDY DESIGN/MATERIALS AND METHODS: The study was conducted using our recently developed multi-wavelength PMI system, which operates at four laser wavelengths (760, 808, 860, and 980 nm). An agar tissue-simulating phantom containing water, lipid, and ink was illuminated using these wavelengths, and the slight internal laser-induced temperature rise was measured using magnetic resonance thermometry (MRT). The phantom optical absorption was recovered at the used wavelengths using our dedicated PMI image reconstruction algorithm. These absorption maps were then used to resolve the concentration of the tissue chromophores, and thus deduce its optical absorption spectrum in the NIR region based on the Beer-Lambert law. RESULTS: The optical absorption of the phantom was successfully recovered at the used four wavelengths with an average error of ~1.9%. The recovered absorption coefficient was then used to simulate temperature variations inside the phantom. A comparison between the modeled temperature maps and the MRT measured ones showed that these maps are in a good agreement with an average pseudo R2 statistic of 0.992. These absorption values were used to successfully recover the concentration of the used chromophores. Finally, these concentrations are used to accurately calculate the total absorption spectrum of the phantom in the NIR spectral window with an average error as low as ~2.3%. CONCLUSIONS: Multi-wavelength PMI demonstrated a great ability to assess the distribution of tissue chromophores, thus providing its total absorption at any wavelength within the NIR spectral range. Therefore, applications of photothermal therapy applied at NIR wavelengths can benefit from the absorption spectrum recovered by PMI to determine important parameters such as laser power as well as the laser exposure time needed to attain a specific increase in temperature prior to treatment. Lasers Surg. Med. 00:00-00, 2020. © 2020 Wiley Periodicals LLC.

Effects of the time dependence of a bioluminescent source on the tomographic reconstruction.

There are two goals in this simulation study: (1) to show that the time variation of the bioluminescence source can cause artifacts in the tomographic images such that quantification and localization becomes impossible; and (2) to show that the a priori knowledge of the light kinetics can be used to eliminate these artifacts. These goals are motivated by the fact that the half-life of luciferase has been reported as 30 min to 2 h in vivo. We perform two-dimensional simulations. We consider a 40 mm diameter circular region with an inclusion of 6 mm diameter located 10 mm away from the center. The measurement data is simulated using a finite-element-based forward solver. We model the noncontact measurements such that four-wavelength data is collected from four 90 degrees apart views. The results show that the ratio of the total imaging time to the half-life of the exponentially decaying bioluminescent source is the deciding factor in the reconstruction of the source. It is also demonstrated that a priori knowledge of the source kinetics is required to perform tomographic bioluminescence imaging of short half-life bioluminescent sources and the use of spatial a priori information alone is not adequate.

Diffuse optical tomographic reconstruction using multifrequency data.

We investigated the use of multifrequency diffuse optical tomography (MF-DOT) data for the reconstruction of the optical parameters. The experiments were performed in a 63 mm diameter cylindrical phantom containing a 15 mm diameter cylindrical object. Modulation frequencies ranging from 110 MHz to 280 MHz were used in the phantom experiments changing the contrast in absorption of the object with respect to the phantom while keeping the scattering value the same. The diffusion equation was solved using the finite element method. The sensitivity information from each frequency was combined to form a single Jacobian. The inverse problem was solved iteratively by minimizing the difference between the measurements and forward problem using single and multiple modulation frequency data. A multiparameter Tikhonov scheme was used for regularization. The phantom results show that the peak absorption coefficient in a region of interest was obtained with an error less then 5% using two-frequency reconstruction for absorption contrast values up to 2.2 times higher than background and 10% for the absorption contrast values larger than 2.2. The use of two-frequency data is sufficient to improve the quantitative accuracy compared with the single frequency reconstruction with appropriate selection of these frequencies.