Genauso wie Strahlaufhärtungseffekte können auch Streueffekte unerwünschte Bildartefakte im rekonstruierten Tomographiebild einer durchstrahlten Schicht eines Patienten oder eines be- liebigen anderen Untersuchungsobjekts hervorrufen. Diese Bildartefakte spiegeln Strukturen vor, die keiner realen Vor- lage des Untersuchungsobjekts entsprechen und deshalb zu Fehlinterpretationen des Tomographiebilds verleiten. Beson- ders im medizinischen Sektor können solche Fehlinterpre- tationen gravierende Folgen bis hin zur Gefährdung des Lebens des Patienten haben.

Zur Unterdrückung des Streustrahlungsanteils in den mit einem Detektor gemessenen Strahlungsintensitätsmesswerten ist es bekannt, die das Untersuchungsobjekt durchstrahlende Röntgen- strahlung detektorseitig zu kollimieren. Kollimatoren werden in der Regel aus Wolfram hergestellt, das sich wegen seiner hohen Schwächung hierfür sehr gut eignet. Wolfram hat aller- dings den Nachteil, dass es sehr teuer ist. Dieser Kosten- nachteil fällt besonders dann ins Gewicht, wenn als Detektor eine Detektormatrix mit einer Vielzahl in mehreren überein- anderliegenden Detektorzeilen nebeneinander angeordneter De- tektorelemente verwendet wird. Bei solchen Detektoren muss mit zunehmender Zeilenzahl die Schachttiefe der für jedes einzelne Detektorelement vorgesehenen Kollimatorschächte ver- größert werden. Der konstruktive und der Materialaufwand wür- den ab einer gewissen Zeilenzahl des Detektors als nicht mehr akzeptabel angesehen.

Aufgabe der Erfindung ist es daher, mit geringerem Aufwand die Vermeidung streustrahlungsbedingter Bildartefakte bei mehrzeiligen Detektoren zu ermöglichen.

Bei der Lösung dieser Aufgabe geht die Erfindung nach einem ersten Aspekt von einem Verfahren zur Streustrahlungskorrek- tur von Strahlungsintensitätsmesswerten aus, welche in einer Röntgen-Computertomographieeinrichtung mittels einer in einem Tomographiemessfeld der Computertomographieeinrichtung lie- genden Detektormatrix erhalten werden, die eine Vielzahl in mehreren übereinanderliegenden Detektorzeilen nebeneinander angeordneter Detektorelemente aufweist.

Erfindungsgemäß ist dabei vorgesehen, dass zunächst mindes- tens eine Referenzverteilung der Streustrahlungsintensität in Zeilenrichtung der Detektormatrix ermittelt wird und dass sodann ausgehend von dieser mindestens einen Referenzvertei- lung ein Streustrahlungsanteil jedes Intensitätsmesswerts ermittelt wird und die Intensitätsmesswerte in Abhängigkeit von ihrem jeweiligen Streustrahlungsanteil korrigiert werden, wobei der Streustrahlungsanteil der Intensitätsmesswerte min- destens einer Teilanzahl der Detektorzeilen durch Rekursion in folgender Weise ermittelt wird : a) der Streustrahlungsanteil der Intensitätsmesswerte einer aktuellen Detektorzeile der Rekursion wird aus den Inten- sitätsmesswerten dieser aktuellen Detektorzeile und einem Primärstrahlüngsanteil der Intensitätsmesswerte einer vorher- gehenden Detektorzeile der Rekursion ermittelt, b) der Primärstrahlungsanteil der Intensitätsmesswerte der vorhergehenden Detektorzeile wird aus den Intensitätsmess- werten dieser vorhergehenden Detektorzeile und deren Streu- strahlungsanteil ermittelt, und c) als Streustrahlungsanteil der Intensitätsmesswerte einer ersten Detektorzeile der Rekursion werden Intensitätswerte aus der Referenzverteilung der Streustrahlungsintensität ver- wendet.

Wenn hier von Primärstrahlung die Rede ist, so wird darunter derjenige Strahlungsanteil der auf die Detektorelemente ein- fallenden Gesamtstrahlung verstanden, der ungestreut, also auf direktem Weg von der Strahlungsquelle der Computertomo- graphieeinrichtung zu der Detektormatrix gelangt. Als Tomo- graphiemessfeld wird dann ein mit Detektorelementen bestück- ter Messbereich verstanden, in dem die gemessene Gesamtstrah- lung einen Primärstrahlungsanteil enthält. In der Regel wird das Tomographiemessfeld durch eine quellenseitige Blendenan- ordnung festgelegt.

