### Abstract

This paper describes an iterative perturbation approach for imaging the absorption properties of a dense scattering medium. This method iteratively adjusts the current estimate until the calculated photon fluxes for the estimated medium match the detected readings. The inverse update in each iteration is accomplished by solving a linear perturbation equation. It is similar to the compensation theory method used in electrical impedance tomography. A comparison was made between the methods of conjugate gradient descent and projection onto convex sets for the solution of the perturbation equation. The former converges more rapidly, but can yield an inaccurate solution where the problem is underdetermined. The latter can incorporate many types of a priori information to reach a correct solution, but progresses very slowly. A multigrid progressive reconstruction technique is proposed which computes the fine details with the help of the coarse structure. It is quite effective in forcing the correct solution and reducing computation time. These methods have been used to reconstruct several inhomogeneous media containing simple structures, from steady-state reflectance data. Two sets of data are tested: one calculated according to the perturbation model, and the other using Monte Carlo methods. When the difference between the absorption distributions of the test medium and the initial estimate is localized, a single step of the perturbation approach can resolve the absorption distribution reasonably well to within five transport mean free path lengths from the surface. At greater depths, the reconstruction is less reliable.

Original language | English (US) |
---|---|

Title of host publication | Proceedings of SPIE - The International Society for Optical Engineering |

Publisher | Publ by Int Soc for Optical Engineering |

Pages | 58-71 |

Number of pages | 14 |

Volume | 1641 |

ISBN (Print) | 0819407879 |

State | Published - 1992 |

Event | Physiological Monitoring and Early Detection Diagnostic Methods - Los Angeles, CA, USA Duration: Jan 22 1992 → Jan 23 1992 |

### Other

Other | Physiological Monitoring and Early Detection Diagnostic Methods |
---|---|

City | Los Angeles, CA, USA |

Period | 1/22/92 → 1/23/92 |

### Fingerprint

### ASJC Scopus subject areas

- Electrical and Electronic Engineering
- Condensed Matter Physics

### Cite this

*Proceedings of SPIE - The International Society for Optical Engineering*(Vol. 1641, pp. 58-71). Publ by Int Soc for Optical Engineering.

**Imaging of scattering media by diffusion tomography : an iterative perturbation approach.** / Wang, Yao; Chang, Jeng Hua; Aronson, Raphael; Barbour, Randall L.; Graber, Harry L.; Lubowsky, Jack.

Research output: Chapter in Book/Report/Conference proceeding › Conference contribution

*Proceedings of SPIE - The International Society for Optical Engineering.*vol. 1641, Publ by Int Soc for Optical Engineering, pp. 58-71, Physiological Monitoring and Early Detection Diagnostic Methods, Los Angeles, CA, USA, 1/22/92.

}

TY - GEN

T1 - Imaging of scattering media by diffusion tomography

T2 - an iterative perturbation approach

AU - Wang, Yao

AU - Chang, Jeng Hua

AU - Aronson, Raphael

AU - Barbour, Randall L.

AU - Graber, Harry L.

AU - Lubowsky, Jack

PY - 1992

Y1 - 1992

N2 - This paper describes an iterative perturbation approach for imaging the absorption properties of a dense scattering medium. This method iteratively adjusts the current estimate until the calculated photon fluxes for the estimated medium match the detected readings. The inverse update in each iteration is accomplished by solving a linear perturbation equation. It is similar to the compensation theory method used in electrical impedance tomography. A comparison was made between the methods of conjugate gradient descent and projection onto convex sets for the solution of the perturbation equation. The former converges more rapidly, but can yield an inaccurate solution where the problem is underdetermined. The latter can incorporate many types of a priori information to reach a correct solution, but progresses very slowly. A multigrid progressive reconstruction technique is proposed which computes the fine details with the help of the coarse structure. It is quite effective in forcing the correct solution and reducing computation time. These methods have been used to reconstruct several inhomogeneous media containing simple structures, from steady-state reflectance data. Two sets of data are tested: one calculated according to the perturbation model, and the other using Monte Carlo methods. When the difference between the absorption distributions of the test medium and the initial estimate is localized, a single step of the perturbation approach can resolve the absorption distribution reasonably well to within five transport mean free path lengths from the surface. At greater depths, the reconstruction is less reliable.

AB - This paper describes an iterative perturbation approach for imaging the absorption properties of a dense scattering medium. This method iteratively adjusts the current estimate until the calculated photon fluxes for the estimated medium match the detected readings. The inverse update in each iteration is accomplished by solving a linear perturbation equation. It is similar to the compensation theory method used in electrical impedance tomography. A comparison was made between the methods of conjugate gradient descent and projection onto convex sets for the solution of the perturbation equation. The former converges more rapidly, but can yield an inaccurate solution where the problem is underdetermined. The latter can incorporate many types of a priori information to reach a correct solution, but progresses very slowly. A multigrid progressive reconstruction technique is proposed which computes the fine details with the help of the coarse structure. It is quite effective in forcing the correct solution and reducing computation time. These methods have been used to reconstruct several inhomogeneous media containing simple structures, from steady-state reflectance data. Two sets of data are tested: one calculated according to the perturbation model, and the other using Monte Carlo methods. When the difference between the absorption distributions of the test medium and the initial estimate is localized, a single step of the perturbation approach can resolve the absorption distribution reasonably well to within five transport mean free path lengths from the surface. At greater depths, the reconstruction is less reliable.

UR - http://www.scopus.com/inward/record.url?scp=0026439503&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0026439503&partnerID=8YFLogxK

M3 - Conference contribution

AN - SCOPUS:0026439503

SN - 0819407879

VL - 1641

SP - 58

EP - 71

BT - Proceedings of SPIE - The International Society for Optical Engineering

PB - Publ by Int Soc for Optical Engineering

ER -