Search

Significance: The Cherenkov emission spectrum overlaps with that of ambient room light sources. Choice of room lighting devices dramatically affects the efficient detection of Cherenkov emission during patient treatment. Aim: To determine optimal room light sources allowing Cherenkov emission imaging in normally lit radiotherapy treatment delivery rooms. Approach: A variety of commercial light sources and long-pass (LP) filters were surveyed for spectral band separation from the red to near-infrared Cherenkov light emitted by tissue. Their effects on signal-to-noise ratio (SNR), Cherenkov to background signal ratio, and image artifacts were quantified by imaging irradiated tissue equivalent phantoms with an intensified time-gated CMOS camera. Results: Because Cherenkov emission from tissue lies largely in the near-infrared spectrum, a controlled choice of ambient light that avoids this spectral band is ideal, along with a camera that is maximally sensitive to it. An RGB LED light source produced the best SNR out of all sources that mimic room light temperature. A 675-nm LP filter on the camera input further reduced ambient light detected (optical density > 3), achieving maximal SNR for Cherenkov emission near 40. Reduction of the room light signal reduced artifacts from specular reflection on the tissue surface and also minimized spurious Cherenkov signals from non-tissue features such as bolus. Conclusions: LP filtering during image acquisition for near-infrared light in tandem with narrow band LED illuminated rooms improves image quality, trading off the loss of red wavelengths for better removal of room light in the image. This spectral filtering is also critically important to remove specular reflection in the images and allow for imaging of Cherenkov emission through clear bolus. Beyond time-gated external beam therapy systems, the spectral separation methods can be utilized for background removal for continuous treatment delivery methods including proton pencil beam scanning systems and brachytherapy.

Significance: As linear array transducers are widely used in clinical ultrasound imaging, photoacoustic tomography (PAT) with linear arrays is similarly suitable for clinical applications. However, due to the limited-view problem, a linear array has limited performance and leads to artifacts and blurring, which has hindered its broader application. There is a need to address the limited-view problem in PAT imaging with linear arrays. Aim: We investigate potential approaches for improving PAT reconstruction from linear array, by optimizing the detection geometry and implementing iterative reconstruction. Approach: PAT imaging with a single-array, dual-probe configurations in parallel-shape and L-shape, and square-shape configuration are compared in simulations and phantom experiments. An iterative model-based algorithm based on the variance-reduced stochastic gradient descent (VR-SGD) method is implemented. The optimum configuration found in simulation is validated on phantom experiments. Results: PAT imaging with dual-probe detection and VR-SGD algorithm is found to improve the limited-view problem compared to a single probe and provide comparable performance as full-view geometry in simulation. This configuration is validated in experiments where more complete structure is obtained with reduced artifacts compared with a single array. Conclusions: PAT with dual-probe detection and iterative reconstruction is a promising solution to the limited-view problem of linear arrays.

Significance: Quantitative stiffness information can be a powerful aid for tumor or fibrosis diagnosis. Currently, very promising elastography approaches developed for non-contact biomedical imaging are based on transient shear-waves imaging. Transient elastography offers quantitative stiffness information by tracking the propagation of a wave front. The most common method used to compute stiffness from the acquired propagation movie is based on shear-wave time-of-flight calculations. Aim: We introduce an approach to transient shear-wave elastography with spatially coherent sources, able to yield full-field quantitative stiffness maps with reduced artifacts compared to typical artifacts observed in time-of-flight. Approach: A noise-correlation algorithm developed for passive elastography is adapted to spatially coherent narrow or any band sources. This noise-correlation-inspired (NCi) method is employed in parallel with a classic time-of-flight approach. Testing is done on simulation images, experimental validation is conducted with a digital holography setup on controlled homogeneous samples, and full-field quantitative stiffness maps are presented for heterogeneous samples and ex-vivo biological tissues. Results: The NCi approach is first validated on simulations images. Stiffness images processed by the NCi approach on simulated inclusions display significantly less artifacts than with a time-of-flight reconstruction. The adaptability of the NCi algorithm to narrow or any band shear-wave sources was tested successfully. Experimental testing on homogeneous samples demonstrates similar values for both the time-of-flight and the NCi approach. Soft inclusions in agarose sample could be resolved using the NCi method and feasibility on ex-vivo biological tissues is presented. Conclusions: The presented NCi approach was successful in computing quantitative full-field stiffness maps with narrow and broadband source signals on simulation and experimental images from a digital holography setup. Results in heterogeneous media show that the NCi approach could provide stiffness maps with less artifacts than with time-of-flight, demonstrating that a NCi algorithm is a promising approach for shear-wave transient elastography with spatially coherent sources.

Significance: Speckle noise limits the diagnostic capabilities of optical coherence tomography (OCT) images, causing both a reduction in contrast and a less accurate assessment of the microstructural morphology of the tissue. Aim: We present a speckle-noise reduction method for OCT volumes that exploits the advantages of adaptive-noise wavelet thresholding with a wavelet compounding method applied to several frames acquired from consecutive positions. The method takes advantage of the wavelet representation of the speckle statistics, calculated properly from a homogeneous sample or a region of the noisy volume. Approach: The proposed method was first compared quantitatively with different state-of-the-art approaches by being applied to three different clinical dermatological OCT volumes with three different OCT settings. The method was also applied to a public retinal spectral-domain OCT dataset to demonstrate its applicability to different imaging modalities. Results: The results based on four different metrics demonstrate that the proposed method achieved the best performance among the tested techniques in suppressing noise and preserving structural information. Conclusions: The proposed OCT denoising technique has the potential to adapt to different image OCT settings and noise environments and to improve image quality prior to clinical diagnosis based on visual assessment.

Negative ion mode paper spray mass spectrometry (PS-MS) suffers from intense background noise and unstable MS signal. For the first time, we reported fluorinated boron nitride nanosheet (h-FBN) assisted negative ion PS-MS for the detection of a series of molecules. We demonstrated that the introduction of h-FBN can greatly improve the detection sensitivity and signal stability in the negative ion mode.

Significance: High density diffuse optical tomography (HD-DOT) as applied in functional near-infrared spectroscopy (fNIRS) is largely limited to continuous wave (CW) data. Using a single modulation frequency, frequency domain (FD) HD-DOT has recently demonstrated better localization of focal activation as compared to CW data. We show that combining CW and FD measurements and multiple modulation frequencies increases imaging performance in fNIRS. Aim: We evaluate the benefits of multiple modulation frequencies, combining different frequencies as well as CW data in fNIRS HD-DOT. Approach: A layered model was used, with activation occurring within a cortex layer. CW and FD measurements were simulated at 78, 141, and 203 MHz with and without noise. The localization error, full width half maximum, and effective resolution were evaluated. Results: Across the average of the three metrics, at 141 MHz, FD performed 8.4% better than CW, and the combination of CW and FD was 21.7% better than CW. FD measurements at 203 MHz performed 5% better than 78 MHz. Moreover, the three combined modulation frequencies of FD and CW performed up to 3.92% better than 141 MHz alone. Conclusions: We show that combining CW and FD measurements offers better performance than FD alone, with higher modulation frequencies increasing accuracy. Combining CW and FD measurements at multiple modulation frequencies yields the best overall performance.

