Novel needle-probe single-fiber reflectance spectroscopy to quantify sub-surface myoglobin forms in beef psoas major steaks during retail display

Meat discoloration starts at the interface between the bright red oxymyoglobin layer and the interior deoxymyoglobin layer. Currently, limited tools are available to characterize myoglobin forms formed within the sub-surface of meat. The objective was to demonstrate a needle-probe based single-fiber reflectance (SfR) spectroscopy approach for characterizing sub-surface myoglobin forms of beef psoas major muscles during retail storage. A 400-μm fiber was placed in a 17-gauge needle, and the assembly was inserted into the muscle at five depths of 1 mm increment and 1 cm lateral shift. Metmyoglobin content increased at all depths during display and content at 1 mm was greater compared to that of 2 to 5 mm depth. The a* values decreased (P < 0.05) during retail display aligning with the sub-surface formation of metmyoglobin. In summary, the results suggest that needle-probe SfR spectroscopy can determine interior myoglobin forms and characterize meat discoloration.

Diffuse photon remission associated with the center-illuminated-area-detection geometry. II. Approach to the time-domain model

This part proposes a model of time-dependent diffuse photon remission for the center-illuminated-area-detection (CIAD) geometry, by virtue of area integration of the radially resolved time-dependent diffuse photon remission formulated with the master-slave dual-source scheme demonstrated in Part I for steady-state measurements. The time-domain model is assessed against Monte Carlo (MC) simulations limiting to only the Heyney-Greenstein scattering phase function for CIAD of physical scales and medium properties relevant to single-fiber reflectance (SfR) and over a 2 ns duration, in compliance with the timespan of the only experimental report of SfR demonstrated with a 50 µm gradient index fiber. The time-domain model-MC assessments are carried out for an absorption coefficient ranging three orders of magnitude over [0.001,0.01,0.1,1]m m -1 at a fixed scattering, and a reduced scattering coefficient ranging three orders of magnitude over [0.01,0.1,1,10]m m -1 at a fixed absorption, among others. Photons of shorter and longer propagation times, relative to the diameter of the area of collection, respond differently to the scattering and absorption changes. Limited comparisons of MC between CIAD and a top-hat geometry as the idealization of SfR reveal that the time-domain photon remissions of the two geometries differ appreciably in only the early arriving photons.

Diffuse photon remission associated with the center-illuminated-area-detection geometry: Part I, an approach to the steady-state model

Diffuse photon remission associated with the center-illuminated-area-detection (CIAD) geometry has been useful for non-contact sensing and may inform single-fiber reflectance (SfR). This series of work advances model approaches that help enrich the understanding and applicability of the photon remission by CIAD. The general approach is to derive the diffuse photon remission by the area integration of the radially resolved diffuse reflectance while limiting the analysis to a medium exhibiting only the Heyney-Greenstein (HG) scattering phase function. Part I assesses the steady-state photon remission in CIAD over a reduced scattering scaled diameter of μ s ' d _{a r e a} ∈[0.5×10^-3,10³] that covers the range extensively modeled for SfR. The corresponding radially resolved diffuse reflectance is obtained by concatenating an empirical expression for the semi-ballistic region near the point-of-illumination and a formula utilizing a master-slave dual-source scheme over the semi-diffusive to a diffusive regime while being constrained by an extrapolated zero-boundary condition. The terminal algebraic photon remission is examined against Monte Carlo simulations for an absorption coefficient over [0.001,1]m m ^-1, a reduced scattering coefficient over [0.01,1000]m m ^-1, a HG scattering anisotropy factor within [0.5,0.95], and a diameter of the area of collection ranging [50,1000]µm. The algebraic model is also applied to phantom data acquired over a ∼2c m non-contact CIAD configuration and with a 200 µm SfR probe. The model approach will be extended in a subsequent work towards the time-of-flight characteristics of CIAD.

Daily Quantification of Myoglobin Forms on Beef Longissimus Lumborum Steaks Over 7 Days of Display by Near-infrared Diffuse Reflectance Spectroscopy

Near-infrared diffuse reflectance spectroscopy (NIR-DRS) was utilized to develop an algorithm using approximately 18 wavelengths spanning 480 to 650 nm to determine oxymyoglobin (OxyMb), deoxymyoglobin (DeoxyMb), and metmyoglobin (MetMb) contents on beef longissimus lumborum muscles. Daily changes in subsurface myoglobin redox forms were evaluated for 7 d using NIR-DRS and compared with the surface color as assessed by a HunterLab MiniScan spectrophotometer as a reference modality. Both measurements revealed that MetMb increased steadily over the duration of display, showing high correlation (R2 = 0.91) between the 2 methods. Comparatively, whereas NIR-DRS revealed OxyMb to have decreased steadily over the period of display, the HunterLab MiniScan spectrophotometer indicated a much later onset of the apparent decrease of OxyMb, resulting in a moderate correlation (R2 = 0.64) between the 2 methods. No correlation was found between the 2 methods regarding the changes of DeoxyMb over the duration of display. The newly developed NIR-DRS algorithm has potential as an alternative method of color assessment in postrigor skeletal muscle.

Diffuse photon-remission associated with single-fiber geometry may be a simple scaling of that collected over the same area when under centered-illumination

Robust models for single-fiber reflectance (SFR) are relatively complex [Opt. Lett.45, 2078 (2020)OPLEDP0146-959210.1364/OL.385845] due to overlapping of the illumination and collection areas that entails probability weighting of the spatial integration of photon-remission. We demonstrate, via analytical means for limiting cases and Monte Carlo simulation of broader conditions, that diffuse photon-remission collected by single-fiber geometry may be scaled over the center-illuminated photon-remission. We specify for a medium revealing Henyey-Greenstein (HG) scattering anisotropy that the diffuse photon-remission from a sub-diffusive area of a top-hat illumination is ∼84.9% of that collected over the same area when under a centered-illumination. This ratio remains consistent over a reduced-scattering fiber-size product of μs'd_fib=[10^-5,10⁰], for absorption varying 3 orders of magnitude. When applied to hemoglobin oxygenation changes induced in an aqueous phantom using a 200 µm single-fiber probe, the center-illumination-scaled model of SFR produced fitting results agreeing with reference measurements.

Diffuse photon remission from thick opaque media of the high absorption/scattering ratio beyond what is accountable by the Kubelka-Munk function

The Kubelka-Munk (KM) theory of diffuse photon remission from opaque media is widely applied to quality-control processes. Recent works based on radiative transfer revealed that the KM function as the backbone parameter of the method may saturate at strong absorption to cause the KM approach to be unfit to predict the change of diffuse reflectance from the medium at strong absorption. We demonstrate by empirical means based on Monte Carlo results that diffuse photon remission from a strong-absorbing medium depends simply upon the absorption/scattering ratio when evaluated over a large area centered at the point of illumination differing in geometry from those convenient for the KM approach. Our empirical prediction gives ∼11% mean errors of the diffuse photon remission from thick opaque medium having an absorption coefficient ranging 0.001 to up to 1000 times stronger than the reduced-scattering coefficient. A slight modification to the KM function in terms of the absorption weighting and absorption-scattering coupling for use within the KM approach also noticeably improves the prediction of diffuse photon remission from thick opaque medium of strong absorption. Our empirical model and the KM approach using the modified KM function were compared against measurements from a thick opaque medium, of which the absorption coefficient was changed over four orders of magnitude.

