zum Seiteninhalt

Physikalisch-Technische Bundesanstalt

FachabteilungenAbt. 8 Medizinphysik und metrologische Informationstechnik 8.4 Mathematische Modellierung und Datenanalyse8.41 Modellierung und Simulation > Simulation und Rekonstruktion zur optischen Tomographie
Simulation und Rekonstruktion zur optischen Tomographie
Arbeitsgruppe 8.41



Die optische Tomographie beruht auf der Durchleuchtung von Gewebe im Infrarotbereich mit potentieller Anwendung u.a. zur Brustkrebserkennung. Grundannahme ist hier, dass Tumore Änderungen in den optischen Eigenschaften gegenüber dem gesunden Gewebe verursachen können. Die Streulichtausbreitung kann unter bestimmten Bedingungen als Diffusionsprozeß beschrieben werden, welche für biologisches Gewebe erfüllt sind. In der resultierenden parabolischen Differentialgleichung für die Photonendichte treten die verteilten, optischen Parameter (Lichtabsorption µa und –transportstreuung µs´) als Koeffizienten auf. In lokalen Gebieten mit erhöhter Absorption oder Streuungen können Tumore erwartet werden. Mathematisch führt die Rekonstruktion der optischen Parameter auf die Lösung eines inversen Problems.
Mit der Lösung des Vorwärtsproblems der Lichtausbreitung, d.h. der FEM-Simulation der Laufzeitverteilungen an den Detektorpositionen  kann ein beliebiges Experiment simuliert werden, d.h. ein virtuelles Experiment durchgeführt werden.

Beispiel 1

Simulationsrechnungen zur Lichtausbreitung in einem Phantom (linkes Bild oben) in umgebender Flüssigkeit (linkes Bild unten). Der rote Pfeil weist auf die Einfallstelle des Laserimpulses, in deren Umgebung das Finite-Element-Gitter aufgrund der großen Gradienten der Photonendichte in einem Anfangsintervall stark verfeinert werden muss (Kooperation mit Arbeitsgruppe 8.31).

 

Nach oben

Beispiel 2

2D FEM-Simulation zur Lichtausbreitung in einem bis auf einen mittigen Absorber homogenen Gebiet. Im Absorber wurde der Absorptionskoeffizient auf das 10-fache des umgebenden Gebiets erhöht.

 

Nach oben

Beispiel 3

Im Rahmen eines abgeschlossenen, BMBF-geförderten Forschungsprojekt wurden zwei Algorithmen zur Rekonstruktion der optischen Eigenschaften entwickelt. Der erste beruht auf einem gewichteten “output least-squares” Problem mit Regularisierung. Für den zweiten wurde ein Störungsansatz ohne Regularisierung verwendet.

 

Nach oben




Literatur

  • H. Wabnitz, A. Liebert, M. Möller, D. Grosenick, R. Model, H. Rinneberg (2002). Scanning laser-pulse mammography: matching fluid and off-axis measurements. OSA Technical Digest (Optical Society of America, Washington DC), 686-688.
  • R. Model, M. Orlt, M. Walzel, R. Hünlich (1998). Optical imaging: three-dimensional approximation and perturbation approaches for time-domain data.   Appl. Optics 37, 7968-7976.
  • R. Model, M. Orlt, M. Walzel, R. Hünlich (1998).  Reconstruction algorithm for near  infrared imaging in turbid media by means of time-domain data. In  O. Minet, G. Müller, and J. Beuthan (eds), ” Optical Tomography, Fundamentals and Applications in Medicine”, SPIE MS 147, Optical Engeneering Press, Bellingham, 537-548.
  • R. Model, R. Hünlich, M. Orlt, M. Walzel (1997). Mathematische Behandlung der Streulichttomographie von dicken Gewebeschichten und Phantomen. In VDI Technologiezentrum, ”Lasermedizinforschung” , 130-132.
  • M. Orlt, M. Walzel, R. Model (1997). Influence of different kinds of errors on imaging results in optical tomography.   In P. Ciarlini,  M.G. Cox, F. Pavese, D.  Richter (eds) „Advanced mathematical tools in metrology III“ , World Scientific Singapore, New York, London, Hong Kong, 277-279.
  • R. Model, R. Hünlich, M. Orlt, M. Walzel (1997). Reconstruction algorithm for near infrared imaging in turbid media by means of time-domain data.  J. Opt. Soc. Am. A   14, 313-324.
  • R. Model, R. Hünlich, M. Orlt, M. Walzel (1996). NIR imaging in random media using time domain data. Proc. SPIE 2925,  77-88.    
  • R. Model, R. Hünlich (1996). Optical imaging of highly scattering media.   Z. angew. Math.  Mech. 76, 483-484.
  • R. Model (1996). Reconstruction algorithms for optical tomography. In P. Ciarlini, M.G. Cox,  F.  Pavese, D.  Richter (eds) „Advanced mathematical tools in metrology II“ , World Scientific Singapore, New York, London, Hong Kong, 172-181.
  • M. Orlt, M. Walzel, R. Model (1995). Transillumination imaging performance using time domain data. Proc. SPIE 2626, 346-357.
  • R. Model, R. Hünlich (1995). Parameter sensitivity in near infrared imaging.   Proc. SPIE 2626, 56-65.
  • R. Model, R. Hünlich, M. Orlt, M. Walzel, (1995). I mage reconstruction for random media by diffusion tomography.  Proc. SPIE 2389, 400-410.
  • R. Model, R. Hünlich, D. Richter, H. Rinneberg, H. Wabnitz, M. Walzel (1995). Imaging in random media: Simulating light transport by numerical integration of the diffusion equation Proc. SPIE  2326, 11-22.
  • D. Richter, R. Model, H. Rinneberg (1993). Simulation in der Streulichttomographie. In: 4. Fachgespräch der Fachgruppe Physik ,  Informatik, Informationstechnik (PII),  Tagungssammelband ,  25-34.

Nach oben




Kontakt

Physikalisch-Technische Bundesanstalt
Arbeitsgruppe 8.41 Modellierung und Simulation
Abbestr. 2-12
10587 Berlin

Nach oben


© Physikalisch-Technische Bundesanstalt, letzte Änderung: 2012-05-24, Webmaster Abteilung 8 Seite drucken DruckansichtPDF-Export PDF