Bei der erfindungsgemäßen Lösung erfolgt für alle Detektor- zeilen eine rechnerische Abschätzung des Streustrahlungs- anteils anhand mindestens einer Referenzverteilung. Hierdurch kann auf teure Kollimatorschächte verzichtet werden. Die Re- kursion, die zumindest für eine Teilanzahl der Detektor- zeilen angewendet wird, bietet die Grundlage, um ein sich über die Detektorzeilen hinweg änderndes Profil des Streu- strahlungsanteils zu berücksichtigen.

Bei einer ersten Ausgestaltung des erfindungsgemäßen Ver- fahrens wird die mindestens eine Referenzverteilung der Streustrahlungsintensität aus Referenzintensitätsmesswerten gewonnen, welche durch Strahlungsintensitätsmessung außerhalb des Tomographiemessfelds erhalten werden. Hierbei wird aus- genutzt, dass außerhalb des Tomographiemessfelds keine Pri- märstrahlung auftritt. Dort angeordnete Messelemente detek- tieren deshalb ausschließlich Streustrahlung. Hieraus kann leicht eine Verteilung der Streustrahlung in Zeilenrichtung ermittelt werden, die dann als Referenzverteilung verwendet wird.

Zweckmäßigerweise wird die Strahlungsintensitätsmessung ober- halb einer ersten Detektorzeile der Detektormatrix oder/und unterhalb einer letzten Detektorzeile der Detektormatrix durchgeführt werden.

Im allgemeinen kann das Ortsprofil der Streustrahlung durch eine vergleichsweise niederfrequente Funktion dargestellt werden. Deshalb genügt es, wenn in Zeilenrichtung nur in ei- nem relativ groben Raster Messwerte für die Streustrahlung aufgenommen werden. Mit anderen Worten werden die Referenz- intensitätsmesswerte bevorzugt an in Zeilenrichtung der De- tektormatrix im Abstand voneinander liegenden Messorten ge- wonnen, deren Anzahl kleiner, insbesondere wesentlich kleiner als die Anzahl der Detektorelemente pro Detektorzeile ist.

Die Referenzverteilung der Streustrahlungsintensität kann dann leicht durch Interpolation der Referenzintensitätsmess- werte gewonnen werden.

Es kann sogar eine Referenzverteilung durch Strahlungsinten- sitätsmessung oberhalb der ersten Detektorzeile der Detektor- matrix und eine weitere Referenzverteilung durch Strahlungs- intensitätsmessung unterhalb der letzten Detektorzeile der Detektormatrix gewonnen werden.

Die Rekursion sollte zweckmäßigerweise zumindest in einer randseitigen Detektorzeile der Detektormatrix begonnen wer- den. Hier wird die Annahme, dass sich die Streustrahlungs- intensität außerhalb des Tomographiemessfelds und die Streu- strahlungsintensität in einer randseitigen Detektorzeile- wenn überhaupt-nur unwesentlich unterscheiden, in der Regel zutreffen. Deshalb wird der Fehler vernachlässigbar sein, der sich ergibt, wenn als Streustrahlungsanteil der Intensitäts- messwerte der randseitigen Detektorzeile Intensitätswerte aus der Referenzverteilung verwendet werden.

Bei einer zweiten Ausgestaltung des erfindungsgemäßen Ver- fahrens wird die mindestens eine Referenzverteilung der Streustrahlungsintensität unter Verwendung der Intensitäts- messwerte mindestens einer Detektorzeile der Detektormatrix berechnet. Insbesondere kann dabei die Referenzverteilung auf Grundlage eines mathematischen Faltungsmodells berechnet wer- den. Für eine Computertomographieeinrichtung mit einzeilig

angeordneten Detektorelementen ist ein solches Faltungsmodell beispielsweise aus B. Ohnesorge :"Untersuchungen der Scatter- Korrektur in Elektronenstrahl-Computertomographen", Lehrstuhl für Nachrichtentechnik der Universität Erlangen-Nürnberg, Diplomarbeit, 1994, bekannt. Durch Anpassung dieses Faltungs- modells an eine mehrzeilige Detektormatrix kann die Streu- strahlungsverteilung für eine Detektorzeile der Matrix rech- nerisch aus den erhaltenen Intensitätsmesswerten dieser De- tektorzeile abgeschätzt werden.

Es könnte eingewendet werden, dass grundsätzlich in allen Detektorzeilen die Streustrahlungsverteilung jeweils mit Hil- fe des obigen Faltungsmodells berechnet werden könnte und dass dann eine Rekursion überflüssig sei. Faltungsoperationen können jedoch sehr rechenintensiv sein. Die Anwendung der Rekursion für zumindest eine Teilanzahl der Detektorzeilen erlaubt es hingegen, den Rechenaufwand in vertretbaren Gren- zen zu halten und gleichzeitig mögliche Änderungen der Streu- strahlungsverteilung von Detektorzeile zu Detektorzeile be- rücksichtigen zu können.

Zweckmäßigerweise wird die Referenzverteilung unter Verwen- dung der Intensitätsmesswerte einer mittleren Detektor-zeile der Detektormatrix berechnet werden und die Rekursion zumin- dest in dieser mittleren Detektorzeile zu oberen und unteren Detektorzeilen hin begonnen werden. Es versteht sich jedoch, dass auch anhand der Intensitätsmesswerte einer anderen De- tektorzeile, insbesondere sogar einer randseitigen Detektor- zeile, die Referenzverteilung berechnet werden kann.

Um die Qualität der erhaltenen Ergebnisse für den Streustrah- lungsanteil der Intensitätsmesswerte zu verbessern, kann die Rekursion nach einer Teilanzahl von Detektorzeilen beendet werden und in einer nachfolgenden Detektorzeile eine weitere Rekursion gestartet werden. Dabei kann die weitere Rekursion auf Grundlage derselben oder einer anderen Referenzverteilung der Streustrahlungsintensität gestartet werden.

Wenn das Untersuchungsobjekt vergleichsweise kontrastreiche Strukturen enthält, können sich die Intensitätsmesswerte von Detektorzeile zu Detektorzeile oder/und innerhalb einer De- tektorzeile von Detektorelement zu Detektorelement relativ stark ändern, und zwar nicht aufgrund einer rapiden Änderung der Streustrahlung (die sich-wie bereits gesagt-örtlich im Regelfall nur vergleichsweise langsam ändert) sondern auf- grund örtlich wechselnder Schwächungseigenschaften des durch- strahlten Materials. Damit derartige Unstetigkeiten in der gemessenen Gesamtintensität die letztendlich zur Korrektur der Intensitätsmesswerte herangezogenen Streustrahlungsantei- le nicht wesentlich verfälschen, werden die nach Durchführung der Rekursion ermittelten Streustrahlungsanteile in Spalten- und gewünschtenfalls auch in Zeilenrichtung der Detektormat- rix vorzugsweise tiefpassgefiltert. Durch die Tiefpassfilte- rung werden aus den rekursiv ermittelten Streustrahlungsan- teilen solche Intensitätsänderungen herausgefiltert, die eine vergleichsweise hohe Ortsfrequenz besitzen. Diese beruhen regelmäßig auf Änderungen der Schwächungseigenschaften. Die gefilterten Streustrahlungsanteile geben so sehr gut das nie- derfrequente Profil der Streustrahlung wieder. Die Intensi- tätsmesswerte werden dann in Abhängigkeit von ihrem jeweili- gen gefilterten Streustrahlungsanteil korrigiert.

Eine Verfeinerung der erhaltenen Abschätzung für den Streu- strahlungsanteil der Intensitätsmesswerte ist möglich, wenn, für jeden Intensitätsmesswert ausgehend von zwei verschiede- nen Referenzverteilungen zwei Werte des Streustrahlungsan- teils ermittelt werden und die Intensitätsmesswerte entspre- chend einem jeweiligen gemittelten Streustrahlungsanteil kor- rigiert werden.

Unabhängig von der rekursiven Ermittlung der Streustrahlungs- anteile soll auch der Gedanke, die Referenzverteilung durch Strahlungsintensitätsmessung außerhalb des Tomographiemess- felds zu ermitteln, selbständigen Schutz im Rahmen der Erfin-