Proton imaging is a promising technology for proton radiotherapy as it can be used for: (1) direct sampling of the tissue stopping power, (2) input information for multi-modality RSP reconstruction, (3) gold-standard calibration against concurrent techniques, (4) tracking motion and (5) pre-treatment positioning. However, no end-to-end characterization of the image quality (signal-to-noise ratio and spatial resolution, blurring uncertainty) against the dose has been done. This work aims to establish a model relating these characteristics and to describe their relationship with proton energy and object size. The imaging noise originates from two processes: the Coulomb scattering with the nucleus, producing a path deviation, and the energy loss straggling with electrons. The noise is found to increases with thickness crossed and, independently, decreases with decreasing energy. The scattering noise is dominant around high-gradient edge whereas the straggling noise is maximal in homogeneous regions. Image quality metrics are found to behave oppositely against energy: lower energy minimizes both the noise and the spatial resolution, with the optimal energy choice depending on the application and location in the imaged object. In conclusion, the model presented will help define an optimal usage of proton imaging to reach the promised application of this technology and establish a fair comparison with other imaging techniques.

Significance: Speckle noise is an inherent limitation of optical coherence tomography (OCT) images that makes clinical interpretation challenging. The recent emergence of deep learning could offer a reliable method to reduce noise in OCT images. Aim: We sought to investigate the use of deep features (VGG) to limit the effect of blurriness and increase perceptual sharpness and to evaluate its impact on the performance of OCT image denoising (DnCNN). Approach: Fifty-one macula-centered OCT pairs were used in training of the network. Another set of 20 OCT pair was used for testing. The DnCNN model was cascaded with a VGG network that acted as a perceptual loss function instead of the traditional losses of L1 and L2. The VGG network remains fixed during the training process. We focused on the individual layers of the VGG-16 network to decipher the contribution of each distinctive layer as a loss function to produce denoised OCT images that were perceptually sharp and that preserved the faint features (retinal layer boundaries) essential for interpretation. The peak signal-to-noise ratio (PSNR), edge-preserving index, and no-reference image sharpness/blurriness [perceptual sharpness index (PSI), just noticeable blur (JNB), and spectral and spatial sharpness measure (S3)] metrics were used to compare deep feature losses with the traditional losses. Results: The deep feature loss produced images with high perceptual sharpness measures at the cost of less smoothness (PSNR) in OCT images. The deep feature loss outperformed the traditional losses (L1 and L2) for all of the evaluation metrics except for PSNR. The PSI, S3, and JNB estimates of deep feature loss performance were 0.31, 0.30, and 16.53, respectively. For L1 and L2 losses performance, the PSI, S3, and JNB were 0.21 and 0.21, 0.17 and 0.16, and 14.46 and 14.34, respectively. Conclusions: We demonstrate the potential of deep feature loss in denoising OCT images. Our preliminary findings suggest research directions for further investigation.

Significance: Coherence, a fundamental property of waves and fields, plays a key role in photoacoustic image reconstruction. Previously, techniques such as short-lag spatial coherence (SLSC) and filtered delay, multiply, and sum (FDMAS) have utilized spatial coherence to improve the reconstructed resolution and contrast with respect to delay-and-sum (DAS). While SLSC uses spatial coherence directly as the imaging contrast, FDMAS employs spatial coherence implicitly. Despite being more robust against noise, both techniques have their own drawbacks: SLSC does not preserve a relative signal magnitude, and FDMAS shows a reduced contrast-to-noise ratio. Aim: To overcome these limitations, our aim is to develop a beamforming algorithm—generalized spatial coherence (GSC)—that unifies SLSC and FDMAS into a single equation and outperforms both beamformers. Approach: We demonstrated the application of GSC in photoacoustic computed tomography (PACT) through simulation and experiments and compared it to previous beamformers: DAS, FDMAS, and SLSC. Results: GSC outperforms the imaging metrics of previous state-of-the-art coherence-based beamformers in both simulation and experiments. Conclusions: GSC is an innovative reconstruction algorithm for PACT, which combines the strengths of FDMAS and SLSC expanding PACT’s applications.

Significance: Photoacoustic imaging (PAI) can be used to infer molecular information about myocardial health non-invasively in vivo using optical excitation at ultrasonic spatial resolution. For clinical and preclinical linear array imaging systems, conventional delay-and-sum (DAS) beamforming is typically used. However, DAS cardiac PA images are prone to artifacts such as diffuse quasi-static clutter with temporally varying noise-reducing myocardial signal specificity. Typically, multiple frame averaging schemes are utilized to improve the quality of cardiac PAI, which affects the spatial and temporal resolution and reduces sensitivity to subtle PA signal variation. Furthermore, frame averaging might corrupt myocardial oxygen saturation quantification due to the presence of natural cardiac wall motion. In this paper, a spatiotemporal singular value decomposition (SVD) processing algorithm is proposed to reduce DAS PAI artifacts and subsequent enhancement of myocardial signal specificity. Aim: Demonstrate enhancement of PA signals from myocardial tissue compared to surrounding tissues and blood inside the left-ventricular (LV) chamber using spatiotemporal SVD processing with electrocardiogram (ECG) and respiratory signal (ECG-R) gated in vivo murine cardiac PAI. Approach: In vivo murine cardiac PAI was performed by collecting single wavelength (850 nm) photoacoustic channel data on eight healthy mice. A three-dimensional (3D) volume of complex PAI data over a cardiac cycle was reconstructed using a custom ECG-R gating algorithm and DAS beamforming. Spatiotemporal SVD was applied on a two-dimensional Casorati matrix generated using the 3D volume of PAI data. The singular value spectrum (SVS) was then filtered to remove contributions from diffuse quasi-static clutter and random noise. Finally, SVD processed beamformed images were derived using filtered SVS and inverse SVD computations. Results: Qualitative comparison with DAS and minimum variance (MV) beamforming shows that SVD processed images had better myocardial signal specificity, contrast, and target detectability. DAS, MV, and SVD images were quantitatively evaluated by calculating contrast ratio (CR), generalized contrast-to-noise ratio (gCNR), and signal-to-noise ratio (SNR). Quantitative evaluations were done at three cardiac time points (during systole, at end-systole (ES), and during diastole) identified from co-registered ultrasound M-Mode image. Mean CR, gCNR, and SNR values of SVD images at ES were 245, 115.15, and 258.17 times higher than DAS images with statistical significance evaluated with one-way analysis of variance. Conclusions: Our results suggest that significantly better-quality images can be realized using spatiotemporal SVD processing for in vivo murine cardiac PAI.

Significance: Frequency domain diffuse optical spectroscopy (FD-DOS) uses intensity modulated light to measure the absorption and reduced scattering coefficients of turbid media such as biological tissue. Some FD-DOS instruments utilize a single modulation frequency, whereas others use hundreds of frequencies. The effect of modulation frequency choice and measurement bandwidth on optical property (OP) extraction accuracy has not yet been fully characterized. Aim: We aim to assess the effect of modulation frequency selection on OP extraction error and develop a high-speed look-up table (LUT) approach for OP estimation. Approach: We first used noise-free simulations of light transport in homogeneous media to determine optimized iterative inversion model parameters and developed a new multi-frequency LUT method to increase the speed of inversion. We then used experimentally derived noise models for two FD-DOS instruments to generate realistic simulated data for a broad range of OPs and modulation frequencies to test OP extraction accuracy. Results: We found that repeated measurements at a single low-frequency (110 MHz) yielded essentially identical OP errors as a broadband frequency sweep (35 evenly spaced frequencies between 50 and 253 MHz) for these noise models. The inclusion of modulation frequencies >300 MHz diminished overall performance for one of the instruments. Additionally, we developed a LUT inversion algorithm capable of increasing inversion speeds by up to 6 × , with 1000 inversions / s and ∼1 % error when a single modulation frequency was used. Conclusion: These results suggest that simpler single-frequency systems are likely sufficient for many applications and pave the way for a new generation of simpler digital FD-DOS systems capable of rapid, large-volume measurements with real-time feedback.

