DE69106954T2 - Verfahren zur Ermittlung der Reziprozitätsabweichungsfunktion für einen Röntgenfilm. - Google Patents

Description

• Die Erfindung betrifft die radiologischen Systeme, die einen Röntgenf ilm verwenden, und insbesondere in derartigen Systemen ein Verfahren, das die Wirkung der Nichtreziprozität des Röntgenfilms zu bestimmen ermöglicht.
• Ein radiologisches System enthält im wesentlichen eine Röntgenröhre und einen Empfänger für eine derartige Strahlung, zwischen denen ein zu untersuchender Gegenstand wie etwa ein Körperteil eines Patienten angeordnet ist. Der Bildempfänger, der beispielsweise ein Film-Bildschirm-Paar ist, liefert nach einer geeigneten Belichtungsdauer und der Entwicklung des Films ein Bild des Gegenstandes. Damit das Bild des Gegenstandes unter besseren Bedingungen ausgewertet werden kann, müssen die verschiedenen Punkte, die dieses Bild ergeben, untereinander einen ausreichenden Kontrast besitzen, d.h., daß die Schwärzung des Röntgenfilms korrekt sein muß, und dies trotz der Opazitätsunterschiede, die der geröntgte Gegenstand aufweisen kann, von einer Röntgenaufnahme zur folgenden.
• Die Schwärzung steht mit der Energiemenge der auf das Film-Bildschirm-Paar auftreffenden Strahlung in Beziehung, d.h. mit dem Produkt der Intensität der Strahlung, der der Röntgenfilm unterworfen ist, oder der "Film"-Dosisrate, mit der Zeit, in der der Film dieser Strahlung ausgesetzt ist. Um daher eine konstante Schwärzung des Films von einer Röntgenaufnahme zur folgenden zu erhalten, ist es bekannt, während der Untersuchung die auf den Film auftreffende Energie mittels einer Detektorzelle zu messen, die im allgemeinen vor dem Empfänger angeordnet ist, der für die Röntgenstrahlung empfindlich ist und der einen zur "Film"-Dosisrate proportionalen Strom liefert. Dieser Strom wird ab dem Beginn der Aufnahme in einer Integratorschaltung integriert, die einen Wert liefert, der im Verlauf der Aufnahme ansteigt. Dieser ansteigende Wert wird während der Aufnahmezeit mit einem festen Sollwert verglichen, der im voraus in Abhängigkeit von den Eigenschaften des Films erstellt worden ist. Das Ende der Aufnahmezeit ist durch den Zeitpunkt bestimmt, in welchem der Vergleich ergibt, daß der die auf den Film auftreffende Energie darstellende Wert gleich dem Sollwert ist.
• In dem Fall, in dem der Röntgenfilm direkt der Röntgenstrahlung unterworfen ist und die Veränderung der Aufnahmezeit von einer Untersuchung zur nächsten ausreichend gering ist, wird von einer Aufnahme zur nächsten eine konstante Schwärzung des Films unabhängig von der Dauer der Belichtungszeit S erhalten, vorausgesetzt, daß das Produkt aus der Aufnahmezeit S mit der Dosisrate F konstant ist, d.h., daß der aus der Integration sich ergebende Wert konstant bleiben muß. Dies ist nur dann richtig, falls die Eigenschaften des Films dem Reziprozitätsgesetz gehorchen, das angibt, daß die optische Dichte des Films zum Produkt F × S Proportional ist, und falls das Ansprechverhalten des Films von der Qualität des auftreffenden Röntgenstrahlbündels unabhängig ist.
• Dieses Reziprozitätsgesetz ist in dem Fall, in dem die Veränderung der Aufnahmezeit groß ist, nicht mehr gültig.
• Außerdem hängt in dem Fall, in dem der Röntgenfilm einem Verstärkerbildschirm zugeordnet ist, die Schwärzung des Films von der Qualität des Spektrums ab. Das Ansprechverhalten des Bildschirms hängt nämlich von der energetischen Verteilung des Spektrums der empfangenen Strahlung ab, was bedeutet, daß es gegenüber der spektralen Härtung und gegenüber der Spannungsänderung der Röntgenröhre empfindlich ist.
• Schließlich gibt es bestimmte Anwendungen, für die es ungünstig ist, wenn die Detektorzelle vor dem Film angeordnet ist, beispielsweise bei der Mammographie, weil die Strahlungsenergie derart ist, daß die Detektorzelle dann auf dem Film sichtbar wäre. In diesem Fall ist sie hinter dem Bildempfänger angeordnet, dies ruft jedoch eine zusätzliche Schwierigkeit hervor, weil das von der Detektorzelle aufgenommene Signal dasjenige ist, das nicht zur Schwärzung des Films beigetragen hat. Daraus folgt, daß die von der Detektorzelle ausgeführte Messung im allgemeinen nicht die auf den Röntgenfilm auftreffende Belichtung repräsentiert.
• Die Abweichung vom Reziprozitätsgesetz, die sich je nach Filmtyp verändert, repräsentiert die relative Veränderung der Belichtung, die notwendig ist, um eine konstante optische Dichte zu erhalten, wenn sich die Aufnahmezeit S verändert, während das Röntgenstrahlungsspektrum konstant ist. Dies äußert sich in der Tatsache, daß für den Erhalt derselben optischen Dichte des Films die Belichtung für eine Aufnahmezeit S = 0,1 Sekunden beispielsweise 1 sein muß, für S = 1 Sekunde 1,3 sein muß und für S = 4 Sekunden 2 sein muß.
• Diese Abweichung vom Reziprozitätsgesetz ist durch das unter dem Namen Schwarzschild-Effekt bekannte Phänomen bedingt. Dieser Effekt ist insbesondere in dem Buch mit dem Titel "CHIMIE ET PHYSIQUE PHOTOGRAPHIQUES" von Pierre GLAFKIDES - 4. Ausgabe, Seiten 234 bis 238, und herausgegeben von PUBLICATIONS PHOTO-CINÉMA Paul MONTEL.
• Um diese Abweichung vom Reziprozitätsgesetz zu berücksichtigen, sind verschiedene Lösungen vorgeschlagen worden, wobei eine von ihnen in dem französischen Patent 2 584 504 (EP-A-0 208 607) beschrieben ist. In diesem Patent ist vorgesehen, den integrierten Wert des von der Detektorzelle gelieferten Signals mit einem Sollwert zu vergleichen, der sich während der Aufnahme gemäß einem vorgegebenen Gesetz verändert. Genauer wird ausgehend vom Beginn jeder Aufnahmezeit zur Differenz der Werte des integrierten Signals und des Sollwerts ein zusätzlicher Wert hinzugefügt, der in Abhängigkeit von der Zeit gemäß einem im voraus bestimmten Gesetz, beispielsweise exponentiel wächst.
• Dieses im voraus bestimmte Gesetz, sei es exponentiel oder anders, berücksichtigt die Abweichung vom Reziprozitätsgesetz nur unvollkommen, insbesondere berücksichtigt sie die Veränderungen der vom Film effektiv empfangenen Lichtintensität nicht. Außerdem berücksichtigt diese Korrektur die Auswirkungen anderer Phänomene wie etwa die Härtung der Röntgenstrahlung wegen der Dicke des durchquerten Gegenstands sowie die Modifikation des Spektrums aufgrund der Spannung der Röntgenröhre nicht.
• Darüber hinaus ist in diesem Verfahren die Detektorzelle vor dem Bildempfänger angeordnet.
• Eine andere Lösung, die die verschiedenen auftretenden Effekte berücksichtigt, insbesondere die Veränderungen des Röhrenstroms, die spektrale Härtung, die durch die Dicke des durchquerten Gegenstandes bedingt ist, und die Modifikation des Spektrums, die durch die Röntgenspannung bedingt ist, sowie bei Vorhandensein eines Verstärkerbildschirms das Absorptions-Ansprechverhalten desselben, ist in der heute eingereichten Patentanmeldung (EP-A-0 465 360) mit dem Titel "PROCEDE DE DETERMINATION AUTOMATIQUE DE LA DUREE D'EXPOSITION D'UN FILM RADIOGRAPHIQUE ET SYSTEME DE MISE EN OEUVRE" beschrieben.
• In diesem Verfahren wird, um die Veränderungen des Stroms der Röntgenröhre und allgemeiner die Veränderungen der Photonenrate auf dem Röntgenfilm zu berücksichtigen, ein Nichtreziprozitäts-Koeffizient CNRD verwendet, der durch die auf den Film auftreffende Photonenrate durch besondere Messungen und Rechnungen ausgedrückt wird.
• Das Ziel der vorliegenden Erfindung ist daher, ein Verfahren zur Bestimmung der die Wirkung der Nichtreziprozität eines Röntgenfilms repräsentierenden Funktion, zu schaffen, wobei die Nichtreziprozität als Funktion der Photonenrate auf dem Röntgenfilm dargestellt wird.
• Die Erfindung betrifft ein Verfahren zur Bestimmung der die Wirkung der Nichtreziprozität eines Röntgenfilms repräsentierenden Funktion in einem radiologischen System gemäß Anspruch 1.
• Die Koeffizienten CNRD(d) können gemäß einer Funktion wie etwa
• modellhaft dargestellt werden, wobei die Parameter A'&sub0;, A'&sub1; und A'&sub2; ausgehend von den Werten di und CNRD(di) geschätzt werden.
• Die Operation (a) des Verfahrens gemäß Anspruch 1 kann die folgenden Operationen enthalten:
• (a1) Modifizieren des Heizstroms der Röhre in der Weise, daß verschiedene Werte des Stroms erhalten werden;
• (a2) Ablesen der Werte M(ti), die von der Integratorschaltung (16) für verschiedene Aufnahmezeiten (ti) ausgegeben werden, um so eine feste optische Dichte DOrefo des Films zu erhalten;
• (a3) Berechnen des Verhältnisses
• M(ti)/M(tref) (29)
• das den Koeffizienten CNRT(ti) angibt, wobei M(tref)
• der Wert M(ti) für ti = tref ist.
• Die Operation (a) des Verfahrens kann außerdem durch die folgenden Operationen ausgeführt werden:
• (g1) Verwirklichen eines ersten Sensitogramms Srefo mit Hilfe eines Sensitometers mit veränderlicher Zeit, wenn die Aufnahmezeit für eine Referenzzeit trefo gesteuert wird;
• (g&sub2;) Verwirklichen von q Sensitogrammen S&sub1; bis Sq mit Hilfe desselben Sensitometers mit veränderlicher Zeit für q verschiedene Aufnahmezeiten ti;
• (g&sub3;) Wählen einer optischen Referenzdichte DOrefo,
• (g&sub4;) Messen der Beleuchtungsstufe Echrefo, Ech&sub1;,.. .,Echi, ...,Echq, die der optischen Dichte DOrefo entsprechen, für jedes Sensitogramm,
• (g&sub5;) Berechnen der Koeffizienten CNRT(ti) durch die Gleichung:
• mit K = 2/log&sub1;&sub0;(2).
• Die Koeffizienten CNRT(ti) können in Form eines analytischen Modells modellhaft dargestellt werden:
• CNRD(t) = A&sub0; + A&sub1; log t + A&sub2;[log t]², (18)
• wobei die Parameter A&sub0;, A&sub1; und A&sub2; anhand der Werte ti und CNRT(ti) geschätzt werden.
• Für die Operation (b) des Verfahrens gemäß Anspruch 1 kann die Referenzbelichtung Lref durch verschiedene Eichungen, beispielsweise mit Hilfe einer Detektorzelle bestimmt werden, welche das vom Schirm emittierte Licht mißt.
• Weitere Ziele, Merkmale und Vorteile der vorliegenden Erfindung werden deutlich beim Lesen der folgenden Beschreibung des Verfahrens gemäß der Erfindung und eines besonderen Ausführungsbeispiels des radiologischen Systems für dessen Ausfünrung, wobei die Beschreibung mit Bezug auf die beigefügten Zeichnungen gegeben wird, in denen:
• - Fig. 1 ein Funktionsschema eines radiologischen Systems ist, das die Ausführung des Verfahrens gemäß der Erfindung ermöglicht;
• - Fig. 2 ein Diagramm ist, das die Kurven zeigt, welche bei der Ausführung des Eichverfahrens erhalten werden, das im Verfahren gemäß der Erfindung verwendet wird,
• - Fig. 3 ein Diagramm ist, das eine Kurve der Veränderung der Nichtreziprozitäts-Koeffizienten CNRT als Funktion der Aufnahmezeit t zeigt,
• - Fig. 4 ein Diagramm ist, das eine Kurve der Veränderung der Nichtreziprozitäts-Koeffizienten CNRD als Funktion der inversen Rate d zeigt, und
• - Fig. 5 ein Diagramm ist, das die Kurven der Veränderung der optischen Dichte eines radiologischen Films als Funktion der Belichtung zeigt.
• Ein radiologisches System, auf das das Verfahren zur automatischen Bestimmung der Belichtungszeit eines zu röntgenden Gegenstandes 13 gemäß der Erfindung angewandt wird, enthält eine Quelle 11 für Röntgenstrahlung wie etwa eine Röntgenröhre, die ein Bündel 14 von Röntgenstrahlen liefert, die diesen Gegenstand 13 belichten, sowie einen Bildempfänger 17 wie etwa ein Film-Bildschirm-Paar, der so angeordnet ist, daß er die Röntgenstrahlen empfängt, die den Gegenstand durchquert haben, und nach einer geeigneten Aufnahmedauer S und der Entwicklung des Films ein Bild des Gegenstandes 13 liefert.
• Für die Ausführung des Verfahrens der Erfindung enthält das System außerdem eine Detektorzelle 12, die im Fall eines Röntgenfilms mit Verstärkerbildschirm hinter dem Bildempfänger 17 angeordnet ist. Diese Zelle kann im Falle eines Films ohne Verstärkerbildschirm vor dem Empfänger angeordnet sein. Die Detektorzelle 12 ermöglicht die Umwandlung einer physikalischen Größe, die für die Röntgenstrahlung charakteristisch ist, die den Gegenstand und den Bildempfänger durchguert hat, etwa das KERMA oder die Energiefluenz, in ein Meßsignal L, das beispielsweise elektrischer Art ist. Das Signal L, das von der Detektorzelle 12 geliefert wird, wird in eine Schaltung 16 eingegeben, die eine Integration des elektrischen Signals während der Aufnahmedauer S ausführt. Das Signal M, das sich aus der Integration ergibt, ist ein Maß für die Strahlung, die den Gegenstand 13 während der Aufnahmedauer S durchquert hat.
• Der Quelle 11 der Röntgenstrahlung ist eine Versorgungsvorrichtung 15 zugeordnet, die eine veränderliche Versorgungshochspannung Vm der Röntgenröhre liefert und ein Gerät zum Messen des Anodenstroms I der Röhre enthält. Um die Dauer der Aufnahmezeit S abzuwandeln, enthalten die Versorgungsvorrichtung 15 und die Röntgenröhre Mittel zum Anregen der Röntgenemission in einem genauen Zeitpunkt und zum Beenden derselben nach einer veränderlichen Dauer S, die in Abhängigkeit von dem von der Schaltung 16 gelieferten Signal M und den Werten von I, S und Vm und genauer vom Verhältnis M/I × S, das Wirkungsgrad D genannt wird und von der Einrichtung 18 berechnet wird, bestimmt wird. Die Werte des Wirkungsgrades D werden von einem Rechner oder Mikroprozessor 19 verarbeitet.
• Das Verfahren gemäß der Erfindung wird nun im Rahmen eines allgemeineren Verfahrens beschrieben, welches dasjenige der automatischen Bestimmung der Belichtungsdauer eines Röntgenfilms ist, wobei das Verfahren den Gegenstand einer Patentanmeldung (EP-A-0 465 360) bildet, die am gleichen Tag wie die vorliegende Anmeldung eingereicht worden ist und den Titel "PROCEDE DE DETERMINATION AUTOMATIQUE DE LA DUREE D'EXPOSITION D'UN FILM RADIOGRAPHIQUE ET SYSTEME DE MISE EN OEUVRE" hat.
• Die erste Operation des Verfahrens zur automatischen Bestimmung der Belichtungsdauer besteht darin, eine Eichung des radio1ogischen Systems von Fig. 1 auszuführen, die zu einer Funktion der Schätzung der vom Röntgenfilm gesehenen Belichtung führt. Diese Eichung und die Schätzfunktion sind in einer Patentanmeldung (EP-A-0 465 362) beschrieben, die am gleichen Tag wie die vorliegende Anmeldung eingereicht worden ist und den Titel "PROCEDE D'ESTIMATION ET D'ETALONNAGE DE LA LUMINATION RECUE PAR UN FILM RADIOGRAPHIQUE" hat.
• Für das Verständnis der folgenden Beschreibung wird daran erinnert, daß das Verfahren zum Schätzen der von einem Röntgenfilm empfangenen Belichtung auf Eichungen, die zur Definition einer zur Photonenrate auf den Film proportionalen Funktion, die Filmrate genannt wird, sowie auf einer Eichung basiert, die die Verbindung zwischen der Filmraten-Funktion und der vom Film unter festen Referenzbedingungen empfangenen Belichtung, die zu einer gegebenen Schwärzung des Films führt, herzustellen ermöglicht. Diese letztere Eichung wird in der folgenden Beschreibung genauer dargelegt.
• Die Eichungen, die die Definition der Filmrate ermöglichen, sind aus einem Eichungsverfahren abgeleitet, das in der Patentanmeldung Nr. 89 07686 (die der EP-A-0 402 244 entspricht), die am 9. Juni 1939 eingereicht worden ist und den Titel "PROCEDE D'ETALONNAGE D'UN SYSTEME RADIOLOGIQUE ET DE MESURE DE L'EPAISSEUR EQUIVALENTE D'UN OBJET" hat. Dieses Verfahren besteht darin, den Wirkungsgrad D der Zelle für jedes Normal der gewählten Versorgungsspannungen Vm zu messen. Genauer wird mit einem ersten Normal der Dicke E&sub1; eine Messung des Wirkungsgrades D1m für jeden Wert Vm, die eine bestimmte Gruppe bilden, ausgeführt. Diese Werte D1m als Funktion der Spannung Vm können in einem Diagramm aufgezeichnet werden, um die Punkte 21' von Fig. 2 zu erhalten.
• Die Messungen des Wirkungsgrades D werden für ein anderes Normal der Dicke E&sub2; ausgeführt, wobei die Werte D2m erhalten werden, die den Punkten 22' von Fig. 2 entsprechen, usw. um die anderen Serien von Punkten 23', 24' und 25' zu erhalten, die den Wirkungsgraden D3m, D4m bzw. D5m sowie den Dicken E&sub3;, E&sub4; bzw. E&sub5; entsprechen.
• Es ist zu bemerken, daß in Fig. 2 die Wirkungsgrade Dpm als logarithmische Ordinaten aufgetragen sind, während die Versorgungsspannungen als Abszissen von 20 Kilovolt bis 44 Kilovolt aufgetragen sind.
• Diese Serien von Punkten 21' bis 25' dienen dazu, die Parameter eines analyitschen Modells zu definieren, das das Verhalten des Wirkungsgrades D als Funktion der Parameter Vm und Et für eine gegebene Konfiguration des radiologischen Systems beschreibt. Dieses analytische Modell wird folgendermaßen geschrieben:
• D = f(Vm,Ep) (1)
• Die Parameter des analytischen Modells können mit Hilfe von herkömmlichen Schätzwerkzeugen wie etwa der Methode der Minimierung des quadratischen Fehlers eingestellt werden.
• Die Kurven 21 bis 25 zeigen den Wert des Wirkungsgrades D, der durch das analytische Modell gegeben ist, das durch den folgenden Ausdruck repräsentiert ist:
• D = f(Vm, Ep) = exp[f&sub1;(Vm) + Ep × f&sub2;(Vm)], (2)
• in dem f&sub1;(Vm) und f&sub2;(Vm) Polynome zweiten Grades sind, deren Ausdruck gegeben ist durch
• f&sub1;(Vm) =A&sub0; + A&sub1;Vm + A&sub2;Vm²
• f&sub2;(Vm) =B&sub0; + B&sub1;Vm + B&sub2;Vm²
• Die inverse Funktion der durch die Formel (2) ausgedrückten Funktion ermöglicht die Berechnung von Ep, wenn D und Vm bekannt sind, indem die folgende Formel (3) verwendet wird:
• Die Funktion f&sub2;(Vm) kann nämlich für übliche Werte von Vm nicht Null werden, weil bei den betrachteten Spannungen Vm der Wirkungsgrad D stets von der Dicke Ep abhängt.
• Mit anderen Worten, ein Wertepaar (Ep,Vm) entspricht einer Messung des Wirkungsgrades D' welche die Bestimmung von Ep in Abhängigkeit von Vm und von D ermöglicht. Während einer radiologischen Untersuchung ermöglicht eine Messung des Wirkungsgrades D, die mit einer gegebenen Versorgungsspannung Vm ausgeführt wird, eine äquivalente Dicke zu bestimmen, die in den für Ep verwendeten Einheiten ausgedrückt wird.
• Diese Eichung wird zweimal mit bezüglich des Empfängers 17 unterschiedlichen Konfigurationen des Röntgensystems ausgeführt. Die erste dieser Eichungen wird mit dem Empfänger 17 ohne Verstärkerbildschirm ausgeführt. Gemäß Gleichung (1) wird eine Funktion f' bestimmt, die Werte des Wirkungsgrades der Zelle 12 ergibt, die mit Dse bezeichnet werden, etwa:
• Dse= f'(Vm,Ep) (4)
• sowie die inverse Funktion
• Ep = g'(Vm,Dse) (5)
• Die zweite Operation des Verfahrens besteht darin, eine zweite Eichung mit Verstärkerbildschirm und Empfänger 17 auszuführen, wobei dann eine Serie von Werten des Wirkungsgrades Dc erhalten wird und wie oben die Funktion f" wie etwa:
• Dc = f"(Vm,Ep) (6)
• und die inverse Funktion
• Ep = g"(Vm,Dc) (7)
• bestimmt werden.
• Aus den beiden vorangehenden Eichungen wird eine Funktion Df abgeleitet, die den Wirkungsgrad auf dem Film repräsentiert, etwa:
• Df = Dse -Dc,
• also Df = f'(Vm,Ep) - f"(Vm,Ep). (8)
• Diese Funktion Df berücksichtigt die durch die zusätzliche Filterung zwischen dem Verstärkerbildschirm und der Detektorzelle 12 bedingte Modifikation des Spektrums der Röntgenstrahlung, die beispielsweise von der Ausgangsseite der das Film-Bildschirm-Paar enthaltenden Kassette ausgeht, nicht. Um sie zu berücksichtigen, wird in Gleichung (8) Ep durch (Ep - (zusätzliches Filter)) ersetzt, wobei (zusätzliches Filter) die für den geröntgten Gegenstand äquivalente Dicke ist, die dieser Filterung entspricht.
• Diese äquivalente Dicke wird erhalten, indem beispielsweise in das Bündel 14 ein dieser Filterung äquivalenter Gegenstand eingesetzt wird und indem je nach Konfiguration der Maschine die Eichfunktion g' oder g" verwendet wird, die die äquivalente Dicke bestimmt.
• Da das Produkt Df × I × t zu der im Verstärkerbildschirm während einer Zeit t und bei einem Anodenstrom I absorbierten Energie proportional ist, ist die Größe Df × I, die mit Filmrate bezeichnet wird, zur Rate der auf den Film auftreffenden Photonen proportional und wird in Meßeinheiten des Signals der Detektorzelle 12 ausgedrückt.
• Diese Proportionalitätsbeziehung wird umso besser erfüllt, je besser die Anzahl der vom Verstärkerbildschirm emittierten Photonen ihrerseits zur absorbierten Energie proportional ist. Wenn die Anzahl der vom Verstärker emittierten Photonen einem anderen Abhängigkeitsgesetz von der absorbierten Energie gehorcht, muß dieses andere Gesetz auf Df × I angewandt werden, um die "Filmrate" zu erhalten.
• Eine letzte Eichung besteht darin, die obenbeschriebenen elektrischen Funktionen mit einem Wert der Schwärzung des Films, d.h. einer optischen Dichte, die am Ende der Belichtung erhalten werden soll, in Beziehung zu setzen. Die Wahl dieses Wertes wird durch den Anwender in Abhängigkeit vom Film-Bildschirm-Paar, vom Typ der Diagnose, vom zu untersuchenden Körperteil des Patienten und von seinen Gewohnheiten bei Röntgenuntersuchungen ausgeführt. Die Wahl ermöglicht die Bestimmung der Referenzbelichtung, die mit Lref bezeichnet wird, d.h. der Belichtung, die der Film unter festen Referenzbedingungen empfangen muß, um eine solche Schwärzung zu erreichen. Das Verfahren der Bestimmung von Lref wird später beschrieben.
• Die Eichoperationen werden nicht bei jeder Röntgenuntersuchung eines Gegenstandes oder eines Patienten ausgeführt, sondern nur von Zeit zu Zeit, um die Veränderungen der Eigenschaften des radiologischen Systems im Laufe der Zeit zu berücksichtigen, insbesondere die Alterung der Röntgenröhre. Die Ergebnisse dieser Operationen werden im Speicher des Mikroprozessors 19 in Form von Funktionen gespeichert, die durch die Gleichungen 4 bis 8 repräsentiert sind, was bedeutet, daß der Mikroprozessor 19 Ep berechnen kann, wenn er Dc kennt, woraufhin er Df berechnen kann.
• Während der radiologischen Untersuchung des Patienten besteht das Verfahren außerdem darin, die folgenden Hauptoperationen auszuführen:
• (e1) Anordnen des zu röntgenden Gegenstandes oder Patienten,
• (e2) Auslösen des Beginns der Aufnahme durch den Anwender,
• (e3) Messen des Wirkungsgrades Dc nach einer bestimmten Zeit t' seit dem Beginn der Aufnahme,
• (e4) Berechnen der äquivalenten Dicke anhand der Messung des Wirkungsgrades Dc,
• (e5) Berechnen des Wirkungsgrades Df auf Höhe des Films,
• (e6) Schätzen der vom Film empfangenen Belichtung seit dem Beginn der Aufnahme,
• (e7) Berechnen der Belichtung, die noch erfaßt werden muß, um die gewählte Schwärzung zu erhalten,
• (e8) vorausgehende Berechnung der mA s, die noch in die Röntgenröhre geliefert werden müssen, um die gewählte Schwärzung zu erhalten,
• (e9) Messen der mA s, bezeichnet mit mAsmes, die je nach Fall seit dem Beginn der Aufnahme oder der vorangehenden Messung geliefert wurden,
• (e10) Beenden der Röntgenstrahlung, wenn die mAsmes größer oder gleich den berechneten mAs sind, oder Rückkehr zur Operation (e3) im entgegengesetzten Fall.
• Es ist zu beachten, daß mit Belichtung das Produkt der empfangenen Lichtmenge, beispielsweise der Beleuchtung EC der empfindlichen Fläche, mit der Dauer der Bestrahlung oder Aufnahme bezeichnet wird.
• Die Operation (e3) besteht darin, den integrierten Wert D zu messen, der von der Vorrichtung 18 nach einer bestimmten Zeit t' seit dem Beginn der Aufnahme geliefert wird, da der Integrierer 16 je nach Fall entweder am Beginn der Aufnahme oder nach der letzten Messung auf Null zurückgesetzt worden ist. Die Zeit t' der Integration entspricht je nach Fall der seit dem Beginn der Aufnahme verstrichenen Zeit oder der seit der letzten Messung verstrichenen Zeit.
• Die Operation (e4) wird vom Mikroprozessor 19 anhand der ersten Eichung des radiologischen Systems wie oben beschrieben ausgeführt: Sie ist durch die Gleichung (7) bestimmt; dann wird ein Wert E&sub1; der äquivalenten Dicke erhalten.
• Es ist zu bemerken, daß es für die zweite Iteration des Verfahrens und die folgenden nicht notwendig ist, die Operation (e4) auszuführen, sofern die Schätzung der äquivalenten Dicke bei der ersten Iteration ausreichend genau gewesen ist. Die Operation (e5) besteht darin, den Wirkungsgrad des Films Df1, der der Dicke E&sub1; entspricht, zu berechnen, indem die durch die Gleichung (8) definierte Funktion verwendet wird, was ermöglicht, insbesondere den Einfluß des Empfängerbildschirms zu berücksichtigen. Diese Operation ist oben kurz beschrieben worden.
• Die Operation (e6) besteht darin, die vom Film seit dem Beginn der Aufnahme empfangene Belichtung Lf zu schätzen, indem die folgende Gleichung angewandt wird:
• Lf = Lam + Df1 × δmA s, (9)
• in welcher Lam die vom Film vor der Operation (e3) empfangene Belichtung ist und δmA s die Anzahl der mA s ist, die in die Röhre während der Zeit t' geliefert werden und die durch das Produkt des Stroms I der Röhre mit der Integrationszeit S definiert ist.
• Die Operation (e7) besteht darin, die Belichtung Lra zu berechnen, die aufzunehmen bleibt, um die gewählte Schwärzung zu erhalten; sie ist durch die folgende Gleichung bestimmt:
• Lra = Lref - Lf. (10)
• Die Operation (es) besteht darin, die mA s zu berechnen, die noch geliefert werden müssen, um die gewählte Schwärzung zu erhalten, und die durch die folgende Gleichung gegeben sind:
• mAsr = Lra/Df1. (11)
• Es ist dann möglich, die Anzahl der während der Berechnungen gelieferten mA s abzuleiten, nämlich mASc. Die tatsächlich noch auf zunehmenden mA s, die mit mAsra bezeichnet sind, sind folgendermaßen definiert:
• mAsra = mAsr - mAsc (12)
• oder mASc = I × tc, (13)
• wobei tc die Rechenzeit ist.
• Die Operation (e10) besteht in der Ausführung einer Wahl: entweder die Aufnahme beenden oder sie je nach dem Wert der noch zu liefernden mA s oder der noch abzuwartenden Aufnahmezeit fortzusetzen oder die Schätzung des vorläufigen Wertes des Aufnahmeendzeitpunktes erneut berechnen.
• Das Kriterium für das Aufnahmeende kann folgendermaßen lauten: Wenn der Wert
• Dif(mA s) = mAsra - mAsmes (15)
• null oder kleiner als ein fester Wert Val&sub0; ist, hält der Mikroprozessor 19 die Röntgenstrahlung durch Einwirken auf die Versorgungsvorrichtung 15 an. Im entgegengesetzten Fall wird zur Operation (e3) zurückgekehrt.
• Es ist möglich, eine zusätzliche Prüfung hinsichtlich des Wertes der verbleibenden Aufnahmezeit trc, welche durch die folgende Beziehung definiert ist, in Betracht zu ziehen:
• trc = mAsra/I (14)
• Diese zusätzliche Prüfung besteht darin, den Schätzwert mAsra in dem Fall, in dem trc kleiner als ein Wert t&sub0; ist, nicht zu modifizieren. Dann endet das Aufnahmeende in einer offenen Schleife, indem nur die Aufnahmeende-Operationen fortgeführt werden, d.h. Erniedrigen der Anzahl der gelieferten mA s und Beenden der Aufnahme, wenn diese Anzahl kleiner oder gleich Null wird. Ein möglicher Wert von t&sub0; ist ein Wert, der im wesentlichen gleich dem Zeitintervall ist, der zwei der Operation (e3) entprechende Messungen trennt. Somit enthält die Operation (e10) in diesem Fall zwei Prüfungen:
• - Eine erste Prüfung von mAsra, die entscheidet, ob die Aufnahme endet oder nicht,
• - dann eine Prüfung von trc, um zu entscheiden, ob eine neue Schätzung der noch zu liefernden mAs unternommen wird oder ob der Wert mAsra bis zum Aufnahmeende fest bleibt. In diesem letzteren Fall erfolgt die Aufnahmeende-Prüfung periodisch mit dem Wert mAsra.
• Außerdem können die Operationen der Schätzung der verbleibenden Zeit und diejenige des Unterbrechens der Aufnahme entkoppelt sein, um die Genauigkeit des Belichters zu verfeinern. Somit wird das Verfahren in der folgenden Weise zerlegt: in einem Arbeitsvorgang TE, der dazu bestimmt ist, die vor dem Aufnahmeende noch zu liefernden mA s zu schätzen, und in einem Arbeitsvorgang TC der Unterbrechung der Aufnahme. Dies sind zwei unabhängige Arbeitsvorgänge, die parallel ablaufen.
• Der Schätz-Arbeitsvorgang TE der noch zu liefernden mA s ist durch die Operationen (e3) bis (e8) gebildet, denen eine Operation (e'9) der Umwandlung der mA s in ein Signal in Einheiten der Zelle 12 hinzugefügt ist, etwa:
• CEziel = mAsra × Dc. (16)
• Dieser Schätz-Arbeitsvorgang TE wird während der Aufnahme periodisch erneuert, beispielsweise in Zeitpunkten t&sub1;,t&sub2;, ... ,tn, welche Meßzeitpunkte sind, die durch eine Dauer getrennt sind, die wenigstens gleich der Rechenzeit tc ist. Am Ende des Schätz-Arbeitsvorgangs TE wird der Zielwert CEziel des Unterbrechungs-Arbeitsvorgangs TC aktualisiert. Diese Aktualisierung muß das von der Detektorzelle 12 zwischen dem Meßzeitpunkt am Beginn der Operation (e3) und dem Zeitpunkt der Aktualisierung des Wertes CEziel am Ende der Operation TC empf angene Signal berücksichtigen.
• Der Arbeitsvorgang TC der Unterbrechung der Aufnahme ist ein Arbeitsvorgang, der darin besteht, einen gegebenen Wert oder ein Ziel in Abhängigkeit von dem von dieser Zelle tatsächlich empfangenen Signal zu erniedrigen. Dieser Arbeitsvorgang unterbricht die Aufnahme, sobald der Wert CEziel kleiner oder gleich Val&sub0; wird, der beispielsweise gleich null ist.
• Somit kann die Funktionsweise des Arbeitsvorgangs TC durch die folgenden Operationen zusammengefaßt werden:
• (f1) Messen des von der Zelle 12 integrierten Signals Mm nach einer bestimmten Zeit tTC;
• (f2) Erniedrigen des Zielwerts um diesen Wert: (CEziel-Mm);
• (f3) Beenden der Aufnahme, wenn (CEziel-Mm) kleiner als Val&sub0; ist, andernfalls Rückkehr zu (f1).
• Das eben beschriebene Verfahren arbeitet in dem Maß richtig, in dem für den Empfänger 17 und die Detektorzelle 12 keine Abweichung vom Reziprozitätsgesetz vorliegt. Wenn dies nicht der Fall ist, müssen die operationen (e6) und (e8) dahingehend vervollständigt werden, daß hierfür ein Korrekturkoeffizient berücksichtigt und durch besondere Messungen und Berechnungen bestimmt wird. Dieser Korrekturkoeffizient wird in die Gleichungen (9) und (11) eingeführt, in welche die Belichtung und der Wirkungsgrad des Films eingehen.
• Daher werden die Formeln (9) und (11) zu:
• mit Filmrate = Df1 × I. (17)
• CNRD ist die Funktion, die die Wirkung der Nichtreziprozität repräsentiert, die in Abhängigkeit von der Photonenrate auf den Film ausgedrückt wird.
• Die Funktion CNRD wird durch das Verfahren erhalten, die den Gegenstand der vorliegenden Patentanmeldung bildet. Dieses Verfahren besteht zu allererst darin, die Nichtreziprozitäts- Koeffizienten des Films in Abhängigkeit von der Aufnahmedauer ti zu bestimmen, die mit CNRT(ti) bezeichnet werden. Diese Funktion CNRT wird experimentell bestimmt und kann durch eine analytische Funktion dargestellt werden.
• Genauer besteht das Verfahren darin, für verschiedene Werte IRi der Intensität der Strahlung den Wert ti der Aufnahmedauer zu bestimmen, welche erforderlich ist, um eine feste optische Dichte DOrefo des Films zu erhalten, beispielsweise DOrefo = 1, und die von der Integratorschaltung 16 gelieferten Werte für die verschiedenen Aufnahmezeiten ti aufzunehmen, wobei diese Werte mit M(ti) bezeichnet werden.
• Diese Werte werden mit einem Referenzwert M(tref) verglichen, der beispielsweise derjenige ist, der einer Aufnahmezeit von einer Sekunde entspricht, indem das folgende Verhältnis berechnet wird:
• M(ti)/M(tref). (29)
• Dieses Verhältnis bestimmt den Koeffizienten der zeitlichen Nichtreziprozität CNRT(ti) für die Aufnahmezeit ti. Eine andere Weise der Bestimmung dieser Koeffizienten CNRT(ti) wird später beschrieben.
• Die Koeffizienten CNRT(ti) stehen im Fall beispielsweise einer optischen Dichte DOrefo = 1 und einer Referenz-Aufnahmezeit ref = 1 Sekunde als Funktion der Zeit durch die Kurve von Fig. 3 miteinander in Beziehung. Diese Kurve zeigt, daß die Belichtung, die erforderlich ist, um die gewünschte optische Dichte zu erhalten, mit der Aufnahmezeit wächst. Deshalb liegt in diesem Beispiel das Verhältnis zwischen den Energien für die beiden Aufnahmezeiten von 50 ms und 6,5 s in der Größenordnung von 1,6. Die Kurve von Fig. 3 kann mit Hilfe einer Funktion modellhaft dargestellt werden, welche die Form besitzt:
• CNRD(ti) = A&sub0; + A&sub1; log t + A&sub2;[log t]² , (18)
• deren Parameter A&sub0;, A&sub1; und A&sub2; ausgehend von Meßpunkten durch eine Schätzprozedur der kleinsten Quadrate geschätzt werden.
• Im Prinzip könnte der Schwarzschild-Effekt, der in den Gleichungen (9') und (11') berücksichtigt wird, durch die Funktion CNRT modellhaft dargestellt werden. Der Nutzen der Verwendung der mit der Rate indexierten Funktion CNRD besteht darin, daß die Veränderungen des Anodenstroms berücksichtigt werden können. Daher hat ein automatischer Belichter, der die Funktion CNRD gemäß den Gleichungen (9') und (11') benutzt, beispielsweise den Vorteil, daß die Röhre mit abnehmender Last arbeiten kann.
• Um von den von der Zeit abhängigen Koeffizienten CNRT(t) zu den von der Rate abhängigen Koeffizienten CNRD(d) überzugehen, muß die Tatsache berücksichtigt werden, daß die Koeffizienten CNRT(t) durch Messungen der veränderlichen Aufnahmezeit bestimmt worden sind, ohne die Werte der Photonenrate auf dem Film zu kennnen. Wenn für jede Aufnahmezeit ti die Filmrate di gemessen wird, ist der Wert des Koeffizienten CNRD(di) für di gleich demjenigen des Koeffizienten CNRT(ti) für die Aufnahmezeit ti entsprechend der Beziehung:
• CNRD(di) = CNRT(ti) (19)
• Die verschiedenen Werte von CNRD(di) stehen als Funktion der inversen Rate 1/d durch die Kurve (Fig. 4) miteinander in Beziehung. Diese Kurve kann mit Hilfe einer Funktion der folgenden Form modellhaft dargestellt werden:
• Die Werte di können nicht durch Eichung gegeben sein, vor allem weil sie in Meßeinheiten der Zelle 12 ausgedrückt sind, die nicht notwendigerweise diejenige ist, die bei der Eichung verwendet wird. Somit müssen im allgemeinen die Werte di mit den bekannten Werten ti durch die folgende Beziehung verbunden werden:
• Lref × CNRT(ti) = di × ti (21)
• oder aber
• Es wird hier daran erinnert, daß Lref die Belichtung ist, die der Film unter festen und bekannten radiologischen Bedingungen empf ängt, wenn der Film eine gegebene Schwärzung erreicht hat, und bei der die Wirkung der Nichtreziprozität korrigiert ist.
• Um die Definition der Funktion CNRD abzuschließen, bleibt das Verfahren der Auswertung der Referenzbelichtung zu betrachten.
• Dieses Verfahren ist in der obenerwähnten Patentanmeldung (EP-A-0 465 362) mit dem Titel "PROCEDE D'ESTIMATION ET D'ETALONNAGE DE LA LUMINATION RECUE PAR UN FILM RADIOGRAPHIQUE" beschrieben. Die Referenzbelichtung hängt von der optischen Dichte ab, die auf dem Film erhalten werden soll. Um diese Belichtung zu bestimmen, besteht die erste Stufe darin, ein Sensitogramm des verwendeten Filmtyps zu verwirklichen, anschließend muß unter mit einem bekannten Dickennormal bestimmten radiologischen Bedingungen eine Aufnahme hergestellt werden.
• Die bestimmten radiologischen Bedingungen sind beispielsweise
• - eine optische Referenzdichte DOrefo, die in Abhängigkeit von den Gewohnheiten des Anwenders gewählt wird, beispielsweise DOrefo = 1
• - ein Dickennormal Eo,
• - eine Versorgungsspannung Vo,
• - ein Wert der Aufnahmezeit to,
• - ein Wert des Produkts Io × to.
• Für diese Aufnahme werden die optische Dichte DOm sowie die Werte Mo, Io, to gemessen, wodurch es möglich ist, die äquivalente Dicke Ep mit Hilfe der Formel (7) zu berechnen. Anschließend wird der Wirkungsgrad Dfo auf dem Film mit Hilfe der Formel (6) berechnet, wodurch es möglich ist, die vom Film empfangene Belichtung Lfilm mit der folgenden Formel zu berechnen:
• Lfilm = Dfo × Io × to (23)
• Die optische Referenzdichte DOrefo ermöglicht die Berechnung der DOrefo entsprechenden Beleuchtungsstufe Echref auf der sensitometrischen Kurve des verwendeten Films (Fig. 5), wobei diese Kurve mit Hilfe eines Sensitometers und eines Densitometers aufgezeichnet worden ist. Dadurch können die Eigenschaften der zu entwickelnden Maschine, die verwendet wird, berücksichtigt werden. Die Kurve ist beispielsweise in Form einer Funktion im Mikroprozessor 19 (Fig. 1) gespeichert.
• Die gemessene optische Dichte DOm ermöglicht die Berechnung der Meßstufe Echm, die gleich dem Wert der Beleuchtungsstufe ist, die auf der sensitometrischen Kurve (Fig. 5) DOm entspricht.
• Mit den Werten Lfilm der Belichtung auf dem Film, der Referenzstufe Echref und der Meßstufe Echm ist es möglich, die Referenzbelichtung Lref zu berechnen, um die optische Dichte DOrefo zu erhalten, indem die Gleichung verwendet wird, die die Maßstabsänderung zwischen der Belichtung und der Beleuchtungsstufe der Abszissenachse der sensitometrischen Kurve (Fig. 5) definiert:
• Dieser Gleichung (24) wird entnommen:
• mit K = 2/log&sub1;&sub0;(2). (26)
• Die sensitometrische Konstante K entspricht dem beibehaltenen Abstand für die Beleuchtungsstufen.
• Der Wert Lref hängt über Lfilm durch die Gleichungen (22) und (25) von t&sub0; ab. Somit ist der Wert Lref für die Wirkungen der Nichtreziprozität des Films empfindlich. Um den Einfluß der Nichtreziprozität auf den Wert Lref zu korrigieren, ist es ausreichend, in der Gleichung (23) den Wert Lfilm zu verwenden, der definiert ist durch:
• Diese Referenzbelichtung Lref ist diejenige, die in Gleichung (10) verwendet werden muß, um die optische Referenzdichte DOrefo zu erhalten, beispielsweise DOrefo = 1, wobei die Formel (25) zeigt, daß sie insbesondere von der Differenz zwischen der Referenzstufe und der Meßstufe abhängt.
• Die Kenntnis der vom Film empfangenen Belichtung ermöglicht die Kenntnis von di durch Anwenden der Formel (22) und daraus die Ableitung von CNRD(dj) durch die Formel (20).
• Für eine optische Dichte des Röntgenfilms, die vom DOrefo = 1 verschieden ist, ist es notwendig, erneut die obenbeschriebenen Operationen auszuführen, derart, daß die neuen Werte von CNRT(ti) und von Lref bestimmt werden.
• Um diese Operationen zu vereinfachen, können die Koeffizienten CNRT(ti) dadurch erhalten werden, daß die folgenden Operationen ausgeführt werden:
• (g1) Verwirklichen eines ersten Sensitogramms Srefo mit Hilfe eines Sensitometers mit veränderlicher Zeit, wenn die Aufnahmezeit für eine Referenzzeit trefo gesteuert wird;
• (g2) Verwirklichen von q Sensitogrammen S&sub1; bis Sq (Fig. 5) mit Hilfe desselben Sensitometers mit veränderlicher Zeit für q verschiedene Aufnahmezeiten;
• (g3) Wählen einer optischen Referenzdichte DOrefo, z.B. DOrefo = 1;
• (g4) Messen der Beleuchtungsstufe Echrefo, Ech1, ..., Echi, ..., Echq, die der optische Dichte DOrefo (D0refo = 1) entsprechen, für jedes Sensitogramm;
• (g5) Berechnen der Koeffizienten CNRT(ti) durch die Gleichung:
• Wenn der Anwender entscheidet, mit einer unterschiedlichen optischen Dichte zu arbeiten, wird, um die obenbeschriebene Eichung zu vermeiden, vorgeschlagen, die absichtlich korrigierte optische Dichte für die Schwärzung DOcvn zu verwenden. Dann muß die Referenzbelichtung Lref, die in Gleichung (10) verwendet wird, durch die korrigierte Belichtung Lcvn ersetzt werden, die ausgedrückt wird durch:
• wobei
• - CVN die absichtliche Korrektur der Schwärzung ist, die durch eine ganze Zahl beispielsweise von -10 bis +10 ausgedrückt wird,
• - P die Schrittweite der optischen Dichte ist, beispielsweise 0,1,
• - Γ die Steigung des linearen Teils der sensitometrischen Kurve (Fig. 5) ist.
• Das eben beschriebene Verfahren zeigt, daß seine Ausführung eine bestimmte Anzahl von Eichungen erfordert, die zusammengefaßt die folgenden sind:
• (a) die Eichung des radiologischen Systems in der Weise, daß die analytischen Modelle bestimmt werden,
• Dse = f'(Vm,Ep) (4)
• mit Kassette ohne Bildschirm und
• Dc = f"(Vm,Ep) (6)
• Ep = g"(Vm,Dc) (7)
• mit Kassette und Bildschirm.
• Die Differenz Df = (Dse - Dc) (Gleichung (8)) erlaubt die Ableitung des vom Bildschirm absorbierten Wirkungsgrades;
• (b) Die Eichung des Films in der Weise, daß das Gesetz der Nichtreziprozität CNRT(t) bestimmt wird, das als Funktion der Zeit ausgedrückt ist; dieses Gesetz wird dazu verwendet, das Gesetz der Nichtreziprozität CNRD(d) zu bestimmen, das als Funktion der Rate ausgedrückt ist;
• (c) die Eichung der Referenzbelichtung Lref Wenn diese verschiedenen Eichungen ausgeführt worden sind, sind die verschiedenen eigentlichen Operationen des Verfahrens die folgenden:
• (d) Wählen des Wertes der Schwärzung oder des Wertes der absichtlichen Korrektur der Schwärzung durch den Anwender in der Weise, daß die Ziel-Belichtung Lcvn bestimmt wird, die der Film unter den festen Referenzbedingungen empfangen muß, um die gewählte Schwärzung (oder optische Dichte) zu erreichen. Die Belichtung Lcvn wird ausgehend von der Gleichung (27) berechnet oder die Belichtung Lref wird durch Eichung (c) und die Gleichungen (25) und (26) bestimmt.
• (e1) Anordnen des Gegenstandes, der geröntgt werden soll;
• (e2) Auslösen des Beginns der Aufnahme durch den Anwender;
• (e3) Messen des Wirkungsgrades Dc1 auf Höhe der Zelle (12) nach einer Zeit t';
• (e4) Messen der äquivalenten Dicke E&sub1; durch Gleichung (7);
• (e5) Berechnen des Wirkungsgrades Df1 auf Höhe des Films für die Dicke E&sub1; durch die Gleichung (8);
• (e6) Berechnen der vom Film empfangenen Belichtung Lf durch die Gleichung: CNRD(Filmrate
• (e7) Berechnen der Belichtung Lra, die noch aufgenommen werden muß, um die gewählte Schwärzung (oder optische Dichte) zu erhalten, durch die Gleichung:
• Lra = Lcvn - Lf (10')
• (e8) vorausgehende Berechnung der noch zu liefernden mA s (mAsra), um die Schwärzung (oder optische Dichte) zu erhalten, durch die Gleichung: CNRD(Filmrate
• (e9) Messen der mA s, die seit Beginn der Operation (e3) geliefert worden sind;
• (e10) - Beenden der Belichtung, wenn die in (e9) gemessenen mA s gleich oder größer als mAsra sind,
• - oder Rückkehr zur Operation (e3), wenn die bei (e9) gemessenen mA s kleiner als mAsra sind.
• Die eben gegebene Beschreibung des Verfahrens entspricht einer bestimmten Konfiguration des radiologischen Systems. In dem Fall, in dem dieses System mehrere Konfigurationen annehmen kann, beispielsweise durch die Wahl
• - des Anodenmaterials,
• - der Abmessungen des Brennpunkts,
• - des Filters zur Spektrumsmodifikation,
• - der Kollimation,
• - des Vorhandenseins oder der Abwesenheit eines Antistreugitters,
• - des Typs des Bildempfängers,
• - des Typs der Detektorzelle,
• ist es notwendig, für jede dieser Konfigurationen die Eichungen (a), (b) und (c) auszuführen. Die Anzahl dieser Eichungen kann verringert werden, wenn Ähnlichkeiten des Verhaltens von einer Konfiguration zur nächsten berücksichtigt werden, wie dies für die Eichung (a) in der Patentanmeldung Nr. 89 07686, eingereicht am 9. Juni 1989, beschrieben worden ist.
• Bei der Ausführung des Verfahrens durch den Anwender definiert dieser letztere die Konfiguration, deren Eigenschaften an den Mikroprozessor (19) in der Weise übertragen werden, daß dieser letztere die entsprechenden Modelle verwendet.
• Das Verfahren gemäß der Erfindung ist in seiner Anwendung auf einen Empfänger 17 des Typs Film-Bildschirm-Paars beschrieben worden. Es kann auch in dem Fall eines Empfängers 17 ausgeführt werden, der nur einen für die Röntgenstrahlung empfindlichen Film enthält. Mit einem solchen Film
• werden die Operationen (a) und (b) die folgenden:
• (a') Eichung des radiologischen Systems in der Weise, daß das analytische Modell
• D'f = f"'(Vm,Ep) (30)
• bestimmt wird, wobei der Film als Bildempfänger dient.
• Die Modifikationen sind die folgenden: Im Verlauf des Verfahrens
• - wird (e3) zu (e'3): Messen des Wirkungsgrades D'f&sub1; auf Höhe der Zelle 12 nach einer Zeit t';
• - werden (e4) und (e5) unterdrückt und
• - werden die Operationen (e6) und (e8) auf die folgende Weise modifiziert:
• (e'6) Berechnen der vom Film empfangenen Belichtung L'f durch die Gleichung: CNRD(Filmrate)
• (e7) Berechnen der noch aufzunehmenden Belichtung L'ra, um die gewählte Schwärzung (oder optische Dichte) zu erhalten, durch die Gleichung:
• L'ra = Lcvn - L'f; (10")
• vorhergehendes Berechnen der noch zu liefernden mA s (mAs,ra), um die Schwärzung (oder optische Dichte) zu erhalten, durch die Gleichung: CNRD(Filmrate)
• Die Operationen (e9) und die folgenden bleiben unverändert.
• Außerdem ist festzustellen, daß das Sensitometer in diesem Fall vom Typ mit Röntgenemission sein kann.
• Außerdem kann bei einem solchen Empfänger, der keinen Verstärkerbildschirm enthält, die Detektorzelle 12 entweder hinter dem Empfänger 17, wie im Falle des Empfängers des Typs Film- Bildschirm oder vor dem Empfänger 17, falls die Strahlungsenergie dies zuläßt, angeordnet sein.

Claims (11)

1. Verfahren zur Bestimmung der die Wirkung der Nichtreziprozität eines Röntgenfilms repräsentierenden Funktion in einem für die Untersuchung eines Gegenstandes (13) vorgesehenen radiologischen System, das versehen ist mit einer Röntgenröhre (11), deren Versorgungsspannung verschiedene kontinuierlich oder diskret sich verändernde Werte Vm annehmen kann und die ein Bündel (14) von Röntgenstrahlen in Form von Impulsen mit veränderlicher Dauer ti in Richtung des zu untersuchenden Gegenstandes aussendet, einein Bildempfänger (17) für die Röntgenstrahlung, die den Gegenstand (13) durchquert hat, um ein Bild des Gegenstandes zu verwirklichen, einer Detektorzelle (12) für die Röntgenstrahlen, die den zu untersuchenden Gegenstand (13) durchquert haben, mit der eine das Bündel von Röntgenstrahlen kennzeichnende physikalische Größe in ein Meßsignal L umgewandelt werden kann, einer Integratorschaltung (16), die das Meßsignal L während der Dauer ti integriert und ein Signal M ausgibt, und einer Einrichtung (18) zur Berechnung des Wirkungsgrades D, der durch das Verhältnis von M zum Produkt I × ti (oder nA s) des Anodenstroms I der Röhre mit der Dauer ti der Aufnahme gegeben ist, wobei das Verfahren die folgenden Operationen enthält:

(a) Bestimmen der Nichtreziprozitäts-Koeffizienten CNRT(ti) des Films, die für eine optische Dichte DOrefo des Films für verschiedene Werte (IRi) der Strahlungsintensität durch die Aufnahmezeit ti ausgedrückt werden;

(b) Bestimmen der Referenzbelichtung Lref, die von dem Film unter festen und bekannten radiologischen Bedingungen empfangen wird, wenn der Film eine gegebene Schwärzung erreicht und die Wirkung der Nichtreprozität korrigiert ist;

(c) Berechnen der Photonenraten di für den Film durch die Formel:

(d) Bestimmen der Nichtreziprozitäts-Koeffizienten CNRD(di), die durch die Photonenrate di ausgedrückt werden, durch Anwenden der Formel:

CNRD(di) = CNRT(ti). (19)

2. Verfahren gemäß Anspruch 1, dadurch gekennzeichnet, daß es außerdem die folgende Operation enthält:

(e) Bestimmen der Modellfunktion der Koeffizienten CNRD(d), derart, daß

wobei die Parameter A'&sub0;, A'&sub1; und A'&sub2; anhand der Werte di und CNRD(dj) geschätzt werden.

3. Verfahren gemäß Anspruch 1 oder 2, dadurch gekennzeichnet, daß die Operation (a) die folgenden Operationen enthält:

(a&sub1;) Modifizieren des Heizstroms der Röhre (11) in der Weise, daß verschiedene Werte des Stroms erhalten werden;

(a&sub2;) Ablesen der Werte M(ti), die von der Integratorschaltung (16) für verschiedene Aufnahmezeiten (ti) ausgegeben werden, um so eine feste optische Dichte DOrefo des Films zu erhalten;

(a&sub3;) Berechnen des Verhältnisses

M(ti)/M(tref), (29)

das den Koeffizienten CNRT(ti) angibt, wobei M(tref) der Wert M(ti) für ti = tref ist.

4. Verfahren gemäß Anspruch 1 oder 2, dadurch gekennzeichnet, daß die Operation (a) die folgenden Operationen enthält:

(g1) Verwirklichen eines ersten Sensitogramms Srefo mit Hilfe eines Sensitometers mit veränderlicher Zeit, wenn die Aufnahmezeit für eine Referenzzeit trefo gesteuert wird;

(g2) Verwirklichen von q Sensitogrammen S&sub1; bis Sq mit Hilfe desselben Sensitometers mit veränderlicher Zeit für q verschiedene Aufnahmezeiten ti;

(g3) Wählen einer optischen Referenzdichte DOrefo,

(g4) Messen der Beleuchtungsstufe EChrefo, Ech&sub1;,...,Echi, ...,Echq, die der optischen Dichte DOrefo entsprechen, für jedes Sensitogramm,

(g5) Berechnen der Koeffizienten CNRT(ti) durch die Gleichung:

mit K = 2/log&sub1;&sub0;(2).

5. Verfahren gemäß Anspruch 3 oder 4, dadurch gekennzeichnet, daß die Operation (a) außerdem die folgende Operation enthält:

Bestimmen der Modellfunktion der Koeffizienten CNRT (t), derart, daß

CNRD(t)=A&sub0;+A&sub1;logt+A&sub2;[logt]² (18)

wobei die Parameter A&sub0;, A&sub1; und A&sub2; anhand der Werte ti und CNRT(ti) geschätzt werden.

6. Verfahren gemäß einem der Ansprüche 1 bis 5, dadurch gekennzeichnet, daß die Operation (b) in einer Messung der Belichtung Lref mittels einer Detektorzelle besteht.

7. Verfahren gemäß einem der Ansprüche 1 bis 5, dadurch gekennzeichnet, daß die Referenzbelichtung Lref durch ein Eichverfahren bestimmt wird, das in dem Fall, in dem der Bildempfänger (17) aus wenigstens einem Verstärkerschirm und einem für das Licht dieses Schirms empfindlichen Film aufgebaut ist, die folgenden Operationen enthält:

- Verwirklichen einer Aufnahme unter bestimmten radiologischen Bedingungen für eine optische Referenzdichte DOrefo, ein Normal der Dicke Eo, eine Versorgungsspannung Vo, eine Aufnahmezeit to und einen Wert des Produkts Io × to wobei Io der Röhrenstrom ist;

- Messen des Wirkungsgrades Do;

- Berechnen der äquivalenten Dicke Epo durch die Formel:

Epo = g"(Vo, Do) (7)

- Berechnen des Wirkungsgrades Dfo für den Film durch die Formel:

Dfo = f'(Vo, Epo) - f"(Vo, Epo) (8)

wobei f' bzw. f" Funktionen sind, die das Verhalten des Wirkungsgrades D in Abhängigkeit von den Parametern der Versorgungsspannung und der Dicke des Normals für die gegebene Konfiguration des Systems, in der der Empfänger (17) den Verstärkerschirm enthält, bzw. für die Konfiguration des Systems, in der der Empfänger keinen oder keine Verstärkerschirme enthält, beschreiben, und wobei g' die inverse Funktion von f" ist.

- Berechnen der Belichtung Lfilm für den Film durch die Formel:

Lfilm = Dfo × Io × to (23)

- Berechnen der Beleuchtungsstufe Echref, die der optischen Ref erenzdichte DOrefo entspricht, mit Hilfe der sensitoinetrischen Kurve Srefo;

- Messen der optischen Dichte DOm der erhaltenen Aufnahme und Berechnen der Beleuchtungsstufe Echm, die DOm entspricht, mit Hilfe der sensitometrischen Kurve Srefo;

- Berechnen der Referenzbelichtung Lref durch die Formel:

mit K = 2/log&sub1;&sub0;(2) (26)

8. Verfahren gemäß Anspruch 7, dadurch gekennzeichnet, daß für die Berechnung des Wirkungsgrades Dfo für den Film durch die Gleichung (8) die äquivalente Dicke Epo durch (Epo - (zusätzliches Filter)) ersetzt wird, wobei (zusätzliches Filter) die äquivalente Dicke aufgrund der zusätzlichen Filterung ist, die sich aus der Dämpfung zwischen dem Verstärkerschirm und der Detektorzelle (12) ergibt.

9. Verfahren gemäß Anspruch 7 für den Fall, in dem das Empfängerbild (17) aus einem für die Röntgenstrahlung empfindlichen Film gebildet ist und in dem die Detektorzelle (12) vor oder hinter dem Film angeordnet ist, dadurch gekennzeichnet, daß der Wirkungsgrad Dfo für den Film durch die folgende Formel berechnet wird:

D'fo = f"'(Vo, Ep), (30)

in der f"' eine Funktion ist, die das Verhalten des Wirkungsgrades (D) in Abhängigkeit von den Parametern der Versorgungsspannung und der Dicke des Normals für die gegebene Konfiguration des Systems beschreibt.

10. Verfahren gemäß Anspruch 9, dadurch gekennzeichnet, daß in der Gleichung (30) die Dicke Ep durch (Ep - (zusätzliches Filter)) ersetzt ist, wobei (zusätzliches Filter) die äquivalente Dicke aufgrund der zusätzlichen Filterung ist, die sich aus der Dämpfung zwischen dem Empfängerfilm (17) und der Detektorzelle (12) ergibt.

11. Verfahren gemäß einem der Ansprüche 7 bis 10, dadurch gekennzeichnet, daß die Leuchtdichte für den Film durch die Formel berechnet wird: