DE102015226693A1 - CT-Detektionsverfahren und CT-Vorrichtung - Google Patents

CT-Detektionsverfahren und CT-Vorrichtung Download PDF

Info

Publication number
DE102015226693A1
DE102015226693A1 DE102015226693.1A DE102015226693A DE102015226693A1 DE 102015226693 A1 DE102015226693 A1 DE 102015226693A1 DE 102015226693 A DE102015226693 A DE 102015226693A DE 102015226693 A1 DE102015226693 A1 DE 102015226693A1
Authority
DE
Germany
Prior art keywords
virtual sub
focus
projection
focuses
data
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Ceased
Application number
DE102015226693.1A
Other languages
English (en)
Inventor
Yongshun Xiao
Zhiqiang Chen
Ming Chang
Ziran Zhao
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
GRANPECT Co Ltd
Nuctech Co Ltd
Original Assignee
GRANPECT Co Ltd
Nuctech Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by GRANPECT Co Ltd, Nuctech Co Ltd filed Critical GRANPECT Co Ltd
Publication of DE102015226693A1 publication Critical patent/DE102015226693A1/de
Ceased legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G01N23/02Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
    • G01N23/04Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
    • G01N23/046Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/005Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00Investigating materials by wave or particle radiation
    • G01N2223/30Accessories, mechanical or electrical features
    • G01N2223/33Accessories, mechanical or electrical features scanning, i.e. relative motion for measurement of successive object-parts
    • G01N2223/3306Accessories, mechanical or electrical features scanning, i.e. relative motion for measurement of successive object-parts object rotates
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00Image generation
    • G06T2211/40Computed tomography
    • G06T2211/421Filtered back projection [FBP]

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Biochemistry (AREA)
  • Immunology (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pulmonology (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
  • Medical Informatics (AREA)
  • Quality & Reliability (AREA)
  • Computer Vision & Pattern Recognition (AREA)

Abstract

Es wird ein CT-Detektionsverfahren bereitgestellt, das folgende Schritte umfasst: (1) Scannen in Umfangsrichtung eines zu detektierenden Objekts mittels Röntgenstrahlen gemäß einem voreingestellten Winkelabtastwert, der die Anzahl von Abtastpunkten in einem einzelnen Kreis darstellt, um eine Gruppe von Projektionsabtastdaten in verschiedenen Projektionswinkeln zu erhalten, wobei der voreingestellte Winkelabtastwert größer als 1000 ist; (2) Verarbeiten der Projektionsabtastdaten dergestalt, dass Projektionsdaten mehrerer virtueller Sub-Fokusse erhalten werden, die einem großen Fokus einer Strahlungsquelle in einem CT-System entsprechen; und (3) Implementieren einer Bildrekonstruktion gemäß den Projektionsdaten der mehreren virtuellen Sub-Fokusse.

Description

  • QUERVERWEIS AUF VERWANDTE ANMELDUNG
  • Diese Anmeldung beansprucht die Priorität der chinesischen Patentanmeldung Nr. 201410851690.1, eingereicht am 31. Dezember 2014 beim Staatlichen Amt für geistiges Eigentum von China, deren gesamte Offenbarung hiermit durch Bezugnahme in den vorliegenden Text aufgenommen wird.
  • HINTERGRUND DER ERFINDUNG
  • Gebiet der Erfindung
  • Die vorliegende Offenbarung betrifft ein CT-Detektionsverfahren und eine CT-Vorrichtung und betrifft insbesondere ein Bildrekonstruktionsverfahren in einem Röntgen-CT-Bildgabesystem.
  • Beschreibung des Standes der Technik
  • Ein CT-System wird in der Regel in medizinischen Bildgabevorrichtungen und industriellen Detektionsanordnungen verwendet. Es scannt ein zu detektierendes Objekt mittels Röntgenstrahlung ab, bildet das zu detektierende Objekt unter Verwendung von Scan-Daten ab und detektiert dadurch das Innere des zu detektierenden Objekts. Das CT-System muss eine höhere räumliche Auflösung besitzen, falls Messungen von geometrischen Abmessungen des zu detektierenden Objekts und eine Detektion kleinerer Strukturen des zu detektierenden Objekts mittels des CT-Systems implementiert werden sollen.
  • Das Patentdokument CN102525527A offenbart ein Verfahren zur Verarbeitung von Projektionsdaten, bei dem eine Röntgenberechnungs-Tomografievorrichtung für die Bildgabe verwendet wird. Die Vorrichtung umfasst ein Portal, das eine Röntgenstrahlungsquelle und einen zweidimensionales-Array-Typ-Röntgenstrahlungsdetektor trägt. Die Röntgenstrahlungsquelle und der zweidimensionales-Array-Typ-Röntgenstrahlungsdetektor sind an einem rotierenden Ring auf einander gegenüberliegenden Seiten eines zu detektierenden Probanden installiert.
  • Gemäß dem in dem oben genannten Patentdokument offenbarten Verfahren werden die Projektionsdaten unter Verwendung eines Konusstrahls entlang einer kreisförmigen Trajektorie erfasst; ein Gewichtswert für die Projektionsdaten wird gemäß einer normalisierten Electrocardiogram Gated Reconstruction(EGR)-Gewichtungsfunktion WEGR(β, γ) auf der Basis einer EGR-Gewichtungsfunktion UEGR(Φ(β)) bestimmt; die UEGR(Φ(β)) wird durch eine Summierung von UEGR(Φ(βc n))) von n = –NPI zu n = NPI normalisiert, wobei Φ eine Herzphase ist, β ein Betrachtungswinkel ist, βc n ein komplementärer Betrachtungswinkel ist, PI = π; und jede der Projektionsdaten wird durch den Gewichtswert gewichtet, der durch das WEGR(β, γ) bestimmt wurde, wobei γ ein Fächerwinkel des Konusstrahls ist.
  • Jedoch kann die räumliche Auflösung des Bildes, das unter Verwendung des obigen Verfahrens erhalten wird, nicht die Anforderung für die oben angesprochene Anwendung erfüllen.
  • Die räumliche Auflösung ist in der Tat ein wichtiger Parameter für ein Röntgen-CT-System, der durch einige Hardware-Faktoren beschränkt wird, wie zum Beispiel die Fokusgröße einer Strahlenquelle, die Abmessungen eines Detektors, das Vergrößerungsverhältnis eines Scan-Systems und dergleichen.
  • Das Verstärkungsverhältnis des Scan-Systems kann anhand verschiedener Systemparameter sowie der Abmessungen des zu detektierenden Objekts gewählt werden. Infolge dessen bilden die Fokusgröße der Strahlenquelle und die Abmessungen des Detektors Faktoren, welche die räumliche Auflösung des Systems begrenzen. Da Röntgenstrahlen eine relativ hohe Energie besitzen, hat des Weiteren in einem industriellen CT-System der Fokus der Strahlenquelle in der Regel eine Größe im Bereich von 1,0 mm bis 2,0 mm, um die Anforderung an die Wärmeableitung von dem Ziel zu erfüllen. Jedoch liegen die Abmessungen der Detektoreinheit in der Größenordnung des Submillimeterbereichs, und dementsprechend bildet die Fokusgröße der Strahlenquelle einen wesentlichen Faktor bei der Beschränkung der räumlichen Auflösung des CT-Systems.
  • KURZFASSUNG DER ERFINDUNG
  • Um die oben angesprochenen Probleme zu überwinden, wird ein CT-Detektionsverfahren bereitgestellt, das folgende Schritte umfasst: (1) Scannen in Umfangsrichtung eines zu detektierenden Objekts mittels Röntgenstrahlen gemäß einem voreingestellten Winkelabtastwert, der die Anzahl von Abtastpunkten in einem einzelnen Kreis darstellt, um eine Gruppe von Projektionsabtastdaten in verschiedenen Projektionswinkeln zu erhalten, wobei der voreingestellte Winkelabtastwert größer als 1000 ist; (2) Verarbeiten der Projektionsabtastdaten dergestalt, dass Projektionsdaten mehrerer virtueller Sub-Fokusse erhalten werden, die einem großen Fokus einer Strahlungsquelle in einem CT-System entsprechen; und (3) Implementieren einer Bildrekonstruktion gemäß den Projektionsdaten der mehreren virtuellen Sub-Fokusse.
  • In einer Ausführungsform der vorliegenden Erfindung ist der voreingestellte Winkelabtastwert größer als 8000.
  • In einer Ausführungsform der vorliegenden Erfindung wird der voreingestellte Winkelabtastwert N gemäß der folgenden Formel erhalten: N > 2π/tan–1(a/R), wobei a die Fokusgröße einer Röntgenstrahlungsquelle ist und R eine Distanz zwischen der Röntgenstrahlungsquelle und einem Drehmittelpunkt des zu detektierenden Objekts ist.
  • In einer Ausführungsform der vorliegenden Erfindung werden in Schritt (2) die Projektionsdaten der mehreren virtuellen Sub-Fokusse gemäß den Projektionsabtastdaten in verschiedenen Projektionswinkeln und einer Intensität, die den mehreren virtuellen Sub-Fokussen entspricht, erhalten.
  • In einer Ausführungsform der vorliegenden Erfindung werden die Projektionsdaten der mehreren virtuellen Sub-Fokusse unter Verwendung der folgenden Formel erhalten:
    Figure DE102015226693A1_0002
    wobei D eine Distanz zwischen der Röntgenstrahlungsquelle und einem Detektor ist, ai eine Distanz eines virtuellen Sub-Fokus ist, der von einer idealen Punktquelle abweicht, M die Anzahl der virtuellen Sub-Fokusse ist, t die Position des Detektors darstellt, i ein tiefgestellter Index eines virtuellen Sub-Fokus ist, βi die Projektionswinkelposition eines einzelnen virtuellen Sub-Fokus ist und q(k, t) ein Dämpfungskoeffizient des Röntgenstrahls ist, der durch die folgenden Formel erhalten wird: q(k, t) = –ln[I(k, t)/I0(t)]
  • In einer Ausführungsform der vorliegenden Erfindung wird vor dem Schritt (1) der große Fokus des CT-Systems mittels eines Fokusgrößenmessverfahrens in mehrere virtuelle Sub-Fokusse fein eingeteilt, und die Intensität, die den mehreren virtuellen Sub-Fokussen entspricht, wird erhalten.
  • In einer Ausführungsform der vorliegenden Erfindung umfasst der Schritt zum Erhalten der Intensität, die den mehreren virtuellen Sub-Fokussen entspricht, Folgendes:
    Messen der Größe des großen Fokus des CT-Systems mittels eines Kleinloch-Bildgabeverfahrens;
    Feines Einteilen des großen Fokus in die mehreren äquivalenten virtuellen Sub-Fokusse; und
    Bestimmen der Intensität, die den mehreren äquivalenten virtuellen Sub-Fokussen entspricht, gemäß einer Verteilungskurve der Fokusintensität.
  • In einer Ausführungsform der vorliegenden Erfindung umfasst das Verfahren des Weiteren folgenden Schritt:
    Korrigieren der Projektionsabtastdaten vor dem Schritt (2).
  • In einer Ausführungsform der vorliegenden Erfindung umfasst der Schritt zum Korrigieren der Projektionsabtastdaten Folgendes:
    Messen der Intensität I0(t) des durch die Luft gedämpften Strahls;
    Berechnen des Dämpfungskoeffizienten q(k, t) des Röntgenstrahls in verschiedenen Winkeln mittels der folgenden Formel: q(k, t) = –ln[I(k, t)/I0(t)], wobei k ein Winkelabtastzählwert ist, der von einem anfänglichen Wert 1 aus größer wird, bis der Wert k größer ist als ein voreingestellter Winkelabtastzählwert.
  • Es wird außerdem eine CT-Vorrichtung bereitgestellt, die Folgendes umfasst: eine Strahlungsquelle, einen Detektor, welcher der Strahlungsquelle entspricht, eine Voreinstellwinkelabtastwert-Bestimmungseinheit, eine Projektionsabtastdaten-Verarbeitungseinheit und eine Rekonstruktionseinheit, wobei die Voreinstellwinkelabtastwert-Bestimmungseinheit dafür vorgesehen ist, einen voreingestellten Winkelabtastwert zu bestimmen, der die Anzahl von Abtastpunkten in einem einzelnen Kreis darstellt, wobei der voreingestellte Winkelabtastwert größer als 1000 ist; wobei die Strahlungsquelle und der Detektor auf jeder Seite eines zu detektierenden Objekts angeordnet sind und dazu dienen, das zu detektierende Objekt mittels Röntgenstrahlen gemäß dem voreingestellten Winkelabtastwert in Umfangsrichtung zu scannen, um eine Gruppe von Projektionsabtastdaten in verschiedenen Projektionswinkeln zu erhalten; wobei die Projektionsabtastdaten-Verarbeitungseinheit dazu ausgelegt ist, die Projektionsabtastdaten dergestalt zu verarbeiten, dass Projektionsdaten mehrerer virtueller Sub-Fokusse erhalten werden, die einem großen Fokus einer Strahlungsquelle in einem CT-System entsprechen; und die Rekonstruktionseinheit dazu ausgelegt ist, eine Bildrekonstruktion gemäß den Projektionsdaten der mehreren virtuellen Sub-Fokusse zu implementieren.
  • In einer Ausführungsform der vorliegenden Erfindung ist der voreingestellte Winkelabtastwert größer als 8000, und der voreingestellte Winkelabtastwert N wird gemäß der folgenden Formel erhalten: N > 2π/tan–1(a/R), wobei a die Fokusgröße einer Röntgenstrahlungsquelle ist und R eine Distanz zwischen der Röntgenstrahlungsquelle und einem Drehmittelpunkt des zu detektierenden Objekts ist.
  • In einer Ausführungsform der vorliegenden Erfindung erhält die Projektionsabtastdaten-Verarbeitungseinheit die Projektionsdaten der mehreren virtuellen Sub-Fokusse gemäß den Projektionsabtastdaten in verschiedenen Projektionswinkeln und einer Intensität, die den mehreren virtuellen Sub-Fokussen entspricht, und die Projektionsdaten der mehreren virtuellen Sub-Fokusse werden mittels der folgenden Formel erhalten:
    Figure DE102015226693A1_0003
    wobei D eine Distanz zwischen der Röntgenstrahlungsquelle und dem Detektor ist, ai eine Distanz eines einzelnen virtuellen Sub-Fokus ist, der von einer idealen Punktquelle abweicht, M die Anzahl der virtuellen Sub-Fokusse ist, t die Position des Detektors darstellt, i ein tiefgestellter Index eines einzelnen virtuellen Sub-Fokus ist, βi die Projektionswinkelposition eines einzelnen virtuellen Sub-Fokus ist und q(k, t) ein Dämpfungskoeffizient des Röntgenstrahls ist, der durch die folgenden Formel erhalten wird: q(k, t) = –ln[I(k, t)/I0(t)].
  • In einer Ausführungsform der vorliegenden Erfindung umfasst die Vorrichtung des Weiteren eine Fokusgrößenmessungseinheit, eine Fokusteilungseinheit und eine Intensitätsbestimmungseinheit, wobei die Fokusgrößenmessungseinheit dafür vorgesehen ist, die Größe des großen Fokus des CT-Systems mittels eines Kleinloch-Bildgabeverfahrens zu messen; wobei die Fokusteilungseinheit dafür vorgesehen ist, den großen Fokus in die mehreren äquivalenten virtuellen Sub-Fokusse fein einzuteilen; und die Intensitätsbestimmungseinheit dafür vorgesehen ist, eine Intensität, die den mehreren äquivalenten virtuellen Sub-Fokussen entspricht, gemäß einer Verteilungskurve der Fokusintensität zu bestimmen.
  • In einer Ausführungsform der vorliegenden Erfindung umfasst die Vorrichtung des Weiteren eine Datenkorrektureinheit zum Korrigieren der Projektionsabtastdaten, wobei auf der Basis der gemessenen Intensität I0(t) des durch die Luft gedämpften Strahls Dämpfungskoeffizienten q(k, t) des Röntgenstrahls in verschiedenen Winkeln mittels der folgenden Formel berechnet werden: q(k, t) = –ln[I(k, t)/I0(t)], wobei k ein Winkelabtastzählwert ist, der von einem anfänglichen Wert 1 aus größer wird, bis der Wert k größer ist als ein voreingestellter Winkelabtastzählwert.
  • Gemäß dem in der vorliegenden Offenbarung bereitgestellten Verfahren werden die räumliche Auflösung verbessert, die Einschränkung der Fokusgröße der Röntgenstrahlungsquelle beseitigt, die sich auf die räumliche Auflösung auswirkt, und der Scan-Modus des Systems vereinfacht, wodurch eine höher Strahlausnutzung, eine höhere räumliche Auflösung und eine höhere Bildrekonstruktionseffizienz erreicht werden.
  • KURZE BESCHREIBUNG DER ZEICHNUNGEN
  • Um die technischen Lösungen in Ausführungsformen der vorliegenden Erfindung oder im Stand der Technik deutlicher zu veranschaulichen, werden im Folgenden kurz die beiliegenden Zeichnungen zum Beschreiben der Ausführungsformen der vorliegenden Erfindung oder des Standes der Technik vorgestellt. Es versteht sich, dass die beiliegenden Zeichnungen nur einige Ausführungsformen der vorliegenden Erfindung sind und dass der Fachmann andere Zeichnungen aus diesen beiliegenden Zeichnungen erhalten kann, ohne einen erfinderischen Schritt tun zu müssen.
  • 1 ist eine Strukturdarstellung eines Röntgen-CT-Bildgabesystems gemäß einer Ausführungsform der vorliegenden Erfindung;
  • 2 ist eine Ansicht, die eine Fokusteilung des in der vorliegenden Offenbarung bereitgestellten Verfahrens zeigt; und
  • 3 ist eine schematische Ansicht, die eine Bestimmung der Fokusgröße und von virtuellen Fokussen mittels eines Lochblendenverfahrens zeigt, das durch das in der vorliegenden Offenbarung bereitgestellte Verfahren verwendet wird.
  • DETAILLIERTE BESCHREIBUNG BEVORZUGTER AUSFÜHRUNGSFORMEN DER ERFINDUNG
  • Die detaillierten Ausführungsformen der vorliegenden Erfindung werden im Folgenden in Kombination mit den beiliegenden Figuren und den Beispielen ausführlich beschrieben. Die Ausführungsformen sind lediglich beispielhaft und sollen nicht den Schutzumfang der vorliegenden Erfindung einschränken.
  • Es wird ein CT-Detektionsverfahren bereitgestellt, das folgende Schritte umfasst: (1) Scannen in Umfangsrichtung eines zu detektierenden Objekts mittels Röntgenstrahlen gemäß einem voreingestellten Winkelabtastwert, der die Anzahl von Abtastpunkten in einem einzelnen Kreis darstellt, um eine Gruppe von Projektionsabtastdaten in verschiedenen Projektionswinkeln zu erhalten, wobei der voreingestellte Winkelabtastwert größer als 1000 ist; zum Beispiel kann der Wert 2048, 3072, 4096, 5120, 6144, 7168 usw. sein; (2) Verarbeiten der Projektionsabtastdaten dergestalt, dass Projektionsdaten mehrerer virtueller Sub-Fokusse erhalten werden, die einem großen Fokus einer Strahlungsquelle in einem CT-System entsprechen; und (3) Implementieren einer Bildrekonstruktion gemäß den Projektionsdaten der mehreren virtuellen Sub-Fokusse.
  • In einer bevorzugten Ausführungsform kann der voreingestellte Winkelabtastwert größer als 8000 sein. Der voreingestellte Winkelabtastwert N wird gemäß der folgenden Formel erhalten: N > 2π/tan–1(a/R), wobei a die Fokusgröße einer Röntgenstrahlungsquelle ist und R eine Distanz zwischen der Röntgenstrahlungsquelle und einem Drehmittelpunkt des zu detektierenden Objekts ist. Bevorzugt werden in Schritt (2) die Projektionsdaten der mehreren virtuellen Sub-Fokusse gemäß den Projektionsabtastdaten in verschiedenen Projektionswinkeln und einer Intensität, die den mehreren virtuellen Sub-Fokussen entspricht, erhalten. Diese Formel funktioniert auch in dem Fall, dass der voreingestellte Winkelabtastwert in anderen Zahlenbereichen liegt.
  • In einer bevorzugten Ausführungsform werden die Projektionsdaten der mehreren virtuellen Sub-Fokusse unter Verwendung der folgenden Formel erhalten:
    Figure DE102015226693A1_0004
    Figure DE102015226693A1_0005
    wobei D eine Distanz zwischen der Röntgenstrahlungsquelle und einem Detektor ist, ai eine Distanz eines einzelnen virtuellen Sub-Fokus ist, der von einer idealen Punktquelle abweicht, M die Anzahl der virtuellen Sub-Fokusse ist, t die Position des Detektors darstellt, i ein tiefgestellter Index eines einzelnen virtuellen Sub-Fokus ist, βi die Projektionswinkelposition eines einzelnen virtuellen Sub-Fokus ist und q(k, t) ein Dämpfungskoeffizient des Röntgenstrahls ist, der durch die folgenden Formel erhalten wird: q(k, t) = –ln[I(k, t)/I0(t)].
  • Bevorzugt wird vor dem Schritt (1) der große Fokus des CT-Systems mittels eines Fokusgrößenmessverfahrens in mehrere virtuelle Sub-Fokusse fein eingeteilt, und die Intensität, die den mehreren virtuellen Sub-Fokussen entspricht, wird erhalten.
  • Bevorzugt umfasst der Schritt zum Erhalten der Intensität, die den mehreren virtuellen Sub-Fokussen entspricht, Folgendes: Messen der Größe des großen Fokus des CT-Systems mittels eines Kleinloch-Bildgabeverfahrens; Feines Einteilen des großen Fokus in die mehreren äquivalenten virtuellen Sub-Fokusse; und Bestimmen der Intensität, die den mehreren äquivalenten virtuellen Sub-Fokussen entspricht, gemäß einer Verteilungskurve der Fokusintensität.
  • In einer Ausführungsform der vorliegenden Erfindung umfasst das Verfahren des Weiteren folgenden Schritt: Korrigieren der Projektionsabtastdaten vor dem Schritt (2). Insbesondere umfasst der Schritt zum Korrigieren der Projektionsabtastdaten Folgendes: Messen der Intensität I0(t) des durch die Luft gedämpften Strahls; Berechnen des Dämpfungskoeffizienten q(k, t) des Röntgenstrahls in verschiedenen Winkeln mittels der folgenden Formel: q(k, t) = –ln[I(k, t)/I0(t)], wobei k ein Winkelabtastzählwert ist, der von einem anfänglichen Wert 1 aus größer wird, bis der Wert k größer ist als ein voreingestellter Winkelabtastzählwert.
  • In einer weiteren Ausführungsform der vorliegenden Erfindung können alternativ zuerst die Projektionsdaten der mehreren virtuellen Sub-Fokusse wiederhergestellt werden, dann werden die Daten korrigiert, und zum Schluss wird die Bildrekonstruktion implementiert.
  • Gemäß einer bevorzugten Ausführungsform der vorliegenden Erfindung wird eine CT-Vorrichtung bereitgestellt, die Folgendes umfasst: eine Strahlungsquelle, einen Detektor, der der Strahlungsquelle entspricht, eine Voreinstellwinkelabtastwert-Bestimmungseinheit, eine Projektionsabtastdaten-Verarbeitungseinheit und eine Rekonstruktionseinheit, wobei die Voreinstellwinkelabtastwert-Bestimmungseinheit dafür vorgesehen ist, einen voreingestellten Winkelabtastwert zu bestimmen, der die Anzahl von Abtastpunkten in einem einzelnen Kreis darstellt, wobei der voreingestellte Winkelabtastwert größer als 1000 ist; wobei die Strahlungsquelle und der Detektor auf jeder Seite eines zu detektierenden Objekts angeordnet sind und dazu dienen, das zu detektierende Objekt mittels Röntgenstrahlen gemäß dem voreingestellten Winkelabtastwert in Umfangsrichtung zu scannen, um eine Gruppe von Projektionsabtastdaten in verschiedenen Projektionswinkeln zu erhalten; wobei die Projektionsabtastdaten-Verarbeitungseinheit dazu ausgelegt ist, die Projektionsabtastdaten dergestalt zu verarbeiten, dass Projektionsdaten mehrerer virtueller Sub-Fokusse erhalten werden, die einem großen Fokus einer Strahlungsquelle in einem CT-System entsprechen; und die Rekonstruktionseinheit dazu ausgelegt ist, eine Bildrekonstruktion gemäß den Projektionsdaten der mehreren virtuellen Sub-Fokusse zu implementieren.
  • Bevorzugt ist der voreingestellte Winkelabtastwert größer als 8000, und der voreingestellte Winkelabtastwert N wird gemäß der folgenden Formel erhalten: N > 2π/tan–1(a/R), wobei a die Fokusgröße einer Röntgenstrahlungsquelle ist und R eine Distanz zwischen der Röntgenstrahlungsquelle und einem Drehmittelpunkt des zu detektierenden Objekts ist.
  • Bevorzugt erhält die Projektionsabtastdaten-Verarbeitungseinheit die Projektionsdaten der mehreren virtuellen Sub-Fokusse gemäß den Projektionsabtastdaten in verschiedenen Projektionswinkeln und einer Intensität, die den mehreren virtuellen Sub-Fokussen entspricht, und die Projektionsdaten der mehreren virtuellen Sub-Fokusse werden mittels der folgenden Formel erhalten:
    Figure DE102015226693A1_0006
    Figure DE102015226693A1_0007
    wobei D eine Distanz zwischen der Röntgenstrahlungsquelle und dem Detektor ist, ai eine Distanz eines einzelnen virtuellen Sub-Fokus ist, der von einer idealen Punktquelle abweicht, M die Anzahl der virtuellen Sub-Fokusse ist, t die Position des Detektors darstellt, i ein tiefgestellter Index eines einzelnen virtuellen Sub-Fokus ist, βi die Projektionswinkelposition eines einzelnen virtuellen Sub-Fokus ist und q(k, t) ein Dämpfungskoeffizient des Röntgenstrahls ist, der durch die folgenden Formel erhalten wird: q(k, t) = –ln[I(k, t)/I0(t)].
  • In einer bevorzugten Ausführungsform der vorliegenden Erfindung umfasst die Vorrichtung des Weiteren eine Fokusgrößenmessungseinheit, eine Fokusteilungseinheit und eine Intensitätsbestimmungseinheit, wobei die Fokusgrößenmessungseinheit dafür vorgesehen ist, die Größe des großen Fokus des CT-Systems mittels eines Kleinloch-Bildgabeverfahrens zu messen; wobei die Fokusteilungseinheit dafür vorgesehen ist, den großen Fokus in die mehreren äquivalenten virtuellen Sub-Fokusse fein einzuteilen; und die Intensitätsbestimmungseinheit dafür vorgesehen ist, eine Intensität, die den mehreren äquivalenten virtuellen Sub-Fokussen entspricht, gemäß einer Verteilungskurve der Fokusintensität zu bestimmen.
  • In einer bevorzugten Ausführungsform der vorliegenden Erfindung umfasst die Vorrichtung des Weiteren eine Datenkorrektureinheit zum Korrigieren der Projektionsabtastdaten, wobei auf der Basis der gemessenen Intensität I0(t) des durch die Luft gedämpften Strahls Dämpfungskoeffizienten q(k, t) des Röntgenstrahls in verschiedenen Winkeln mittels der folgenden Formel berechnet werden: q(k, t) = –ln[I(k, t)/I0(t)], wobei k ein Winkelabtastzählwert ist, der von einem anfänglichen Wert 1 aus größer wird, bis der Wert k größer ist als ein voreingestellter Winkelabtastzählwert.
  • Die oben beschriebenen Ausführungsformen sollen die vorliegende Erfindung nicht einschränken, und der Fachmann kann die vorliegende Erfindung gewünschtenfalls abändern oder modifizieren.
  • Im Folgenden werden Ausführungsformen der vorliegenden Erfindung vorgestellt, welche die oben beschriebenen Ausführungsformen abändern und in denen ein Bildrekonstruktionsverfahren in einem Röntgen-CT-Bildgabesystem bereitgestellt wird, das folgende Schritte umfasst:
    • (1) Aktivieren eines Strahlgenerators und eines Detektors in dem Röntgen-CT-Bildgabesystem, Messen der Intensität I0(t) des durch die Luft gedämpften Strahls an einer Position t und Übermitteln der gemessenen Daten an einen Datenprozessor in dem Röntgen-CT-Bildgabesystem;
    • (2) Anordnen eines zu detektierenden Objekts in dem Röntgen-CT-Bildgabesystem zum Scannen, Messen der Intensität I(k, t) des Strahls, der durch das zu detektierende Objekt gedämpft wird, in der momentanen Winkelposition an der Position t und Übermitteln der gemessenen Daten an den Datenprozessor in dem Röntgen-CT-Bildgabesystem, wobei k ein Winkelabtastzählwert ist, der von einem anfänglichen Wert 1 aus größer wird, bis der Wert k größer ist als ein voreingestellter Winkelabtastzählwert;
    • (3) Berechnen von Dämpfungskoeffizienten q(k, t) des Röntgenstrahls in verschiedenen Winkeln unter Verwendung der folgenden Formel (1): q(k, t) = –ln[I(k, t)/I0(t)] (1)
    • (4) Rekonstruieren eines Lineardämpfungskoeffizienten-Verteilungsbildes des zu detektierenden Objekts mittels eines gefilterten Rückprojektionsalgorithmus auf der Basis der Daten q(k, t) in Schritt (3).
  • In 1 ist ein CT-Bildgabesystem gezeigt, das in der Ausführungsformen der vorliegenden Erfindung verwendet wird. Das System enthält Folgendes:
    eine Röntgenstrahlungs-Generierungsvorrichtung, die ein Fokusziel mit einer speziellen Größe zum Generieren von Röntgenstrahlung und einen Röntgenstrahlaustritt umfasst;
    eine mechanische Bewegungsvorrichtung und eine Steuereinheit, wobei die mechanische Bewegungsvorrichtung und die Steuereinheit Folgendes umfassen: einen Objektträger, einen Rahmen, der mit einem Röntgenapparat und einem Detektor versehen ist, ein Steuerungssystem und einen Hochpräzisions-Winkelcodierer; wobei der Objektträger in einer kreisförmigen Orbitalrotationsbewegung und/oder einer linearen Translationsbewegung angetrieben wird und der Rahmen in einer linearen Translationsbewegung und/oder einer kreisförmigen Orbitalrotationsbewegung angetrieben wird, wodurch eine in Umfangsrichtung verlaufende Scan-Trajektorie erreicht wird; wobei in der Ausführungsform der vorliegenden Erfindung eine Fächerstrahl-Umfangsscan-Trajektorie unter Verwendung eines ortsfesten Rahmens und eines drehbaren Objektträgers als ein Beispiel zur Beschreibung entwickelt wird und der Hochpräzisions-Winkelcodierer dafür verwendet wird, eine präzise Winkelpositionierung und eine hoch-dichte Winkelabtastung zu erreichen;
    ein Datenerfassungssystem, das eine oder mehrere Reihen von Lineararray-Detektoreinheiten und einen Lesekreis, einen Auslösererfassungssignalkreis und einen Datenübertragungskreis umfasst;
    einen primären Steuerungs- und Datenverarbeitungscomputer zum Steuern des Betriebsprozesses des CT-Systems, einschließlich einer mechanischen Rotation, einer elektrischen Steuerung und einer Sicherheitsverriegelungssteuerung.
  • Wie dies oben beschrieben ist, werden gemäß dem Verfahren der vorliegenden Erfindung in einer bevorzugten Ausführungsform die folgenden Schritte vor dem Implementieren von Schritt (1) implementiert:
    Messen und Erhalten der Fokusgröße einer Röntgenstrahlungsquelle und Messen und Erhalten einer Distanz R zwischen der Röntgenstrahlungsquelle und einem Drehmittelpunkt eines Drehtisches in dem Röntgen-CT-Bildgabesystem und Bestimmen eines voreingestellten Winkelabtastzählwertes N gemäß der folgenden Formel, auf der Basis von a und R: N > 2π/tan–1(a/R).
  • Bei Verwendung der oben beschriebenen Verfahrensweise übersteigt die Anzahl der Abtastungen die traditionelle CT-Vorrichtung um eine Größenordnung, wodurch die räumliche Auflösung des rekonstruierten Bildes signifikant verbessert wird.
  • In einer Ausführungsform der vorliegenden Erfindung umfasst das Bildrekonstruktionsverfahren des Weiteren folgenden Schritt: Bestimmen, ob der Wert k den voreingestellten Winkelabtastzählwert N übersteigt, wenn Schritt (2) implementiert wird, dann:
    wenn k größer ist als N, Implementieren von Schritt (3);
    wenn k kleiner ist als N, Drehen des zu detektierenden Objekts zum nächsten Scan-Winkel und Fortsetzen der Implementierung von Schritt (2).
  • In einer Ausführungsform der vorliegenden Erfindung umfasst das Bildrekonstruktionsverfahren des Weiteren folgende Schritte: Teilen der Röntgenstrahlungsquelle in mehrere virtuelle Sub-Fokusse und Bestimmen der Position ai und der Intensität w(ai) jedes virtuellen Sub-Fokus;
    nach der Implementierung von Schritt (3), Erhalten von Wiederherstellungsdaten pciri, ti) der virtuellen Sub-Fokusse zum Implementieren einer Bildrekonstruktion in dem gefilterten Rückprojektionsalgorithmus mittels der folgenden Formel:
    Figure DE102015226693A1_0008
    wobei D eine Distanz zwischen der Röntgenstrahlungsquelle und dem Detektor ist, ai eine Distanz eines einzelnen virtuellen Sub-Fokus ist, der von einer idealen Punktquelle s abweicht, M die Anzahl der virtuellen Sub-Fokusse ist, t die Position des Detektors darstellt, i ein tiefgestellter Index eines einzelnen virtuellen Sub-Fokus ist, βi die Projektionswinkelposition eines einzelnen virtuellen Sub-Fokus ist; und
    Rekonstruieren eines Lineardämpfungskoeffizienten-Verteilungsbildes des zu detektierenden Objekts mittels eines gefilterten Rückprojektionsalgorithmus auf der Basis der Wiederherstellungsdaten pciri, ti).
  • Gemäß einer Ausführungsform der vorliegenden Erfindung wird zum Verbessern der räumlichen Auflösung des CT-Systems der Fokus bevorzugt in mehrere virtuelle Sub-Fokusse geteilt, und ein äquivalenter idealer Fokus (die ideale Punktquelle) wird auf der Basis der virtuellen Sub-Fokusse erhalten, wie dies in 2 gezeigt ist, welche die Fokusteilung veranschaulicht.
  • Der große Fokus des CT-Bildgabesystems wird gleichmäßig in mehrere virtuelle Sub-Fokusse Si geteilt. Die Projektionsdaten q(k, t) des großen Fokus, der durch die Detektoreinheit Det in einem Projektionswinkel k(2π/N) gemessen wurde, können einer Akkumulation der Projektionsdaten der geteilten virtuellen Sub-Fokusse äquivalent sein. Gemäß der in 2 gezeigten geometrischen Beziehung ist ein Strahlenweg, der durch den virtuellen Sub-Fokus Si und die Detektoreinheit Det bestimmt wird, mit einem Strahlenweg identisch, der durch die ideale Punktquelle in einem Projektionswinkel βi und die Detektoreinheit Deti im Fall eines kreisförmigen Orbit-Scans bestimmt wird. Ein Winkel zwischen dem Projektionswinkel βi und dem momentanen Projektionswinkel k(2π/N) wird als θi bezeichnet und wird durch die Distanz ai des virtuellen Sub-Fokus Si, der von der idealen Punktquelle S abweicht, und einen Scan-Radius R des kreisförmigen Orbits bestimmt. Somit können die Projektionsdaten q(k, t) des großen Fokus des CT-Bildgabesystems im Fall des idealen kreisförmigen Orbit-Scans als eine gewichtete Summe der Projektionsdaten pciri, ti) der virtuellen Sub-Fokusse dargestellt werden, wobei das Gewicht die Intensität w(ai) der virtuellen Sub-Fokusse ist. Durch die obige Vorgehensweise wird die Wiederherstellung der Projektionsdaten der virtuellen Sub-Fokusse erreicht, wodurch eine hochauflösende CT-Bildrekonstruktion erreicht wird.
  • Wie dies oben beschrieben ist, wird in dem in der vorliegenden Erfindung bereitgestellten Verfahren die Bildrekonstruktion durch den gefilterten Rückprojektionsalgorithmus implementiert.
  • Die detaillierten Implementierungsschritte der Rekonstruktion unter Verwendung des gefilterten Rückprojektionsalgorithmus sind folgende:
    erstens, Gewichten der Projektionsdaten:
    Figure DE102015226693A1_0009
    dann Implementieren einer eindimensionalen Rampenfilterung in der Richtung des Detektors:
    Figure DE102015226693A1_0010
    wobei t ein Abtastintervall in der Richtung des Detektors ist, was in der Regel eine Detektoreinheit-Beabstandung ist;
    Implementieren einer gewichteten Rückprojektion an den gefilterten Daten (in der Winkelrichtung):
    Figure DE102015226693A1_0011
    wobei ƒ(r, ϕ) eine polare Koordinatendarstellung des zu rekonstruierenden Bildes ist und Δβ = 2π/N ein Abtastintervall des Projektionswinkels ist.
  • In der vorliegenden Offenbarung ist es möglich, das in 3 veranschaulichte Verfahren zu verwenden, um die Fokusgröße erhalten, wie dies insbesondere in 3 gezeigt ist:
    Verwendung eines Materials mit einem hohen Absorptionskoeffizienten, um eine Ablenkplatte zu bilden, und Ausbilden eines Lochs darin, dann Erhalten eines Bildes des Fokus mittels einer Kleinloch-Bildgabetheorie, dadurch Bestimmen der Fokusgröße a, der Position ai und der Intensität w(ai) des virtuellen Sub-Fokus durch Messen der Graustufenverteilung des Fokusbildes; wobei bi ein Abtastintervall des virtuellen Sub-Fokus ist, das in der Regel als bi = b/M definiert werden kann, wobei M die Anzahl der virtuellen Sub-Fokusse ist und Fi ein Grauwert des Fokusbildes bei bi ist.
  • Die oben angesprochenen Parameter erfüllen die folgende Beziehung:
    Figure DE102015226693A1_0012
  • Im Folgenden wird die ganz besonders bevorzugte Ausführungsform der vorliegenden Erfindung beschrieben. Es ist anzumerken, dass die zu beschreibende Ausführungsform nicht die vorliegende Erfindung einschränken soll, und jede der oben beschriebenen Ausführungsformen kann die Aufgabe des Verbesserns der räumlichen Auflösung auf verschiedenen technischen Wegen erfüllen. Die unten besprochene ganz besonders bevorzugte Ausführungsform kombiniert alle oben angesprochenen Mittel, um den optimalen technischen Effekt zu erreichen. Dadurch soll die vorliegende Erfindung allerdings nicht eingeschränkt werden; das heißt, jede der oben beschriebenen Ausführungsformen kann die Probleme lösen, die mit der vorliegenden Erfindung angegangen werden.
  • In der ganz besonders bevorzugte Ausführungsform:
    wird das in 1 gezeigte CT-Bildgabesystem verwendet. Um die beste räumliche Auflösung zu erhalten, umfasst das Bildrekonstruktionsverfahren folgende Schritte:
    • (1) Aktivieren der Röntgenstrahlungsquelle und des Detektors, Implementieren einer Fokusteilung und einer Fokusgrößenmessung unter Verwendung des in den 2 und 3 gezeigten Verfahrens dergestalt, dass eine Intensitätsverteilung des Fokus des Systems erhalten wird, dann Teilen der Röntgenstrahlungsquelle in mehrere virtuelle Sub-Fokusse; zum Beispiel kann der große Fokus in 5 virtuelle Sub-Fokusse geteilt werden; dadurch Bestimmen der Position ai und der Intensität w(ai) jedes virtuellen Sub-Fokus;
    • (2) Bestimmen der Winkelabtastungsanzahl N unter Verwendung der folgenden Formel, auf der Basis der gemessenen Fokusgröße a und der Scan-Geometrie-Parameter R: N > 2π/tan–1(a/R)
    • (3) Aussenden des Röntgenstrahls von der Röntgenstrahlungsquelle, Messen der Intensität des durch die Luft gedämpften Strahls, Umwandeln der Intensität des durch die Detektoreinheit empfangenen Röntgenstrahls an der Position t in einen Wert I0(t) und Senden dieser Daten an den Datenverarbeitungscomputer zum Speichern;
    • (4) Aussenden des Röntgenstrahls von der Röntgenstrahlungsquelle zum Implementieren eines Winkel-Scans des zu detektierenden Objekts, Umwandeln der Intensität des durch die Detektoreinheit empfangenen Röntgenstrahls an der Position t in einen Wert I(k, t) und Senden dieser Daten an den Datenverarbeitungscomputer, wobei k ein Winkelabtastzählwert mit einem anfänglichen Wert 1 ist;
    • (5) wenn der Wert K um 1 erhöht wird, dann Bestimmen, ob der Wert k den voreingestellten Winkelabtastzählwert N übersteigt; wenn nicht, dann Drehen des zu detektierenden Objekts zum nächsten Scan-Winkel und Zurückkehren zu (4), wenn ja, Gehen zu (6);
    • (6) vorübergehendes Korrigieren der Projektionsdaten gemäß der folgenden Formel, dergestalt, dass der Dämpfungskoeffizient q(k, t) des Röntgenstrahls, der durch jede Detektorseinheit in verschiedenen Projektionswinkeln empfangen wird, erhalten wird: q(k, t) = –ln[I(k, t)/I0(t)]
    • (7) die äquivalente Beziehung zwischen den Projektionsdaten q(k, t) des großen Fokus, die in Schritt (6) erhalten wurden, und den Projektionsdaten pciri, ti) der äquivalenten Sub-Fokusse mit den gleichen Geometrie-Scan-Parametern wird durch die folgende Formel dargestellt, wobei die gemessenen Daten von den verschiedenen Detektoreinheiten in verschiedenen Winkeln lineare Gleichungen bilden, wobei q(k, t) und w(ai) bekannt sind und pciri, ti) unbekannt ist; der Prozessor kann die Wiederherstellung der Projektionsdaten pciri, ti) der Sub-Fokusse durch Lösen dieser linearen Gleichungen erreichen:
      Figure DE102015226693A1_0013
      wobei R eine Distanz zwischen der Röntgenstrahlungsquelle und dem Drehmittelpunkt O des Drehtisches ist, D eine Distanz zwischen der Strahlenquelle und dem Detektor ist, ai eine Distanz von einer virtuellen Punktquelle Si ist, die von der idealen Punktquelle S abweicht, t die Position der Detektoreinheit darstellt und i ein tiefgestellter Index eines einzelnen virtuellen Sub-Fokus ist.
    • (8) der Prozessor rekonstruiert das Lineardämpfungskoeffizienten-Verteilungsbild des zu detektierenden Objekts unter Verwendung des gefilterten Rückprojektionsalgorithmus (FBP) auf der Basis der wiederhergestellten Projektionsdaten pciri, ti) der Sub-Fokusse, die in Schritt (7) erhalten wurde.
  • Die technischen Lösungen der vorliegenden Erfindung verbessern die räumliche Auflösung des CT-Systems ohne Verändern des existierenden CT-Systems und bewirken eine hohe Effizienz der Datenverarbeitung und Bildrekonstruktion.
  • Im Stand der Technik sind die Kanten des Bildes verschwommen, und die räumliche Auflösung ist schlecht, nur 2,4 LP/mm. Im Gegensatz dazu sind in den Ausführungsformen der vorliegenden Erfindung die Kanten Bildes scharf, und die räumliche Auflösung ist viel besser als 3,0 LP/mm.
  • Die obigen Beschreibungen sind lediglich beispielhafte Implementierungsformen der vorliegenden Erfindung. Der Schutzumfang der vorliegenden Erfindung ist nicht darauf beschränkt. Variationen oder Ersetzungen, die dem Fachmann sofort in den Sinn kommen, ohne vom technischen Geltungsbereich der vorliegenden Erfindung abzuweichen, sollen unter den Schutzumfang der vorliegenden Erfindung fallen. Darum unterliegt der Schutzumfang der vorliegenden Erfindung dem Schutzumfang der Ansprüche.
  • ZITATE ENTHALTEN IN DER BESCHREIBUNG
  • Diese Liste der vom Anmelder aufgeführten Dokumente wurde automatisiert erzeugt und ist ausschließlich zur besseren Information des Lesers aufgenommen. Die Liste ist nicht Bestandteil der deutschen Patent- bzw. Gebrauchsmusteranmeldung. Das DPMA übernimmt keinerlei Haftung für etwaige Fehler oder Auslassungen.
  • Zitierte Patentliteratur
    • CN 102525527 A [0004]