Significance: The measurement of human vital signs based on photoplethysmography imaging (PPGI) can be severely affected by the interference of various factors in the measurement process; therefore, a lot of complex signal processing techniques are used to remove the influence of the interference. Aim: We comprehensively analyze several methods for color channel combination in the color spaces currently used in PPGI and determine the combination method that can improve the quality of the pulse signal, which results in a modified plane-orthogonal-to-skin based method (POS). Approach: Based on the analysis of the previous studies, 13 methods for color channel combination in the different color spaces, which can be seen as having potential abilities in measuring vital signs, were compared by employing the average value of signal-to-noise ratio (SNR) and the box-plot in the public databases UBFC-RPPG and PURE. In addition, the pulse signal was extracted through the dual-color space transformation (sRGB → intensity normalized RGB → YCbCr) and fine-tuning on the CbCr plane. Results: Among the 13 methods for color channel combination, the signal extracted by the Cb+Cr combination in the YCbCr color space includes the most pulse information. Furthermore, the average SNR of the modified POS for all the used databases is improved by 69.3% compared to POS. Conclusions: The methods using prior knowledge are not only simple to calculate but can significantly increase the SNR, which will provide a great help in the practical use of vital sign measurements based on PPGI.

Significance: Non-contact, camera-based heart rate variability estimation is desirable in numerous applications, including medical, automotive, and entertainment. Unfortunately, camera-based HRV accuracy and reliability suffer due to two challenges: (a) darker skin tones result in lower SNR and (b) relative motion induces measurement artifacts. Aim: We propose an algorithm HRVCam that provides sufficient robustness to low SNR and motion-induced artifacts commonly present in imaging photoplethysmography (iPPG) signals. Approach: HRVCam computes camera-based HRV from the instantaneous frequency of the iPPG signal. HRVCam uses automatic adaptive bandwidth filtering along with discrete energy separation to estimate the instantaneous frequency. The parameters of HRVCam use the observed characteristics of HRV and iPPG signals. Results: We capture a new dataset containing 16 participants with diverse skin tones. We demonstrate that HRVCam reduces the error in camera-based HRV metrics significantly (more than 50% reduction) for videos with dark skin and face motion. Conclusion: HRVCam can be used on top of iPPG estimation algorithms to provide robust HRV measurements making camera-based HRV practical.

Significance: Particle field holography is a versatile technique to determine the size and distribution of moving or stationary particles in air or in a liquid without significant disturbance of the sample volume. Although this technique is applied in biological sample analysis, it is limited to small sample volumes, thus increasing the number of measurements per sample. In this work, we characterize the maximum achievable volume limit based on the specification of a given sensor to realize the development of a potentially low-cost, single-shot, large-volume holographic microscope.Aim: We present mathematical formulas that will aid in the design and development and improve the focusing speed for the numerical reconstruction of registered holograms in particle field holographic microscopes. Our proposed methodology has potential application in the detection of Schistosoma haematobium eggs in human urine samples.Approach: Using the Fraunhofer holography theory for opaque objects, we derived an exact formula for the maximum diffraction-limited volume for an in-line holographic setup. The proof-of-concept device built based on the derived formulas was experimentally validated with urine spiked with cultured Schistosoma haematobium eggs.Results: Results obtained show that for urine spiked with Schistosoma haematobium eggs, the volume thickness is limited to several millimeters due to scattering properties of the sample. The distances of the target particles could be estimated directly from the hologram fringes.Conclusion: The methodology proposed will aid in the development of large-volume holographic microscopes. (C) The Authors. Published by SPIE under a Creative Commons Attribution 4.0 Unported License.

Significance: Photoacoustic tomography (PAT) is a promising emergent modality for the screening and staging of breast cancer. To minimize barriers to clinical translation, it is common to develop PAT systems based upon existing ultrasound hardware, which can entail significant design challenges in terms of light delivery. This often results in inherently non-uniform fluence within the tissue and should be accounted for during image reconstruction. Aim: We aim to integrate PAT into an automated breast ultrasound scanner with minimal change to the existing system. Approach: We designed and implemented an illuminator that directs spatially non-uniform light to the tissue near the acquisition plane of the imaging array. We developed a graphics processing unit-accelerated reconstruction method, which accounts for this illumination geometry by modeling the structure of the light in the sample. We quantified the performance of this system using a custom, modular photoacoustic phantom and graphite rods embedded in chicken breast tissue. Results: Our illuminator provides a fluence of 2.5 mJ cm−2 at the tissue surface, which was sufficient to attain a signal-to-noise ratio (SNR) of 8 dB at 2 cm in chicken breast tissue and image 0.25-mm features at depths of up to 3 cm in a medium with moderate optical scattering. Our reconstruction scheme is 200× faster than a CPU implementation; it provides a 25% increase in SNR at 2 cm in chicken breast tissue and lowers image error by an average of 31% at imaging depths >1.5 cm compared with a method that does not account for the inhomogeneity of the illumination or the transducer directivity. Conclusions: A fan-shaped illumination geometry is feasible for PAT; however, it is important to account for non-uniform fluence in illumination scenarios such as this. Future work will focus on increasing fluence and further optimizing the ultrasound hardware to improve SNR and overall image quality.

Reducing the bit-depth is an effective approach to lower the cost of optical coherence tomography (OCT) systems and increase the transmission efficiency in data acquisition and telemedicine. However, a low bit-depth will lead to the degeneration of the detection sensitivity thus reduce the signal-to-noise ratio (SNR) of OCT images. In this paper, we propose to use deep learning for the reconstruction of the high SNR OCT images from the low bit-depth acquisition. Its feasibility was preliminarily evaluated by applying the proposed method to the quantized $3\sim8$-bit data from native 12-bit interference fringes. We employed a pixel-to-pixel generative adversarial network architecture in the low to high bit-depth OCT image transition. Retinal OCT data of a healthy subject from a homemade spectral-domain OCT system was used in the study. Extensively qualitative and quantitative results show this deep-learning-based approach could significantly improve the SNR of the low bit-depth OCT images especially at the choroidal region. Superior similarity and SNR between the reconstructed images and the original 12-bit OCT images could be derived when the bit-depth $\geq 5$. This work demonstrates the proper integration of OCT and deep learning could benefit the development of healthcare in low-resource settings., Comment: submitted

Significance: Photoacoustic imaging (PAI) has been greatly developed in a broad range of diagnostic applications. The efficiency of light to sound conversion in PAI is limited by the ubiquitous noise arising from the tissue background, leading to a low signal-to-noise ratio (SNR), and thus a poor quality of images. Frame averaging has been widely used to reduce the noise; however, it compromises the temporal resolution of PAI. Aim: We propose an approach for photoacoustic (PA) signal denoising based on a combination of low-pass filtering and sparse coding (LPFSC). Approach: LPFSC method is based on the fact that PA signal can be modeled as the sum of low frequency and sparse components, which allows for the reduction of noise levels using a hybrid alternating direction method of multipliers in an optimization process. Results: LPFSC method was evaluated using in-silico and experimental phantoms. The results show a 26% improvement in the peak SNR of PA signal compared to the averaging method for in-silico data. On average, LPFSC method offers a 63% improvement in the image contrast-to-noise ratio and a 33% improvement in the structural similarity index compared to the averaging method for objects located at three different depths, ranging from 10 to 20 mm, in a porcine tissue phantom. Conclusions: The proposed method is an effective tool for PA signal denoising, whereas it ultimately improves the quality of reconstructed images, especially at higher depths, without limiting the image acquisition speed.