dung genießen. Die Erfindung sieht nach einem zweiten Aspekt daher ferner ein Verfahren zur Streustrahlungskorrektur von Strahlungsintensitätsmesswerten vor, welche in einer Röntgen- Computertomographieeinrichtung mittels einer in einem Tomo- graphiemessfeld der Computertomographieeinrichtung liegenden Detektormatrix erhalten werden, die eine Vielzahl in mehreren übereinanderliegenden Detektorzeilen nebeneinander angeordne- ter Detektorelemente aufweist. Erfindungsgemäß ist dabei vor- gesehen, dass zunächst mindestens eine Referenzverteilung der Streustrahlungsintensität in Zeilenrichtung der Detektormat- rix aus Referenzintensitätsmesswerten gewonnen wird, welche durch Strahlungsintensitätsmessung außerhalb des Tomographie- messfelds erhalten werden, und dass sodann ausgehend von die- ser mindestens einen Referenzverteilung ein Streustrahlungs- anteil jedes Intensitätsmesswerts ermittelt wird und die In- tensitätsmesswerte in Abhängigkeit von ihrem jeweiligen Streustrahlungsanteil korrigiert werden.

Um den Streustrahlenanteil der Intensitätsmesswerte abzu- schätzen, kann auch hier die früher erläuterte Rekursion mit den Schritten a) bis c) zur Anwendung kommen. Es ist aber auch denkbar, für jeden Intensitätsmesswert einen Intensi- tätswert aus der Referenzverteilung der Streustrahlungsinten- sität als Streustrahlungsanteil zu verwenden. Dabei wird für jede Detektorzeile einfach die Referenzverteilung als Streu- strahlungsverteilung unmittelbar übernommen. In Fällen, wo sich die Streustrahlungsintensität über die Detektormatrix hinweg tatsächlich nur wenig ändert, können so bereits recht gute Ergebnisse erzielt werden. Ist dagegen mit merklichen Änderungen der Streustrahlungsintensität zu rechnen, wird die rekursive Vorgehensweise vorgezogen werden.

Es versteht sich, dass das Verfahren nach dem zweiten Aspekt durch weitere Merkmale des Verfahrens nach dem ersten Aspekt ausgestaltet werden kann.

Gegenstand der Erfindung ist außerdem eine Röntgen-Computer- tomographieeinrichtung, welche zur Durchführung des Verfah- rens nach dem ersten oder/und zweiten Aspekt ausgelegt ist.

Insbesondere kann bei dieser Computertomographieeinrichtung eine außerhalb des Tomographiemessfelds angeordnete Zusatzde- tektoranordnung zur Gewinnung der Referenzintensitätsmesswer- te vorgesehen sein. Die Zusatzdetektoranordnung kann oberhalb einer ersten Detektorzeile der Detektormatrix oder/und unter- halb einer letzten Detektorzeile der Detektormatrix mehrere in Zeilenrichtung der Detektormatrix im Abstand voneinander angeordnete Zusatzdetektorelemente aufweisen, deren jedes einen der Referenzintensitätsmesswerte liefert. Die Anzahl der Zusatzdetektorelemente in Zeilenrichtung der Detektormat- rix ist dabei vorzugsweise kleiner, insbesondere wesentlich kleiner als die Anzahl der Detektorelemente pro Detektorzei- le.

Die Erfindung wird im Folgenden anhand der beigefügten Zeich- nungen näher erläutert. Es stellen dar : Fig. 1 schematisch eine Ausführungsform eines erfindungsge- mäßen CT-Scanners mit mehrzeiliger Detektormatrix, und Fig. 2 schematisch eine Draufsicht auf die Detektormatrix bei Betrachtung in Pfeilrichtung II in Fig. 1.

Der in den Figuren gezeigte CT-Scanner umfasst eine Röntgen- quelle 10 sowie eine Detektoranordnung 12. Die Röntgenquelle 10 sendet Röntgenstrahlung fächerförmig aus, wie bei 14 ange- deutet. Ein im Strahlengang zwischen der Röntgenquelle 10 und der Detektoranordnung 12 angeordnetes Untersuchungsobjekt 16 wird von der Röntgenstrahlung durchdrungen. Die Detektoran- ordnung 12 detektiert die hinter dem Untersuchungsobjekt 16 vorhandene Röntgenstrahlung. Speziell umfasst die Detektoran- ordnung 12 eine Detektormatrix 18 aus einer Vielzahl von De- tektorelementen 20, welche auf mehrere übereinanderliegende