Claims (18)

  1. CT-Detektionsverfahren, das folgende Schritte umfasst: (1) Scannen in Umfangsrichtung eines zu detektierenden Objekts mittels Röntgenstrahlen gemäß einem voreingestellten Winkelabtastwert, der die Anzahl von Abtastpunkten in einem einzelnen Kreis darstellt, um eine Gruppe von Projektionsabtastdaten in verschiedenen Projektionswinkeln zu erhalten, wobei der voreingestellte Winkelabtastwert größer als 1000 ist; (2) Verarbeiten der Projektionsabtastdaten dergestalt, dass Projektionsdaten mehrerer virtueller Sub-Fokusse erhalten werden, die einem großen Fokus einer Strahlungsquelle in einem CT-System entsprechen; und (3) Implementieren einer Bildrekonstruktion gemäß den Projektionsdaten der mehreren virtuellen Sub-Fokusse.
  2. Verfahren nach Anspruch 1, wobei der voreingestellte Winkelabtastwert größer ist als 8000.
  3. Verfahren nach Anspruch 1 oder 2, wobei der voreingestellte Winkelabtastwert N gemäß der folgenden Formel erhalten wird: N > 2π/tan–1(a/R), wobei a die Fokusgröße einer Röntgenstrahlungsquelle ist und R eine Distanz zwischen der Röntgenstrahlungsquelle und einem Drehmittelpunkt des zu detektierenden Objekts ist.
  4. Verfahren nach Anspruch 1, wobei in Schritt (2) die Projektionsdaten der mehreren virtuellen Sub-Fokusse gemäß den Projektionsabtastdaten in verschiedenen Projektionswinkeln und einer Intensität, die den mehreren virtuellen Sub-Fokussen entspricht, erhalten werden.
  5. Verfahren nach Anspruch 4, wobei die Projektionsdaten der mehreren virtuellen Sub-Fokusse unter Verwendung der folgenden Formel erhalten werden:
    Figure DE102015226693A1_0014
    wobei D eine Distanz zwischen der Röntgenstrahlungsquelle und einem Detektor ist, ai eine Distanz eines einzelnen virtuellen Sub-Fokus ist, der von einer idealen Punktquelle abweicht, M die Anzahl der virtuellen Sub-Fokusse ist, t die Position des Detektors darstellt, i ein tiefgestellter Index eines einzelnen virtuellen Sub-Fokus ist, βi die Projektionswinkelposition eines einzelnen virtuellen Sub-Fokus ist und q(k, t) ein Dämpfungskoeffizient des Röntgenstrahls ist, der durch die folgenden Formel erhalten wird: q(k, t) = –ln[I(k, t)/I0(t)].
  6. Verfahren nach Anspruch 1, wobei vor dem Schritt (1) der große Fokus des CT-Systems mittels eines Fokusgrößenmessverfahrens in mehrere virtuelle Sub-Fokusse fein eingeteilt wird und die Intensität, die den mehreren virtuellen Sub-Fokussen entspricht, erhalten wird.
  7. Verfahren nach Anspruch 6, wobei der Schritt zum Erhalten der Intensität, die den mehreren virtuellen Sub-Fokussen entspricht, Folgendes umfasst: Messen der Größe des großen Fokus des CT-Systems mittels eines Kleinloch-Bildgabeverfahrens; Feines Einteilen des großen Fokus in die mehreren äquivalenten virtuellen Sub-Fokusse; und Bestimmen der Intensität, die den mehreren äquivalenten virtuellen Sub-Fokussen entspricht, gemäß einer Verteilungskurve der Fokusintensität.
  8. Verfahren nach Anspruch 1 oder 5, das des Weiteren folgenden Schritt umfasst: Korrigieren der Projektionsabtastdaten vor dem Schritt (2).
  9. Verfahren nach Anspruch 8, wobei der Schritt zum Korrigieren der Projektionsabtastdaten Folgendes umfasst: Messen der Intensität I0(t) des durch die Luft gedämpften Strahls; Berechnen des Dämpfungskoeffizienten q(k, t) des Röntgenstrahls in verschiedenen Winkeln mittels der folgenden Formel: q(k, t) = –ln[I(k, t)/I0(t)], wobei k ein Winkelabtastzählwert ist, der von einem anfänglichen Wert 1 aus größer wird, bis der Wert k größer ist als ein voreingestellter Winkelabtastzählwert.
  10. CT-Vorrichtung, die Folgendes umfasst: eine Strahlungsquelle, einen Detektor, der der Strahlungsquelle entspricht, eine Voreinstellwinkelabtastwert-Bestimmungseinheit, eine Projektionsabtastdaten-Verarbeitungseinheit und eine Rekonstruktionseinheit, wobei die Voreinstellwinkelabtastwert-Bestimmungseinheit dafür vorgesehen ist, einen voreingestellten Winkelabtastwert zu bestimmen, der die Anzahl von Abtastpunkten in einem einzelnen Kreis darstellt, wobei der voreingestellte Winkelabtastwert größer als 1000 ist; die Strahlungsquelle und der Detektor auf jeder Seite eines zu detektierenden Objekts angeordnet sind und dazu dienen, das zu detektierende Objekt mittels Röntgenstrahlen gemäß dem voreingestellten Winkelabtastwert in Umfangsrichtung zu scannen, um eine Gruppe von Projektionsabtastdaten in verschiedenen Projektionswinkeln zu erhalten; die Projektionsabtastdaten-Verarbeitungseinheit dazu ausgelegt ist, die Projektionsabtastdaten dergestalt zu verarbeiten, dass Projektionsdaten mehrerer virtueller Sub-Fokusse erhalten werden, die einem großen Fokus einer Strahlungsquelle in einem CT-System entsprechen; und die Rekonstruktionseinheit dazu ausgelegt ist, eine Bildrekonstruktion gemäß den Projektionsdaten der mehreren virtuellen Sub-Fokusse zu implementieren.
  11. CT-Vorrichtung nach Anspruch 10, wobei der voreingestellte Winkelabtastwert größer ist als 8000 und der voreingestellte Winkelabtastwert N gemäß der folgenden Formel erhalten wird: N > 2π/tan–1(a/R), wobei a die Fokusgröße einer Röntgenstrahlungsquelle ist und R eine Distanz zwischen der Röntgenstrahlungsquelle und einem Drehmittelpunkt des zu detektierenden Objekts ist.
  12. CT-Vorrichtung nach Anspruch 10, wobei die Projektionsabtastdaten-Verarbeitungseinheit die Projektionsdaten der mehreren virtuellen Sub-Fokusse gemäß den Projektionsabtastdaten in verschiedenen Projektionswinkeln und einer Intensität, die den mehreren virtuellen Sub-Fokussen entspricht, erhält, und die Projektionsdaten der mehreren virtuellen Sub-Fokusse mittels der folgenden Formel erhalten werden:
    Figure DE102015226693A1_0015
    wobei D eine Distanz zwischen der Röntgenstrahlungsquelle und dem Detektor ist, ai eine Distanz eines einzelnen virtuellen Sub-Fokus ist, der von einer idealen Punktquelle abweicht, M die Anzahl der virtuellen Sub-Fokusse ist, t die Position des Detektors darstellt, i ein tiefgestellter Index eines einzelnen virtuellen Sub-Fokus ist, βi die Projektionswinkelposition eines einzelnen virtuellen Sub-Fokus ist und q(k, t) ein Dämpfungskoeffizient des Röntgenstrahls ist, der durch die folgenden Formel erhalten wird: q(k, t) = –ln[I(k, t)/I0(t)].
  13. CT-Vorrichtung nach Anspruch 10, die des Weiteren eine Fokusgrößenmessungseinheit, eine Fokusteilungseinheit und eine Intensitätsbestimmungseinheit umfasst, wobei die Fokusgrößenmessungseinheit dafür vorgesehen ist, die Größe des großen Fokus des CT-Systems mittels eines Kleinloch-Bildgabeverfahrens zu messen; die Fokusteilungseinheit dafür vorgesehen ist, den großen Fokus in die mehreren äquivalenten virtuellen Sub-Fokusse fein einzuteilen; und die Intensitätsbestimmungseinheit dafür vorgesehen ist, eine Intensität, die den mehreren äquivalenten virtuellen Sub-Fokussen entspricht, gemäß einer Verteilungskurve der Fokusintensität zu bestimmen.
  14. CT-Vorrichtung nach Anspruch 10, die des Weiteren eine Datenkorrektureinheit zum Korrigieren der Projektionsabtastdaten umfasst, wobei auf der Basis der gemessenen Intensität I0(t) des durch die Luft gedämpften Strahls Dämpfungskoeffizienten q(k, t) des Röntgenstrahls in verschiedenen Winkeln mittels der folgenden Formel berechnet werden: q(k, t) = –ln[I(k, t)/I0(t)], wobei k ein Winkelabtastzählwert ist, der von einem anfänglichen Wert 1 aus größer wird, bis der Wert k größer ist als ein voreingestellter Winkelabtastzählwert.
  15. Bildrekonstruktionsverfahren in einem Röntgen-CT-Bildgabesystem, das folgende Schritte umfasst: (1) Aktivieren eines Strahlgenerators und eines Detektors in dem Röntgen-CT-Bildgabesystem, Messen der Intensität I0(t) des durch die Luft gedämpften Strahls an einer Position t und Übermitteln der gemessenen Daten an einen Datenprozessor in dem Röntgen-CT-Bildgabesystem; (2) Anordnen eines zu detektierenden Objekts in dem Röntgen-CT-Bildgabesystem zum Scannen, Messen der Intensität I(k, t) des Strahls, der durch das zu detektierende Objekt gedämpft wird, in der momentanen Winkelposition an der Position t und Übermitteln der gemessenen Daten an den Datenprozessor in dem Röntgen-CT-Bildgabesystem, wobei k ein Winkelabtastzählwert ist, der von einem anfänglichen Wert 1 aus größer wird, bis der Wert k größer ist als ein voreingestellter Winkelabtastzählwert; (3) Berechnen von Dämpfungskoeffizienten q(k, t) des Röntgenstrahls in verschiedenen Winkeln unter Verwendung der folgenden Formel (1): q(k, t) = –ln[I(k, t)/I0(t)] (1) (4) Rekonstruieren eines Lineardämpfungskoeffizienten-Verteilungsbildes des zu detektierenden Objekts mittels eines gefilterten Rückprojektionsalgorithmus auf der Basis der Daten q(k, t) in Schritt (3).
  16. Bildrekonstruktionsverfahren nach Anspruch 15, wobei die folgenden Schritte vor dem Implementieren von Schritt (1) implementiert werden: Messen und Erhalten der Fokusgröße einer Röntgenstrahlungsquelle und Messen und Erhalten einer Distanz R zwischen der Röntgenstrahlungsquelle und einem Drehmittelpunkt eines Drehtisches in dem Röntgen-CT-Bildgabesystem und Bestimmen eines voreingestellten Winkelabtastzählwertes N gemäß der folgenden Formel, auf der Basis von a und R: N > 2π/tan–1(a/R).
  17. Bildrekonstruktionsverfahren nach Anspruch 16, das des Weiteren folgenden Schritt umfasst: Bestimmen, ob der Wert k den voreingestellten Winkelabtastzählwert N übersteigt, wenn Schritt (2) implementiert wird, dann: wenn k größer ist als N, Implementieren von Schritt (3); wenn k kleiner ist als N, Drehen des zu detektierenden Objekts zum nächsten Scan-Winkel und Fortsetzen der Implementierung von Schritt (2).
  18. Bildrekonstruktionsverfahren nach Anspruch 16, das des Weiteren folgende Schritte umfasst: Teilen der Röntgenstrahlungsquelle in mehrere virtuelle Sub-Fokusse und Bestimmen der Position ai und der Intensität w(ai) jedes virtuellen Sub-Fokus; nach der Implementierung von Schritt (3), Erhalten von Wiederherstellungsdaten pciri, ti) der virtuellen Sub-Fokusse zum Implementieren einer Bildrekonstruktion in dem gefilterten Rückprojektionsalgorithmus mittels der folgenden Formel:
    Figure DE102015226693A1_0016
    Figure DE102015226693A1_0017
    wobei D eine Distanz zwischen der Röntgenstrahlungsquelle und dem Detektor ist, ai eine Distanz eines einzelnen virtuellen Sub-Fokus ist, der von einer idealen Punktquelle s abweicht, M die Anzahl der virtuellen Sub-Fokusse ist, t die Position des Detektors darstellt, i ein tiefgestellter Index eines einzelnen virtuellen Sub-Fokus ist, βi die Projektionswinkelposition eines einzelnen virtuellen Sub-Fokus ist; und Rekonstruieren eines Lineardämpfungskoeffizienten-Verteilungsbildes des zu detektierenden Objekts mittels eines gefilterten Rückprojektionsalgorithmus auf der Basis der Wiederherstellungsdaten pciri, ti).