Significance: Researchers have made great progress in single-image super-resolution (SISR) using deep convolutional neural networks. However, in the field of leukocyte imaging, the performance of existing SISR methods is still limited as it fails to thoroughly explore the geometry and structural consistency of leukocytes. The inaccurate super-resolution (SR) results will hinder the pathological study of leukocytes, since the structure and cell lineage determine the types of leukocyte and will significantly affect the subsequent inspection. Aim: We propose a deep network that takes full use of the geometry prior and structural consistency of the leukocyte images. We establish and annotate a leukocyte dataset, which contains five main types of leukocytes (basophil, eosinophil, monocyte, lymphocyte, and neutrophil), for learning the structure and geometry information. Approach: Our model is composed of two modules: prior network and SR network. The prior network estimates the parsing map of the low-resolution (LR) image, and then the SR network takes both the estimated parsing map and LR image as input to predict the final high-resolution image. Result: Experiments show that the geometry prior and structural consistency in use obviously improves the SR performance of leukocyte images, enhancing the peak-signal-to-noise ratio (PSNR) by about 0.4 dB in our benchmark. Conclusion: As proved by our leukocyte SR benchmark, the proposed method significantly outperforms state-of-the-art SR methods. Our method not only improves the PSNR and structural similarity indices, but also accurately preserves the structural details of leukocytes. The proposed method is believed to have potential use in the wide-field cell prescreening by simply using a low-power objective.

Significance: Diffuse correlation spectroscopy (DCS) is an established optical modality that enables noninvasive measurements of blood flow in deep tissue by quantifying the temporal light intensity fluctuations generated by dynamic scattering of moving red blood cells. Compared with near-infrared spectroscopy, DCS is hampered by a limited signal-to-noise ratio (SNR) due to the need to use small detection apertures to preserve speckle contrast. However, DCS is a dynamic light scattering technique and does not rely on hemoglobin contrast; thus, there are significant SNR advantages to using longer wavelengths (>1000 nm) for the DCS measurement due to a variety of biophysical and regulatory factors. Aim: We offer a quantitative assessment of the benefits and challenges of operating DCS at 1064 nm versus the typical 765 to 850 nm wavelength through simulations and experimental demonstrations. Approach: We evaluate the photon budget, depth sensitivity, and SNR for detecting blood flow changes using numerical simulations. We discuss continuous wave (CW) and time-domain (TD) DCS hardware considerations for 1064 nm operation. We report proof-of-concept measurements in tissue-like phantoms and healthy adult volunteers. Results: DCS at 1064 nm offers higher intrinsic sensitivity to deep tissue compared with DCS measurements at the typically used wavelength range (765 to 850 nm) due to increased photon counts and a slower autocorrelation decay. These advantages are explored using simulations and are demonstrated using phantom and in vivo measurements. We show the first high-speed (cardiac pulsation-resolved), high-SNR measurements at large source–detector separation (3 cm) for CW-DCS and late temporal gates (1 ns) for TD-DCS. Conclusions: DCS at 1064 nm offers a leap forward in the ability to monitor deep tissue blood flow and could be especially useful in increasing the reliability of cerebral blood flow monitoring in adults.

Significance: The use of diffuse correlation spectroscopy (DCS) has shown efficacy in research studies as a technique capable of noninvasively monitoring blood flow in tissue with applications in neuromonitoring, exercise science, and breast cancer management. The ability of DCS to resolve blood flow in these tissues is related to the optical sensitivity and signal-to-noise ratio (SNR) of the measurements, which in some cases, particularly adult cerebral blood flow measurements, is inadequate in a significant portion of the population. Improvements to DCS sensitivity and SNR could allow for greater clinical translation of this technique. Aim: Interferometric diffuse correlation spectroscopy (iDCS) was characterized and compared to traditional homodyne DCS to determine possible benefits of utilizing heterodyne detection. Approach: An iDCS system was constructed by modifying a homodyne DCS system with fused fiber couplers to create a Mach–Zehnder interferometer. Comparisons between homodyne and heterodyne detection were performed using an intralipid phantom characterized at two extended source–detector separations (2.4, 3.6 cm), different photon count rates, and a range of reference arm power levels. Characterization of the iDCS signal mixing was compared to theory. Precision of the estimation of the diffusion coefficient and SNR of the autocorrelation curve were compared between different measurement conditions that mimicked what would be seen in vivo. Results: The mixture of signals present in the heterodyne autocorrelation function was found to agree with the derived theory and resulted in accurate measurement of the diffusion coefficient of the phantom. Improvement of the SNR of the autocorrelation curve up to ∼2× and up to 80% reduction in the variability of the diffusion coefficient fit were observed for all measurement cases as a function of increased reference arm power. Conclusions: iDCS has the potential to improve characterization of blood flow in tissue at extended source–detector separations, enhancing depth sensitivity and SNR.

Significance: Photoacoustic-based visual servoing is a promising technique for surgical tool tip tracking and automated visualization of photoacoustic targets during interventional procedures. However, one outstanding challenge has been the reliability of obtaining segmentations using low-energy light sources that operate within existing laser safety limits. Aim: We developed the first known graphical processing unit (GPU)-based real-time implementation of short-lag spatial coherence (SLSC) beamforming for photoacoustic imaging and applied this real-time algorithm to improve signal segmentation during photoacoustic-based visual servoing with low-energy lasers. Approach: A 1-mm-core-diameter optical fiber was inserted into ex vivo bovine tissue. Photoacoustic-based visual servoing was implemented as the fiber was manually displaced by a translation stage, which provided ground truth measurements of the fiber displacement. GPU-SLSC results were compared with a central processing unit (CPU)-SLSC approach and an amplitude-based delay-and-sum (DAS) beamforming approach. Performance was additionally evaluated with in vivo cardiac data. Results: The GPU-SLSC implementation achieved frame rates up to 41.2 Hz, representing a factor of 348 speedup when compared with offline CPU-SLSC. In addition, GPU-SLSC successfully recovered low-energy signals (i.e., ≤268 μJ) with mean ± standard deviation of signal-to-noise ratios of 11.2±2.4 (compared with 3.5±0.8 with conventional DAS beamforming). When energies were lower than the safety limit for skin (i.e., 394.6 μJ for 900-nm wavelength laser light), the median and interquartile range (IQR) of visual servoing tracking errors obtained with GPU-SLSC were 0.64 and 0.52 mm, respectively (which were lower than the median and IQR obtained with DAS by 1.39 and 8.45 mm, respectively). GPU-SLSC additionally reduced the percentage of failed segmentations when applied to in vivo cardiac data. Conclusions: Results are promising for the use of low-energy, miniaturized lasers to perform GPU-SLSC photoacoustic-based visual servoing in the operating room with laser pulse repetition frequencies as high as 41.2 Hz.

Rapid and sensitive detection of metabolites and chemical residues in human tears is highly beneficial for understanding eye health. In this study, Schirmer paper was used for noninvasive microsampling of human tears, and then paper spray mass spectrometry (PSMS) was performed for direct analysis of human tears. Schirmer PSMS was successfully used for rapid diagnosis of dry-eye syndrome by detecting the volume and metabolites of human tears. Drugs of abuse, therapeutic drugs, and pharmacodynamics in human tears were also investigated by Schirmer PSMS. Furthermore, specific markers of environmental exposures in the air to human eyes, including volatile organic compounds, aerosol, and smoke, were unambiguously sampled and detected in human tears using Schirmer PSMS. Excellent analytical performances were achieved, including single-use, low-sample consumption (1.0 μL), rapid analysis (the whole analytical procedure completed within 3 min), high sensitivity (absolute limit of detection less than or equal to 0.5 pg, signal-to-noise ratio greater than or equal to 3), good reproducibility (relative standard deviation less than 10%, n = 3), and accurate quantitation (average deviation less than 3%, n = 3). Overall, our results showed that Schirmer PSMS is a highly effective method for direct tear analysis and is expected to be a convenient tool for human tear analysis in significant clinical applications.