Laparoscopic diffuse reflectance spectroscopy of an underlying tubular inclusion: a phantom study

We demonstrate diffuse reflectance spectroscopy (DRS) of a subsurface tubular inclusion by using a fiber probe having a single source-detector pair attached to a laparoscopic bipolar device. A forward model was also developed for DRS sensing of an underlying long absorbing tubular inclusion set in parallel to the tissue surface, normal to the line of sight of the source-detector pair, and equidistant from the source and the detector. The model agreed with measurements performed at 500 nm and using a 10 mm source-detector separation (SDS) on an aqueous tissue phantom embedding a tubing of 2 or 4 mm inner diameter that contained 9.1% to 33.3% red dye at a depth of up to 11.5 mm. When tested on solid phantoms using the 10 mm SDS, a tubular inclusion of $ \ge 3\;{\rm mm}$≥3mm inner diameter containing 0.05% red dye at a background absorption coefficient of $ 0.021\;{\rm mm}^{-1} $0.021mm^-1 caused $ \ge 8\% $≥8% change of the signal at 500 nm versus the baseline when the inclusion was shallower than 5 mm. When assessed on avian muscle tissue having a 4 mm tubular inclusion embedded at an edge depth of 2 mm, DRS with the 10 mm SDS differentiated the following contents of the inclusion: 33.3% red dye (mimicking blood), 33.3% green dye, 33.3% yellow dye (mimicking bile), water (mimicking urine), and air.

Simple analytical total diffuse reflectance over a reduced-scattering-pathlength scaled dimension of [10-5, 10-1] from a medium with HG scattering anisotropy

