Multiwavelength photo-magnetic imaging algorithm improved for direct chromophore concentration recovery using spectral constraints.

Multiwavelength photo-magnetic imaging (PMI) is a novel combination of diffuse optics and magnetic resonance imaging, to the best of our knowledge, that yields tissue chromophore concentration maps with high resolution and quantitative accuracy. Here, we present the first experimental results, to the best of our knowledge, obtained using a spectrally constrained PMI image reconstruction method, where chromophore concentration maps are directly recovered, unlike the conventional two-step approach that requires an intermediate step of reconstructing wavelength-dependent absorption coefficient maps. The imposition of the prior spectral information into the PMI inverse problem improves the reconstructed image quality and allows recovery of highly quantitative concentration maps, which are crucial for effective cancer detection and characterization. The obtained results demonstrate the higher performance of the direct reconstruction method. Indeed, the reconstructed concentration maps are not only of higher quality but also obtained approximately 2 times faster than the conventional method.

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.

T1rho and T2 mapping of ankle cartilage of female and male ballet dancers.

BACKGROUND: Since ballet dancers begin their training before skeletal maturity, accurate and non-invasive identification of cartilage diseases is clinically important. Angle-dependent analysis of T1rho and T2 sequences can be useful for quantification of the composition of cartilage. PURPOSE: To investigate the angle-dependent T1rho and T2 profiles of ankle cartilage in non-dancers and dancers. MATERIAL AND METHODS: Ten female non-dancers, ten female dancers, and 9 male dancers were evaluated using T1rho and T2 mapping sequences. Manual segmentation of talar and tibial cartilage on these images was performed by two radiologists. Inter- and intra-rater reliabilities were calculated using intraclass correlation coefficients (ICCs) and Bland-Altman analysis. Mean thickness and volume of cartilage were estimated. Angle-dependent relaxation time profiles of talar and tibial cartilage were created. RESULTS: ICCs of the number of segmented pixels were poor to excellent. Bland-Altman plots indicated that differences were associated with segment sizes. Segmented cartilage on T1rho demonstrated larger thickness and volume than those on T2 in all populations. Male dancers showed larger cartilage thickness and volume than female dancers and non-dancers. Each cartilage demonstrated angular-dependent T1rho and T2 profiles. Minimal T1rho and T2 values were observed at approximately 180°-200°; higher values were seen at the angle closer to the magic angle. Minimal T2 value of talar cartilage of dancers was larger than that of non-dancers. CONCLUSION: In this small cohort study, regional and sex variations of ankle cartilage T1rho and T2 values in dancers and non-dancers were demonstrated using an angle-dependent approach.

Ex vivo validation of photo-magnetic imaging.

We recently introduced a new high-resolution diffuse optical imaging technique termed photo-magnetic imaging (PMI), which utilizes magnetic resonance thermometry (MRT) to monitor the 3D temperature distribution induced in a medium illuminated with a near-infrared light. The spatiotemporal temperature distribution due to light absorption can be accurately estimated using a combined photon propagation and heat diffusion model. High-resolution optical absorption images are then obtained by iteratively minimizing the error between the measured and modeled temperature distributions. We have previously demonstrated the feasibility of PMI with experimental studies using tissue simulating agarose phantoms. In this Letter, we present the preliminary ex vivo PMI results obtained with a chicken breast sample. Similarly to the results obtained on phantoms, the reconstructed images reveal that PMI can quantitatively resolve an inclusion with a 3 mm diameter embedded deep in a biological tissue sample with only 10% error. These encouraging results demonstrate the high performance of PMI in ex vivo biological tissue and its potential for in vivo imaging.

Comprehensive analytical model for CW laser induced heat in turbid media.

In this work, we present a new analytical approach to model continuous wave laser induced temperature in highly homogeneous turbid media. First, the diffusion equation is used to model light transport and a comprehensive solution is derived analytically by obtaining a special Greens' function. Next, the time-dependent bio-heat equation is used to describe the induced heat increase and propagation within the medium. The bio-heat equation is solved analytically utilizing the separation of variables technique. Our theoretical model is successfully validated using numerical simulations and experimental studies with agarose phantoms and ex-vivo chicken breast samples. The encouraging results show that our method can be implemented as a simulation tool to determine important laser parameters that govern the magnitude of temperature rise within homogenous biological tissue or organs.