Purpose: Photon-counting silicon strip detectors are attracting interest for use in next-generation CT scanners. For CT detectors in a clinical environment, it is desirable to have a low power consumption. However, decreasing the power consumption leads to higher noise. This is particularly detrimental for silicon detectors, which require a low noise floor to obtain a good dose efficiency. The increase in noise can be mitigated using a longer shaping time in the readout electronics. This also results in longer pulses, which requires an increased deadtime, thereby degrading the count-rate performance. However, as the photon flux varies greatly during a typical CT scan, not all projection lines require a high count-rate capability. We propose adjusting the shaping time to counteract the increased noise that results from decreasing the power consumption. Approach: To show the potential of increasing the shaping time to decrease the noise level, synchrotron measurements were performed using a detector prototype with two shaping time settings. From the measurements, a simulation model was developed and used to predict the performance of a future channel design. Results: Based on the synchrotron measurements, we show that increasing the shaping time from 28.1 to 39.4 ns decreases the noise and increases the signal-to-noise ratio with 6.5% at low count rates. With the developed simulation model, we predict that a 50% decrease in power can be attained in a proposed future detector design by increasing the shaping time with a factor of 1.875. Conclusion: Our results show that the shaping time can be an important tool to adapt the pulse length and noise level to the photon flux and thereby optimize the dose efficiency of photon-counting silicon detectors. (C) The Authors. Published by SPIE under a Creative Commons Attribution 4.0 Unported License. Funding Agencies|Erling-Persson Family Foundation; MedTechLabs; European Unions Horizon 2020 Research and Innovation Programme [830294]; project CALIPSO plus from the EU Framework Programme for Research and Innovation HORIZON 2020 [730872]

Significance: The large background, narrow dynamic range, and detector saturation have been the common limiting factors in stimulated emission (SE)-based pump-probe microscopy, attributed to the very small signal overriding the very intense laser probe beam. To better differentiate the signal of interest from the background, lock-in detection is used to measure the fluorescence quenching, which is termed spontaneous loss (SL). The advantages are manifold. The spontaneous fluorescence signal can be well separated from both the pump and the probe beams with filters, thus eliminating the background, enlarging the dynamic range, and avoiding the saturation of the detector. Aim: We propose and demonstrate an integrated pump-probe microscopy technique based on lock-in detection for background removal and dynamic range enhancement through SL detection. Approach: The experimental setup is configured with a pulsed diode laser at a wavelength λpu=635 nm, acting as a pump (excitation) and a mode-locked Ti:sapphire laser at a central wavelength λpr=780 nm, serving as the probe beam (stimulation). Both pulse trains are temporally synchronized through high precision delay control by adjusting the length of the triggering cables. The pump and probe beams are alternatively modulated at different frequencies f1 and f2 to extract the stimulated gain (SG) and SL signal. Results: SG signal shows saturation due to the irradiation of the intense probe beam onto the photodetector. However, the detector saturation does not occur at high probe beam power for SL detection. The fluorescence lifetime images are acquired with reduced background. The theoretical signal-to-noise ratios for SG and SL are also estimated by photon statistics. Conclusion: We have confirmed that the detection of SL allows the elimination of the large background without photodetector saturation, which commonly exists in SG configuration. This modality would allow unprecedented manipulation and investigation of fluorophores in fluorescence imaging.

SIGNIFICANCE: Colorectal cancer is one of the major causes of cancer-related deaths worldwide. Surgical removal of the cancerous growth is the primary treatment for this disease. A colorectal cancer surgery, however, is often unsuccessful due to the anastomotic failure that may occur following the surgical incision. Prevention of an anastomotic failure requires continuous monitoring of intestinal tissue viability during and after colorectal surgery. To date, no clinical technology exists for the dynamic and continuous monitoring of the intestinal perfusion.\ud \ud AIM: A dual-wavelength indwelling bowel photoplethysmography (PPG) sensor for the continuous monitoring of intestinal viability was proposed and characterized through a set of in silico and in vivo investigations.\ud \ud APPROACH: The in silico investigation was based on a Monte Carlo model that was executed to quantify the variables such as penetration depth and detected intensity with respect to the sensor-tissue separations and tissue perfusion. Utilizing the simulated information, an indwelling reflectance PPG sensor was designed and tested on 20 healthy volunteers. Two sets of in vivo studies were performed using the driving current intensities 20 and 40 mA for a comparative analysis, using buccal tissue as a proxy tissue-site.\ud \ud RESULTS: Both simulated and experimental results showed the efficacy of the sensor to acquire good signals through the "contact" to a "noncontact" separation of 5 mm. A very slow wavelength-dependent variation was shown in the detected intensity at the normal and hypoxic states of the tissue, whereas a decay in the intensity was found with the increasing submucosal-blood volume. The simulated detected-to-incident-photon-ratio and the experimental signal-to-noise ratio exhibited strong positive correlations, with the Pearson product-moment correlation coefficient R ranging between 0.65 and 0.87.\ud \ud CONCLUSIONS: The detailed feasibility analysis presented will lead to clinical trials utilizing the proposed sensor.

Significance: It is commonly assumed that using the objective lens to create a tightly focused light spot for illumination provides a twofold resolution improvement over the Rayleigh resolution limit and that resolution improvement is independent of object properties. Nevertheless, such an assumption has not been carefully examined. We examine this assumption by analyzing the performance of two super-resolution methods, known as image scanning microscopy (ISM) and illumination-enhanced sparsity (IES). Aim: We aim to identify the fundamental differences between the two methods, and to provide examples that help researchers determine which method to utilize for different imaging conditions. Approach: We input the same image datasets into the two methods and analyze their restorations. In numerical simulations, we design objects of distinct brightness and sparsity levels for imaging. We use biological imaging experiments to verify the simulation results. Results: The resolution of IES often exceeds twice the Rayleigh resolution limit when imaging sparse objects. A decrease in object sparsity negatively affects the resolution improvement in both methods. Conclusions: The IES method is superior for imaging sparse objects with its main features being bright and small against a dark, large background. For objects that are largely bright with small dark features, the ISM method is favorable.