DE102015226693.1A 2014-12-31 2015-12-23 CT-Detektionsverfahren und CT-Vorrichtung Ceased DE102015226693A1 (de)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201410851690.1A CN105806858B (zh) 2014-12-31 2014-12-31 Ct检测方法和ct设备
CN201410851690.1 2014-12-31

Publications (1)

Publication Number Publication Date
DE102015226693A1 true DE102015226693A1 (de) 2016-06-30

Family

ID=56116856

Family Applications (1)

Application Number Title Priority Date Filing Date
DE102015226693.1A Ceased DE102015226693A1 (de) 2014-12-31 2015-12-23 CT-Detektionsverfahren und CT-Vorrichtung

Country Status (5)

Country Link
US (1) US10175182B2 (de)
JP (1) JP6120940B2 (de)
CN (1) CN105806858B (de)
DE (1) DE102015226693A1 (de)
RU (1) RU2619839C1 (de)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110264541A (zh) * 2019-07-12 2019-09-20 四川明峰医疗科技有限公司 一种z向飞焦点扫描方式和图像重建方法
CN110389139A (zh) * 2018-04-20 2019-10-29 Fei 公司 用于感兴趣区域的断层摄影的扫描轨迹
DE102018132025A1 (de) * 2018-12-13 2020-06-18 Intom GmbH Industrielles Röntgen-Werkstückvermessungssystem und Verfahren zum Betreiben eines solchen

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10031091B2 (en) * 2016-07-27 2018-07-24 Avonix Imaging, LLC Shifting turn table for x-ray imaging manipulator
CN108811488B (zh) * 2017-02-27 2022-05-03 纳米维景(上海)医疗科技有限公司 适应大视野要求的静态实时ct成像系统及其成像方法
CN108195854A (zh) * 2017-12-26 2018-06-22 中国计量科学研究院 一种x射线空气衰减系数检测方法
CN108226195B (zh) * 2017-12-28 2023-10-13 清华大学 Ct检查系统和ct成像方法
JP7143567B2 (ja) * 2018-09-14 2022-09-29 株式会社島津テクノリサーチ 材料試験機および放射線ct装置
CN109685867B (zh) * 2018-12-29 2023-05-30 上海联影医疗科技股份有限公司 一种ct系统、图像重建方法、装置、设备及介质
CN111265231B (zh) * 2019-04-15 2021-08-31 清华大学 分布式光源ct图像重建方法与系统
EP3764325B1 (de) * 2019-07-12 2024-05-15 Canon Medical Systems Corporation System und programm zur vermeidung von brennfleck-unschärfe
CN110766686A (zh) * 2019-10-31 2020-02-07 上海联影医疗科技有限公司 Ct投影数据处理方法、系统、可读存储介质和设备
CN111157558B (zh) * 2020-03-03 2022-08-09 国网安徽省电力有限公司电力科学研究院 一种在役盆式绝缘子的x射线数字成像检测方法
CN115105108B (zh) * 2022-06-30 2023-06-09 赛诺威盛科技(北京)股份有限公司 散焦校正方法、装置、散焦校正模体及电子装置
CN115511831B (zh) * 2022-09-27 2023-04-04 佳木斯大学 一种组织胚胎病理切片的数据分析处理系统及方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102525527A (zh) 2010-10-13 2012-07-04 株式会社东芝 Ct成像的投影数据加权方法