Zeilen verteilt sind und in jeder Zeile in Richtung eines Fächerwinkels ß nebeneinander angeordnet sind. In Fig. 2 sind beispielhaft vier solcher Detektorzeilen gezeigt ; es versteht sich jedoch, dass die Anzahl der Detektorzeilen hiervon be- liebig abweichen kann und beispielsweise stattdessen 8,16 oder 24 betragen kann. Die Größe des Strahlenfächers 14 in Richtung des Fächerwinkels ß ist mittels einer Blendenanord- nung 22 einstellbar, welche zwischen der Röntgenquelle 10 und dem Untersuchungsobjekt 16 angeordnet ist. In Spaltenrichtung der Detektormatrix 18, d. h. in einer Richtung z in Fig. 2, erfolgt ebenfalls eine Begrenzung der von der Röntgenquelle 10 ausgesandten Strahlung durch eine vergleichbare Blenden- anordnung (nicht gezeigt). Die Blendenanordnung 22 und die soeben angesprochene z-Blendenanordnung definieren im Bereich der Detektoranordnung 12 ein Tomographiemessfeld, innerhalb dessen Primärstrahlung detektierbar ist, die auf geradem Weg von der Röntgenquelle 10 auf die Detektoranordnung 12 trifft, ohne in dem Untersuchungsobjekt 16 gestreut zu werden. Die Detektormatrix 18 liegt vollständig innerhalb dieses Tomogra- phiemessfelds. Jede Position in Richtung des Fächerwinkels ß, an der sich ein Detektorelement 20 befindet, entspricht einem Projektionskanal.

Jedes Detektorelement 20 detektiert die in seinem Raumbereich einfallende Strahlung und liefert ein entsprechendes Intensi- tätsmesssignal IG (n, k) an eine elektronische Auswerte-und Rekonstruktionseinheit 24. Dabei steht der Index n für die Nummer der Zeile der Detektormatrix 18, in der sich das betreffende Detektorelement 20 befindet, k repräsentiert die Kanalnummer. Die Auswerte-und Rekonstruktionseinheit 24 führt an den eingehenden Intensitätsmesssignalen Is (n, k) zunächst eine Streustrahlungskorrektur durch, indem sie von den Intensitätsmesssignalen IG (n, k) einen Streustrahlungsan- teil Is (n, k) subtrahiert. Es verbleibt ein Primärstrahlungs- anteil Ip (n, k), der für die Intensität der auf das jeweilige Detektorelement 20 einfallenden Primärstrahlung repräsentativ ist. Sodann ermittelt die Auswerte-und Rekonstruktionsein-

heit 24 aus den Intensitätswerten Ip (n, k) Schwächungswerte, die sie zur Rekonstruktion eines auf einem Monitor 26 ange- zeigten Tomographiebilds der durchstrahlten Schicht des Un- tersuchungsobjekts 16 verwendet. Es versteht sich, dass der CT-Scanner zur Rekonstruktion des Tomographiebilds Projektio- nen aus einer Vielzahl unterschiedlicher Richtungen benötigt.

Zu diesem Zweck ist die Röntgenquelle 10 in Pfeilrichtung 28 um das Untersuchungsobjekt 16 bewegbar.

Um die Streustrahlungskorrektur durchführen zu können, ist der CT-Scanner dazu ausgelegt, zunächst eine Referenzvertei- lung der Streustrahlungsintensität in Zeilenrichtung zu er- mitteln. Diese Referenzverteilung gibt für jeden Kanal k ei- nen Referenzwert Isref (k) für die Streustrahlungsintensität an. Zur Ermittlung der Referenzverteilung umfasst die Detek- toranordnung 12 zusätzlich zur Detektormatrix 18 eine Mehr- zahl von Zusatzdetektorelementen 30 (siehe Fig. 2). Diese liegen außerhalb des Tomographiemessfelds und werden demge- mäss nicht von Primärstrahlung sondern von ausschließlich von Streustrahlung getroffen. Die Zusatzdetektorelemente 30 er- lauben es demnach, messtechnisch Informationen über die In- tensität der Streustrahlung zu erhalten. Auch die Zusatzde- tektorelemente 30 sind mit der Auswerte-und Rekonstruktions- einheit 24 verbunden und liefern ihre Messsignale an selbige.

Die Zusatzdetektorelemente 30 sind oberhalb der in z-Richtung obersten oder/und unterhalb der in z-Richtung untersten Zeile der Detektormatrix 18 angeordnet. Da sich die örtliche Ver- teilung der Streustrahlung im allgemeinen durch eine ver- gleichsweise niederfrequente Funktion beschreiben lässt, ge- nügt in Zeilenrichtung ein grobes Raster der Zusatzdetektor- elemente 30, weshalb im Vergleich zu der pro Zeile vorhande- nen Anzahl von Detektorelementen 20 vorzugsweise nur eine wesentlich geringere, beispielsweise um eine Größenordnung kleinere Anzahl von Zusatzdetektorelementen 30 in Zeilenrich- tung vorgesehen ist. Durch Interpolation ermittelt die Aus- werte-und Rekonstruktionseinheit 24 dann aus den von den

