For a more thorough explanation of the PTOCT instrumentation, image collection protocols, and image corrections algorithms, please see our previously published work.
16 Briefly, a spectral-domain OCT system (Leica Microsystems, Inc.) was modified to allow for PTOCT imaging. The OCT system's broadband light source was centered at 860 nm (40 nm bandwidth), while the additional laser used for photothermal imaging was centered at 750 nm with a power of 8 mW and amplitude-modulated by an acousto-optics modulator (Brimrose, Sparks, MD) at f
0 = 500 Hz (50% duty-cycle, square wave). Scan patterns were 700 repeated A-scans (M-scans) at each sample location, and 400 A-scans/B-scans. Collected data were processed using custom MATLAB code.
27 The data were first resampled from wavelength to wavenumber, dispersion compensated, and background subtracted. A Chirp Z transform was used to obtain the OCT magnitude and phase signal as a function of depth. The temporal derivative of the phase signal was then taken over time. A Fourier transform of the phase signal was performed to go from the time domain to the frequency domain. The amplitude of the peak at f
0 = 500 Hz minus the average amplitude of the signal at surrounding frequencies was taken as the PTOCT signal. The PTOCT signal is expressed as the change in the optical path length in the tissue due to the photothermal excitation in units of nanometers. Images were corrected for breathing artifacts associated with retinal imaging in live animals following an algorithm described in Guizar-Sicairos et al.
16,47 Once processed, data analysis was performed by hand-selecting the lesion of anatomic laser damage on individual B-scans that have been imaged in a given retina and calculating both the lesion volume and the lesion-associated PTOCT signal inside the selected volume by collating available B-scan measurements into volumes. Lesions were excluded from measurement if there were gross morphological anomalies on OCT, most commonly retinal detachment, and individual B-scans were excluded if the image quality was degraded to a point where photothermal signal measurement was infeasible. These volume measurements were then used to compute a photothermal signal density, defined as the PT signal divided by lesion volume. Selections were made to include the lesions but exclude the adjacent RPE layer, since the RPE layer contains melanin, which also produces a photothermal signal. Selections were made based on the OCT image while blinded to the PTOCT signal present in the lesion. Additional nonphotothermal OCT scans of lesions were used to estimate the extent of anatomic laser damage by measuring the size of the lesion on the B-scan with the maximum area of anatomic laser damage for a given set of scans, which is the lesion diameter.