Bearbeiter: R. Hünlich, R. Model (PTB)
Kooperation: Physikalisch-Technische Bundesanstalt Berlin-Charlottenburg (PTB)
Förderung: BMBF/PTB
Beschreibung der Forschungsarbeit:
Die Lichtausbreitung in optisch dichten Medien kann in guter Näherung durch
eine Diffusions-Absorptionsgleichung für die Photonendichte
beschrieben
werden, die als optische Parameter des Mediums den Transportstreukoeffizienten
und den Absorptionskoeffizienten enthält. Für Anwendungen in der
optischen Tomographie ist die Modellierung der Situation am Rande des Objektes
von großer Bedeutung, da diese die Laufzeitkurven des transmittierten Lichtes
in Detektorpositionen am Rande des Objekts, die die Eingangsinformation für
die Bilderkennung liefern, stark beeinflussen.
Das Testobjekt (eine Küvette, mit einem trüben Medium gefüllt und zusätzlich einen dünnen beweglichen Absorber enthaltend), die Durchführung der Messungen (zeitaufgelöste Messung der Reaktion auf eingestrahlte Laserimpulse im ps-Bereich), die Meßdaten und erste Ergebnisse zur Parameteridentifikation sind in [4] beschrieben. Dort sind sowohl auf dem Rand der Küvette als auch am Rand des Absorbers Randbedingungen 1. Art benutzt worden. Inzwischen wurde in Analogie zum Algorithmus zur Bilderkennung (siehe [1,3]) ein Algorithmus zur Parameteridentifikation entwickelt, der es gestattet, bei Randbedingungen 3. Art neben den optischen Parametern auch die Randparameter (auf dem Rand der Küvette bzw. des Absorbers) zu bestimmen, wobei aus meßtechnischen Gründen noch zwei weitere unbekannte Parameter (eine Zeitverschiebung und ein Skalierungsfaktor) zu berücksichtigen sind. Der Algorithmus beruht auf einer FEM-Lösung des Vorwärtsproblems und einer Minimierungsstrategie des Fehlerfunktionals.
Die Abbildungen (entnommen aus [1,2]) zeigen gemessene und simulierte Laufzeitkurven für verschiedene Absorberpositionen, wobei alle simulierten Laufzeitkurven der oberen bzw. unteren Abbildung jeweils mit dem einheitlichen Parametersatz berechnet worden sind, der sich nach Identifikation ergeben hat. Diese Abbildungen zeigen deutlich, daß nur Randbedingungen 3. Art die Meßdaten gut beschreiben, zumal sich physikalisch plausible Werte für die Randparameter ergeben haben.
Die gleichzeitige Identifizierung der optischen Parameter des Mediums und der Parameter aus den Randbedingungen führte auf ein extrem schlecht konditioniertes Problem. Für Anwendungen in der optischen Tomographie sollten die Randparameter vorab bestimmt werden, um bei der Lösung des inversen Problems zur Bilderkennung nicht noch zusätzliche Konvergenzschwierigkeiten zu erhalten.
Projektliteratur: