Virtual Photonics Technology Initiative

2017 Laboratory G: Recovery of Tissue Properties from Time-Resolved and Temporal Frequency-Domain Reflectance Measurements

From VirtualPhotonics

Jump to: navigation, search

GOAL: This GUI Interaction aims to examine (a) the impact of optical absorption and scattering on temporally-resolved and temporal frequency-domain reflectance signals; and (b) the impact of optical properties and measurement selection on the tissue region probed by the detected photons.

Bring up the Silverlight GUI

I. Sensitivity of Time-Resolved Reflectance R(t) to Optical Properties

  1. Select the Forward Solver/Analysis Panel.
  2. Click Clear All and return the Y Axis Spacing back to Linear.
  3. Uncheck use spectral panel inputs.
  4. In the Fwd Solver drop down menu select "Scaled Monte Carlo - NURBS (g=0.8, n=1.4)".
  5. In the Solution Domain select Time-Domain, R(ρ,t).
  6. For the Independent Axis choose t and set ρ = 10 mm.
  7. In Detection Times choose Begin = 0 ns and End = 0.5 ns with Number = 201 time points (1 point every 5 ps).
  8. In Optical Properties: enter μa = 0.01mm-1, μ's = 1mm-1.
  9. Click the Plot Reflectance button.
  10. Confirm that the Hold On checkbox is selected.
  11. Repeat the above for μa = 0.03, 0.1 and 0.3 mm-1.
    Question: Note the difference in the magnitude and shape of these plots. What do you believe is responsible for this? Hint: It may helpful to view the results under both linear and log y-axis spacing.
  12. Click the Clear All button.
  13. Start again with Optical Properties: μa = 0.01mm-1, μ's = 1mm-1.
  14. Click the Plot Reflectance button.
  15. Confirm that the Hold On checkbox is selected.
  16. Repeat the above for μ's = 0.5 and 1.5 mm-1.

Questions:

  1. Note that no photons are detected before a finite time in the time-resolved reflectance signal. Can you independently calculate the minimal delay time?
  2. Note that the peak reflectance values are different and not located at the same time point. Can you speculate as to the origin of these features? Hint: It may help to use both linear and log y-axis spacing.
  3. You are designing a time-resolved optical imaging system to detect early formation of a fibrous tumor. What is the approximate time resolution and source detector separation necessary to differentiate normal breast with μ's=0.6 mm-1 from a fibroid tumor with μ's=1.2 mm-1 using such a system?

II. Optical Property Recovery using Temporally-Resolved Reflectance Measurements: Impact of Noise and Initial Guess

  1. Select the Inverse Solver Panel.
  2. For Fwd Solver: select "Scaled Monte Carlo - NURBS (g=0.8, n=1.4)", for Inv Solver: select "Standard Diffusion (Analytic - Isotropic Point Source)".
  3. In Solution Domain select "Time-Domain(R(ρ,t))".
  4. Select t as the Independent Axis and ρ=10mm.
  5. In Detection Times choose Begin = 0 ns and End = 1.0 ns with Number = 51 time points (1 point every 20 ps).
  6. Set Optimization Parameters to: μa and μ's.
  7. Simulate measured data: set Forward Simulation Optical Properties: to: μa = 0.01 mm-1, μ's = 1 mm-1, g = 0.8 and n = 1.4 and 2% noise.
  8. Confirm the Hold On checkbox is checked.
  9. Click the Plot Measured Data button.
  10. Set Initial Guess Optical Properties: to: μa = 0.02 mm-1, μ's = 1.2 mm-1, g = 0.8 and n = 1.4.
  11. Click the Plot Initial Guess button.
  12. Click the Run Inverse Solver button.

Questions:

  1. To what optical property values did the inverse solver converge? (Scroll to the bottom of the page to see the output).
  2. Why are the converged values not exactly the forward simulation optical properties?
  3. Perform the same analysis with 0% noise added to the simulated measured data. How accurate are the converged properties now?
  4. Perform the same analysis with initial guess μa = 0.05 mm-1, μ's = 0.7 mm-1, g = 0.8 and n=1.4. How accurate are the converged properties now?
  5. How would you modify the Detection Times to improve the inverse solution? Run the inverse solution with this new time window and check your results.


III. Sensitivity of Temporal Frequency Domain Reflectance to Optical Properties

