DE10217543B4

DE10217543B4 - Verfahren zur spektrometrischen Bestimmung der Sauerstoffsättigung von Blut in Gegenwart optischer Störgrößen

Info

Publication number: DE10217543B4
Application number: DE10217543A
Authority: DE
Inventors: Dr.-Ing.habil. Schweitzer Dietrich; Dr.rer.nat. Hammer Martin; Dr.rer.nat. Thamm Eike
Original assignee: Carl Zeiss Meditec AG
Current assignee: Carl Zeiss Meditec AG
Priority date: 2002-04-17
Filing date: 2002-04-17
Publication date: 2012-12-06
Anticipated expiration: 2022-04-18
Also published as: US20060063994A1; DE10217543A1; WO2003086193A1; US7333842B2; US20080194931A1

Abstract

Verfahren zur Bestimmung der Sauerstoffsättigung von Blut in Gegenwart optischer Störgrößen, insbesondere durch ein das Blutgefäß umgebendes biologisches Gewebe und/oder des Blutes bzw. des Blutgefäßes selbst, bei dem durch Transmissions- bzw. Reflexionsmessung in einem Mess-Spektrum spektrale Messwerte (Mi) bei für Hämoglobin (Hb) und Oxihämoglobin (HbO2) isosbestischen Wellenlängen und mindestens ein anderer Messwert (Ma) bei einer Wellenlänge, bei welcher sich die Referenzwerte von Hämoglobin und Oxihämoglobin unterscheiden, erzeugt sowie mit bekannten Referenzwerten der Referenz-Spektren von Hämoglobin und Oxihämoglobin verglichen werden, dadurch gekennzeichnet, a) dass in dem Mess-Spektrum zumindest zwei besagte spektrale Messwerte (Mi1, Mi2), bei für Hämoglobin und Oxihämoglobin isosbestischen Wellenlängen (λi1, λi2) und zumindest der andere Messwert (Ma) bei einer Wellenlänge (λa), bei welcher sich die Referenzwerte von Hämoglobin und Oxihämoglobin in den Referenz-Spektren möglichst weitgehend unterscheiden, logarithmisch erfasst werden, wobei zumindest aus zwei der Messwerte (Mi1, Mi2) für isosbestische Wellenlängen (λi1, λi2) eine lineare Hilfsfunktion (FH) generiert wird,...

Description