RNA sequencing (RNA-seq) has become a standard procedure to investigate transcriptional changes between conditions and is routinely used in research and clinics. While standard differential expression (DE) analysis between two conditions has been extensively studied, and improved over the past decades, RNA-seq time course (TC) DE analysis algorithms are still in their early stages. In this study, we compare, for the first time, existing TC RNA-seq tools on an extensive simulation data set and validated the best performing tools on published data. Surprisingly, TC tools were outperformed by the classical pairwise comparison approach on short time series (, Briefings in Bioinformatics, 20 (1), ISSN:1467-5463, ISSN:1477-4054

As photoacoustic imaging (PAI) technology matures, computational modeling will increasingly represent a critical tool for facilitating clinical translation through predictive simulation of real-world performance under a wide range of device and biological conditions. While modeling currently offers a rapid, inexpensive tool for device development and prediction of fundamental image quality metrics (e.g., spatial resolution and contrast ratio), rigorous verification and validation will be required of models used to provide regulatory-grade data that effectively complements and/or replaces in vivo testing. To address methods for establishing model credibility, we developed an integrated computational model of PAI by coupling a previously developed three-dimensional Monte Carlo model of tissue light transport with a two-dimensional (2D) acoustic wave propagation model implemented in the well-known k-Wave toolbox. We then evaluated ability of the model to predict basic image quality metrics by applying standardized verification and validation principles for computational models. The model was verified against published simulation data and validated against phantom experiments using a custom PAI system. Furthermore, we used the model to conduct a parametric study of optical and acoustic design parameters. Results suggest that computationally economical 2D acoustic models can adequately predict spatial resolution, but metrics such as signal-to-noise ratio and penetration depth were difficult to replicate due to challenges in modeling strong clutter observed in experimental images. Parametric studies provided quantitative insight into complex relationships between transducer characteristics and image quality as well as optimal selection of optical beam geometry to ensure adequate image uniformity. Multidomain PAI simulation tools provide high-quality tools to aid device development and prediction of real-world performance, but further work is needed to improve model fidelity, especially in reproducing image noise and clutter.

A significant number of patients receiving breast-conserving surgery (BCS) for invasive carcinoma and ductal carcinoma in situ (DCIS) may need reoperation following tumor-positive margins from final histopathology tests. All current intraoperative margin assessment modalities have specific limitations. As a first step towards the development of a compact system for intraoperative specimen imaging based on edge illumination x-ray phase contrast, we prove that the system’s dimensions can be reduced without affecting imaging performance. We analysed the variation in noise and contrast to noise ratio (CNR) with decreasing system length using the edge illumination x-ray phase contrast imaging setup. Two-(planar) and three-(computed tomography (CT)) dimensional imaging acquisitions of custom phantoms and a breast tissue specimen were made. Dedicated phase retrieval algorithms were used to separate refraction and absorption signals. A ‘single-shot’ retrieval method was also used, to retrieve thickness map images, due to its simple acquisition procedure and reduced acquisition times. Experimental results were compared to numerical simulations where appropriate. The relative contribution of dark noise signal in integrating detectors is significant for low photon count statistics acquisitions. Under constant exposure factors and magnification, a more compact system provides an increase in CNR. Superior CNR results were obtained for refraction and thickness map images when compared to absorption images. Results indicate that the ‘single-shot’ acquisition method is preferable for a compact CT intraoperative specimen scanner; it allows for shorter acquisition times and its combination of the absorption and refraction signals ultimately leads to a higher contrast. The first CT images of a breast specimen acquired with the compact system provided promising results when compared to those of the longer length system.

Since it was first demonstrated more than a decade ago, the single-pixel camera concept has been used in numerous applications in which it is necessary or advantageous to reduce the channel count, cost, or data volume. Here, three-dimensional (3-D), compressed-sensing photoacoustic tomography (PAT) is demonstrated experimentally using a single-pixel camera. A large area collimated laser beam is reflected from a planar Fabry–Pérot ultrasound sensor onto a digital micromirror device, which patterns the light using a scrambled Hadamard basis before it is collected into a single photodetector. In this way, inner products of the Hadamard patterns and the distribution of thickness changes of the FP sensor—induced by the photoacoustic waves—are recorded. The initial distribution of acoustic pressure giving rise to those photoacoustic waves is recovered directly from the measured signals using an accelerated proximal gradient-type algorithm to solve a model-based minimization with total variation regularization. Using this approach, it is shown that 3-D PAT of imaging phantoms can be obtained with compression rates as low as 10%. Compressed sensing approaches to photoacoustic imaging, such as this, have the potential to reduce the data acquisition time as well as the volume of data it is necessary to acquire, both of which are becoming increasingly important in the drive for faster imaging systems giving higher resolution images with larger fields of view.

Fluorescence microscopy plays a vital role in understanding the subcellular structures of living cells. However, it requires considerable effort in sample preparation related to chemical fixation, staining, cost, and time. To reduce those factors, we present a virtual fluorescence staining method based on deep neural networks (VirFluoNet) to transform fluorescence images of molecular labels into other molecular fluorescence labels in the same field-of-view. To achieve this goal, we develop and train a conditional generative adversarial network (cGAN) to perform digital fluorescence imaging demonstrated on human osteosarcoma U2OS cell fluorescence images captured under Cell Painting staining protocol. A detailed comparative analysis is also conducted on the performance of the cGAN network between predicting fluorescence channels based on phase contrast or based on another fluorescence channel using human breast cancer MDA-MB-231 cell line as a test case. In addition, we implement a deep learning model to perform autofocusing on another human U2OS fluorescence dataset as a preprocessing step to defocus an out-focus channel in U2OS dataset. A quantitative index of image prediction error is introduced based on signal pixel-wise spatial and intensity differences with ground truth to evaluate the performance of prediction to high-complex and throughput fluorescence. This index provides a rational way to perform image segmentation on error signals and to understand the likelihood of mis-interpreting biology from the predicted image. In total, these findings contribute to the utility of deep learning image regression for fluorescence microscopy datasets of biological cells, balanced against savings of cost, time, and experimental effort. Furthermore, the approach introduced here holds promise for modeling the internal relationships between organelles and biomolecules within living cells., 20 pages, 9 figures

Multiphoton microscopy provides a suitable technique for imaging biological tissues with submicrometer resolution. Usually a Gaussian beam (GB) is used for illumination, leading to a reduced power efficiency in the multiphoton response and vignetting for a square-shaped imaging area. A flat-top beam (FTB) provides a uniform spatial intensity distribution that equalizes the probability of a multiphoton effect across the imaging area. We employ a customized widefield multiphoton microscope to compare the performance of a square-shaped FTB illumination with that based on using a GB, for both two-photon fluorescence (TPF) and second-harmonic generation (SHG) imaging. The variation in signal-to-noise ratio across TPF images of fluorescent dyes spans ∼5.6 dB for the GB and ∼1.2 dB for the FTB illumination, respectively. For the GB modality, TPF images of mouse colon and Convallaria root, and SHG images of chicken tendon and human breast biopsy tissue showcase ∼20% area that are not imaged due to either insufficient or lack of illumination. For quantitative analysis that depends on the illuminated area, this effect can potentially lead to inaccuracies. This work emphasizes the applicability of FTB illumination to multiphoton applications.

Near-infrared spectroscopy (NIRS) is widely used in biomedical optics with applications ranging from basic science, such as in functional neuroimaging, to clinical, as in pulse oximetry. Despite the relatively low absorption of tissue in the near-infrared, there is still a significant amount of optical attenuation produced by the highly scattering nature of tissue. Because of this, designers of NIRS systems have to balance source optical power and source–detector separation to maximize the signal-to-noise ratio (SNR). However, theoretical estimations of SNR neglect the effects of speckle. Speckle manifests as fluctuations of the optical power received at the detector. These fluctuations are caused by interference of the multiple random paths taken by photons in tissue. We present a model for the NIRS SNR that includes the effects of speckle. We performed experimental validations with a NIRS system to show that it agrees with our model. Additionally, we performed computer simulations based on the model to estimate the contribution of speckle noise for different collection areas and source–detector separations. We show that at short source–detector separation, speckle contributes most of the noise when using long coherence length sources. Considering this additional noise is especially important for hybrid applications that use NIRS and speckle contrast simultaneously, such as in diffuse correlation spectroscopy.

To achieve fast imaging and large field of view (FOV), we improved our multimodal imaging system, which integrated optical resolution photoacoustic microscopy, optical coherence tomography (OCT), and confocal fluorescence microscopy in one platform, by combining optical scanning with mechanical scanning. To ensure good focusing of the objective lens over all the imaged area, we employed OCT-guided dynamic focusing. Different from our previous point-by-point dynamic focusing, we employed an area-by-area focusing adjustment strategy, in which each fast optical scanning area has a fixed focusing depth. We have demonstrated the performance of the system by imaging biological samples ex vivo (plant leaf) and in vivo (mouse ear). The system has achieved uniform resolution in an FOV of 10 mm×10 mm with an imaging time of about 5 min.

The silicon photomultiplier (SIPM) is an emerging detector technology that enables both high sensitivity and high dynamic range detection of visible and near-infrared light at a fraction of the cost of conventional vacuum tube photomultiplier tubes (PMTs). A low-cost detection circuit is presented and the performance of a commercial SIPM is evaluated for high-speed laser scanning microscopy applications. For moderate-to-high-speed fluorescent imaging applications, the measurements and imaging results indicate that the SIPM exceeds the sensitivity of GaAsP PMTs, while providing higher dynamic range and better saturation behavior. For low speed or applications requiring large detector areas, the GaAsP PMT retains a sensitivity advantage due to large area and lower dark counts. The calculations presented show that, above a critical detection bandwidth, the SIPM sensitivity exceeds that of a GaAsP PMT.

Despite the utility of tumour characterisation using quantitative parameter maps from multi-b-value diffusion-weighted MRI (DWI), clinicians often prefer the use of the image with highest diffusion-weighting (b-value), for instance for defining regions of interest (ROIs). However, these images are typically degraded by noise, as they do not utilize the information from the full acquisition. We present a principal component analysis (PCA) approach for model-free denoising of DWI data. PCA-denoising was compared to synthetic MRI, where a diffusion model is fitted for each voxel and a denoised image at a given b-value is generated from the model fit. A quantitative comparison of systematic and random errors was performed on data simulated using several diffusion models (mono-exponential, bi-exponential, stretched-exponential and kurtosis). A qualitative visual comparison was also performed for in vivo images in six healthy volunteers and three pancreatic cancer patients. In simulations, the reduction in random errors from PCA-denoising was substantial (up to 55%) and similar to synthetic MRI (up to 53%). Model-based synthetic MRI denoising resulted in substantial (up to 29% of signal) systematic errors, whereas PCA-denoising was able to denoise without introducing systematic errors (less than 2%). In vivo, the signal-to-noise ratio (SNR) and sharpness of PCA-denoised images were superior to synthetic MRI, resulting in clearer tumour boundaries. In the presence of motion, PCA-denoising did not cause image blurring, unlike image averaging or synthetic MRI. Multi-b-value MRI can be denoised model-free with our PCA-denoising strategy that reduces noise to a level similar to synthetic MRI, but without introducing systematic errors associated with the synthetic MRI method.

Blockage of healthy blood vessels by blood clots can lead to serious or even life-threatening complications. The use of a combined ultrasound (US) and photoacoustic (PA) imaging was explored for blood clot monitoring during microbubble-assisted sonothrombolysis. PA imaging is an emerging hybrid imaging modality that has garnered the attention of the biomedical imaging community in recent years. It enables the study of the composition of a blood clot due to its sensitivity toward optical absorption. Here, in vitro imaging of the side of a blood clot facing the microbubbles was done over time. The US and PA signal-to-noise (SNR) ratio value changes during microbubble-assisted sonothrombolysis were studied for two different local environments: blood clot in deionized water and blood clot in blood. In the first case, US and PA SNR values increased by 4.6% and reduced by 20.8%, respectively after 30 min of sonothrombolysis treatment. After 10 min of sonothrombolysis treatment of the blood clot in blood, the US and PA SNR values increased by 7.7% and 38.3%, respectively. The US and PA SNR value changes were recorded in response to its local environment. This technique can be used to determine the final composition of the blood clot which may, in turn, help in the administration of clot-dissolving drugs. MOE (Min. of Education, S’pore) Published version

Significance: Nanosecond-pulsed laser has proven to be used to obtain the velocity of blood using the speckle contrast method. Without the scanning time, it has potential for achieving fast two-dimensional blood flow images in a photoacoustic imaging system with the same pulsed laser. Aim: Our study aimed to evaluate the qualities of regional cerebral blood flow (rCBF) obtained in a laser speckle contrast imaging (LSCI) system using continuous wave (cw) and nanosecond pulse laser sources. Approach: First, a LSCI system consisting of a cw laser with a wavelength of 632.8 nm and a cw laser/nanosecond pulse laser with a wavelength of 532 nm was developed. This system was used to obtain rCBF images of mouse in vivo with two different laser sources. Results: Continuous wave lasers (532 and 632.8 nm) show different imaging characteristics for rCBF imaging. The rCBF images obtained using 532-nm nanosecond pulse laser showed higher resolution than those using 532-nm cw laser. There was no significant difference in the results using nanosecond pulse laser among various pulse widths or repetition rates. Conclusions: It is proved that a nanosecond pulse laser could be used for LSCI.

We analyze transducer-matched multipulse excitation as a method for improving of the signal-to-noise ratio (SNR) for diode laser-based photoacoustic systems. We discuss the principle of the technique, its advantages, and potential drawbacks and perform measurements to analyze the obtainable SNR increase. We show in experiment and computationally that a lower boundary estimate of 1.2 to 1.8 fold SNR improvement can be provided using transducer-matched pulse bursts, depending on the transducer and particular arrangement. Finally, we analyze implications that the transducer resonance effects may have on the recently introduced advanced photoacoustic techniques. The findings are of immediate interest to modalities utilizing dense pulse sequences and systems possessing limited pulse energy. In particular, transducer-matched multipulse excitation may be beneficial for diode-based photoacoustic systems operated with transducers in the range of 1 to 5 MHz since the required hardware is readily available.

Live-subject microscopies, including microendoscopy and other related technologies, offer promise for basic biology research as well as the optical biopsy of disease in the clinic. However, cellular resolution generally comes with the trade-off of a microscopic field-of-view. Microimage mosaicking enables stitching many small scenes together to aid visualization, quantitative interpretation, and mapping of microscale features, for example, to guide surgical intervention. The development of hyperspectral and multispectral systems for biomedical applications provides motivation for adapting mosaicking algorithms to process a number of simultaneous spectral channels. We present an algorithm that mosaics multichannel video by correlating channels of consecutive frames as a basis for efficiently calculating image alignments. We characterize the noise tolerance of the algorithm by using simulated video with known ground-truth alignments to quantify mosaicking accuracy and speed, showing that multiplexed molecular imaging enhances mosaic accuracy by leveraging observations of distinct molecular constituents to inform frame alignment. A simple mathematical model is introduced to characterize the noise suppression provided by a given group of spectral channels, thus predicting the performance of selected subsets of data channels in order to balance mosaic computation accuracy and speed. The characteristic noise tolerance of a given number of channels is shown to improve through selection of an optimal subset of channels that maximizes this model. We also demonstrate that the multichannel algorithm produces higher quality mosaics than the analogous single-channel methods in an empirical test case. To compensate for the increased data rate of hyperspectral video compared to single-channel systems, we employ parallel processing via GPUs to alleviate computational bottlenecks and to achieve real-time mosaicking even for video-rate multichannel systems anticipated in the future. This implementation paves the way for real-time multichannel mosaicking to accompany next-generation hyperspectral and multispectral video microscopy.

The compromise between lateral resolution and usable imaging depth range is a bottleneck for optical coherence tomography (OCT). Existing solutions for optical coherence microscopy (OCM) suffer from either large data size and long acquisition time or a nonideal point spread function. We present volumetric OCM of mouse brain ex vivo with a large depth coverage by leveraging computational adaptive optics (CAO) to significantly reduce the number of OCM volumes that need to be acquired with a Gaussian beam focused at different depths. We demonstrate volumetric reconstruction of ex-vivo mouse brain with lateral resolution of 2.2 μm, axial resolution of 4.7 μm, and depth range of ∼1.2 mm optical path length, using only 11 OCT data volumes acquired on a spectral-domain OCM system. Compared to focus scanning with step size equal to the Rayleigh length of the beam, this is a factor of 4 fewer datasets required for volumetric imaging. Coregistered two-photon microscopy confirmed that CAO-OCM reconstructions can visualize various tissue microstructures in the brain. Our results also highlight the limitations of CAO in highly scattering media, particularly when attempting to reconstruct far from the focal plane or when imaging deep within the sample.

A fiber-based endoscopic imaging system combining narrowband red-green-blue (RGB) reflectance with optical coherence tomography (OCT) and autofluorescence imaging (AFI) has been developed. The system uses a submillimeter diameter rotary-pullback double-clad fiber imaging catheter for sample illumination and detection. The imaging capabilities of each modality are presented and demonstrated with images of a multicolored card, fingerprints, and tongue mucosa. Broadband imaging, which was done to compare with narrowband sources, revealed better contrast but worse color consistency compared with narrowband RGB reflectance. The measured resolution of the endoscopic system is 25 μm in both the rotary direction and the pullback direction. OCT can be performed simultaneously with either narrowband RGB reflectance imaging or AFI.

Angularly resolved elastic light scattering is an established technique for probing the average size of organelles in biological tissue and cellular ensembles. Focusing of the incident light to illuminate no more than one cell at a time restricts the minimum forward-scattering angle θmin that can be detected. Series of simulated single-cell angular-scattering patterns have been generated to explore how size estimates vary as a function of θmin. At a setting of θmin=20 deg, the size estimates hop unstably between multiple minima in the solution space as simulated noise (mimicking experimentally observed levels) is varied. As θmin is reduced from 20 deg to 10 deg, the instability vanishes, and the variance of estimates near the correct answer also decreases. The simulations thus suggest that robust Mie theory fits to single-cell scattering at 785 nm excitation require measurements down to at least 15 deg. Notably, no such instability was observed at θmin=20 deg for narrow bead distributions. Accurate sizing of traditional calibration beads is, therefore, insufficient proof that an angular-scattering system is capable of robust analysis of single cells. Experimental support for the simulation results is also presented using measurements on cells fixed with formaldehyde.

Intravascular photoacoustics (IV-PA) is an emerging atherosclerosis imaging modality that provides chemical-specific optical information of arterial walls with acoustic depth penetration and resolution. As lipid composition of atherosclerotic plaques is considered to be one of the primary indicators for plaque vulnerability, many IV-PA applications are calibrated so as to target plaque necrotic cores. Based on the mode of optical excitation and the corresponding signal processing technique, IV-PA is categorized into two different modalities. The pulse-based IV-PA has been the universal IV-PA imaging mode with its high peak power and straightforward time-domain signal processing technique. As an alternative, the low power continuous-wave (CW)-based IV-PA has been under intense development as a radar-like frequency-domain signal processing modality. The two state-of-the-art types of IV-PA are reviewed in terms of their physics and imaging capabilities, with major emphasis on frequency-swept CW-based IV-PA that has been recently introduced in the field.

A multi-spectral laser speckle contrast imaging (MS-LSCI) system is proposed using only a single wavelength-swept laser, which provides both highly coherent and multi-spectral outputs to simultaneously generate laser speckle contrast images and multi-spectral images, respectively. Using a laser light swept from 770 to 821 nm at a repetition rate of 5 Hz and a CCD camera of 335 fps, 67 multi-spectral frame images are acquired in 0.76 nm wavebands over 51 nm spectral range. The spectral sub-windowing method of single wavelength-swept laser source is used to solve the lack of spectral information from a few individual light sources, which is a limitation of conventional MS-LSCI systems. In addition to the speckle flow index from the LSCI frames, the multi-spectrally encoded images can generate additional images of spectral absorbance. To further examine the performance of the MS-LSCI system, an in vivo cuff-induced ischemia experiment was conducted to show the real-time imaging of hemodynamic and blood oxygen saturation changes simultaneously over the entire 2.5 cm×4.5 cm field of view.

Optical coherence tomography angiography (OCTA) has been widely used for en face visualization of the microvasculature, but is challenged for real three-dimensional (3-D) topologic imaging due to the "tail" artifacts that appear below large vessels. Further, OCTA is generally incapable of differentiating descending arterioles from ascending venules. We introduce a normalized field autocorrelation function-based OCTA (g1-OCTA), which minimizes the tail artifacts and is capable of distinguishing penetrating arterioles from venules in the 3-D image. g1 ( τ ) is calculated from repeated optical coherence tomography (OCT) acquisitions for each spatial location. The decay amplitude of g1 ( τ ) is retrieved to represent the dynamics for each voxel. To account for the small g1 ( τ ) decay in capillaries where red blood cells are flowing slowly and discontinuously, Intralipid is injected to enhance the OCT signal. We demonstrate that the proposed technique realizes 3-D OCTA with negligible tail projections and the penetrating arteries are readily identified. In addition, compared to regular OCTA, the proposed g1-OCTA largely increased the depth-of-field. This technique provides a more accurate rendering of the vascular 3-D anatomy and has the potential for more quantitative characterization of vascular networks.

The optical attenuation coefficient (OAC) reflects the optical properties of various tissues or tissues of the same type under different physiological conditions. Quantitative measurement of OAC from optical coherence tomography (OCT) signals can provide additional information and can increase the potential for OCT applications. We present an optimized depth-resolved estimation (ODRE) method that derives a precise mapping between the measured OCT signal and the OAC. In contrast to previous depth-resolved estimation (DRE) methods, the optimized method can estimate the OAC in any depth range and ignore whether the light is completely attenuated. Numerical simulations and phantom experiments are used to verify its validity, and this method is applied to detect cerebral damage. In combination with OCT angiography, real-time observation of the change of blood perfusion and the degree of cerebral damage in mice with focal cerebral ischemia provides important information to help us understand the temporal relationship between brain damage and ischemia.

Physiological fluctuations are expected to be a dominant source of noise in blood oxygenation level-dependent (BOLD) magnetic resonance imaging (MRI) experiments to assess tumour oxygenation and angiogenesis. This work investigates the impact of various physiological noise regressors: retrospective image correction (RETROICOR), heart rate (HR) and respiratory volume per unit time (RVT), on signal variance and the detection of BOLD contrast in the breast in response to a modulated respiratory stimulus. BOLD MRI was performed at 3 T in ten volunteers at rest and during cycles of oxygen and carbogen gas breathing. RETROICOR was optimized using F-tests to determine which cardiac and respiratory phase terms accounted for a significant amount of signal variance. A nested regression analysis was performed to assess the effect of RETROICOR, HR and RVT on the model fit residuals, temporal signal-to-noise ratio, and BOLD activation parameters. The optimized RETROICOR model accounted for the largest amount of signal variance (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$ \Delta R_{\text{adj}}^{2}$ \end{document}ΔRadj2 = 3.3 ± 2.1%) and improved the detection of BOLD activation (P = 0.002). Inclusion of HR and RVT regressors explained additional signal variance, but had a negative impact on activation parameter estimation (P

Physics in medicine and biology 61(9), 3500-3526 (2016). doi:10.1088/0031-9155/61/9/3500, Published by IOP Publ., Bristol

Raman spectroscopy has proven to be an indispensable technique for the identification of various types of analytes due to the fingerprint vibration spectrum obtained. Paper microfluidics has also emerged as a low cost, easy to fabricate and portable approach for point of care testing. However, due to inherent background fluorescence, combining Raman spectroscopy with paper microfluidics is to date an unmet challenge in the absence of using surface enhanced mechanisms. We describe the first use of wavelength modulated Raman spectroscopy (WMRS) for analysis on a paper microfluidics platform. This study demonstrates the ability to suppress the background fluorescence of the paper using WMRS and the subsequent implementation of this technique for pharmaceutical analysis. The results of this study demonstrate that it is possible to discriminate between both paracetamol and ibuprofen, whilst, also being able to detect the presence of each analyte quantitatively at nanomolar concentrations. Publisher PDF