Zusatzdetektorelementen 30 gelieferten Messsignalen die Refe- renzverteilung Isref (k). Die Zusatzdetektorelemente 30 sind zweckmäßigerweise in Zeilenrichtung in gleichmäßigen Abstän- den verteilt ; dies ist jedoch nicht zwingend. Selbstverständ- lich ist es nicht ausgeschlossen, eine der Anzahl der Detek- torelemente 20 gleiche Anzahl von Zusatzdetektorelementen 30 in Zeilenrichtung vorzusehen.

Für die in der Detektorzeile n (n=l,.., L) gemessene Gesamtin- tensität IG (n, k) im Kanal k (k=l,..., N) gilt : Is (n, k) = Ip (n, k) + Is (n, k) (1) Ziel der in der Auswerte-und Rekonstruktionseinheit 24 durchgeführten Streustrahlungskorrektur ist es zunächst, den Streustrahlungsanteil Is (n, k) möglichst genau abzuschätzen, um anschließend möglichst genaue Werte für den Primärstrah- lungsanteil Ip (n, k) zur Verfügung zu haben, die der Bildre- konstruktion zugeführt werden können.

Die Streustrahlungsabschätzung beginnt in der obersten oder der untersten Detektorzeile, je nachdem, ob die Referenzver- teilung Isref (k) aus den Messsignalen oberhalb oder unterhalb der Detektormatrix 18 liegender Zusatzdetektorelemente 30 gewonnen wurde. Im folgenden wird davon ausgegangen, dass in der obersten Detektorzeile begonnen wird. Zur vereinfachten Notation wird dabei die Kanalnummer nicht mehr explizit ange- geben. Die folgenden Überlegungen gelten jedoch für beliebige Winkelpositionen im Strahlenfächer und demnach für beliebige Kanalnummern. Für die oberste Detektorzeile gilt dann : IG (1)-IP (l) + Is (l) (2) Zur Ermittlung der Primärstrahlungsintensität Ip (l) in der obersten (ersten) Detektorzeile wird die Annahme gemacht, dass sich Isref und der Streustrahlungsanteil Is (l) der ersten Detektorzeile-wenn überhaupt-so nur vernachlässigbar un-

terscheiden. Die Primärstrahlungsintensität Ip (l) kann daher auf einfache Weise wie folgt berechnet werden : Ip (l) = IG (1)-ISreE (3) Die Primärstrahlungsintensitäten in allen weiteren Detektor- zeilen können nun analog unter der Annahme ermittelt werden, dass die Primärstrahlungsintensität Ip (n-l) der n-l-ten De- tektorzeile näherungsweise der Primärstrahlungsintensität Ip (n) der n-ten Zeile entspricht. Unter dieser Annahme kann die Streustrahlungsintensität Is (n) in der n-ten Zeile aus der tatsächlich gemessenen Gesamtintensität IG (n) in dieser Zeile und der Primärstrahlungsintensität Ip (n-l) in der vor- hergehenden Zeile n-1 rekursiv wie folgt berechnet werden : Is (n) = IG (n)-Ip (n-l) (4) Sodann kann die Primärstrahlungsintensität Ip (n) der n-ten Zeile abgeschätzt werden gemäß : Ip (n) = IG (n)-Is (n) (5) Die Annahme Ip (n-l) == Ip (n) ist bei kontrastarmen Strukturen im allgemeinen gerechtfertigt. Enthält das Untersuchungsob- jekt 16 jedoch kontrastreiche Strukturen, wie beispielsweise Knochen, so können zwischen aufeinander-folgenden Zeilen o- der/und Kanälen signifikante Änderungen der gemessenen Ge- samtintensität auftreten. Damit sich solche Signalunstetig- keiten beim Übergang von Zeile n-1 zu Zeile n in der obigen Rekursion nicht auf die Berechnung der Streustrahlungsinten- sitäten übertragen und daher die Ip (n)-Werte verfälschen, werden die geschätzten Werte Is (n) der Streustrahlungsinten- sität beispielsweise mit einem Medianfilter wählbarer Länge tiefpassgefiltert. Durch die Tiefpassfilterung werden die angesprochenen Unstetigkeiten beseitigt. Die gefilterten Is (n)-Werte geben dann eine sehr gute Abschätzung der tat- ., sächlichen Streustrahlungsintensität wieder. Aus den gefil-