First let us examine the sensitivity of Temporal Frequency Domain Reflectance to optical absorption

  1. Go to the Forward Solver/Analysis Panel
  2. For Fwd Solver: select "Scaled Monte Carlo - NURBS (g=0.8, n=1.4)"
  3. In Solution Domain select "Frequency Domain(R(ρ,ft))".
  4. Select Independent Axis ft and set ρ=10mm.
  5. In Temporal Frequencies choose Begin = 0 GHz and End = 2.0 GHz with Number = 101 frequency points (1 point every 20 MHz).
  6. In Optical Properties: enter μa = 0.01mm-1, μ's=1mm-1, n=1.4.
  7. Click the Plot Reflectance button.
  8. The Plot Toggle radio buttons toggle the plot from real/imag results to phase and amplitude results. The phase is shown in units of degrees.
  9. Confirm the Hold On checkbox is checked.
  10. Fix μ's=1mm-1. Repeat the above steps for μa = 0.03 and 0.1 mm-1.
  11. Note the trend of decreasing reflectance with increasing absorption.
  12. Note what temporal frequency regime shows the most sensitivity to μa changes.
  13. Is the temporal frequency that shows the most sensitivity to μa the same for real/imag, phase and amplitude?

Now let us examine the sensitivity of Temporal Frequency Domain Reflectance to optical scattering

  1. Click the Clear All button and toggle back to Linear y-axis spacing.
  2. In Temporal Frequencies choose Begin = 0 GHz and End = 2.0 GHz with Number = 101 time points (1 point every 20 GHz).
  3. In Optical Properties: enter μa = 0.01mm-1, μ's = 0.5 mm-1, n=1.4.
  4. Click the Plot Reflectance button.
  5. Confirm the Hold On checkbox is checked.
  6. Fix μa=0.01mm-1. Repeat the above steps for μ's = 1.0 and 1.5 mm-1.
  7. Note what temporal frequency regime shows the most sensitivity to μ's changes.
  8. Is the temporal frequency domain that shows the most sensitivity to μ's the same for real/imag, phase and amplitude?

IV. Optical Property Recovery using Temporal Frequency Domain Reflectance Measurements: Impact of Noise and Initial Guess

  1. Click Clear All and set Normalization to None.
  2. Select the Inverse Solver Panel.
  3. For Fwd Solver: select "Scaled Monte Carlo - NURBS (g=0.8, n=1.4)", for Inv Solver: select "Standard Diffusion (Analytic - Isotropic Point Source)".
  4. In Solution Domain select "Temporal Frequency(R(ρ,ft))".
  5. Select Independent Axis ft and set ρ=10mm.
  6. In Temporal Frequencies choose Begin = 0 GHz and End = 0.5 GHz with Number = 51 time points (1 point every 20 GHz).
  7. Set Optimization Parameters to: μa and μ's.
  8. Simulate measured data: set Forward Simulation Optical Properties: to: μa = 0.05 mm-1, μ's = 1 mm-1, g = 0.8 and n = 1.4 and 2% noise.
  9. Confirm the Hold On checkbox is checked.
  10. Click the Plot Measured Data button.
  11. Set Initial Guess Optical Properties: to: μa = 0.01 mm-1, μ's = 1.5 mm-1, g = 0.8 and n = 1.4.
  12. Click the Plot Initial Guess button.
  13. Click the Run Inverse Solver button.

Questions:

  1. To what optical property values did the inverse solver converge? (Scroll to the bottom of the page to see the output).
  2. Perform the same analysis with 0% noise added to the simulated measured data. How accurate are the converged properties now?
  3. Perform the same analysis with initial guess μa = 0.001 mm-1, μ's = 0.5 mm-1, g = 0.8 and n=1.4. How accurate are the converged properties now?

V. Effect of Measurement Range on Sensitivity to Optical Absorption and Scattering in Temporal Frequency Domain Reflectance

  1. Go to the Inverse Solver Panel.
  2. Follow the instructions provided in Section IV except modify the Temporal Frequency begin and end values to those obtained in Section III. "Sensitivity of Temporal Frequency Domain Reflectance to Optical Properties". Note that the inverse solution fits the real/imag measurements.
  3. Rerun the inverse solver.

Questions

  1. Were you able to improve the μa and μ's converged properties?
  2. In what temporal frequency domain is reflectance most sensitive to μa? Why is this?
  3. In what temporal frequency domain is reflectance most sensitive to μ's? Why is this?
Personal tools