Quantitative MRI Analysis of the Talocrural and Talonavicular Joints in Ballet Dancers.

The ankles of ballet dancers are routinely under heavy loading that may lead to osteoarthritic changes. It would be clinically useful to identify such pathology as early as possible in a dancer's career. Therefore, the purpose of this study was to compare quantitative measurements in magnetic resonance (MR) images of the talocrural and talonavicular joints in ballet dancers and healthy non-dancers for use in formulating prediction of chronic injury and degenerative joint disease in these locations. Quantitative measurements in MR images of the talocrural and talonavicular joints were compared in 10 female ballet dancers, 10 healthy female non-dancers, and nine male ballet dancers. Fat-suppressed density-weighted proton, T1rho, and T2 mapping images were acquired with a 3.0 T MR scanner. Medial and lateral subchondral bone distance between the tibia and talus (MSBD and LSBD), axial navicular-talus axis angle (ANT angle), sagittal talar neck angle against the posterior talocalcaneal joint (TN angle), and curvature of navicular surface at the talonavicular joint were measured on sagittal images. The medial subchondral bone distance was found to be significantly larger in female dancers than female non-dancers (4.05 mm vs. 2.75 mm, p < 0.05), whereas there were no significant differences in LSBD (2.63 mm vs. 2.63 mm, p = 0.87). Axial navicular talus angles in female dancers were significantly larger than those in female non-dancers (38.9° vs. 24.3°, p < 0.05). There was a tendency for the TN angle to be smaller and navicular curvature (NC) to be larger in female dancers compared to female non-dancers, though the differences were not significant (TN angle: 16.6° vs. 22.3°, p = 0.09, and NC: 0.186 vs. 0.165, p = 0.28). There were no significant differences in T1rho or T2 values of talonavicular joint cartilage. These results show that the bony anatomy of dancers' ankles may adapt to the stresses placed on them by ballet.

Resolving tissue chromophore concentration at MRI resolution using multi-wavelength photo-magnetic imaging.

Photo-magnetic imaging (PMI) is an emerging optical imaging modality that showed great performance on providing absorption maps with high resolution and quantitative accuracy. As a multi-modality technology, PMI warms up the imaged object using a near infrared laser while temperature variation is measured using magnetic resonance imaging. By probing tissue at multiple wavelengths, concentration of the main tissue chromophores such as oxy- and deoxy-hemoglobin, lipid, and water are obtained then used to derive functional parameters such as total hemoglobin concentration and relative oxygen saturation. In this paper, we present a multi-wavelength PMI system that was custom-built to host five different laser wavelengths. After recovering the high-resolution absorption maps, a least-squared minimization process was used to resolve the different chromophore concentration. The performance of the system was experimentally tested on a phantom with two different dyes. Their concentrations were successfully assessed with high spatial resolution and average accuracy of nearly 80%.

Automatic Breast and Fibroglandular Tissue Segmentation in Breast MRI Using Deep Learning by a Fully-Convolutional Residual Neural Network U-Net.

RATIONALE AND OBJECTIVES: Breast segmentation using the U-net architecture was implemented and tested in independent validation datasets to quantify fibroglandular tissue volume in breast MRI. MATERIALS AND METHODS: Two datasets were used. The training set was MRI of 286 patients with unilateral breast cancer. The segmentation was done on the contralateral normal breasts. The ground truth for the breast and fibroglandular tissue (FGT) was obtained by using a template-based segmentation method. The U-net deep learning algorithm was implemented to analyze the training set, and the final model was obtained using 10-fold cross-validation. The independent validation set was MRI of 28 normal volunteers acquired using four different MR scanners. Dice Similarity Coefficient (DSC), voxel-based accuracy, and Pearson's correlation were used to evaluate the performance. RESULTS: For the 10-fold cross-validation in the initial training set of 286 patients, the DSC range was 0.83-0.98 (mean 0.95 ± 0.02) for breast and 0.73-0.97 (mean 0.91 ± 0.03) for FGT; and the accuracy range was 0.92-0.99 (mean 0.98 ± 0.01) for breast and 0.87-0.99 (mean 0.97 ± 0.01) for FGT. For the entire 224 testing breasts of the 28 normal volunteers in the validation datasets, the mean DSC was 0.86 ± 0.05 for breast, 0.83 ± 0.06 for FGT; and the mean accuracy was 0.94 ± 0.03 for breast and 0.93 ± 0.04 for FGT. The testing results for MRI acquired using four different scanners were comparable. CONCLUSION: Deep learning based on the U-net algorithm can achieve accurate segmentation results for the breast and FGT on MRI. It may provide a reliable and efficient method to process large number of MR images for quantitative analysis of breast density.