terten Is (n)-Werten werden anschließend durch Einsetzen in obige Gleichung (5) neue Ip (n)-Werte berechnet, die für die Bildrekonstruktion herangezogen werden.

Die Tiefpassfilterung kann als eindimensionale Filterung in z-Richtung oder auch als zweidimensionale Filterung in z-und Zeilenrichtung durchgeführt werden.

Sind Zusatzdetektorelemente 30 ober-und unterhalb der Detek- tormatrix 18 vorgesehen, so können zwei Referenzverteilungen ISreftl und ISreft2 ermittelt werden, nämlich eine (ISreftl) aus den Messsignalen der oberhalb der Detektormatrix 18 liegenden Zusatzdetektorelemente 30 und die andere (ISreft2) aus den Messsignalen der unterhalb der Detektormatrix 18 liegenden Zusatzdetektorelemente 30. Es kann dann obiges Verfahren zur rekursiven Abschätzung der Primärstrahlungsintensitäten zwei- mal durchgeführt werden, nämlich einmal beginnend in der o- bersten Detektorzeile auf Grundlage der Referenzverteilung ISreftl und einmal beginnend in der untersten Detektorzeile auf Grundlage der Referenzverteilung ISrefB2. Es werden so für jedes Detektorelement 20 zwei Werte Ip, l und Ip, 2 der Primär- strahlungsintensität erhalten, die anschließend gemittelt werden. Für die Bildrekonstruktion werden dann die gemittel- ten Intensitätswerte verwendet.

In manchen Fällen kann es bereits genügen, die anhand der Zusatzdetektorelemente 30 gewonnene Referenzverteilung Isref als Modell für die Streustrahlungsverteilung aller Detektor- zeilen der Detektormatrix 18 zu verwenden. Die Primärstrah- lungsintensitäten Ip (n) lassen sich dann einfach wie folgt berechnen : Ip(n) = IG (n)-Isref (6) Es ist auch denkbar, die Rekursion nicht über alle Detektor- zeilen hinweg fortzuführen, sondern sie nach einer Teilanzahl der Detektorzeilen, beispielsweise nach jeder zweiten, drit-

ten oder vierten Detektorzeile oder nach der Hälfte der De- tektorzeilen, abzubrechen und sodann in einer neuen Detektor- zeile eine neue Rekursion zu starten. In dieser neuen Detek- torzeile wird dann analog zu Gleichung (3) wieder die Annahme gemacht, dass die Streustrahlungsverteilung dieser Zeile der Referenzverteilung Isref entspricht. Es ist sogar vorstellbar, beim Neustart der Rekursion von einer anderen Referenzvertei- lung auszugehen. Bei obigem Beispiel mit Zusatzdetektorele- menten 30 ober-und unterhalb der Detektormatrix könnte es etwa sinnvoll sein, für die obere Hälfte der Detektorzeilen eine Rekursion auf Grundlage der Referenzverteilung Isrefil durchzuführen und für die untere Hälfte der Detektorzeilen eine Rekursion auf Grundlage der Referenzverteilung ISref, 2 durchzuführen, insbesondere dann, wenn die Detektormatrix 18 eine große Anzahl von Zeilen aufweist, beispielsweise 16,24 oder 32.

Die Referenzverteilung Isref kann auch auf andere Weise als mit Hilfe der Zusatzdetektorelemente 30 ermittelt werden. So ist es beispielsweise möglich, aus den Intensitätsmesswerten IG (n) einer Zeile der Detektormatrix 18 die zugehörige Streu- strahlungsverteilung Is (n) dieser Zeile rechnerisch zu bestimmen. Diese Streustrahlungsverteilung Is (n) kann dann als Referenzverteilung Isref verwendet werden, um für die üb- rigen Zeilen der Detektormatrix 18 den Streustrahlungsanteil der Intensitätsmesswerte dieser Zeilen mittels obigen Rekur- sionsverfahrens abzuschätzen.

Zur rechnerischen Ermittlung einer Streustrahlungsverteilung aus Intensitätsmesswerten ist aus der weiter oben zitierten Literatur von B. Ohnesorge ein Faltungsmodell für ein einzei- liges Detektorsystem bekannt. Diesem Modell liegt die Idee zugrunde, dass die prinzipiellen Streuwinkelabhängigkeiten der differentiellen Wirkungsquerschnitte und Streuenergien von Compton-und Raleigh-Streuung die Annahme rechtfertigen, dass die Streubeiträge in einem Detektorkanal k, der zu einem Fächerwinkel ßk zugehörig ist, mit dem Winkelabstand im Fä-