Model approximation is necessary for reflectance assessment of tissue at sub-diffusive to non-diffusive scale. For tissue probing over a sub-diffusive circular area centered on the point of incidence, we demonstrate simple analytical steady-state total diffuse reflectance from a semi-infinite medium with the Henyey-Greenstein (HG) scattering anisotropy (factor $g$g). Two physical constraints are abided to: (1) the total diffuse reflectance is the integration of the radial diffuse reflectance; (2) the radial and total diffuse reflectance at $g \gt {0}$g>0 analytically must resort to their respective forms corresponding to isotropic scattering as $g$g becomes zero. Steady-state radial diffuse reflectance near the point of incidence from a semi-infinite medium of $g \approx 0$g≈0 is developed based on the radiative transfer for isotropic scattering, then integrated to find the total diffuse reflectance for $g \approx 0$g≈0. The radial diffuse reflectance for $g \ge 0.5$g≥0.5 is semi-empirically formulated by comparing to Monte Carlo simulation results and abiding to the second constraint. Its integration leads to a total diffuse reflectance for $g \ge 0.5$g≥0.5 that is also bounded by the second constraint. Over a collection diameter of the reduced-scattering pathlength ($1/\mu _s^{ \prime}$1/μs') scaled size of [${{10}^{ - 5}}$10^-5, ${{10}^{ - 1}}$10^-1] for $g = [{0.5},{0.95}]$g=[0.5,0.95] and the absorption coefficient as strong as the reduced scattering coefficient, the simple analytical total diffuse reflectance is found to be accurate, with an average error of 16.1%.

Erratum: Magneto-thermo-acoustics from magnetic nanoparticles by short bursting or frequency chirped alternating magnetic field: a theoretical feasibility analysis. Med. Phys. 40(6): p. 063301 (2013)

We correct one typographical error that has occurred in four equations in Med Phys, 2013. 40(6): p. 063301.

Transcutaneous transmission of photobiomodulation light to the spinal canal of dog as measured from cadaver dogs using a multi-channel intra-spinal probe

The target level photobiomodulation (PBM) irradiances along the thoracic to lumbar segment of the interior spinal canal in six cadaver dogs resulting from surface illumination at 980 nm were measured. Following a lateral hemi-laminectomy, a flexible probe fabricated on a plastic tubular substrate of 6.325 mm diameter incorporating nine miniature photodetectors was embedded in the thoracic to lumbar segment of the spinal canal. Intra-spinal irradiances at the nine photodetector sites, spanning an approximate 8 cm length caudal to T13, were measured for various applied powers of continuous wave (CW) surface illumination at 980 nm with a maximal power of 10 W corresponding to a surface irradiance of 3.14 W/cm². The surface illumination conditions differed in skin transmission when the probe was off-contact with tissue and probe-skin contact when the skin was in place. For each condition of surface illumination, the beam was directed to respectively T13 (surface site 1), a spinal column site 4 cm caudal to T13 (surface site 5), and a spinal column site 8 cm caudal to T13 (surface site 9). Off-contact surface irradiation of 3.14 W/cm² at surface sites 1, 5, and 9 transmitted respectively 234.0 ± 120.7 μW/cm², 230.7 ± 178.3 μW/cm², and 130.2 ± 169.6 μW/cm² to the spinal canal without the skin, and respectively 35.7 ± 33.2 μW/cm², 50.9 ± 75.3 μW/cm², and 15.7 ± 16.3 μW/cm² with the skin. Transmission with skin was as low as 12% of the transmission without the skin. On-contact surface irradiation of 3.14 W/cm² at surface sites 1, 5, and 9 transmitted respectively 44.6 ± 43.1 μW/cm², 85.4 ± 139.1 μW/cm², and 22.0 ± 23.6 μW/cm² to the spinal canal. On-contact application increased transmission by a maximum of 67% comparing to off-contact application. The information gathered highlights the need to clinically consider the impact of skin transmission and on-contact application technique when attempting to treat spinal cord disease with PBM.

Simple empirical master-slave dual-source configuration within the diffusion approximation enhances modeling of spatially resolved diffuse reflectance at short-path and with low-scattering from a semi-infinite homogeneous medium: erratum

We correct one typographical error of three equations in Appl. Opt.56, 1447 (2017)APOPAI0003-693510.1364/AO.56.001447.

Simple empirical master-slave dual-source configuration within the diffusion approximation enhances modeling of spatially resolved diffuse reflectance at short-path and with low-scattering from a semi-infinite homogeneous medium: erratum

We correct one typographical error of three equations in Appl. Opt.56, 1447 (2017)APOPAI0003-693510.1364/AO.56.001447.

Identification of critical microRNAs in gastrointestinal stromal tumor patients treated with Imatinib

Gastrointestinal stromal tumors (GISTs) are the most common mesenchymal neoplasms of the gastrointestinal tract. Imatinib mesylate was considered to be a breakthrough drug in clinical treatment of GIST, but GIST patients showed resistance against it. We aimed to identify critical microRNAs (miRNAs) related to imatinib resistance in imatinib-treated GIST patients. Microarray datasets under the accession number of GSE63159 and GSE45901 were downloaded from the Gene Expression Omnibus (GEO) database. The differentially expressed miRNAs (DEMs) that are related to imatinib resistance were identified. GO function and KEGG pathway enrichment analyses were performed, and lncRNA-miRNA-target gene regulatory networks were constructed. Finally, the critical miRNAs and their target genes that are related to imatinib resistance or sensitivity were identified. In total, 20 DEMs in the GSE63159 dataset (7 significantly up-regulated and 13 down-regulated) and 23 DEMs in the GSE45901 dataset (8 up-regulated and 15 down-regulated) were identified. In lncRNA-miRNA-target gene regulatory networks, five critical miRNAs and 109 target genes were identified. GO function and KEGG pathway enrichment analysis showed that the target genes of DEMs were mainly involved in several signaling pathways, such as focal adhesion and the GnRH signaling pathway. Among the five miRNAs, the overexpression of hsa-miR-28-5p and hsa-miR-125a-5p had significant correlation to imatinib resistance or imatinib sensitivity in GIST patients. Hsa-miR-28-5p and hsa-miR-125a-5p may be involved in the development and progression of GIST, and they may be able to serve as prognostic markers for imatinib-response in GIST patients.

Perspective review on solid-organ transplant: needs in point-of-care optical biomarkers

Solid-organ transplant is one of the most complex areas of modern medicine involving surgery. There are challenging opportunities in solid-organ transplant, specifically regarding the deficiencies in pathology workflow or gaps in pathology support, which may await alleviations or even de novo solutions, by means of point-of-care, or point-of-procedure optical biomarkers. Focusing the discussions of pathology workflow on donor liver assessment, we analyze the undermet need for intraoperative, real-time, and nondestructive assessment of the donor injuries (such as fibrosis, steatosis, and necrosis) that are the most significant predictors of post-transplant viability. We also identify an unmet need for real-time and nondestructive characterization of ischemia or irreversible injuries to the donor liver, earlier than appearing on morphological histology examined with light microscopy. Point-of-procedure laparoscopic optical biomarkers of liver injuries and tissue ischemia may also facilitate post-transplant management that is currently difficult for or devoid of pathological consultation due to lack of tools. The potential and pitfalls of point-of-procedure optical biomarkers for liver assessment are exemplified in breadth for steatosis. The more general and overarching challenges of point-of-procedure optical biomarkers for liver transplant pathology, including the shielding effect of the liver capsule that was quantitated only recently, are projected. The technological and presentational benchmarks that a candidate technology of point-of-procedure optical biomarkers for transplant pathology must demonstrate to motivate clinical translation are also foreseen.

Effects of capsule on surface diffuse reflectance spectroscopy of the subcapsular parenchyma of a solid organ

We hypothesize that the capsular optical properties and thickness combined affect how accurate the diffuse reflectance on the surface of a capsular solid organ represents that on the subcapsular parenchyma. Monte Carlo simulations on two-layer geometries evaluated how a thin superficial layer with the thickness from 10 to 1000  μm affected the surface diffuse reflectance over a source-detector separation spanning 0.01 to 10 mm. The simulations represented the superficial layer presenting various contrasts concerning refractive index, anisotropy factor, absorption coefficient, and reduced scattering coefficient, versus those of the subsurface main medium. An analytical approach modeled the effects of the superficial layer of various thicknesses and optical properties on diffuse reflectance. Diffuse reflectance spectroscopy was performed ex vivo on 10 fresh human livers and 9 fresh human kidneys using a surface probe with a 3-mm source-detector separation. The difference of the device-specific diffuse reflectance on the organ between with the capsule and without the capsule has significantly greater spectral variation in the kidney than in the liver. The significantly greater spectral deviation of surface diffuse reflectance between with and without the capsule in the kidney than in the liver was analytically accountable by considering the much thicker capsule of the kidney than of the liver.

Flexible nine-channel photodetector probe facilitated intraspinal multisite transcutaneous photobiomodulation therapy dosimetry in cadaver dogs

Noninvasive photobiomodulation therapy (PBMT) of spinal cord disease remains speculative due to the lack of evidence for whether photobiomodulatory irradiances can be transcutaneously delivered to the spinal cord under a clinically acceptable PBMT surface irradiation protocol. We developed a flexible nine-channel photodetection probe for deployment within the spinal canal of a cadaver dog after hemilaminectomy to measure transcutaneously transmitted PBMT irradiance at nine sites over an eight-cm spinal canal length. The probe was built upon a 6.325-mm tubular stem, to the surface of which nine photodiodes were epoxied at approximately 1 cm apart. The photodiode has a form factor of 4.80  mm×2.10  mm×1.15  mm (length×width×height). Each photodiode was individually calibrated to deliver 1 V per 7.58  μW/cm2 continuous irradiance at 850 nm. The outputs of eight photodiodes were logged concurrently using a data acquisition module interfacing eight channels of differential analog signals, while the output of the ninth photodiode was measured by a precision multimeter. This flexible probe rendered simultaneous intraspinal (nine-site) measurements of transcutaneous PBMT irradiations at 980 nm in a pilot cadaver dog model. At a surface continuous irradiance of 3.14  W/cm2 applied off-contact between L1 and L2, intraspinal irradiances picked up by nine photodiodes had a maximum of 327.48  μW/cm2 without the skin and 5.68  μW/cm2 with the skin.

Freehand diffuse optical spectroscopy imaging for intraoperative identification of major venous and arterial vessels underlying peritoneal fat: an in vivo demonstration in a pig model

Inadvertent injury to important anatomic structures is a significant risk in minimally invasive surgery (MIS) that potentially requires conversion to an open procedure, which results in increased morbidity and mortality. Surgeons operating minimal-invasively currently do not have an easy-to-use, real-time device to aid in intraoperative identification of important anatomic structures that underlie tissue planes. We demonstrate freehand diffuse optical spectroscopy (DOS) imaging for intraoperatively identifying major underlying veins and arteries. An applicator probe that can be affixed to and detached from an 8-mm laparoscopic instrument has been developed. The 10-mm DOS source-detector separation renders sampling of tissue heterogeneities a few millimeters deep. DOS spectra acquired consecutively during freehand movement of the applicator probe on the tissue surface are displayed as a temporal and spectral image to assist in spatially resolved identification of the underlying structures. Open surgery identifications of the vena cava and aorta underlying peritoneal fat of ∼4  mm in thickness using the applicator probe under room light were demonstrated repeatedly in multiple pigs in vivo.

In vivo assessment of diet-induced rat hepatic steatosis development by percutaneous single-fiber spectroscopy detects scattering spectral changes due to fatty infiltration

This study explores percutaneous single-fiber spectroscopy (SfS) of rat livers undergoing fatty infiltration. Eight test rats were fed a methionine-choline-deficient (MCD) diet, and four control rats were fed a normal diet. Two test rats and one control rat were euthanized on days 12, 28, 49, and 77 following initiation of the diet, after percutaneous SfS of the liver under transabdominal ultrasound guidance. Histology of each set of the two euthanized test rats showed mild and mild hepatic lipid accumulations on day 12, moderate and severe on day 28, severe and mild on day 49, and moderate and mild on day 77. Livers with moderate or higher lipid accumulation generally presented higher spectral reflectance intensity when compared to lean livers. Livers of the eight test rats on day 12, two of which had mild lipid accumulation, revealed an average scattering power of 0.37±0.14 in comparison to 0.07±0.14 for the four control rats (p<0.01 ). When livers of the test rats with various levels of fatty infiltration were combined, the average scattering power was 0.36±0.15 0.36±0.15 in comparison to 0.14±0.24 of the control rats (0.05<p<0.1). Increasing lipid accumulation in concentration and size seemed to cause an increase of the scattering power prior to increasing total spectral reflectance.

On the geometry dependence of differential pathlength factor for near-infrared spectroscopy. I. Steady-state with homogeneous medium

This work analytically examines some dependences of the differential pathlength factor (DPF) for steady-state photon diffusion in a homogeneous medium on the shape, dimension, and absorption and reduced scattering coefficients of the medium. The medium geometries considered include a semi-infinite geometry, an infinite-length cylinder evaluated along the azimuthal direction, and a sphere. Steady-state photon fluence rate in the cylinder and sphere geometries is represented by a form involving the physical source, its image with respect to the associated extrapolated half-plane, and a radius-dependent term, leading to simplified formula for estimating the DPFs. With the source-detector distance and medium optical properties held fixed across all three geometries, and equal radii for the cylinder and sphere, the DPF is the greatest in the semi-infinite and the smallest in the sphere geometry. When compared to the results from finite-element method, the DPFs analytically estimated for 10 to 25 mm source–detector separations on a sphere of 50 mm radius with μa=0.01  mm(−1) and μ′s=1.0  mm(−1) are on average less than 5% different. The approximation for sphere, generally valid for a diameter≥20 times of the effective attenuation pathlength, may be useful for rapid estimation of DPFs in near-infrared spectroscopy of an infant head and for short source–detector separation.

Photon diffusion in a homogeneous medium bounded externally or internally by an infinitely long circular cylindrical applicator. VI. Time-domain analysis

Part VI analytically examines time-domain (TD) photon diffusion in a homogeneous medium enclosed by a "concave" circular cylindrical applicator or enclosing a "convex" circular cylindrical applicator, both geometries being infinite in the longitudinal dimension. The aim is to assess characteristics of TD photon diffusion, in response to a spatially and temporally impulsive source, versus the line-of-sight source-detector distance along the azimuthal or longitudinal direction on the concave or convex medium-applicator interface. By comparing to their counterparts evaluated along a straight line on a semi-infinite medium-applicator interface versus the same source-detector distance, the following patterns are indicated: (1) the peak photon fluence rate is always reached sooner in concave and later in convex geometry; (2) the peak photon fluence rate decreases slower along the azimuthal and faster along the longitudinal direction on the concave interface, and conversely on the convex interface; (3) the total photon fluence decreases slower along the azimuthal and faster along the longitudinal direction on the concave interface, and conversely on the convex interface; (4) the ratio between the peak photon fluence rate and the total fluence is always greater in concave geometry and smaller in convex geometry. The total fluence is equivalent to the steady-state photon fluence analyzed in Part I [J. Opt. Soc. Am. A27, 648 (2010)10.1364/JOSAA.27.000648JOAOD61084-7529]. The patterns of peak fluence rate, time to reaching peak fluence rate, and the ratio of these two, correspond to those of AC amplitude, phase, and modulation depth of frequency-domain results demonstrated in Part IV [J. Opt. Soc. Am. A29, 1445 (2012)10.1364/JOSAA.29.001445JOAOD61084-7529].

Percutaneous single-fiber reflectance spectroscopy of canine intervertebral disc: is there a potential for in situ probing of mineral degeneration?

BACKGROUND AND OBJECTIVES

Intervertebral disc herniation is a common disease in chondrodystrophic dogs, and a similar neurologic condition also occurs in humans. Percutaneous laser disc ablation (PLDA) is a minimally invasive procedure used increasingly for prevention of disc herniation. Currently, PLDA is performed on thoracolumbar discs with the same laser energy applied regardless of the differing extent of degeneration among mineralized discs. In a previous study performed on 15 normal and 6 degenerated intervertebral discs in chondrodystrophoid canine species, it was demonstrated that percutaneous single-fiber reflectance spectroscopy (SfRS) detected increased light scattering from mineralized intervertebral discs when comparing to normal discs. The objective of this study is to evaluate how SfRS evaluation of mineralized discs in situ fairs with X-ray radiography and computed tomography (CT) diagnoses and if SfRS sensing of the scattering changes correlates with the level of mineral degeneration in nucleus pulposus.

MATERIALS AND METHODS

Percutaneous SfRS was performed on a total of 28 intervertebral discs of three dogs post-mortem, through a 20 gauge spinal needle standard to PLDA. The raw SfRS measurement was normalized to extract a dimension-less spectral intensity profile, from which the average over 600-900 nm was used as the SfRS intensity index to compare among the measured discs. The discs were imaged prior to percutaneous SfRS by radiography and CT, and harvested after percutaneous SfRS for histopathologic examinations.

RESULTS

Five among 10 discs of dog #1, six among 9 discs of dog #2, and nine out of 9 discs of dog #3 were determined by histopathology to have central focal or multi-focal areas of mineralization occupying 5-75% of the examined area of nucleus pulposus. The overall numbers of discs with detectable and undetectable central mineralization were 20 and 8, respectively. CT resulted in one false positive (FP) and four false negative (FN) diagnoses for dog #1, three FP and zero FN diagnoses for dog #2, and zero FP and one FN diagnosis for dog #3. Of the total 28 discs the CT had an overall positive predictive value (PPV) of 78.8% and an overall negative predictive value (NPV) of 44.4%. X-ray radiography gave five FN diagnoses for dog #1, two FN diagnoses for dog #2, and eight FN diagnoses for dog #3. Of the total 28 discs the radiography had an overall PPV of 100% and an overall NPV of 30.4%. The receiver-operating-characteristic analysis of the SfRS measurement was performed on 24 discs that had a central mineralization not greater than 50%. An area-under-curve of 0.6758 infers that the SfRS intensity weakly indicates the level of mineralization.

CONCLUSIONS

Percutaneous SfRS may be useful as an in situ sensing tool for assessing the level of mineral degeneration in intervertebral discs for the prospect of disc-specific dosage adjustment in PLDA.

Magnetothermoacoustics from magnetic nanoparticles by short bursting or frequency chirped alternating magnetic field: a theoretical feasibility analysis

PURPOSE

To propose an alternative method of thermoacoustic wave generation based on heating of magnetic nanoparticles (MNPs) using alternating magnetic field (AMF).

METHODS

The feasibility of thermoacoustic wave generation from MNPs by applying a short-burst of AMF or a frequency-modulated AMF is theoretically analyzed. As the relaxation of MNPs is strongly dependent upon the amplitude and frequency of AMF, either an amplitude modulated, fixed frequency AMF (termed time-domain AMF) or a frequency modulated, constant amplitude AMF (termed frequency-domain AMF) will result in time-varying heat dissipation from MNPs, which has the potential to generate thermoacoustic waves. Following Rosensweig's model of specific power loss of MNPs in a steady-state AMF, the time-resolved heat dissipations of MNPs of superparamagnetic size when exposed to a short bursting of AMF and∕or to a linearly frequency chirped AMF are derived, and the resulted acoustic propagation is presented. Based on experimentally measured temperature-rise characteristics of a superparamagnetic iron-oxide nanoparticle (SPION) matrix in a steady-state AMF of various frequencies, the heat dissipations of the SPION under time-domain and frequency-domain AMF configurations that could have practical utility for thermoacoustic wave generation are estimated.

RESULTS

The initial rates of the temperature-rise of the SPION matrix were measured at an iron-weight concentration of 0.8 mg∕ml and an AMF frequency of 88.8 kHz to 1.105 MHz. The measured initial rates of temperature-rise were modeled by Rosensweig's theory, and projected to 10 MHz AMF frequency, at which a 1 μs bursting corresponding to a 1.55 mm axial resolution of acoustic detection could contain 10 complete cycles of AMF oscillation and the power dissipation is approximately 84 times of that at 1 MHz. Exposing the SPION matrix to a 1 μs bursting of AMF at 10 MHz frequency and 100 Oe field intensity would produce a volumetric heat dissipation of 7.7 μJ∕cm(3) over the microsecond duration of the AMF burst. If the SPION matrix is exposed to a 1 ms long AMF train at 100 Oe field intensity that chirps linearly from 1 to 10 MHz, the volumetric heat dissipation produced over each 2π phase change of the AMF oscillation is estimated to increase from 0.15 to 1.1 μJ∕cm(3) within the millisecond duration of the chirping of AMF.

CONCLUSIONS

The heat dissipations upon SPION (∼1 mg∕ml iron-weight concentration) by a 1 μs bursting of 100 Oe AMF at 10 MHz and a 1 ms train of 100 Oe AMF that chirps linearly from 1 to 10 MHz were estimated to determine the potential of thermal-acoustic wave generation. Although thermoacoustic wave generation from MNPs by time- or frequency-domain AMF applications is predicted, the experimental generation of such a wave remains challenging.

A geometric-sensitivity-difference method improves object depth-localization for continuous-wave fluorescence diffuse optical tomography: an in silico study in an axial outward-imaging geometry

A geometric-sensitivity-difference (GSD) based reconstruction method is demonstrated in fluorescence diffuse optical tomography (FDOT) for improving the depth-localization of objects. The GSD method optimizes the data-model fit based on paired-measurements between source-detector pairs sharing either the source or the detector channel, as comparing to conventional methods that optimize the data-model fit based on un-paired measurements of individual source-detector pairs. This in silico study is limited to continuous-wave and 2-dimension, for a circular-array outward-imaging geometry of which the native sensitivity of measurement varies strongly with respect to the depth of the object. The outcomes of GSD method are compared to that of two conventional methods: one is the baseline method which does not involve any scheme to compensate the variation of native sensitivity; the other applies a depth-adapted weight to counteract the depth-variance of the native sensitivity. These three methods were evaluated using synthetic data corresponding to the following conditions of the object: (1) Single object with a 3-folds of positive contrast of fluorescence over the background was set at edge-depths of 0.5 mm, 5 mm and 10 mm; (2) Two objects with identical 3-folds of positive or 1/3-folds of negative contrast of fluorescence over the background were set at a fixed edge-depth of 10 mm and different azimuthal separations of 45 degree, 135 degree, and 180 degree; (3) Two objects with identical 3-folds of positive or 1/3-folds of negative contrast of fluorescence over the background were set at a fixed azimuthal separation of 90° and at edge-depths of 0.5 mm, 5 mm and 10 mm. The GSD method outperforms the other two methods in localizing a single anomaly and resolving two anomalies, for the anomaly possessing either the 3 folds positive or 1/3-folds negative contrast of fluorescence over the background. The case of objects with negative contrast over the background has specific implications to imaging zinc-specific fluorophore uptake in prostate.

Diffuse optical tomography to monitor the photocoagulation front during interstitial photothermal therapy: Numerical simulations and measurements in tissue-simulating phantoms

Near-infrared interstitial photothermal therapy (PTT) is currently undergoing clinical trials as an alternative to watchful waiting or radical treatments in patients with low/intermediate-risk focal prostate cancer. Currently, magnetic resonance imaging (MRI)-based thermography is used to monitor thermal energy delivery and determine indirectly the completeness of the target tumor destruction while avoiding damage to adjacent normal tissues, particularly the rectal wall. As an alternative, transrectal diffuse optical tomography (TRDOT) is being developed to image directly the photocoagulation boundary based on the changes in tissue optical properties, particularly scattering. An established diffusion-theory finite-element software platform was used to perform forward simulations to determine the sensitivity of changes in the optical signal resulting from a growing coagulated lesion with optical scattering contrast, for varying light source-detector separations in both longitudinal and transverse imaging geometries. The simulations were validated experimentally in tissue-simulating phantoms using an existing continuous-wave TRDOT system, in a configuration that is representative of one potential intended clinical use. This provides critical guidance for the optimum design of the transrectal applicator probe, in terms of achieving maximum sensitivity to the presence of the coagulation boundary and, consequently, the highest accuracy in determining the boundary location relative to the rectal wall.

Toward transduodenal diffuse optical tomography of proximal pancreas

We demonstrate the feasibility of diffuse optical tomography (DOT) of the proximal pancreas by using optical applicator channels deployed longitudinally along the exterior surface of a duodenoscope. As the duodenum that nearly encircles the proximal pancreas forms a natural "C-loop" that is approximately three-quarters of a circle of 5-6 cm in diameter, a multichannel optical applicator attached to a duodenoscope has the potential to perform transduodenal DOT sampling of the bulk proximal pancreas wherein most cancers and many cystic lesions occur. The feasibility of transduodenal DOT is demonstrated on normal porcine pancreas tissues containing an introduced gelatinous inclusion of approximately 3 cm in diameter, by using nine source channels and six detector channels attached to a duodenoscope. Concurrent ultrasonography of the gelatinous inclusion in the porcine pancreas parenchyma provided a coarse, albeit indispensable, anatomic prior to transduodenal DOT in reconstructing a contrast of optical properties in the pancreas.

Photon diffusion in a homogeneous medium bounded externally or internally by an infinitely long circular cylindrical applicator. V. Steady-state fluorescence

As Part V in our series, this paper examines steady-state fluorescence photon diffusion in a homogenous medium that contains a homogenous distribution of fluorophores, and is enclosed by a "concave" circular cylindrical applicator or is enclosing a "convex" circular cylindrical applicator, both geometries being infinite in the longitudinal dimension. The aim is to predict by analytics and examine with the finite-element method the changing characteristics of the fluorescence-wavelength photon-fluence rate and the ratio (sometimes called the Born ratio) of it versus the excitation-wavelength photon-fluence rate, with respect to the source-detector distance. The analysis is performed for a source and a detector located on the medium-applicator interface and aligned either azimuthally or longitudinally in both concave and convex geometries. When compared to its steady-state counterparts on a semi-infinite medium-applicator interface with the same line-of-sight source-detector distance, the fluorescence-wavelength photon-fluence rate reduces faster along the longitudinal direction and slower along the azimuthal direction in the concave geometry, and conversely in the convex geometry. However, the Born ratio increases slower in both azimuthal and longitudinal directions in the concave geometry and faster in both directions in the convex geometry, respectively, when compared to that in the semi-infinite geometry.

Photon diffusion in a homogeneous medium bounded externally or internally by an infinitely long circular cylindrical applicator. IV. Frequency-domain analysis

Part IV examines frequency-domain photon diffusion in a homogeneous medium enclosed by a "concave" circular cylindrical applicator or enclosing a "convex" circular cylindrical applicator, both geometries being infinite in the longitudinal dimension. The aim is to assess by analogical and finite-element methods the changes of AC amplitude, modulation depth, and phase with respect to the line-of-sight source-detector distance for a source and a detector located along the azimuthal or longitudinal direction on the concave or convex medium-applicator interface. By comparing to their counterparts along a straight line on a semi-infinite medium-applicator interface, for the same line-of-sight source-detector distance, it is found that: (1) the decay-rate of AC photon fluence is smaller along the azimuthal direction and greater along the longitudinal direction on the concave interface, (2) the decay-rate of AC photon fluence is greater along the azimuthal direction and smaller along the longitudinal direction on the convex interface, (3) the modulation depth along both azimuthal and longitudinal directions decays more slowly on the concave interface and faster on the convex interface, and (4) the phase along both azimuthal and longitudinal directions increases more slowly on the concave interface and faster on the convex interface.

Photon diffusion in a homogeneous medium bounded externally or internally by an infinitely long circular cylindrical applicator. III. Synthetic study of continuous-wave photon fluence rate along unique spiral paths

This is Part III of the work that examines photon diffusion in a scattering-dominant medium enclosed by a "concave" circular cylindrical applicator or enclosing a "convex" circular cylindrical applicator. In Part II of this work Zhang et al. [J. Opt. Soc. Am. A, 66 (2011)] predicted that, on the tissue-applicator interface of either "concave" or "convex" geometry, there exists a unique set of spiral paths, along which the steady-state photon fluence rate decays at a rate equal to that along a straight line on a planar semi-infinite interface, for the same line-of-sight source-detector distance. This phenomenon of steady-state photon diffusion is referred to as "straight-line-resembling-spiral paths" (abbreviated as "spiral paths"). This Part III study develops analytic approaches to the spiral paths associated with geometry of a large radial dimension and presents spiral paths found numerically for geometry of a small radial dimension. This Part III study also examines whether the spiral paths associated with a homogeneous medium are a good approximation for the medium containing heterogeneity. The heterogeneity is limited to an anomaly that is aligned azimuthally with the spiral paths and has either positive or negative contrast of the absorption or scattering coefficient over the background medium. For a weak-contrast anomaly the perturbation by it to the photon fluence rate along the spiral paths is found by applying a well-established perturbation analysis in cylindrical coordinates. For a strong-contrast anomaly the change by it to the photon fluence rate along the spiral paths is computed using the finite-element method. For the investigated heterogeneous-medium cases the photon fluence rate along the homogeneous-medium associated spiral paths is macroscopically indistinguishable from, and microscopically close to, that along a straight line on a planar semi-infinite interface.

A low-cost needle-based single-fiber reflectance spectroscopy method to probe scattering changes associated with mineralization in intervertebral discs in chondrodystrophoid canine species – A pilot study

Intervertebral disc herniation is a common disease in chondrodystrophic dogs, and a similar neurologic condition also occurs in humans. Percutaneous laser disc ablation (PLDA) is a minimally invasive procedure used increasingly for prevention of disc herniation. PLDA is performed on thoracolumbar discs to which the same laser energy is applied regardless of their mineral content. Knowledge of individual disc mineral composition would allow laser energy dosage adjustments and more accurate treatment of degenerative discs. Usually, PLDA is guided by radiography/fluoroscopy, which has a limited sensitivity of approximately 60% for identification of mineralized discs. An imaging or sensing technology that provides a more accurate pre-operative

A pilot study was performed on a total of 21 intervertebral discs from two cadaveric dogs (“Dog A” and “Dog B”). The discs were imaged by computed tomography (CT), radiography, and SFR spectroscopy, before histopathologic examination. SFR spectroscopy in the visible/near-infrared band was performed on the nucleus pulposus of the intervertebral disc through a 20-gauge spinal needle placed percutaneously for PLDA. A normalization method was applied to the raw remission spectra to extract a dimension-less and wavelength-dependent intensity profile in the 500–950 nm spectral range.

In total, six discs were determined to be degenerative on histopathology, five discs of “Dog A” and one disc of “Dog B”. CT diagnosed all six degenerated discs, whereas radiography missed two of the five degenerated discs of “Dog A”. The wavelength-dependent mean scattering intensity profiles of the six degenerated discs were noticeably higher than the mean scattering intensity profiles of the 15 “normal” or insignificantly mineralized discs over the entire spectral range. The mean scattering intensities, averaged over each of the entire profiles, were 2.79±0.58 (mean±SD) for the six degenerated discs and 1.48±0.37 for the 15 “normal” or insignificantly mineralized discs. A two-sample

SFR spectroscopy measurements indicate that the increase of light scattering intensity across the entire 500–950 nm spectral range is associated with the mineralization in canine intervertebral discs. However, the scattering characteristics of the nucleus pulposus measured in this study may not necessarily represent the optical properties of the nucleus pulposus at the laser wavelength used for PLDA (2100 nm). More studies on cadaveric and eventually

Transrectal ultrasound-integrated spectral optical tomography of hypoxic progression of a regressing tumor in a canine prostate

The objective of this study was to evaluate if transrectal optical tomography implemented at three wavelength bands for spectral detection could monitor changes of the hemoglobin oxygen saturation (StO2) in addition to those of the total hemoglobin concentration ([HbT]) in lesions of a canine prostate, including an induced tumor modeling canine prostate cancer. Near-infrared (NIR) optical tomography was integrated with ultrasound (US) for transrectal imaging. Multi-spectral detection at 705_nm, 785_nm and 808_nm rendered measurements of [HbT] and StO2. Canine transmissible venereal tumor (TVT) cells were injected into the right lobe of a dog's prostate gland, which had a pre-existing cyst in the left lobe. Longitudinal assessments of the prostate were performed weekly over a 63-day duration by NIR imaging concurrent with grey-scale and Doppler US. Ultrasonography revealed a bi-lobular tumor-mass regressing from day-49 to day-63. At day-49 this tumor-mass developed a hypoxic core that became larger and more intense by day-56 and expanded further by day-63. The tumor-mass presented a strong hyper-[HbT] feature on day-56 that was inconsistent with US-visualized blood flow. Histology confirmed two necrotic TVT foci within this tumor-mass. The cyst appeared to have a large anoxic-like interior that was greater in size than its ultrasonographically delineated lesion, and a weak lesional elevation of [HbT]. On day-56, the cyst presented a strong hyper-[HbT] feature consistent with US-resolved blood flow. Histology revealed acute and chronic hemorrhage in the periphery of the cyst. The NIR imaging features of two other TVT nodules and a metastatic lymph node were evaluated retrospectively. Transrectal US-integrated spectral optical tomography seems to enable longitudinal monitoring of intra-lesional oxygenation dynamics in addition to the hemoglobin content of lesions in the canine prostate.

The utility of direct-current as compared to frequency domain measurements in spectrally-constrained diffuse optical tomography toward cancer imaging

This work investigates, by means of analytical and simulation studies, the performance of spectrally-constrained image reconstruction in Continuous-Wave or Direct-Current (DC) and Frequency Domain (FD) near-infrared optical tomography. A recent analytic approach for estimating the accuracy of target recovery and the level of background artifact for optical tomography at single wavelength, based on the analysis of parametric reconstruction uncertainty level (PRUL), is extended to spectrally-constrained optical tomography. The analytical model is implemented to rank three sets of wavelengths that had been used as spectral prior in an independent experimental study. Subsequent simulation appraises the recovery of oxygenated hemoglobin (HbO), deoxygenated hemoglobin (Hb), water (H2O), scattering amplitude (A), and scattering power (b) using DC-only, DC-excluded FD, and DC-included FD, based on the three sets of wavelengths as the spectral prior. The simulation results support the analytic ranking of the performance of the three sets of spectral priors, and generally agree with the performance outcome of DC-only versus that of DC-excluded FD and DC-included FD. Specifically, this study indicate that: 1) the rank of overall quality of chromophore recovery is Hb, H2O, and HbO from the highest to lowest; and in the scattering part the A is always better recovered than b. This outcome does suggest that the DC-only information gives rise to unique solution to the image reconstruction routine under the given spectral prior. 2) DC-information is not-redundant in FD-reconstruction, as the artifact levels of DC-included FD reconstruction are always lower than those of DC-excluded FD. 3) The artifact level as represented by the noise-to-contrast-ratio is almost always the lowest in DC-only, leading to generally better resolution of multiple targets of identical contrasts over the background than in FD. However, the FD could outperform DC in the recovery of scattering properties including both A and b when the spectral prior is less optimal, implying the benefit of phase-information in scattering recovery in the context of spectrally-constrained optical tomography.