Family Cites Families (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4472822A (en) * 1980-05-19 1984-09-18 American Science And Engineering, Inc. X-Ray computed tomography using flying spot mechanical scanning mechanism
JP2002523740A (ja) * 1998-08-18 2002-07-30 ロッキード、マーティン、コーパレイシャン 放射線透過式デジタル型溶接部検査システム
US6504892B1 (en) * 2000-10-13 2003-01-07 University Of Rochester System and method for cone beam volume computed tomography using circle-plus-multiple-arc orbit
US6477221B1 (en) * 2001-02-16 2002-11-05 University Of Rochester System and method for fast parallel cone-beam reconstruction using one or more microprocessors
JP3963086B2 (ja) * 2001-07-13 2007-08-22 株式会社島津製作所 断層再構成ソフトウエアとその記録媒体と断層撮影装置
DE10222702A1 (de) * 2002-05-22 2003-12-24 Siemens Ag Verfahren zur Korrektur der Extrafokalstrahlung einer Röntgenröhre in der Computertomographie
JP4163991B2 (ja) * 2003-04-30 2008-10-08 株式会社モリタ製作所 X線ct撮影装置及び撮影方法
JP2005040328A (ja) * 2003-07-22 2005-02-17 Canon Inc 放射線画像処理装置
JP2006034800A (ja) * 2004-07-29 2006-02-09 Toyota Motor Corp X線源焦点サイズの測定装置および方法
US20110188626A1 (en) * 2005-05-24 2011-08-04 National University Corporation Hokkaido University Damage Evaluation Apparatus, and Damage Evaluation Method
CN101071111B (zh) * 2006-05-08 2011-05-11 清华大学 一种多视角航空集装箱安全检查系统及方法
CN101382508A (zh) * 2007-09-05 2009-03-11 同方威视技术股份有限公司 一种检查航空货运集装箱中违禁物品的装置和方法
CN101396270B (zh) * 2007-09-29 2010-09-01 上海西门子医疗器械有限公司 一种降低检查床对ct图像影响的方法
GB0917950D0 (en) * 2009-10-13 2009-11-25 Shawcor Ltd X-ray inspection method and apparatus for pipeline girth weld inspection
KR20110135136A (ko) * 2010-06-10 2011-12-16 주식회사 하이닉스반도체 반도체 장치의 극미세 패턴 형성을 위한 방법
US8842804B2 (en) * 2011-01-12 2014-09-23 Kabushiki Kaisha Toshiba X-ray computed tomography apparatus
JP2013081527A (ja) * 2011-10-06 2013-05-09 Toshiba Corp X線診断装置
JP6009178B2 (ja) * 2012-03-01 2016-10-19 株式会社神戸製鋼所 アルミニウム材内の介在物可視化方法
US9414790B2 (en) * 2012-10-15 2016-08-16 Elekta Ab (Publ) Dual-energy cone-beam CT scanning
CN103913779B (zh) * 2012-12-31 2017-07-07 清华大学 多能ct成像系统以及成像方法
CN104068885B (zh) * 2013-03-28 2016-07-06 上海西门子医疗器械有限公司 Ct机的扫描触发控制方法、装置及ct机
CN103366389A (zh) * 2013-04-27 2013-10-23 中国人民解放军北京军区总医院 Ct图像重建方法
CN104142353B (zh) * 2014-06-06 2017-01-11 清华大学 可见光成像辅助定位局部ct扫描方法及系统

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102525527A (zh) 2010-10-13 2012-07-04 株式会社东芝 Ct成像的投影数据加权方法

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110389139A (zh) * 2018-04-20 2019-10-29 Fei 公司 用于感兴趣区域的断层摄影的扫描轨迹
DE102018132025A1 (de) * 2018-12-13 2020-06-18 Intom GmbH Industrielles Röntgen-Werkstückvermessungssystem und Verfahren zum Betreiben eines solchen
CN110264541A (zh) * 2019-07-12 2019-09-20 四川明峰医疗科技有限公司 一种z向飞焦点扫描方式和图像重建方法
CN110264541B (zh) * 2019-07-12 2020-10-02 明峰医疗系统股份有限公司 一种z向飞焦点扫描方式和图像重建方法

Also Published As

Publication number Publication date
US20160189403A1 (en) 2016-06-30
JP6120940B2 (ja) 2017-04-26
JP2016126010A (ja) 2016-07-11
CN105806858A (zh) 2016-07-27
CN105806858B (zh) 2019-05-17
RU2619839C1 (ru) 2017-05-18
US10175182B2 (en) 2019-01-08

Similar Documents

Publication Publication Date Title
DE102015226693A1 (de) CT-Detektionsverfahren und CT-Vorrichtung
DE60026681T2 (de) Verfahren und Gerät zur Mehrschichttomografie mit partieller Abtastung
EP2268204B1 (de) Vorrichtung und computer implementiertes verfahren zur rotationsfreien computertomographie
DE102013217852B3 (de) Verfahren und Röntgensystem zur Zwei-Energiespektren-CT-Abtastung und Bildrekonstruktion
DE102009028104A1 (de) Bildrekonstruktionsverfahren für energiereiches Doppelenergie-CT-System
DE102005009263A1 (de) Verfahren und System zur Bildgebung unter Verwendung mehrerer versetzter Röntgenstrahlen-Emissionspunkte
DE102010060989A1 (de) System und Verfahren zur Kompensation von Daten schwacher Signale für Dualenergie-CT
DE112006003039T5 (de) Verfahren und Anordnung zur Röntgenbildgebung
DE10252662A1 (de) Computertomographie-Verfahren mit kohärenten Streustrahlen und Computertomograph
DE102011056347A1 (de) Integrierte Röntgendetektoranordnung und Verfahren zur Herstellung derselben
DE112007001159T5 (de) Verfahren zur Frachtsicherheitskontrolle mit mehreren Betrachtungswinkeln
DE102013205406A1 (de) Röntgenaufnahmesystem zur Röntgenbildgebung bei hohen Bildfrequenzen eines Untersuchungsobjekts mittels direkter Messung des Interferenzmusters
DE10309887A1 (de) Kegelstrahl-Computertomographiesystem
DE102012104599A1 (de) Verfahren und System zur Rekonstruktion tomografischer Bilder
DE102012203086A1 (de) Bildverarbeitungsvorrichtung eines computertomographie-systems
Pelt et al. Ring artifact reduction in synchrotron x-ray tomography through helical acquisition
WO2004070661A1 (de) Voxel-getriebene spiralrekonstruktion für kegelstrahl-ct
DE102012217163A1 (de) Verfahren zur Rekonstruktion von CT-Bilddaten mit gewichteter Rückprojektion, einschließlich Recheneinheit und CT-System für dieses Verfahren
DE102007024409A1 (de) Verfahren und Röntgen-CT-System zur Erzeugung computertomographischer Darstellungen
EP1779327A1 (de) 3d-rekonstruktion mit schräger geometrie
DE10337935A1 (de) Vorrichtung für die Aufnahme von Strukturdaten eines Objekts
DE102013200329B4 (de) Verfahren und Vorrichtung zur Dejustagekorrektur für Bildgebungsverfahren
DE60214022T2 (de) Verfahren zur verringerung von artefakten in objektbildern
EP2644095B1 (de) Computertomographiesystem und Verfahren zur Datenermittlung für eine Störeinfluss-korrigierte Computertomographieaufnahme eines Untersuchungsobjekts
DE102019202452A1 (de) Röntgenuntersuchungsanordnung und Verfahren zum Betreiben einer Röntgenuntersuchungsanordnung

Legal Events

Date Code Title Description
R012 Request for examination validly filed
R002 Refusal decision in examination/registration proceedings
R003 Refusal decision now final