cher (ß-ßk) abnehmen. (In der Herleitung werden nur Ein- fachstreuprozesse berücksichtigt.) Eine zur Beschreibung ver- wendbare"Abstandsfunktion"G (ß) hat dann ein Maximum bei ß=ßk und überstreicht den Winkelbereich (-ßmax+ßk, ßmax+ßk)- Eine vom.. Fächerwinkel ß abhängige Streustrahlungsverteilung Isc (ß) ergibt sich dann wie folgt : Dabei bezeichnen Cm eine Maschinenkonstante und f (Azsl) eine schichtdickenabhängige Gewichtung. Isc, forw (ß) ist eine im Mo- dell der Einfachstreuung berechnete Vorwärts-Streustrah- lungsintensität mit Ksc. forw ist eine Proportionalitätskonstante, Io die Intensität der ungeschwächten Strahlung und I (ß) die im Fächerwinkel ß vom Detektorsystem gemessene Strahlungsintensität. In der Faltungsgleichung (7) werden durch den Faltungskern G (ß) die Streubeiträge aller Strahlen im Fächer zu allen Detektorele- menten berücksichtigt. Er wird üblicherweise als Abstandskern angegeben : A ist ein Parameter, mit dem die Breite gesteuert werden kann. Er kann empirisch aus Bildoptimierungen oder aus einem Vergleich von im Faltungsmodell berechneten und simulierten Streustrahlungsverteilungen bestimmt werden.

Für die Funktion R (ß) gilt :

Weitere Informationen zu obigem Faltungsmodell für Einzeilen- Detektorsysteme können der Literatur von B. Ohnesorge entnom- men werden.

Dieses bekannte Modell kann nun im Rahmen der Erfindung zur Anpassung an mehrzeilige bzw. flächenartige Detektoren, wie beispielsweise in Fig. 2 gezeigt, modifiziert werden. Wegen der Rotationssymmetrie der differentiellen Wirkungsquer- schnitte hinsichtlich der Fächerkoordinate ß und der zeilen- koordinate Zn (Zn = (L/2-n) Az ; (n=l,..., L)) kann Gleichung (7) leicht erweitert werden. (Az stellt die Zeilenhöhe dar.) Die Streustrahlungsintensität Isc (ß, zn) ist dann gegeben durch : Dabei haben CM und f (asz) die gleiche Bedeutung wie oben.

Isc, foa (ß, Zn) ist wiederum die im Modell der Einfachstreuung berechnete Vorwärts-Streustrahlungsintensität mit I (ß, Zn) bezeichnet die im Fächerwinkel ß der n-ten Detektor- zeile gemessene Strahlungsintensität. Für R (ß, Zn) gilt : R (ß,zn)=1, falls ß#[-ßmax,ßmax] und 1#n#L ; 0 sonst (13) Der Abstandskern lautet nun:

Darin bezeichnet A'wiederum den Breitenparameter, ß2 + (Zn/Rfd) 2 misst den Abstand vom Detektorursprung zum Detektor- element im Fächerwinkel ß der n-ten Detektorzeile und Rfd be- zeichnet den Abstand zwischen Fokus und Detektor des CT- Scanners.

Die Streustrahlungsverteilung Isc (ß, zn) kann in vorstehender Weise für eine beliebige Zeile der Detektormatrix 18 berech- net werden. Es wird empfohlen, sie für eine mittlere Detek- torzeile zu berechnen. Die so berechnete Streustrahlungsver- teilung wird sodann als Referenzverteilung Isref für die Re- kursion verwendet. Die Rekursion wird in der Zeile gestartet, aus deren Intensitätsmesswerten die Referenzverteilung be- rechnet wurde. Im Fall einer mittleren Detektorzeile wird sowohl eine Rekursion zu oberen Detektorzeilen als auch eine Rekursion zu unteren Detektorzeilen hin gestartet. Wird nach einer Teilanzahl von Detektorzeilen die Rekursion unterbro- chen, so wird in einer neuen Zeile eine neue Rekursion vor- zugsweise mit einer neuen Referenzverteilung begonnen, die mit Hilfe obigen Faltungsmodells aus den Intensitätsmesswer- ten dieser neuen Zeile berechnet wurde. Auf diese Weise lässt sich eine hohe Qualität bei der Abschätzung der genauen Streustrahlungsanteile erzielen.