Diffuse photon remission along unique spiral paths on a cylindrical interface is modeled by photon remission along a straight line on a semi-infinite interface

We demonstrate that, for a long cylindrical applicator that interfaces concavely or convexly with a scattering-dominant medium, a unique set of spiral-shaped directions exist on the tissue-applicator interface, along which the diffuse photon remission is essentially modeled by the photon remission along a straight line on a semi-infinite interface. This interesting phenomenon, which is validated in steady state in this work by finite-element and Monte Carlo methods, may be particularly useful for simplifying deeper-tissue sensing in endoscopic imaging geometry.

Photon diffusion in a homogeneous medium bounded externally or internally by an infinitely long circular cylindrical applicator. II. Quantitative examinations of the steady-state theory

This is Part II of the work that examines photon diffusion in a homogenous medium enclosed by a concave circular cylindrical applicator or enclosing a convex circular cylindrical applicator. Part I of this work [J. Opt. Soc. Am. A 27, 648 (2010)] analytically examined the steady-state photon diffusion between a source and a detector for two specific cases: (1) the detector is placed only azimuthally with respect to the source, and (2) the detector is placed only longitudinally with respect to the source, in the infinitely long concave and convex applicator geometries. For the first case, it was predicted that the decay rate of photon fluence would become smaller in the concave geometry and greater in the convex geometry than that in the semi-infinite geometry for the same source-detector distance. For the second case, it was projected that the decay rate of photon fluence would be greater in the concave geometry and smaller in the convex geometry than that in the semi-infinite geometry for the same source-detector distance. This Part II of the work quantitatively examines these predictions from Part I through several approaches, including (a) the finite-element method, (b) the Monte Carlo simulation, and (c) experimental measurement. Despite that the quantitative examinations have to be conducted for finite cylinder applicators with large length-to-radius ratio to approximate the infinite-length condition modeled in Part I, the results obtained by these quantitative methods for two concave and three convex applicator dimensions validated the qualitative trend predicted by Part I and verified the quantitative accuracy of the analytic treatment of Part I in the diffusion regime of the measurement, at a given set of absorption and reduced scattering coefficients of the medium.