Die Erfindung betrifft ein Verfahren zur spektrometrischen Bestimmung der Sauerstoffsättigung von Blut in Gegenwart optischer Störgrößen, wie sie beispielsweise ebenfalls pigmentiertes sowie Licht streuendes Umgebungsgewebe und/oder auch die Blutgefäßwandung selbst bilden. Dieses Problem, Sauerstoffsättigung von Blut möglichst ohne den Einfluss dieser die Messgenauigkeit beeinflussenden Faktoren zu bestimmen tritt insbesondere bei nichtinvasiven in vivo- oder in vitro-Untersuchungen von Blutgefäßen auf, die vor, hinter bzw. in dem besagten pigmentierten und streuenden Gewebe liegen, beispielsweise bei der Untersuchung von Blutgefäßen des Augenhintergrundes oder anderer Gewebebereiche des Körpers, wie der Haut und endoskopisch zugänglicher Organe.
Es ist allgemein bekannt, dass sich das Absorptionsspektrum des roten Blutfarbstoffes Hämoglobin mit der Sauerstoffsättigung ändert (beispielsweise van Assendelft O. W., Spectrophotometry of heamoglobin derivatives, Assen: Royal Vangorcum, 1970). Damit ist die Bestimmung der Sauerstoffsättigung einer Hämoglobinprobe durch den Vergleich des Spektrums der Probe mit den Spektren von vollständig oxygeniertem und vollständig reduziertem Hämoglobin möglich.
Neuere Arbeiten zur Oximetrie am Augenhintergrund nach dem Lambert – Beerschen Gesetz, d. h. unter Berücksichtigung lediglich der Absorption, sind u. a. von Smith et al (Smith M. H., Denninghoff K. R., Lompado A., Hillman L. W., Effect of multiple light paths in retinal vessel oximetry, Appl. Opt. 39, 2000, 1183–1193) veröffentlicht worden. Zahlreiche patentrechtlich geschützte Anordnungen und Verfahren beruhen auf diesem Prinzip (beispielsweise US 4,485,820 ; US 5,119,814 ; US 5,308,919 ; US 4,253,744 ; US 4,305,398 ; US 5,776,060 ; US 5,935,076 DE 199 20 157 A1 ; US 5,318,022 ). Bei der in vivo-Messung liegt das Hämoglobin jedoch nicht isoliert vor, sondern es ist eingeschlossen in den Erythrozyten. Die Streuung des Lichtes an den Erythrozyten beeinflusst das Extinktionsspektrum des Blutes erheblich. Nach Ergebnissen der Multiple Scattering Theorie von Twersky (Twersky V., Absorption and multiple scattering by biological suspensions, J. Opt. Soc. Amer. 60, 1970, 1084–1093) lassen sich aber die Einflüsse von Streuung und Absorption separieren. Auf dieser Grundlage haben Pittman und Duling ein Verfahren beschrieben, das die Sauerstoffsättigung in Vollblut aus in Transmission ausgeführten Messungen bei der Wellenlänge 555 nm sowie an den isosbestischen Punkten bei 522 nm und 546 nm bestimmt (Pittman RN, Duling BR. A new method for the measurement of percent oxyhemoglobin. J. Appl. Physiol. 38, 1975, 315–320). Dieses Verfahren ist von Delori zur Bestimmung der Sauerstoffsättigung in retinalen Gefäßen verwandt worden (Delori F. C., Noninvasive technique for oximetry of blood in retinal vessels. Appl. Opt, 27, 1988, 113–1125).
Untersuchungen von Hammer et al (Hammer M., Leistritz S., Leistritz L., Schweitzer D., Light paths in retinal vessel oxymetry, IEEE Trans Biomed Eng 48 (5), 2001, 592–8) zeigen jedoch, dass die auf retinalen Gefäßen gemessenen Reflexionsspektren nicht nur von der Absorption des Hämoglobins und der Streuung im Blut und dem die Gefäße umgebenden Gewebe beeinflusst werden, sondern auch von dem im retinalen Pigmentepithel und in der Aderhaut lokalisierten Melanin. Gleiches gilt auch für Gefäße in der Haut oder anderen pigmentierten Organen.
Die Korrektur der Verfälschung der Hämoglobinspektren durch andere Chromophore und deren Korrektur für die spektroskopische Oximetrie ist in der bisherigen Literatur durch Normierung der auf einem Gefäß gemessenen Spektren auf Messungen neben dem Gefäß versucht worden (z. B. DE 199 20 157 A1 ; US 5,935,076 ; Delori F. C., Noninvasive technique for oximetry of blood in retinal vessels. Appl. Opt. 27, 1988, 113–1125; Schweitzer D., Hammer M., Kraft J., Thamm E., Königsdörffer E., Strobel J., In Vivo Measurement of the Oxygen Saturationm at the Normal Human Eye, IEEE Trans. Biomed. Eng. 46, 1999, 1454–1465). Dieser Ansatz lässt jedoch die überaus komplizierten Verhältnisse (Hammer M., Leistritz S., Leistritz L., Schweitzer D., Light paths in retinal vessel oxymetry, IEEE Trans Biomed Eng 48 (5), 2001, 592–8) der Strahlungsausbreitung im Blutgefäß und dem dieses umgebenden Gewebe außer Acht.
Die exakte Lichtausbreitung im biologischen Gewebe lässt sich noch immer nicht vollständig physikalisch beschreiben. Selbst Anstrengungen, diese Vorgänge zur Störgrößeneliminierung möglichst umfassend und realitätsnah nachzuempfinden ( DE 199 20 157 A1 ) sowie die Optik des die Blutgefäße umgebenden lebenden oder toten biologischen Gewebes zu modellieren ( DE 44 33 827 A1 ), haben nicht zu genaueren Messergebnissen als die vorgenannten Verfahren geführt, welche ohnehin schon relativ zeitaufwendig und rechenintensiv sind. Mit diesem Aufwand sind die Verfahren speziell für Routine- und Vorsorgeuntersuchungen lediglich bedingt geeignet. Insbesondere die für die klinische Praxis interessante Bestimmung der Sauerstoffsättigung an jedem Punkt einer zweidimensionalen, bildhaften Aufnahme erfordert ein Verfahren, das einerseits schnell ist, aber andererseits optische und spektrometrische Störungen durch die Gefäßumgebung kompensiert.
Der Erfindung liegt deshalb die Aufgabe zu Grunde, die Sauerstoffsättigung in einem möglichst einfachen und schnellen Verfahren dennoch mit hoher Genauigkeit zu bestimmen.
Gelöst wird diese Aufgabe durch ein Verfahren zur Bestimmung der Sauerstoffsättigung von Blut mit den im Patentanspruch 1 angegebenen Merkmalen.
Vorteilhaft ist es, wenn die spektralen Messwerte und Referenzdaten logarithmisch erfasst werden und die Hilfs- sowie Referenzfunktion jeweils durch eine Gerade gebildet werden, auf welcher zwei der Mess- bzw. Referenzwerte für isosbestische Wellenlängen liegen. Die aus linearer Hilfs- sowie Referenzfunktion gebildete Korrekturfunktion lässt eine ebenfalls lineare korrigierte Hilfsfunktion des korrigierten Mess-Spektrums entstehen. Die übrigen spektralen Messwerte, d. h. der spektrale Messwert für die dritte isosbestische Wellenlänge und der andere Messwert, bei einer Wellenlänge, bei welcher sich die Referenzwerte von Hämoglobin und Oxihämoglobin in den Referenz-Spektren möglichst weitgehend unterscheiden, werden mit einem konstanten Multiplikator beaufschlagt, welcher so bestimmt wird, dass der durch diese Streckung korrigierte dritte isosbestische Messwert des korrigierten Mess-Spektrums mit dem korrespondierenden Referenzwert der Referenz-Spektren übereinstimmt. In diesem speziellen Fall kann die Differenz der Referenzwerte für Hämoglobin und Oxihämoglobin linear zwischen 0 und 1 skaliert werden. Die Sauerstoffsättigung des Blutes wird aus dem auf die korrigierte Hilfsfunktion des korrigierten Mess-Spektrums umgerechneten anderen Messwert in Relation zu dieser Skala ermittelt.
Für eine übersichtliche zweidimensionale Darstellung der Sauerstoffsättigung des Blutes werden vier monochromatische Einzelbilder bei den besagten Wellenlängen erzeugt und die Sauerstoffsättigung für jeden Bildpunkt ermittelt.
Überraschend hat sich gezeigt, dass sich mit dem vorgeschlagenen Verfahren im Vergleich zu den eingangs genannten Untersuchungsmethoden die Sauerstoffsättigung des Blutes mit gleicher Genauigkeit, jedoch mit wesentlich geringerem Aufwand (minimal vier Messwerte erforderlich) für Messung, Berechnung und Auswertung, bestimmen lässt. Das Verfahren ermöglicht eine zweidimensionale ortsabhängige Darstellung der Messergebnisse zur übersichtlichen und schnellen Auswertbarkeit. Es benötigt nur wenige Messwerte und verwendet lediglich lineare Transformationen. Mit diesen Vorzügen des geringen Aufwandes und der zeiteffektiven Messauswertung ist die vorgeschlagene Methode auch zum Screening sowie für Routine- bzw. Vorsorgeuntersuchungen geeignet.
Die Erfindung soll nachstehend anhand eines in der Zeichnung dargestellten Ausführungsbeispiels näher erläutert werden.
Es zeigen:
1: Diagramm mit spektralen Mess- und Referenzwerten in logarithmischer Darstellung im Wellenlängenbereich zwischen 400 nm und 700 nm, einschließlich dreier isosbestischer Wellenlängen λ_i1 = 522 nm, λ_i2 = 586 nm, λ_i3 = 569 nm sowie der anderen Wellenlänge λ_a = 555 nm
2: Diagramm mit spektralen Referenzwerten gemäß 1 mit eingezeichneter linearer Referenz-Funktion F_R sowie mit korrigierten Messwerten M'
3: Diagramm mit spektralen Referenzwerten gemäß 1 mit eingezeichneter linearer Referenz-Funktion F_R, mit korrigierten Messwerten M'' sowie mit Skalierung zur Ablesung der Sauerstoffsättigung
Die Reflexion oder die Transmission von Gewebe an einem Ort oder mit örtlicher Auflösung in einem Bild wird bei drei isosbestischen Wellenlängen λ_i1, λ_i2 und λ_i3 (λ_i1 = 522 nm, λ_i2 = 586 nm, λ_i3 = 569 nm) als Messdaten M_i1, M_i2 und M_i3 sowie bei einer anderen Wellenlänge λ_a (555 nm), bei der sich die Absorptionskoeffizienten von oxigeniertem und reduziertem Hämoglobin unterscheiden, als Messwert M_a gemessen und mit der Reflexion oder Transmission von Hämoglobin oder Vollblut mit Sauerstoffsättigungen von jeweils 0% und 100% bei diesen Wellenlängen (Referenzdaten R₁, R₂, R₃, R 0% / α und R 100% / α) nach dem folgendem Verfahren verglichen:

1. Alle Meß- und Referenzdaten werden logarithmiert. 1 zeigt die Mess- und die Referenzdaten in logarithmischer Darstellung, in diesem Beispiel die Reflexion einer retinalen Vene (Messdaten) sowie die Transmission einer 0,1 mm dicken Schicht Vollblut (Referenzdaten). Der Anschaulichkeit halber sind die vollständigen Spektren zwischen 400 nm und 700 nm dargestellt. Die in diesem Beispiel verwendeten Wellenlängen λ_i1 = 522 nm, λ_i2 = 586 nm, λ_i3 = 569 nm und λ_a = 555 nm sind eingezeichnet.
2. Eine lineare Hilfsfunktion F_H der Wellenlänge wird im Mess-Spektrum so berechnet, dass ihre Werte bei den isosbestischen Wellenlängen λ_i1 und λ_i1 mit den Messdaten M_i1 und M_i2 bei diesen Wellenlängen übereinstimmen.
3. Eine lineare Referenzfunktion F_R der Wellenlänge wird in den Referenzspektren so berechnet, dass ihre Werte bei den isosbestischen Wellenlängen λ_i1 und λ_i2 mit den Referenzdaten R_i1 und R_i2 bei diesen Wellenlängen übereinstimmen.
4. Die Messdaten werden bei jeder Wellenlänge um die Differenz der linearen Funktionen F_H und F_R additiv so korrigiert, dass sie bei den isosbestischen Wellenlängen λ_i1 und λ_i2 mit den Referenzdaten übereinstimmen: M'_λ = M_λ + F_R – F_H. 2 zeigt wiederum die Referenzdaten, die lineare Referenzfunktion F_R der Wellenlänge sowie die korrigierten Messdaten M'. Diese Korrektur kompensiert zusätzlich zur Absorption des Hämoglobins bestehende Extinktionen, deren Spektren im Wellenlängenbereich 522 nm bis 586 nm als linear im logarithmischen Maßstab angenommen bzw. approximiert werden können. Im hier betrachteten Ausführungsbeispiel sind dies die Absorptionen des Melanin und der vorderen Augenmedien sowie die Streuung im Blut und im umgebenden Gewebe.
5. Die korrigierten Messdaten M' werden so um die lineare Referenzfunktion F_R gestreckt oder gestaucht, dass sie bei der isosbestischen Wellenlänge λ_i3 mit dem Referenzwert R_i3 übereinstimmen:
3 zeigt die in M'' resultierende Spreizung (ggf. auch Stauchung) der korrigierten Messdaten M' um die lineare Referenzfunktion F_R, die so vorgenommen wird, dass korrigierte Messdaten und Referenzdaten bei der isosbestischen Wellenlänge λ_i3 (569 nm) übereinstimmen. Diese Korrektur kompensiert unterschiedliche Absolutwerte der Mess- und Referenzdaten, die durch unterschiedliche Beleuchtungs- und Meßbedingungen entstehen.
6. Die Lage von M''_a auf einer linear zwischen R 0% / α und R 100% / α aufgespannten Skala gibt die Sauerstoffsättigung OS an:

Die skalierte Ablesung der Sauerstoffsättigung (OS) zwischen den Werten 0 und 1 ist ebenfalls in 3 dargestellt. Die abgelesene Sauerstoffsättigung OS beträgt im Ausführungsbeispiel 0,69.
Aufstellung der verwendeten Bezugszeichen

M

– Messdaten

M'

– durch Addition der korrigierten Hilfsfunktion korrigierte Messdaten

M''

– durch Beaufschlagung mit einem Faktor korrigierte Messdaten M'

F_R

– Referenzfunktion

OS

– Sauerstoffsättigung

λ

– Wellenlänge

Claims

Verfahren zur Bestimmung der Sauerstoffsättigung von Blut in Gegenwart optischer Störgrößen, insbesondere durch ein das Blutgefäß umgebendes biologisches Gewebe und/oder des Blutes bzw. des Blutgefäßes selbst, bei dem durch Transmissions- bzw. Reflexionsmessung in einem Mess-Spektrum spektrale Messwerte (M_i) bei für Hämoglobin (Hb) und Oxihämoglobin (HbO₂) isosbestischen Wellenlängen und mindestens ein anderer Messwert (M_a) bei einer Wellenlänge, bei welcher sich die Referenzwerte von Hämoglobin und Oxihämoglobin unterscheiden, erzeugt sowie mit bekannten Referenzwerten der Referenz-Spektren von Hämoglobin und Oxihämoglobin verglichen werden, dadurch gekennzeichnet, a) dass in dem Mess-Spektrum zumindest zwei besagte spektrale Messwerte (M_i1, M_i2), bei für Hämoglobin und Oxihämoglobin isosbestischen Wellenlängen (λ_i1, λ_i2) und zumindest der andere Messwert (M_a) bei einer Wellenlänge (λ_a), bei welcher sich die Referenzwerte von Hämoglobin und Oxihämoglobin in den Referenz-Spektren möglichst weitgehend unterscheiden, logarithmisch erfasst werden, wobei zumindest aus zwei der Messwerte (M_i1, M_i2) für isosbestische Wellenlängen (λ_i1, λ_i2) eine lineare Hilfsfunktion (F_H) generiert wird, b) dass in den Referenz-Spektren aus den mit den für dieselben isosbestischen Wellenlängen (λ_i1, λ_i2) des Hämoglobins und Oxihämoglobins der im Mess-Spektrum ermittelten zumindest zwei Messwerte (M_i1, M_i2) korrespondierenden Referenzwerten (R_i1, R_i2) eine lineare Referenzfunktion (F_R) generiert wird, c) dass aus der Hilfsfunktion (F_H) des Mess-Spektrums, auf welcher die besagten zumindest zwei Messwerte (M_i1, M_i2) für isosbestische Wellenlängen (λ_i1, λ_i2) liegen und aus der Referenzfunktion (F_R) der Referenz-Spektren, auf welcher die zu den zumindest zwei Messwerte (M_i1, M_i2) korrespondierenden zumindest zwei Referenzwerte (R_i1, R_i2) liegen eine lineare Korrekturfunktion (F_K) gebildet wird, mit welcher in einem korrigierten Mess-Spektrum eine ebenfalls lineare korrigierte Hilfsfunktion (F_Hk) erzeugt wird, die mit der Referenzfunktion (F_R) in den Referenz-Spektren identisch ist und d) dass die Sauerstoffsättigung des Blutes aus dem auf die korrigierte Hilfsfunktion (F_Hk) des korrigierten Mess-Spektrums umgerechneten anderen Messwert (M_a'') in Relation zu den Referenzwerten für Hämoglobin und Oxihämoglobin bei dieser Wellenlänge (λ_a) ermittelt wird.
Verfahren nach Anspruch 1, dadurch gekennzeichnet, a) dass drei besagte spektrale Messwerte (M_i1, M_i2, M_i3) bei für Hämoglobin und Oxihämoglabin isosbestischen Wellenlängen (λ_i1, λ_i2, λ_i3) und ein anderer Messwert (M_a) bei einer Wellenlänge (λ_a), bei welcher sich die Referenzwerte von Hämoglobin und Oxihämoglobin in den Referenz-Spektren möglichst weitgehend unterscheiden, logarithmisch erfasst werden, wobei aus zwei logarithmischen Messwerten (M_i1, M_i2) für isosbestische Wellenlängen (λ_i1, λ_i2) eine lineare Hilfsfunktion (F_H) generiert wird, b) dass in den Referenz-Spektren aus den mit den für dieselben isosbestischen Wellenlängen (λ_i1, λ_i2) des Hämoglobins und Oxihämoglobins der im Mess-Spektrum ermittelten Messwerte (M_i1, M_i2) korrespondierenden Referenzwerten (R_i1, R_i2) eine ebenfalls lineare Referenzfunktion (F_R) generiert wird, c) dass aus der Hilfsfunktion (F_H) des Mess-Spektrums und aus der Referenzfunktion (F_R) der Referenz-Spektren eine lineare Korrekturfunktion (F_K) gebildet wird, mit welcher im korrigierten Mess-Spektrum eine ebenfalls lineare korrigierte Hilfsfunktion (F_Hk) erzeugt wird, die mit der linearen Referenzfunktion (F_R) in den Referenz-Spektren identisch ist, d) dass die übrigen korrigierten spektralen Messwerte, d. h. der dritte spektrale Messwert (M_i3') bei für Hämoglobin und Oxihämoglobin isosbestischer Wellenlänge (λ_i3) und der andere Messwert (M_a') bei einer Wellenlänge (λ_a), bei welcher sich die Referenzwerte von Hämoglobin und Oxihämoglobin in den Referenz-Spektren möglichst weitgehend unterscheiden, mit einem konstanten Multiplikator beaufschlagt werden, welcher so bestimmt wird, dass der hierdurch korrigierte dritte spektrale Messwert (M_i3'') des korrigierten Mess-Spektrums mit dem korrespondierenden Referenzwert der Referenz-Spektren übereinstimmt und e) dass die Sauerstoffsättigung des Blutes an dem auf die korrigierte Hilfsfunktion (F_Hk) des korrigierten Mess-Spektrums umgerechneten anderen Messwert (M_a'') auf einer durch die Referenzwerte für Hämoglobin und Oxihämoglobin bei dieser Wellenlänge (λ_a) aufgespannten Skala von 0 bis 1 abgelesen wird.
Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass zur zweidimensionalen Darstellung der Sauerstoffsättigung des Blutes vier monochromatische Einzelbilder der spektralen Messwerte (M_i, M_a) erzeugt werden und dass die Sauerstoffsättigung gemäß der Schritte a) bis d) für jeden Bildpunkt ermittelt wird.