Real-time photo-magnetic imaging.

We previously introduced a new high resolution diffuse optical imaging modality termed, photo-magnetic imaging (PMI). PMI irradiates the object under investigation with near-infrared light and monitors the variations of temperature using magnetic resonance thermometry (MRT). In this paper, we present a real-time PMI image reconstruction algorithm that uses analytic methods to solve the forward problem and assemble the Jacobian matrix much faster. The new algorithm is validated using real MRT measured temperature maps. In fact, it accelerates the reconstruction process by more than 250 times compared to a single iteration of the FEM-based algorithm, which opens the possibility for the real-time PMI.

A True Multi-modality Approach for High Resolution Optical Imaging: Photo-Magnetic Imaging.

Multi-modality imaging leverages the competitive advantage of different imaging systems to improve the overall resolution and quantitative accuracy. Our new technique, Photo-Magnetic Imaging (PMI) is one of these true multi-modality imaging approaches, which can provide quantitative optical absorption map at MRI spatial resolution. PMI uses laser light to illuminate tissue and elevate its temperature while utilizing MR thermometry to measure the laser-induced temperature variation with high spatial resolution. The high-resolution temperature maps are later converted to tissue absorption maps by a finite element based inverse solver that is based on modeling of photon migration and heat diffusion in tissue. Previously, we have demonstrated the feasibility of PMI with phantom studies. Recently, we have managed to reduce the laser power under ANSI limit for maximum skin exposure therefore, we have well positioned PMI for in vivo imaging. Currently we are expanding our system by adding multi-wavelength imaging capability. This will allow us not only to resolve spatial distribution of tissue chromophores but also exogenous contrast agents. Although we test PMIs feasibility with animal studies, our future goal is to use PMI for breast cancer imaging due to its high translational potential.

Photo-magnetic imaging: resolving optical contrast at MRI resolution.

In this paper, we establish the mathematical framework of a novel imaging technique, namely photo-magnetic imaging (PMI). PMI uses a laser to illuminate biological tissues and measure the induced temperature variations using magnetic resonance imaging (MRI). PMI overcomes the limitation of conventional optical imaging and allows imaging of the optical contrast at MRI spatial resolution. The image reconstruction for PMI, using a finite-element-based algorithm with an iterative approach, is presented in this paper. The quantitative accuracy of PMI is investigated for various inclusion sizes, depths and absorption values. Then, a comparison between conventional diffuse optical tomography (DOT) and PMI is carried out to illustrate the superior performance of PMI. An example is presented showing that two 2 mm diameter inclusions embedded 4.5 mm deep and located side by side in a 25 mm diameter circular geometry medium are recovered as a single 6 mm diameter object with DOT. However, these two objects are not only effectively resolved with PMI, but their true concentrations are also recovered successfully.

Laser-induced photo-thermal magnetic imaging.

Due to the strong scattering nature of biological tissue, optical imaging beyond the diffusion limit suffers from low spatial resolution. In this letter, we present an imaging technique, laser-induced photo-thermal magnetic imaging (PMI), which uses laser illumination to induce temperature increase in a medium and magnetic resonance imaging to map the spatially varying temperature, which is proportional to absorbed energy. This technique can provide high-resolution images of optical absorption and can potentially be used for small animal as well as breast cancer and lymph node imaging. First, we describe the theory of PMI, including the modeling of light propagation and heat transfer in tissue. We also present experimental data with corresponding predictions from theoretical models, which show excellent agreement.