Direct-current-based image reconstruction versus direct-current included or excluded frequency-domain reconstruction in diffuse optical tomography

We study the level of image artifacts in optical tomography associated with measurement uncertainty under three reconstruction configurations, namely, by using only direct-current (DC), DC-excluded frequency-domain, and DC-included frequency-domain data. Analytic and synthetic studies demonstrate that, at the same level of measurement uncertainty typical to optical tomography, the ratio of the standard deviation of mu(a) over mu(a) reconstructed by DC only is at least 1.4 times lower than that by frequency-domain methods. The ratio of standard deviations of D (or mu(s)') over D (or mu(s)') reconstructed by DC only are slightly lower than those by frequency-domain methods. Frequency-domain reconstruction including DC generally outperforms that excluding DC, but as the amount of measurements increases, the difference between the two diminishes. Under the condition of a priori structural information, the performances of three reconstruction configurations are seemingly equivalent.

Photon diffusion in a homogeneous medium bounded externally or internally by an infinitely long circular cylindrical applicator. I. Steady-state theory

This work presents an analytic treatment for photon diffusion in a homogeneous medium bounded externally or internally by an infinitely long circular cylindrical applicator. Focusing initially on the steady-state condition, the photon diffusion in these two geometries is solved in cylindrical coordinates by using modified Bessel functions and by applying the extrapolated boundary condition. For large cylinder diameter, the analytic solutions may be simplified to a format employing the physical source and its image source with respect to a semi-infinite geometry and a radius-dependent term to account for the shape and dimension of the cylinder. The analytic solutions and their approximations are evaluated numerically to demonstrate qualitatively the effect of the applicator curvature--either concave or convex--and the radius on the photon fluence rate as a function of the source-detector distance, in comparison with that in the semi-infinite geometry. This work is subjected to quantitative examination in a coming second part and possible extension to time-resolved analysis.

Singular value decomposition based computationally efficient algorithm for rapid dynamic near-infrared diffuse optical tomography

PURPOSE

A computationally efficient algorithm (linear iterative type) based on singular value decomposition (SVD) of the Jacobian has been developed that can be used in rapid dynamic near-infrared (NIR) diffuse optical tomography.

METHODS

Numerical and experimental studies have been conducted to prove the computational efficacy of this SVD-based algorithm over conventional optical image reconstruction algorithms.

RESULTS

These studies indicate that the performance of linear iterative algorithms in terms of contrast recovery (quantitation of optical images) is better compared to nonlinear iterative (conventional) algorithms, provided the initial guess is close to the actual solution. The nonlinear algorithms can provide better quality images compared to the linear iterative type algorithms. Moreover, the analytical and numerical equivalence of the SVD-based algorithm to linear iterative algorithms was also established as a part of this work. It is also demonstrated that the SVD-based image reconstruction typically requires O(NN2) operations per iteration, as contrasted with linear and nonlinear iterative methods that, respectively, require O(NN3) and O(NN6) operations, with "NN" being the number of unknown parameters in the optical image reconstruction procedure.

CONCLUSIONS

This SVD-based computationally efficient algorithm can make the integration of image reconstruction procedure with the data acquisition feasible, in turn making the rapid dynamic NIR tomography viable in the clinic to continuously monitor hemodynamic changes in the tissue pathophysiology.

In vivo trans-rectal ultrasound-coupled optical tomography of a transmissible venereal tumor model in the canine pelvic canal

In vivo trans-rectal near-infrared (NIR) optical tomography was performed concurrently with, albeit reconstructed without spatial a prior of, trans-rectal ultrasound (US) on transmissible venereal tumor (TVT) developed as a model in the canine pelvic canal. Studies were taken longitudinally at prior to, 14 days after, and 35 days after the TVT injection. As the tumor grew, the nodules became increasingly hyperabsorptive and moderately hyperscattering on NIR. The regions of strong NIR contrast, especially on absorption images, correlated well with those of US hypoechoic masses indicative of tumors. Combining the information of trans-rectal NIR and US detected the tumor more accurately than did the US alone at 14 days postinjection.

Trans-rectal ultrasound-coupled near-infrared optical tomography of the prostate, part II: experimental demonstration

We demonstrate trans-rectal optical tomography of the prostate using an endo-rectal near-infrared (NIR) applicator integrated with a transrectal ultrasound (TRUS) probe. The endo-rectal NIR applicator incorporated a design presented in our previously reported work. A continuous-wave NIR optical tomography system is combined with a commercial US scanner to form the dual-modality imager. Sagittal transrectal imaging is performed concurrently by endo-rectal NIR and TRUS. The TRUS ensures accurate positioning of the NIR applicator as well as guides NIR image reconstruction using the spatial prior of the target. The use of a condom, which is standard for TRUS, is found to have minimal effect on trans-rectal NIR imaging. Tests on avian tissues validates that NIR imaging can recover the absorption contrast of a target, and its accuracy is improved when the TRUS spatial prior is incorporated. Trans-rectal NIR/US imaging of a healthy canine prostate in situ is reported.

Trans-rectal ultrasound-coupled near-infrared optical tomography of the prostate, part I: simulation

We investigate the feasibility of trans-rectal optical tomography of the prostate using an endo-rectal near-infrared (NIR) applicator that is to be integrated with a trans-rectal ultrasound (TRUS) probe. Integration with TRUS ensures accurate endo-rectal positioning of the NIR applicator and the utility of using TRUS spatial prior information to guide NIR image reconstruction. The prostate NIR image reconstruction is challenging even with the use of spatial prior owing to the anatomic complexity of the imaging domain. A hierarchical reconstruction algorithm is developed that implements cascaded initial-guesses for nested domains. This hierarchical image reconstruction method is then applied to evaluating a number of NIR applicator designs for integration with a sagittal TRUS transducer. A NIR applicator configuration feasible for instrumentation development is proposed that contains one linear array of optodes on each lateral side of the sagittal TRUS transducer. The performance of this NIR applicator is characterized for the recovery of single tumor mimicking lesion as well as dual targets in the prostate. The results suggest a strong feasibility of transrectal prostate imaging by use of the endo-rectal NIR/US probe.

Minimization of geometric-beam broadening in a grating-based time-domain delay line for optical coherence tomography application

This paper discusses a dispersion effect in a grating-based time-domain delay line that is different from the second- or higher-order dispersion in a grating-based Fourier-domain delay line. When the lateral broadening of the beam profile after grating dispersion exceeds the collection aperture of the reference fiber, the peripheral spectrum is decoupled by the fiber. The loss of reference spectral bandwidth by this geometric-beam broadening thus degrades the axial resolution. The polarizing-beam reflector used in the Fourier-domain delay line for suppression of lateral beam walk-off is implemented in this grating-based time-domain delay line to minimize geometric-beam broadening. Theoretical analysis and experiments are given to validate the axial resolution improvement after geometric-beam broadening is minimized. In vitro and in vivo imaging results are presented to demonstrate the improvement. It is also shown that geometric-beam broadening may exist in other optical coherence tomography reference arm configurations.

Rapid near-infrared diffuse tomography for hemodynamic imaging using a low-coherence wideband light source

Rapid near-infrared (NIR) diffuse optical tomography is implemented using a low-coherence source. The spectral bandwidth of the low-coherence source is dispersed and coupled to linearly bundled fibers, such that "spread"-spectral encoding among the bundled fibers is formed, and could be used for parallel source illumination onto tissue. In comparison with a previous spectral-encoding technique that employed multiple laser diodes, the use of a low-coherence source for spread-spectral encoding presents a few unique characteristics: (1) it provides shift-free spectral encoding; (2) it reduces the reconstruction uncertainty significantly owing to the minimization of spontaneous channel-to-channel intensity fluctuation; and (3) it enables the implementation of NIR tomography into an endoscopic imaging mode. A 20-mW superluminescent diode centered at 840 nm with a 40-nm bandwidth is used as the source, and a sampling speed of 5 Hz is obtained in a 27-mm imaging array consisting of eight sources and eight detection channels. The principles of using a low-coherence source for spread-spectral encoding are elaborated, the characteristic performances are demonstrated, and the preliminary results of imaging hemoglobin absorption variations during 10 s of voluntary breath-holding are presented.

Endoscopic, rapid near-infrared optical tomography

This is believed to be the first demonstration of near-infrared (NIR) optical tomography employed at the endoscope scale and at a rapid sampling speed that allows translation to in vivo use. A spread-spectral-encoding technique based on a broadband light source and linear-to-circular fiber bundling was used to provide endoscopic probing of many source-detector fibers for tomography as well as parallel sampling of all source-detector pairs for rapid imaging. Endoscopic NIR tomography at an 8 Hz frame rate was achieved in phantoms and tissue specimens with a 12 mm probe housing eight sources and eight detectors. This novel approach provides the key feasibility studies to allow this blood-based contrast imaging technology to be attempted in detection of cancer in internal organs via endoscopic interrogation.

Video-rate near-infrared optical tomography using spectrally encoded parallel light delivery

A novel parallel source implementation approach to near-infrared tomography is demonstrated through spectral encoding of the light delivery. This new technique allows many sources to be input into the tissue at the same time, and a high-resolution spectrometer is used to spatially spread out the signals from each spectrally encoded source. The parallel sampling of all sources at all detection locations renders rapid imaging. Acquisition of complete tomographic data sets at a video rate of 35 frames/s is achieved for imaging of a 6.35 mm diameter inclusion with an absorption coefficient of 0.01 mm(-1) and a reduced scattering coefficient of 1.5 mm(-1) that is moving along a circular path inside a 1% Intralipid solution.

Hybrid positron detection and optical coherence tomography system: design, calibration, and experimental validation with rabbit atherosclerotic models

We evaluate the performance of our novel hybrid optical coherence tomography (OCT) and scintillating probe, demonstrate simultaneous OCT imaging and scintillating detection, and validate the system using an atherosclerotic rabbit model. Preliminary data obtained from the rabbit model suggest that our prototype positron probe detects local uptake of fluorodeoxyglucose (FDG) labeled with 18F positron (beta) radionuclide emitter, and the high-uptake regions correlate with sites of injury and extensive atherosclerosis areas. Preliminary data also suggest that coregistered high-resolution OCT images provide imaging of detailed plaque microstructures, which cannot be resolved by positron detection.