EP2080015A1 - Calibration method for two-spectra or multispectra tomography - Google Patents

Calibration method for two-spectra or multispectra tomography

Info

Publication number
EP2080015A1
EP2080015A1 EP07818508A EP07818508A EP2080015A1 EP 2080015 A1 EP2080015 A1 EP 2080015A1 EP 07818508 A EP07818508 A EP 07818508A EP 07818508 A EP07818508 A EP 07818508A EP 2080015 A1 EP2080015 A1 EP 2080015A1
Authority
EP
European Patent Office
Prior art keywords
calibration
spectra
phantom
decomposition
coefficients
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.)
Withdrawn
Application number
EP07818508A
Other languages
German (de)
French (fr)
Inventor
Marc Kachelriess
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.)
CT Imaging GmbH
Original Assignee
VAMP VERFAHREN und APP DER MED
VAMP Verfahren und Apparate der Medizinischen Physik GmbH
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 VAMP VERFAHREN und APP DER MED, VAMP Verfahren und Apparate der Medizinischen Physik GmbH filed Critical VAMP VERFAHREN und APP DER MED
Publication of EP2080015A1 publication Critical patent/EP2080015A1/en
Withdrawn legal-status Critical Current

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/48—Diagnostic techniques
    • A61B6/482—Diagnostic techniques involving multiple energy imaging
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computerised tomographs
    • A61B6/032—Transmission computed tomography [CT]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/58—Testing, adjusting or calibrating apparatus or devices for radiation diagnosis
    • A61B6/582—Calibration
    • A61B6/583—Calibration using calibration phantoms
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/003—Reconstruction from projections, e.g. tomography
    • G06T11/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/40—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for generating radiation specially adapted for radiation diagnosis
    • A61B6/405—Source units specially adapted to modify characteristics of the beam during the data acquisition process
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/30—Accessories, mechanical or electrical features
    • G01N2223/303—Accessories, mechanical or electrical features calibrating, standardising
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00—Image generation
    • G06T2211/40—Computed tomography
    • G06T2211/408—Dual energy

Definitions

  • the invention relates to a method for calibrating computer tomographs in two- or multi-spectral tomography.
  • X-ray computer tomography provides information about the density or attenuation coefficient distribution of an object. Due to the fact that it is measured with polychromatic spectra, this information is a kind of average over all energy values.
  • Dual or multi-spectral tomography measures the object with two or more different spectra (in the diagnostic energy range below 1 MeV, two scans with two different spectra are sufficient for the base material decomposition and in the following description the representations are reduced to two-spectral methods without restriction of generality) and so on create a base function decomposition of the attenuation coefficient distribution. Material-specific images are provided in this way. This basic function decomposition can be made from analytical assumptions about the spectra and the material weakening or from empirical measurements.
  • the linear attenuation coefficient can be considered as a decomposition of spatial and energy dependence
  • the spatial functions f t may represent the density distribution of material 1 and material 2 and the energy functions ⁇ , whose energy dependence may be in the form of the mass attenuation coefficient.
  • decomposition is the location function as a mass fraction (for example, in% of the base material /) of a material to understand. Frequently, the physical underlying cross sections for the photo and Compton effect are used as energy function so that the reconstructed images / ( directly show the cross section distribution.
  • the index L represents the integration line (ray). Since the corrections discussed below are independent of L, this index is not explicitly noted for the sake of simplicity.
  • the spectrum itself is called a normalized function W j (E) with
  • This energy weighting includes not only the X-ray spectrum but also the detector sensitivity and may also include other effects ("Detected Spectrum").
  • Thicknesses various combinations of materials 1 and 2) and interpolated by analytical calculation of the "monochromatic weakenings", see A.J. Coleman and M. Sinclair, A Beam-Hardening Correction Using Dual-Energy Computed Tomography, Physics in Medicine and
  • Accuracies in the micrometer range can be achieved.
  • such empirical methods often use only a limited beam range (eg, the central beam) and then achieve relatively high calibration errors. After all, those on pure
  • Measurements based on scattering measurements do not detect and compensate for measurement errors based on scattered radiation.
  • calibration methods are also described in some cases in which an object of known shape and composition (often a cylindrical phantom) is scanned and the theoretical attenuation profile is adapted to the measured one by a nonlinear fit method.
  • the Positioning accuracy can be a problem in this case, because the profile functions to be apprehended, for example, can not respond to accidental tilting or rotation of the phantom.
  • scattering effects and the non-linearities of the polychromatic attenuation measurement will distort the fit.
  • An object of the present invention is to obviate the aforementioned errors and to provide a simple calibration method that directly incorporates the material decomposition equations calculated.
  • a calibration phantom consisting of homogeneous regions of both base materials (one base material per spectrum) and possibly also of foreign materials is scanned (CT scan).
  • CT scan a calibration phantom consisting of homogeneous regions of both base materials (one base material per spectrum) and possibly also of foreign materials is scanned (CT scan).
  • the raw data is transformed and reconstructed by suitable basic functions. It is only in the image space that the images thus reconstructed are combined to form an object by means of a linear combination, in such a way that the two homogeneous regions are also homogeneous in the image.
  • the coefficients of the linear combination thus determined are then the decomposition coefficients to be used for later measurements.
  • An essential idea of the invention is not to perform the calibration in the raw data space (attenuation space). Instead, images are first created. The calibration (determination the coefficients) then takes place with the generated image data in the spatial domain. An exact knowledge of the dimensions and the position of the phantom is therefore no longer necessary.
  • the phantom consists of at least two materials. Essential to the invention is that the measurements are carried out with at least two voltages. In the case described above, there is a "dual energy" calibration, in other words a mapping from R 2 to R 2 (two voltages, two materials) Two template images (one for each image) are used.
  • a calibration phantom is scanned.
  • the phantom should have at least one homogeneous (ie easily segmentable) area per base material. The phantom is allowed too
  • the templates for material 1 and material 2 are then given by 1 if r in material 1 0 if r in material 2 or air? if r in undef. Ranges 1 if r in material 2 0 if r in material 1 or air if r in undef. areas
  • a weighting image is additionally generated. Its pixel values reflect the security about the pixel content of the associated CT image of the calibration phantom:
  • edges of the respective associated template image may be weighted away by the weighting image
  • the method described can also be combined with analytical methods. Often the case arises that the correction function depends on additional parameters such as the location of the detector pixel which is about to be corrected. Thus, a host of calibration functions is to be determined.
  • the mentioned location dependency can have many causes:
  • the rays emanating from an X-ray tube will have a slightly different spectrum depending on the angle of emission.
  • a hybrid method is advantageous: an analytical precorrection and first material decomposition followed by the image-based calibration method described above.
  • the invention is not limited to two-spectra tomography. A generalization to any higher dimensions (maps from R N to R N ) is possible. Moreover, the invention is not limited to any particular type of computed tomography.
  • Density e.g., teflon, bone
  • FIGS. 6-9 second images of the calibration phantom (second CT
  • Density e.g., water
  • black outside means air. Inside the phantom means one dark color is a low density material (eg water); a bright color means a high density material (eg Teflon, bone). It can be seen that the two inner materials are not homogeneous. It gets lighter towards the edge (see “Water”, but also applies to “Bones”).
  • a low density material eg water
  • a bright color means a high density material (eg Teflon, bone).
  • the standard reconstruction (FIG. 1) is used for segmentation and for determining the weighting and template images. From this, the material decomposition images are reconstructed.
  • the weighted difference images (FIGS. 5, 9) show only minimal deviations from the template and impressively demonstrate the mode of operation of the method.
  • CT images are taken, once with 8OkV, once with 14OkV tube voltage.
  • Dual-energy CT is a modality that takes one and the same object with two different X-ray spectra. Normally, the data is generated by two different voltages, but other methods such as different pre-filtering, post-filtering or layered detectors are also used.
  • DECT is used for energy or material selective reconstruction, with the reduction of
  • Jet hardening artifacts is another useful effect, cf. R. Alvarez and A. Macovski, "Energy-selective reconstructions in x-ray CT,” Phys. Med. Biol. , vol. 21, no. 5, pp. 733-744, 1976; R. Alvarez and E. Seppi, "A Comparison of noise and dose in conventional and energy selective computed tomography," IEEE Transactions on Nuclear Science, vol. NS-26, no. 2, pp. 2853-2856, Apr. 1979; A. Coleman and M. Sinclair, "A beam-hardening correction using dual-energy computed tomography," Phys. Med. Biol. , vol. 30, no. 11, pp.
  • the two-spectra CT is based on the assumption that the attenuation coefficient ⁇ (F, E), which depends on the location F and the photon energy E, can be decomposed as follows:
  • SR denotes the operator of the 2D radon transformation or the operator of the 3D x-ray transformation.
  • Empirical Dual Energy Calibration proposed here presents a novel empirical calibration algorithm. Unlike other methods, EDEC requires neither knowledge of the spectrum nor the attenuation coefficients. The geometry, size and position of the calibration phantom are not needed. EDEC is based on similar principles as the Empirical Cupping Correction (ECC; see K. Sourbelle, M. Kachelr understand, and WA Calender, "Empirical Water Precorrection for Cone-Beam Computed Tomography", IEEE Medical Imaging Conference Record, pp. 431-145) 1871-1875, Oct. 2005, M. Kachelr understand, K. Sourbelle, and W.208, "Empirical cupping correction: A first order raw data precorrection for cone-beam computed tomography," Med. Phys., Vol. no. 5, pp. 1269-1274, May 2006).
  • ECC Empirical Cupping Correction
  • the total number of basis functions is (K + 1) (L + 1).
  • the total number of coefficients in each coefficient vector C 1 is thus 25 for each material, increasing the total number of unknown coefficients to 50.
  • the task of EDEC is to determine these coefficients C 1 . It should be noted that EDEC is not limited to polynomials, other basic functions may also be chosen.
  • the decomposition according to material 2 is completely analogous.
  • the weighting function w (r) is used to eliminate unwanted structures of the calibration phantom that occurs in the optimization process.
  • the calibration phantom contains homogeneous regions with sufficient quantities of material 1 and 2 and ensures that all meaningful combinations of path lengths in material 1 and 2 are detected. To obtain the basic images, these raw data are transferred to the (K + 1) (L + 1) basis functions and reconstructed. A standard reconstruction of the calibration phantom serves as the basis for determining t (r) and w (r), which are determined by thresholding.
  • the material template / (F) represents the a priori knowledge of the regions containing material 1 and those that certainly do not contain material 1 (ie material 2 or air): 1 for re material 1
  • the material template need not be defined, since these regions are suppressed by the weighting.
  • the image of the material template is set equal to one wherever the content of the voxel r is unique. This is true in regions of material 1 or 2 and in air.
  • w (r) is set equal to zero.
  • the weight and material template images are shown in binary and in black and white.
  • EDEC was also used to calibrate measured data from a clinical two-spectrum CT scanner and a micro CT scanner. The experiments with real data confirmed the results of the simulations.
  • image-based empirical two-spectra calibration is a simple, effective, and accurate method for calibrating two-spectrum CT.

Abstract

The invention relates to a method for calibrating computer tomographs in the case of two-spectra or multispectra tomography. The method is characterized by direct calculation of the material decomposition equations P <SUB>i</SUB>

Description

Beschreibungdescription
Kalibriermethode für Zwei- oder MehrspektrentomographieCalibration method for two- or multi-spectral tomography
Die Erfindung betrifft ein Verfahren zur Kalibrierung von Computertomographen bei Zwei- oder Mehrspektrentomographie.The invention relates to a method for calibrating computer tomographs in two- or multi-spectral tomography.
Einführungintroduction
Die Röntgencomputertomographie gibt Auskunft über die Dichtebzw, die Schwächungskoeffizientenverteilung eines Objekts. Aufgrund der Tatsache, daß mit polychromatischen Spektren gemessen wird, ist diese Information eine Art Mittel über alle Energiewerte. Die Zwei- oder Mehrspektrentomographie hingegen mißt das Objekt mit zwei oder mehreren unterschiedlichen Spektren (im diagnostischen Energiebereich unterhalb 1 MeV reichen zur Basismaterialzerlegungen zwei Scans mit zwei verschiedenen Spektren und in der nachfolgenden Beschreibung werden die Darstellungen ohne Beschränkung der Allgemeinheit auf Zweispektrenverfahren reduziert) und kann so eine Basisfunktionszerlegung der Schwächungskoeffizientenverteilung erstellen. Materialspezifische Bilder werden so bereitgestellt. Diese Basisfunktionszerlegung kann aus analytischen Annahmen über die Spektren und die Materialschwächung oder aus empirischen Messungen erfolgen. Analytische Verfahren haben den Nachteil, daß sie stark von der Genauigkeit der Annahmen über die Spektren und Materialien abhängen und werden deshalb nicht routinemäßig eingesetzt. Empirische Verfahren funktionieren derzeit so, daß unterschiedliche Absorber bekannter Dicke gescannt werden und die theoretischen Linienintegrale mit den tatsächlich gemessenen Werten per Fit in Zusammenhang gebracht werden. Nachteil hierbei ist die Abhängigkeit vonX-ray computer tomography provides information about the density or attenuation coefficient distribution of an object. Due to the fact that it is measured with polychromatic spectra, this information is a kind of average over all energy values. Dual or multi-spectral tomography, on the other hand, measures the object with two or more different spectra (in the diagnostic energy range below 1 MeV, two scans with two different spectra are sufficient for the base material decomposition and in the following description the representations are reduced to two-spectral methods without restriction of generality) and so on create a base function decomposition of the attenuation coefficient distribution. Material-specific images are provided in this way. This basic function decomposition can be made from analytical assumptions about the spectra and the material weakening or from empirical measurements. Analytical methods have the disadvantage that they are highly dependent on the accuracy of the assumptions about the spectra and materials and are therefore not routinely used. Empirical methods currently operate by scanning different absorbers of known thickness and associating the theoretical line integrals with the actual measured values by Fit. Disadvantage here is the dependence on
BESTÄTiOUNGSKOPIE der Geometrie der Absorber und von deren Positioniergenauigkeit .BESTÄTiOUNGSKOPIE the geometry of the absorber and its positioning accuracy.
Problemstellung, Zielsetzung, Stand der TechnikProblem definition, objective, state of the art
Der lineare Schwächungskoeffizient läßt sich als Zerlegung von Orts- und Energieabhängigkeit The linear attenuation coefficient can be considered as a decomposition of spatial and energy dependence
II
darstellen. Ohne Beschränkung der Allgemeinheit wird nachfolgend der in Energiebereichen unter 1 MeV übliche Fall zweier Basisfunktionenrepresent. Without limitation of generality, the case of two basic functions which is usual in energy ranges below 1 MeV will be hereafter
μiT,E) = fx(r)Ψ](E)+f2(r)ψ2(E)μiT, E) = f x (r) Ψ] (E) + f 2 (r) ψ 2 (E)
betrachtet, wobei die Ortsfunktionen ft beispielsweise die Dichteverteilung von Material 1 und Material 2 und die Energiefunktionen ψ, deren Energieabhängigkeit in Form des Massenschwächungskoeffizienten darstellen können. Eine andereFor example, the spatial functions f t may represent the density distribution of material 1 and material 2 and the energy functions ψ, whose energy dependence may be in the form of the mass attenuation coefficient. Another
Interpretationsmöglichkeit obiger Zerlegung ist die Ortsfunktion als Massenanteil (beispielsweise in % des Basismaterials / ) eines Materials aufzufassen. Häufig werden als Energiefunktion auch die physikalisch zugrunde liegenden Wirkungsquerschnitte für den Photo- und Comptoneffekt eingesetzt so daß die rekonstruierten Bilder /( direkt die Wirkungsquerschnittsverteilung zeigen.Interpretation option above decomposition is the location function as a mass fraction (for example, in% of the base material /) of a material to understand. Frequently, the physical underlying cross sections for the photo and Compton effect are used as energy function so that the reconstructed images / ( directly show the cross section distribution.
Um die Ortsfunktionen ft rekonstruieren zu können, werden Linienintegrale der Form In order to be able to reconstruct the spatial functions f t , line integrals of the form become
benötigt. Der Index L stellt dabei die Integrationslinie (Strahl) dar. Da die im Folgenden besprochenen Korrekturen unabhängig von L sind, wird dieser Index der Einfachheit halber nicht explizit notiert.needed. The index L represents the integration line (ray). Since the corrections discussed below are independent of L, this index is not explicitly noted for the sake of simplicity.
Die Messung liefert allerdings folgende Daten:However, the measurement provides the following data:
Der Index j zählt hierbei die unterschiedlichen, zur Messung verwendeten Spektren (die hier gemachten Ausführungen beschränken sich auf y=l, 2, die Methode ist jedoch für eine beliebige Anzahl von Spektren und Basisfunktionen anwendbar) . Das Spektrum selbst wird als normierte Funktion Wj(E) mitThe index j here counts the different spectra used for the measurement (the statements made here are limited to y = 1, 2, but the method is applicable to any number of spectra and basis functions). The spectrum itself is called a normalized function W j (E) with
Fläche 1 dargestellt. Diese Energiegewichtung umfaßt nicht nur das Röntgenspektrum, sondern auch noch die Detektorempfindlichkeit und kann auch andere Effekte mit beinhalten („Detected Spectrum") .Area 1 shown. This energy weighting includes not only the X-ray spectrum but also the detector sensitivity and may also include other effects ("Detected Spectrum").
Offensichtlich gilt q} = ^j(P],P2) und man erhält zweiObviously, q } = ^ j (P ] , P 2 ) and we obtain two
Gleichungen mit zwei Unbekannten, nach denen man gerne auflösen würden, um die gewünschte Materialzerlegung zu erhalten. Die Auflösung nach den Unbekannten P1= p,(q\,q2) setzt jedoch die genaue Kenntnis der W7(E) und der ψj{E) voraus. Diese Funktionen sind jedoch selten hinreichend genau bekannt. Zudem können Störeffekte wie beispielsweise Streustrahlung eine solch analytische Inversion fehlschlagen lassen. Stand der Technik sind deshalb so genannte empirische Kalibrierverfahren, bei denen keine oder nur wenige Annahmen gemacht werden müssen. Typischerweise werden die Funktionen Qj = <77(/7 I >/>2) an diskreten Stellen durch Messungen der Schwächung q} (A) verschiedener Absorber A (verschiedeneEquations with two unknowns, which one would like to dissolve to obtain the desired material decomposition. However, the resolution according to the unknowns P 1 = p, (q \ , q 2 ) requires the exact knowledge of W 7 (E) and ψ j {E). However, these functions are rarely known with sufficient accuracy. In addition, parasitics such as stray radiation can cause such an analytical inversion to fail. The prior art are therefore so-called empirical calibration methods in which no or only a few assumptions must be made. Typically, the functions Q j = <7 7 (/ 7 I> /> 2 ) at discrete locations by measurements of the attenuation q } (A) of different absorbers A (various
Dicken, verschiedene Kombinationen der Materialien 1 und 2) und durch analytische Berechnung der „monochromatischen Schwächungen" bestimmt und dazwischen interpoliert, vgl. A. J. Coleman and M. Sinclair, A Beam-Hardening Correction Using Dual-Energy Computed Tomography, Physics in Medicine andThicknesses, various combinations of materials 1 and 2) and interpolated by analytical calculation of the "monochromatic weakenings", see A.J. Coleman and M. Sinclair, A Beam-Hardening Correction Using Dual-Energy Computed Tomography, Physics in Medicine and
Biology, 30 (11) .1251-1256, 1985; sowie US 20050259781 und US 2004184574.Biology, 30 (11) .1251-1256, 1985; and US 20050259781 and US 2004184574.
Dieses Vorgehen hat allerdings einige Nachteile. Einerseits müssen die Abmessungen und die Materialzusammensetzung der Absorber hinreichend genau bekannt sein. Dann setzt das Verfahren eine hohe Positioniergenauigkeit voraus, um auch die Schnittlängen der betrachteten Strahlen mit den Absorbern korrekt zu berücksichtigen. Dies kann beispielsweise den Einsatz bei Mikro-CT hinderlich sein, denn dort müssenHowever, this approach has some disadvantages. On the one hand, the dimensions and the material composition of the absorber must be known with sufficient accuracy. Then, the method requires a high positioning accuracy in order to correctly take into account also the cutting lengths of the considered beams with the absorbers. This can be an obstacle to the use of micro-CT, for example, because there must
Genauigkeiten im Mikrometerbereich erzielt werden. Außerdem nutzen solche empirischen Verfahren oft nur einen eingeschränkten Strahlbereich (beispielsweise den Zentralstrahl) und erzielen dann relativ hohe Kalibrierfehler. Schließlich können die auf reinenAccuracies in the micrometer range can be achieved. In addition, such empirical methods often use only a limited beam range (eg, the central beam) and then achieve relatively high calibration errors. After all, those on pure
Schwächungsmessungen basierenden Verfahren auf Streustrahlung basierende Meßfehler nicht erkennen und ausgleichen. Als Alternative dazu sind teils auch Kalibrierverfahren beschrieben, bei denen ein Objekt bekannter Form und Zusammensetzung (oft ein zylindrisches Phantom) gescannt wird und per nichtlinearen Fitverfahren das theoretische Schwächungsprofil an das gemessene angepaßt wird. Hier treten ganz ähnliche Schwierigkeiten auf . Auch die Positioniergenauigkeit kann dabei ein Problem darstellen, denn die anzufittenden Profilfunktionen können beispielsweise nicht auf versehentliche Kippungen oder Drehungen des Phantoms reagieren. Zudem werden Streustrahleffekte und die Nichtlinearitäten der polychromatischen Schwächungsmessung den Fit verfälschen.Measurements based on scattering measurements do not detect and compensate for measurement errors based on scattered radiation. As an alternative to this, calibration methods are also described in some cases in which an object of known shape and composition (often a cylindrical phantom) is scanned and the theoretical attenuation profile is adapted to the measured one by a nonlinear fit method. Here are quite similar difficulties. Also the Positioning accuracy can be a problem in this case, because the profile functions to be apprehended, for example, can not respond to accidental tilting or rotation of the phantom. In addition, scattering effects and the non-linearities of the polychromatic attenuation measurement will distort the fit.
Eine Aufgabe der vorliegenden Erfindung ist es, die genannten Fehler zu umgehen und eine einfache Kalibriermethode zu schaffen, die direkt die Materialzerlegungsgleichungen errechnet.An object of the present invention is to obviate the aforementioned errors and to provide a simple calibration method that directly incorporates the material decomposition equations calculated.
Lösungsolution
Diese Aufgabe wird durch ein Verfahren nach Anspruch 1 gelöst. Vorteilhafte Ausführungen der Erfindung sind in den Unteransprüchen angegeben.This object is achieved by a method according to claim 1. Advantageous embodiments of the invention are specified in the subclaims.
Statt Absorbern wird ein Kalibrierphantom bestehend aus homogenen Bereichen beider Basismaterialien (je Spektrum ein Basismaterial) und eventuell auch aus Fremdmaterialien gescannt (CT-Scan) . Die Rohdaten werden durch geeignete Basisfunktionen transformiert und rekonstruiert. Erst im Bildraum werden die so rekonstruierten Bilder mittels Linearkombination zum Objekt zusammengesetzt und zwar so, daß die beiden homogenen Bereiche auch im Bild homogen sind. Die so bestimmten Koeffizienten der Linearkombination sind dann die für spätere Messungen anzuwendenden Zerlegungskoeffizienten.Instead of absorbers, a calibration phantom consisting of homogeneous regions of both base materials (one base material per spectrum) and possibly also of foreign materials is scanned (CT scan). The raw data is transformed and reconstructed by suitable basic functions. It is only in the image space that the images thus reconstructed are combined to form an object by means of a linear combination, in such a way that the two homogeneous regions are also homogeneous in the image. The coefficients of the linear combination thus determined are then the decomposition coefficients to be used for later measurements.
Eine Kernidee der Erfindung ist es, die Kalibrierung nicht im Rohdatenraum (Schwächungsraum) durchzuführen. Statt dessen werden zunächst Bilder erzeugt. Die Kalibrierung (Bestimmung der Koeffizienten) erfolgt dann mit den erzeugten Bilddaten im Ortsraum. Eine genaue Kenntnis der Abmessungen und der Lage des Phantoms ist somit nicht mehr erforderlich. Das Phantom besteht dabei aus wenigstens zwei Materialien. Erfindungswesentlich ist es, daß die Messungen mit mindestes zwei Spannungen durchgeführt werden. Im dem oben beschriebenen Fall erfolgt eine „Dual Energy" -Kalibrierung, mit anderen Worten eine Abbildung von R2 nach R2 (zwei Spannungen, zwei Materialien) . Es werden zwei Schablonen- Bilder (für jedes Bild eines) verwendet.An essential idea of the invention is not to perform the calibration in the raw data space (attenuation space). Instead, images are first created. The calibration (determination the coefficients) then takes place with the generated image data in the spatial domain. An exact knowledge of the dimensions and the position of the phantom is therefore no longer necessary. The phantom consists of at least two materials. Essential to the invention is that the measurements are carried out with at least two voltages. In the case described above, there is a "dual energy" calibration, in other words a mapping from R 2 to R 2 (two voltages, two materials) Two template images (one for each image) are used.
Die Idee besteht im einzelnen darin, die gesuchte Inversionsformel direkt zu bestimmen, und zwar durch eine Entwicklung nach Basisfunktionen Bn(qλ,q2) gemäßSpecifically, the idea is to find the inversion formula you are looking for directly determined by a development according to basis functions B n (q λ , q 2 ) according to
n n
Zu bestimmen sind lediglich die Koeffizienten cm . In dem hier beschriebenen Fall wird eine Taylorentwicklung verwendet, so daß die Basisfunktionen die Form mit k = 0,...,K und 1 = 0,...,L haben. Es kann jedoch gleichermaßen jede andere Basisfunktion verwendet werden, sofern die Ergebnisse dies zulassen.Only the coefficients c m are to be determined. In the case described here, Taylor development is used so that the basis functions are the form with k = 0, ..., K and 1 = 0, ..., L have. However, any other basic function can be used equally, as far as the results permit.
Die Bestimmung der Unbekannten erfolgt nun nicht im Rohdatenraum, sondern im Bildraum. Die Linearität derThe determination of the unknown does not take place in the raw data space, but in the image space. The linearity of
Radontransformierten erlaubt nämlich die Basisfunktionen zu rekonstruieren und erst die so entstandenen Basisbilder bn(r) linear zu kombinieren:For radon transformations allows to reconstruct the basis functions and to first linearly combine the resulting base images b n (r):
f,(r) = ∑cιnR-]Bn(gι,q2)= ∑clnbn(r) . Die inverse Radontransformation R ' stellt die Bildrekonstruktion dar und kann für jede beliebige CT- Geometrie mit der auch normale Rekonstruktionen stattfinden können durchgeführt werden.f, (r) = Σc ιn R- ] B n ( gι , q 2 ) = Σc ln b n (r). The inverse radon transformation R 'represents the image reconstruction and can be performed for any CT geometry with which normal reconstructions can also take place.
Um die Koeffizienten cm zu bestimmen, benötigt man noch ein so genanntes Schablonen- oder Templatebild. Dazu wird ein Kalibrierphantom gescannt. Das Phantom sollte je Basismaterial mindestens einen homogenen (d.h. leicht segmentierbaren) Bereich haben. Das Phantom darf auchIn order to determine the coefficients c m , one still needs a so-called stencil or template image. For this a calibration phantom is scanned. The phantom should have at least one homogeneous (ie easily segmentable) area per base material. The phantom is allowed too
Regionen haben, in denen unbekannte Materialien oder Gemische der Basismaterialien vorhanden sind. Ungewünschte Bereiche können durch ein noch zu definierendes Gewichtungsbild ignoriert werden. Die Weglängen der Strahlen durch die beiden homogenen Basismaterialbereiche sollten in etwa den typischen Verhältnissen der später zu scannenden Objekte entsprechen. Eine Standardrekonstruktion des Kalibrierphantoms ist zu erstellen und die Basismaterialbereiche sowie die umgebende Luft sind zu segmentieren (üblicherweise geht dies durch eine einfache schwellwertbasierte Segmentierung) .Have regions where unknown materials or mixtures of base materials are present. Unwanted areas can be ignored by a still-to-be-defined weighting image. The path lengths of the beams through the two homogeneous base material areas should correspond approximately to the typical ratios of the objects to be scanned later. A standard reconstruction of the calibration phantom must be made and the base material areas and surrounding air segmented (usually by simple threshold-based segmentation).
Die Schablonen für Material 1 und Material 2 sind dann gegeben durch 1 falls r in Material 1 0 falls r in Material 2 oder Luft ? falls r in undef. Bereichen 1 falls r in Material 2 0 falls r in Material 1 oder Luft falls r in undef. Bereichen Um Undefinierte oder unbekannte Materialien von der Koeffizientenbestimmung auszuschließen, wird zusätzlich ein Wichtungsbild erzeugt. Dessen Pixelwerte spiegeln die Sicherheit über den Pixelinhalt des zugehörigen CT-Bilds des Kalibrierphantoms wider:The templates for material 1 and material 2 are then given by 1 if r in material 1 0 if r in material 2 or air? if r in undef. Ranges 1 if r in material 2 0 if r in material 1 or air if r in undef. areas In order to exclude undefined or unknown materials from the coefficient determination, a weighting image is additionally generated. Its pixel values reflect the security about the pixel content of the associated CT image of the calibration phantom:
\\ falls r in Material 1 , 2 oder Luft W1 (r) = \\\ if r in material 1, 2 or air W 1 (r) = \
0 falls r in unbekanntem Material oder an Kanten des Schablonenbilds i0 if r in unknown material or on edges of the stencil image i
Die Kanten des jeweils zugehörigen Schablonenbilds können durch das Wichtungsbild weggewichtet werden, umThe edges of the respective associated template image may be weighted away by the weighting image
Punktbildfunktionseffekte aufgrund der begrenztenDot image function effects due to the limited
Ortsauflösung zu vermeiden. Falls nötig, können auch weitereTo avoid spatial resolution. If necessary, others can be added
Effekte, beispielsweise Gradientenverläufe, in die Schablone und die Wichtung mit eingebaut werden.Effects, such as Gradientverläufe, are incorporated in the template and the weighting with.
Die Bestimmung der Koeffizienten erfolgt nun durchThe determination of the coefficients now takes place
Minimierung von mit /(r) = ∑cΛ(r). n Das gesuchte Minimum wird durch Ableitung nach den C1n bestimmt . Dies führt zu dem Gleichungssystem U1 = B1 C1 mitMinimization of with / (r) = ΣcΛ (r). n The searched minimum is determined by derivation after the C 1n . This leads to the system of equations U 1 = B 1 C 1
mit dem Lösungsvektor C1 =B~] a, dessen Einträge die gesuchten Koeffizienten cm sind. Dieses lineare Gleichungssystem ist für jedes Material / zu lösen, im Falle der Dual Energy CT also zwei Mal. Nachdem die Koeffizienten mit Hilfe des Kalibrierscans bestimmt sind, können Folgescans Rohdatenbasiert ohne weiteren Aufwand in ihre Materialanteile zerlegt werden.with the solution vector C 1 = B ~] a, whose entries are the searched coefficients c m . This linear system of equations can be solved for each material, ie twice in the case of Dual Energy CT. After the coefficients with the help of Calibration scans are determined, subsequent scans raw data based can be decomposed into their material shares without further effort.
Die beschriebene Methode kann auch mit analytischen Verfahren kombiniert werden. Oft tritt der Fall auf, daß die Korrekturfunktion von zusätzlichen Parametern wie beispielsweise dem Ort des Detektorpixels, der gerade korrigiert werden soll, abhängt. Somit ist eine Schar von Kalibrierfunktionen zu bestimmen. Die angesprochene Ortsabhängigkeit kann vielerlei Ursachen haben:The method described can also be combined with analytical methods. Often the case arises that the correction function depends on additional parameters such as the location of the detector pixel which is about to be corrected. Thus, a host of calibration functions is to be determined. The mentioned location dependency can have many causes:
Zum Beispiel werden die von einer Röntgenröhre ausgehenden Strahlen je nach Abstrahlwinkel ein leicht unterschiedliches Spektrum aufweisen.For example, the rays emanating from an X-ray tube will have a slightly different spectrum depending on the angle of emission.
Oft setzt man zusätzliche Formfilter ein, die unterschiedliche Strahlen unterschiedlich vorfiltern. - Der unvermeidbare Streustrahlhintergrund ist ortsabhängig und verfälscht somit die Meßwerte, wobei die Größe des Fehlers dann ebenfalls vom Ort abhängt.Often one uses additional form filters which prefilter different rays differently. - The unavoidable scattered background is location-dependent and thus falsifies the measured values, the size of the error then also depends on the location.
In solchen Fällen ist also ein hybrides Verfahren von Vorteil: eine analytische Vorkorrektur und erste Materialzerlegung gefolgt von dem oben beschriebenen bildbasierten Kalibrierverfahren.In such cases, therefore, a hybrid method is advantageous: an analytical precorrection and first material decomposition followed by the image-based calibration method described above.
Unter Verwendung des erfinderischen Grundgedankens und erläuterten Ausführungsformen der Erfindung ergeben sich für einen Fachmann eine Vielzahl weiterer Ausführungsbeispiele, die jedoch an dieser Stelle nicht im Einzelnen beschrieben werden können. Insbesondere ist die Erfindung nicht auf Zweispektren- Tomographie beschränkt. Eine Verallgemeinerung auf beliebige höhere Dimensionen (Abbildungen von RN nach RN) ist möglich. Darüber hinaus ist die Erfindung nicht auf eine bestimmte Art von Computertomographie beschränkt.Using the inventive concept and illustrated embodiments of the invention will be apparent to a person skilled in a variety of other embodiments, which, however, can not be described in detail here. In particular, the invention is not limited to two-spectra tomography. A generalization to any higher dimensions (maps from R N to R N ) is possible. Moreover, the invention is not limited to any particular type of computed tomography.
In diesem Zusammenhang sei darauf hingewiesen, daß alle in der Beschreibung, den Ansprüchen und den Zeichnungen dargestellten Merkmale sowohl einzeln als auch in beliebiger Kombination miteinander erfindungswesentlich sein können.In this context, it should be noted that all features shown in the description, the claims and the drawings can be essential to the invention both individually and in any combination.
Ein Ausführungsbeispiel der Erfindung wird nachfolgend anhand der Zeichnungen näher erläutert. Hierbei zeigen:An embodiment of the invention will be explained in more detail with reference to the drawings. Hereby show:
Fig. 1 eine Standardrekonstruktion,1 a standard reconstruction,
Fig. 2-5 erste Bilder des Kalibrierphantoms (erstes CT-2-5 first images of the calibration phantom (first CT
Spektrum) ; Rekonstruktion für Material mit hoherSpectrum); Reconstruction for material with high
Dichte (z.B. Teflon, Knochen),Density (e.g., teflon, bone),
Fig. 6-9 zweite Bilder des Kalibrierphantoms (zweites CT-FIGS. 6-9 second images of the calibration phantom (second CT
Spektrum) ; Rekonstruktion für Material mit geringerSpectrum); Reconstruction for material with low
Dichte (z.B. Wasser),Density (e.g., water),
Fig. 10 Bilder des Kalibrierphantoms aus Wasser (i=l) und Teflon (i=2) .Fig. 10 Images of the calibration phantom of water (i = l) and Teflon (i = 2).
Eine Beispielapplikation ist in den Abbildungen anhand der Rekonstruktion eines bei 80 kV und bei 140 kV simulierten Yin-Yang-Kalibrierphantoms zu sehen.An example application can be seen in the figures from the reconstruction of a simulated Yin-Yang calibration phantom at 80 kV and 140 kV.
In der Standardrekonstruktion (Fig. 1) bedeutet eine schwarze Färbung (außen) Luft. Im Inneren des Phantoms bedeutet eine dunkle Farbe ein Material mit geringer Dichte (z.B. Wasser); eine helle Farbe bedeutet ein Material mit hoher Dichte (z.B. Teflon, Knochen) . Man erkennt, daß die beiden inneren Materialien nicht homogen sind. Zum Rand hin wird es heller (siehe „Wasser"; gilt aber auch für „Knochen").In the standard reconstruction (Figure 1), black (outside) means air. Inside the phantom means one dark color is a low density material (eg water); a bright color means a high density material (eg Teflon, bone). It can be seen that the two inner materials are not homogeneous. It gets lighter towards the edge (see "Water", but also applies to "Bones").
Die Standardrekonstruktion (Fig. 1) dient zur Segmentierung und zur Bestimmung der Wichtungs- und Templatebilder. Daraus werden die Materialzerlegungsbilder rekonstruiert. Die gewichteten Differenzbilder (Fig. 5, Fig. 9) zeigen nur geringste Abweichungen vom Template und demonstrieren die Funktionsweise der Methode eindrucksvoll.The standard reconstruction (FIG. 1) is used for segmentation and for determining the weighting and template images. From this, the material decomposition images are reconstructed. The weighted difference images (FIGS. 5, 9) show only minimal deviations from the template and impressively demonstrate the mode of operation of the method.
Mit Hilfe der erfindungsgemäßen „Dual-Energy" -Kalibrierung kann also erreicht werden, daß dem Bereich „Wasser", demWith the help of the "dual-energy" calibration according to the invention can thus be achieved that the area "water", the
Bereich „Knochen" und dem Bereich „Luft" jeweils nur ein einziger Graustufenwert zugeordnet ist. In dem beschriebenenOnly one greyscale value is assigned to the "Bones" area and the "Air" area. In the described
Beispiel werden zwei CT-Bilder aufgenommen, einmal mit 8OkV, einmal mit 14OkV Röhrenspannung.For example, two CT images are taken, once with 8OkV, once with 14OkV tube voltage.
Nachfolgend wird eine Ausführung der Erfindung im Detail beschrieben:Hereinafter, an embodiment of the invention will be described in detail:
Für die Qualität materialselektiver Bildgebung mit Zweispektren-CT (Dual Energy CT, DECT) ist eine genaueFor the quality of material-selective imaging with dual-spectral CT (Dual Energy CT, DECT) is an accurate
Kalibrierung der Zerlegungsfunktionen erforderlich. Dafür ist es entscheidend, das detektierte Röntgenspektrum genau zu kennen. Selbst wenn das Spektrum bekannt ist wird die Zuverlässigkeit von DECT durch Streustrahlung geschmälert. Zur Bestimmung einer korrekten Zerlegungsfunktion wird eine empirische Kalibriermethode vorgeschlagen. Im Gegensatz zu anderen Zerlegungsalgorithmen, ist für die hier vorgeschlagene empirische Kalibrierung bei Zweispektren-CT (Empirical Dual Energy Calibration, EDEC) weder die Kenntnis des Spektrums noch der Schwächungskoeffizienten notwendig. Um die gesuchten materialselektiven Rohdaten P1 und p2 als Funktion der gemessenen Schwächungswerte qx und q2 zu erhalten, wird eine polynomiale Abhängigkeit zugrunde gelegt. Die Bestimmung der Koeffizienten des Polynoms geschieht mit der Methode der kleinsten Quadrate, die aufCalibration of the decomposition functions required. For this it is crucial to know the detected X-ray spectrum exactly. Even if the spectrum is known, the reliability of DECT is diminished by stray radiation. To determine a correct decomposition function, an empirical calibration method is proposed. In contrast to other decomposition algorithms, for the empirical calibration proposed here is for two-spectra CT (Empirical Dual Energy Calibration, EDEC) neither the knowledge of the spectrum nor the attenuation coefficients necessary. In order to obtain the desired material-selective raw data P 1 and p 2 as a function of the measured attenuation values q x and q 2 , a polynomial dependence is used. The determination of the coefficients of the polynomial is done with the method of least squares, which
Schwellwertbildern eines Kalibrierphantoms basiert. Es werden keine Voraussetzungen an die genaue Größe oder Position des Phantoms gestellt. Sobald die Zerlegungskoeffizienten bestimmt sind, können die DECT-Rohdaten zerlegt werden, indem sie in das Polynom eingesetzt werden.Threshold images of a calibration phantom based. There are no requirements for the exact size or position of the phantom. Once the decomposition coefficients are determined, the DECT raw data can be decomposed by inserting it into the polynomial.
Um EDEC zu demonstrieren wurde ein Phantom aus Wasser und Teflon simuliert und ein physikalisches Phantom aus Wasser und Hydroxylapatit mit einem In-Vivo Mikro-CT-Scanner aufgenommen. Die Phantome wurden in die zu Grunde liegenden Materialien zerlegt.To demonstrate EDEC, a phantom of water and Teflon was simulated, and a physical phantom of water and hydroxylapatite was recorded with an in-vivo micro-CT scanner. The phantoms were decomposed into the underlying materials.
Zweispektren-CT (Dual Energy CT, DECT) ist eine Modalität, bei der ein und dasselbe Objekt mit zwei verschiedenen Röntgenspektren aufgenommen wird. Normalerweise werden die Daten durch zwei verschiedene Spannungen generiert, es kommen aber auch andere Methoden wie unterschiedliche Vorfilterung, Nachfilterung oder geschichtete Detektoren zur Anwendung.Dual-energy CT (DECT) is a modality that takes one and the same object with two different X-ray spectra. Normally, the data is generated by two different voltages, but other methods such as different pre-filtering, post-filtering or layered detectors are also used.
DECT wird verwendet für die energie- oder materialselektive Rekonstruktion, wobei die Reduzierung vonDECT is used for energy or material selective reconstruction, with the reduction of
Strahlaufhärtungsartefakten ein weiterer nützlicher Effekt ist, vgl. R. Alvarez and A. Macovski, "Energy-selective reconstructions in x-ray CT,"Phys. Med. Biol . , vol. 21, no. 5, pp . 733-744, 1976; R. Alvarez and E. Seppi , "A comparison of noise and dose in conventional and energy selective computed tomography," IEEE Transactions on Nuclear Science, vol. NS-26, no. 2, pp. 2853-2856, Apr. 1979; A. Coleman and M. Sinclair, "A beam-hardening correction using dual-energy computed tomography," Phys . Med. Biol . , vol. 30, no. 11, pp. 1251-1256, 1985; W. A. Kalender, W. Perman, J. Vetter, and E. Klotz, "Evaluation of a prototype dual-energy computed tomographic apparatus . I. Phantom studies," Med. Phys., vol. 13, no. 3, pp. 334-339, May/June 1986; J. Vetter, W. Perman, W. A. Kalender, R. Mazess, and J. Holden, "Evaluation of a prototype dual-energy computed tomographic apparatus . II . Determination of vertebral bone mineral content," Med. Phys., vol. 13, no. 3, pp. 340-343, May/June 1986; W. A. Kalender, E. Klotz, and L. Kostaridou, "An algorithm for noise suppression in dual energy CT material density images," IEEE Transactions on Medical Imaging, vol. 7, no . 3, pp . 218-224, Sept. 1988.Jet hardening artifacts is another useful effect, cf. R. Alvarez and A. Macovski, "Energy-selective reconstructions in x-ray CT," Phys. Med. Biol. , vol. 21, no. 5, pp. 733-744, 1976; R. Alvarez and E. Seppi, "A Comparison of noise and dose in conventional and energy selective computed tomography," IEEE Transactions on Nuclear Science, vol. NS-26, no. 2, pp. 2853-2856, Apr. 1979; A. Coleman and M. Sinclair, "A beam-hardening correction using dual-energy computed tomography," Phys. Med. Biol. , vol. 30, no. 11, pp. 1251-1256, 1985; WA Calender, W. Perman, J. Vetter, and E. Klotz, "Evaluation of a prototype dual-energy computed tomographic apparatus." I. Phantom studies, "Med. Phys., Vol. 13, no. 3, pp. 334-339, May / June 1986; J. Vetter, W. Perman, WA Kalender, R. Mazess, and J. Holden, "Evaluation of a Prototype Dual-Energy Computed Tomographic Apparatus." II. Determination of Vertebral Bone Mineral Content, "Med. Phys., Vol. 13, no. 3, pp. 340-343, May / June 1986; WA Calender, E. Klotz, and L. Kostaridou, "An algorithm for noise suppression in dual energy CT material density images," IEEE Transactions on Medical Imaging, vol. 7, no. 3, pp. 218-224, Sept. 1988.
Obwohl radiographische Zweispektren-Techniken bereits breite Anwendung bei Gepäckkontrollen oder aber auch bei der Knochendichtebestimmung finden, gab es bis dato nur eine klinische Produktimplementation von DECT, vgl. W. A.Although radiographic two-spectra techniques are already widely used in baggage inspection or even in bone density determination, to date only a clinical product implementation of DECT has been found, cf. W.A.
Kalender, W. Perman, J. Vetter, and E. Klotz, "Evaluation of a prototype dual-energy computed tomographic apparatus . I . Phantom studies," Med. Phys., vol. 13, no. 3, pp. 334-339, May/June 1986. Hierbei wurde das duale Spektrum durch ein schnelles Hin- und Herschalten zwischen U1 und U2 nach jederKalender, W. Perman, J. Vetter, and E. Klotz, "Evaluation of a prototype dual-energy computed tomographic apparatus." I. Phantom studies, "Med. Phys., Vol. 13, no. 3, pp. 334-339, May / June 1986. In doing so, the dual spectrum was changed by rapidly switching between U 1 and U 2 after each
Projektion erzeugt. Seit kurzer Zeit steht ein Zweispektren- Spiral-Kegelstrahl-CT mit zwei Röntgenquellen (Somatom Definition, Siemens Medical Solutions, Forchheim) zur Verfügung. Dieses Gerät erlaubt Zweispektren-Aufnahmen und eliminiert das problematische Umschalten zwischen zweiProjection generated. Recently, a two-spectral spiral cone-beam CT with two X-ray sources (Somatom Definition, Siemens Medical Solutions, Forchheim) has been available. This device allows two-spectra recordings and eliminates the problematic switching between two
Spannungen. Mit dieser neuen Technik erscheint eine Renaissance der DECT in naher Zukunft wahrscheinlich. Neben der Anwendung in klinischer CT besteht ein weiteres Ziel darin, das Potential von vorklinischer DECT mit Hilfe eines sehr schnellen In-Vivo Kegelstrahl-Mikro-CT-Scanners (TomoScope 30s, VAMP GmbH, Erlangen) zu untersuchen.Tensions. With this new technology, a renaissance of the DECT seems likely in the near future. In addition to clinical CT, there is another goal to study the potential of preclinical DECT using a very fast in-vivo cone-beam micro-CT scanner (TomoScope 30s, VAMP GmbH, Erlangen, Germany).
Bei der Zweispektren-CT liegt die Annahme zu Grunde, daß sich der Schwächungskoeffizient μ(F,E) , der vom Ort F und der Photonenenergie E abhängt, folgendermaßen zerlegen läßt:The two-spectra CT is based on the assumption that the attenuation coefficient μ (F, E), which depends on the location F and the photon energy E, can be decomposed as follows:
μ(F,E) = fχ(r>,(E)+/2(r>2(E)μ (F, E) = f χ (r>, (E) + / 2 (r> 2 (E)
mit /|(r) und /2(r) als die beiden Materialbilder. Die Εnergieabhängigkeiten ψ,(E) seien a priori bekannt und voneinander unabhängig. Um diese Materialbilder rekonstruieren zu können, müssen die entsprechenden Linienintegrale entlang der Linie Lwith / | (r) and / 2 (r) as the two material images. The energy dependencies ψ, (E) are known a priori and are independent of each other. In order to be able to reconstruct these material images, the corresponding line integrals along the line L have to be reconstructed
P1(L) =Kf1 = ^r/,(F) P1 (L) = Kf 1 = ^ r /, (F)
L bekannt sein. Im Folgenden wird aus Gründen der Übersicht die explizite Angabe der Abhängigkeiten von L vermieden. SR bezeichne den Operator der 2D-Radontransformation bzw. den Operator der 3D-Röntgentransformation. Die Messungen jedoch liefernL be known. In the following, the explicit specification of the dependencies of L is avoided for reasons of clarity. SR denotes the operator of the 2D radon transformation or the operator of the 3D x-ray transformation. However, the measurements deliver
wobei W1 (E) und W2 (E) die beiden detektierten Spektren darstellen, die während des Scans benutzt und auf eins normiert werden. Offensichtlich ist ^7 = q}(px,p2) . Materialzerlegung bedeutet nun, die Bestimmung der inversen P1 = PXg^g2) i wobei die Wj (E) und ψ±(E) genau bekannt sein müssen, was jedoch selten der Fall ist. Um diese Schwierigkeiten zu umgehen, werden empirische Kalibrierverfahren benutzt, die die Messungen bekannter Absorber mit bekannter Position den analytisch hergeleiteten Schnittlängen der Strahlen in den Absorbern gegenüberstellen, vgl. A. Coleman and M. Sinclair, "A beam-hardening correction using dual-energy computed tomography, " Phys . Med. Biol., vol. 30, no. 11, pp. 1251-1256, 1985.where W 1 (E) and W 2 (E) represent the two detected spectra used during the scan and normalized to one. Obviously ^ 7 = q } (p x , p 2 ). Material decomposition now means the determination of the inverse P 1 = PXg ^ g 2 ) i where the Wj (E) and ψ ± (E) must be known exactly, which is seldom the case. Around To avoid difficulties, empirical calibration methods are used which contrast the measurements of known absorbers of known position with the analytically derived cut lengths of the beams in the absorbers, cf. A. Coleman and M. Sinclair, "A beam-hardening correction using dual-energy computed tomography," Phys. Med. Biol., Vol. 30, no. 11, pp. 1251-1256, 1985.
Mit der hier vorgeschlagenen empirischen Kalibrierung (Empirical Dual Energy Calibration, EDEC) wird ein neuartiger empirischer Kalibrierungsalgorithmus präsentiert. Im Gegensatz zu anderen Methoden wird bei EDEC weder die Kenntnis des Spektrums noch der Schwächungskoeffizienten vorausgesetzt. Die Geometrie, Größe und Position des Kalibrierphantoms werden nicht benötigt. EDEC basiert auf ähnlichen Prinzipien wie die Empirische Cupping Korrektur (Empirical Cupping Correction, ECC; vgl. K. Sourbelle, M. Kachelrieß, and W. A. Kalender, "Empirical water precorrection for cone-beam computed tomography," IEEE Medical Imaging Conference Record, pp. 1871-1875, Oct. 2005.; M. Kachelrieß, K. Sourbelle, and W. Kalender, "Empirical cupping correction: A first order rawdata precorrection for cone-beam computed tomography," Med. Phys., vol. 33, no. 5, pp. 1269-1274, May 2006) .The Empirical Dual Energy Calibration (EDEC) proposed here presents a novel empirical calibration algorithm. Unlike other methods, EDEC requires neither knowledge of the spectrum nor the attenuation coefficients. The geometry, size and position of the calibration phantom are not needed. EDEC is based on similar principles as the Empirical Cupping Correction (ECC; see K. Sourbelle, M. Kachelrieß, and WA Calender, "Empirical Water Precorrection for Cone-Beam Computed Tomography", IEEE Medical Imaging Conference Record, pp. 431-145) 1871-1875, Oct. 2005, M. Kachelrieß, K. Sourbelle, and W. Kalender, "Empirical cupping correction: A first order raw data precorrection for cone-beam computed tomography," Med. Phys., Vol. no. 5, pp. 1269-1274, May 2006).
Im Folgenden sollen g^ die Polychromatisehen CT-Rohdatenwerte und pi die gesuchten monochromatischen, materialspezifischen Rohdatenwerte bezeichnen. Der Index J=I, 2 steht für die Röntgenenergien und i=l,2 für die Materialien, nach denen zerlegt wird. Sei P1=DXq^q2), wobei Di eine noch unbekannte Zerlegungsfunktion bezeichnet. Letztere setzt sich aus einer Linearkombination von Basisfunktionen bn(qλ,q2) zusammen: nIn the following, let g ^ denote the polychromatic CT raw data values and pi the monochromatic, material-specific raw data values sought. The index J = I, 2 stands for the X-ray energies and i = L, 2 for the materials to be decomposed. Let P 1 = DXq ^ q 2 ), where Di denotes a still unknown decomposition function. The latter is composed of a linear combination of basis functions b n (q λ , q 2 ): n
Für die erhaltenen Ergebnisse wurden Polynomfunktionen a]-s Bas isf unkt ionen mitFor the results obtained were polynomial functions a] - s Bas isf unkt ions with
k = 0... K, l = 0... L, n = k(L + l) + lk = 0 ... K, l = 0 ... L, n = k (L + l) + l
benutzt. Die Gesamtzahl der Basisfunktionen beträgt (K+l) (L+l) . Für die zu Grunde liegenden Daten wurde K=L=4 verwendet, was in diesem Fall 25 Basisfunktionen entspricht. Die Gesamtzahl der Koeffizienten in jedem Koeffizientenvektor C1 beträgt also 25 für jedes Material, wodurch die Gesamtzahl der unbekannten Koeffizienten auf 50 ansteigt. Die Aufgabe von EDEC ist es diese Koeffizienten C1 zu bestimmen. Es sei noch darauf hingewiesen, daß sich EDEC nicht auf Polynome beschränkt, es können ebenso andere Basisfunktionen gewählt werden.used. The total number of basis functions is (K + 1) (L + 1). For the underlying data K = L = 4 was used, which in this case corresponds to 25 basis functions. The total number of coefficients in each coefficient vector C 1 is thus 25 for each material, increasing the total number of unknown coefficients to 50. The task of EDEC is to determine these coefficients C 1 . It should be noted that EDEC is not limited to polynomials, other basic functions may also be chosen.
Um die Notation zu vereinfachen, wird im Folgenden der Materialindex i weg gelassen und die Zerlegungskoeffizienten c = c, für Material 1 bestimmt . Die Zerlegung nach Material 2 verläuft hierzu völlig analog.In order to simplify the notation, the material index i is left out in the following and the decomposition coefficients c = c, for material 1 are determined. The decomposition according to material 2 is completely analogous.
Aufgrund der Linearität der Röntgentransformation 5H kann eine Menge von (K+l) (L+l) Basisbildern durch Due to the linearity of the X-ray transformation 5H, an amount of (K + 1) (L + 1) basic images can pass through
definiert werden, d.h. fn sind die Rekonstruktionen der Basisfunktionen bn. Gesucht wird nach dem Satz von Koeffizienten c , der die Abweichung der Quadrateare defined, ie f n are the reconstructions of the basis functions b n . We search for the set of coefficients c, which is the deviation of the squares
zwischen den linear kombinierten Basisbildern nbetween the linear combined basic images n
und einer gegebenen Materialvorlage t(r) minimiert. Die Gewichtungsfunktion w(r) wird benutzt um ungewollte Strukturen des Kalibrierphantoms, das in dem Optimierungsprozeß vorkommt, auszuschließen.and a given template t (r) is minimized. The weighting function w (r) is used to eliminate unwanted structures of the calibration phantom that occurs in the optimization process.
Zur Bestimmung der Basisbilder, der Materialvorlage und des Gewichtungsbildes werden zwei Scans eines Kalibrierphantoms bei zwei verschiedenen Energien, und somit zwei verschiedenen Spektren, durchgeführt.To determine the basic images, the material template and the weighting image, two scans of a calibration phantom are performed at two different energies, and thus two different spectra.
Das Kalibrierphantom enthält homogene Bereiche mit ausreichenden Mengen an Material 1 und 2 und stellt sicher, daß alle sinnvollen Kombinationen der Pfadlängen in Material 1 und 2 erfaßt werden. Um die Basisbilder zu erhalten, werden diese Rohdaten den (K+l) (L+l) Basisfunktionen übergeben und rekonstruiert. Eine Standardrekonstruktion des Kalibrierphantoms dient als Grundlage zur Bestimmung von t(r) und w(r) , die durch Schwellwertsetzung bestimmt werden.The calibration phantom contains homogeneous regions with sufficient quantities of material 1 and 2 and ensures that all meaningful combinations of path lengths in material 1 and 2 are detected. To obtain the basic images, these raw data are transferred to the (K + 1) (L + 1) basis functions and reconstructed. A standard reconstruction of the calibration phantom serves as the basis for determining t (r) and w (r), which are determined by thresholding.
Die Materialvorlage /(F) repräsentiert die a priori Kenntnis der Regionen, welche Material 1 enthalten, und jener, die sicher nicht Material 1 beinhalten (d.h. Material 2 oder Luft) : 1 für r e Material 1The material template / (F) represents the a priori knowledge of the regions containing material 1 and those that certainly do not contain material 1 (ie material 2 or air): 1 for re material 1
'(?) = -'(?) = -
O für r G Luft und Material 2O for r G air and material 2
Es sei darauf hingewiesen, daß in Regionen, die eventuell eine Mischung der beiden Basismaterialien 1 und 2 enthalten, die Materialvorlage nicht definiert werden muß, da diese Regionen durch die Gewichtung unterdrückt werden. Das Bild der Materialvorlage wird überall dort gleich eins gesetzt, wo der Inhalt des Voxels r eindeutig ist. Das trifft zu in Regionen des Materials 1 oder 2 und in Luft.It should be noted that in regions which may contain a mixture of the two base materials 1 and 2, the material template need not be defined, since these regions are suppressed by the weighting. The image of the material template is set equal to one wherever the content of the voxel r is unique. This is true in regions of material 1 or 2 and in air.
Für Regionen außerhalb des Meßfeldes, bei Effekten, die der Punktbildfunktion zuzuschreiben sind (diese Bereiche werden mittels Erosion gefunden) , oder bei Materialien, die weder Basismaterial noch Luft sind, wird w(r) gleich Null gesetzt.For regions outside the measurement field, for effects attributable to the dot-image function (those areas are found by erosion), or for materials that are neither base material nor air, w (r) is set equal to zero.
erial erial
Jetzt kann E2 minimiert werden. Leitet man E2 nach Cn ab ergibt daß das lineare System a = B c mitNow E 2 can be minimized. If E 2 is derived from C n , it follows that the linear system a = B c
Die Lösung des Optimierungsproblems ist einfach c = B~x-a. Um das lineare System zu invertieren, sollte man bei einer tatsächlichen Implementierung vermeiden, B~λ explizit auszurechnen und statt dessen eher eine Gauß-Jordan Eliminierung bevorzugen.The solution to the optimization problem is simply c = B ~ x -a. In order to invert the linear system one should avoid, in an actual implementation, B ~ λ explicitly and instead prefer a Gauss-Jordan elimination.
Es wurden Rohdaten eines 45 cm großen Yin-Yang- Kalibrierphantoms bei 80 kV und 140 kV simuliert. Das Phantom bestand aus Wasser und Teflon. Zur Simulation wurden semi- empirische Spektren verwendet, vgl. D. M. Tucker, G. T. Barnes, and D. P. Chakraborty, "Semiempirical model for generating tungsten target x-ray spectra," Med. Phys . , vol. 18, no. 3, pp. 211-218, Mar. /Apr. 1991. Die Ergebnisse der Kalibrierung sind in Fig. 10 dargestellt.Raw data of a 45 cm Yin-Yang calibration phantom at 80 kV and 140 kV were simulated. The phantom was made of water and Teflon. Semi-empirical spectra were used for the simulation, cf. D.M. Tucker, G.T. Barnes, and D.P. Chakraborty, "Semiempirical model for generating tungsten target x-ray spectra," Med. Phys. , vol. 18, no. 3, pp. 211-218, Mar. / Apr. 1991. The results of the calibration are shown in FIG.
Fig. 10 zeigt von oben links nach unten rechts:10 shows from top left to bottom right:
a) Standardrekonstruktion bei 80 kV, b) Wasservorlage ti(r), c) Wassergewichtung wi(r), d) Zerlegung nach Material 1 fχ(r), e) Zerlegung nach Material 2 f2(r) und f) gewichtete Differenz W3. (r) (fi (r) -ti (r) ) .a) standard reconstruction at 80 kV, b) water ti (r), c) water weight wi (r), d) decomposition by material 1 fχ (r), e) decomposition by material 2 f 2 (r) and f) weighted difference W 3. (r) (fi (r) -ti (r)).
Der Standardscan ist mit einem Fenster (C=O HU / W=200 HU) abgebildet. Das Gewichtungs- und das Materialvorlagenbild sind binär und in schwarz/weiß dargestellt. Die Fensterung für die materialselektiven Bilder ist (C_100% / W=20%) und die gewichtete Differenz ist dargestellt mit (C=0% / W=5%) .The standard scan is shown with a window (C = O HU / W = 200 HU). The weight and material template images are shown in binary and in black and white. The windowing for the material-selective images is (C_100% / W = 20%) and the weighted difference is shown as (C = 0% / W = 5%).
Das Differenzbild vom Wasserbild und der Wasservorlage zeigt, daß EDEC sehr zufriedenstellend funktioniert . Zur weiteren Überprüfung des Kalibrierverfahrens wurde außerdem ein ovales CTDI -Testphantom mit den Abmessungen 32 cm x 16 cm simuliert, mit einem Hintergrund aus wasseräquivalentem Material und Einsätzen aus Teflon. Das Testphantom wurde mit den Zerlegungskoeffizienten c, und C2 , die mit dem Yin-Yang- Kalibrierphantom bestimmt wurden, exakt in Wasser und Teflon zerlegt .The difference image of the water image and the water pattern shows that EDEC works very satisfactorily. To further validate the calibration procedure, an oval CTDI test phantom measuring 32 cm x 16 cm was also simulated, with a background of water-equivalent material and Teflon inserts. The test phantom was with the decomposition coefficients c, and C 2 , which were determined with the Yin-Yang calibration phantom, exactly decomposed into water and Teflon.
EDEC wurde außerdem benutzt, um gemessene Daten eines klinischen Zweispektren-CT-Scanners und eines Mikro-CT- Scanners zu kalibrieren. Die Experimente mit realen Daten bestätigten die Ergebnisse der Simulationen.EDEC was also used to calibrate measured data from a clinical two-spectrum CT scanner and a micro CT scanner. The experiments with real data confirmed the results of the simulations.
Zusammenfassend kann festgehalten werden, daß die bildbasierte Empirische Zweispektren-Kalibrierung eine einfache, effektive und akkurate Methode ist, um Zweispektren-CT zu kalibrieren. In summary, it can be stated that image-based empirical two-spectra calibration is a simple, effective, and accurate method for calibrating two-spectrum CT.

Claims

AnsprüchePatentansprüche Claims Patent claims
1. Verfahren zur Kalibrierung von Computertomographen bei1. Method for calibrating computer tomographs
Zwei- oder Mehrspektrentomographie, gekennzeichnet durch ein direktes Berechnen der MaterialzerlegungsgleichungenTwo- or multi-spectral tomography, characterized by a direct calculation of the material decomposition equations
P1 = DJq3)P 1 = DJq 3 )
mit i=l , ... , I1 wobei i die Anzahl der Materialien bezeichnet, und j=l, ..., J1 wobei j die Anzahl der Spektren bezeichnet, wobei die Zerlegungsfunktion D1 als Linearkombination von Basisfunktionen Bn wie folgt darstellbar istwhere i = 1, ..., I 1 where i denotes the number of materials, and j = 1, ..., J 1 where j denotes the number of spectra, the decomposition function D 1 being a linear combination of basis functions B n follows is representable
D1 = ∑ cinBn{q3) nD 1 = Σ c in B n {q 3 ) n
und die Materialzerlegungskoeffizienten cin im Bildraum mittels einer Linearkombination von rekonstruierten Bildern bestimmt werden.and the material decomposition coefficients c in are determined in the image space by means of a linear combination of reconstructed images.
2. Verfahren nach Anspruch 1, wobei zur Bestimmung der Materialzerlegungskoeffizienten C1n ein Kalibrierphantom gescannt wird.2. The method of claim 1, wherein a calibration phantom is scanned to determine the material separation coefficients C 1n .
3. Verfahren nach Anspruch 2, gekennzeichnet durch die Verwendung eines Kalibrierphantoms beliebiger Form und Orientierung, mit entsprechend der Anzahl der Spektren zwei oder mehr segmentierbaren homogenen Bereichen bestehend aus Basismaterialien . 3. The method according to claim 2, characterized by the use of a Kalibrierphantoms any shape and orientation, with corresponding to the number of spectra two or more segmentable homogeneous areas consisting of base materials.
4. Verfahren nach einem der Ansprüche 1 bis 3, gekennzeichnet durch eine Kombination mit einem analytischen Korrekturverfahren .4. The method according to any one of claims 1 to 3, characterized by a combination with an analytical correction method.
5. Verwendung eines Verfahrens nach einem der Ansprüche 1 bis 4 zur Kalibrierung von Computertomographen . 5. Use of a method according to one of claims 1 to 4 for the calibration of computer tomographs.
EP07818508A 2006-10-18 2007-09-27 Calibration method for two-spectra or multispectra tomography Withdrawn EP2080015A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
DE102006049664A DE102006049664A1 (en) 2006-10-18 2006-10-18 Calibration method for two- or multi-spectral tomography
PCT/EP2007/008425 WO2008046498A1 (en) 2006-10-18 2007-09-27 Calibration method for two-spectra or multispectra tomography

Publications (1)

Publication Number Publication Date
EP2080015A1 true EP2080015A1 (en) 2009-07-22

Family

ID=38819960

Family Applications (1)

Application Number Title Priority Date Filing Date
EP07818508A Withdrawn EP2080015A1 (en) 2006-10-18 2007-09-27 Calibration method for two-spectra or multispectra tomography

Country Status (3)

Country Link
EP (1) EP2080015A1 (en)
DE (1) DE102006049664A1 (en)
WO (1) WO2008046498A1 (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102009053664A1 (en) 2009-11-17 2011-05-19 Ziehm Imaging Gmbh Method for the empirical determination of a correction function for the correction of radiation hardening and stray radiation effects in projection radiography and in computed tomography
DE102010040041B3 (en) * 2010-08-31 2012-01-26 Siemens Aktiengesellschaft Method for correcting artifacts due to temporal changes of attenuation values
WO2012176088A1 (en) 2011-06-21 2012-12-27 Koninklijke Philips Electronics N.V. Imaging apparatus
US10758202B2 (en) 2015-09-23 2020-09-01 Koninklijke Philips N.V. Spectral imaging phantom and method
EP3552547B1 (en) 2018-04-13 2022-12-07 Siemens Healthcare GmbH Method for providing conversion information about an imaging data set, x-ray device, computer program and electronically readable data carrier

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6173034B1 (en) * 1999-01-25 2001-01-09 Advanced Optical Technologies, Inc. Method for improved breast x-ray imaging

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
GLEASON S. ET AL: "Reconstruction of Multi-Energy X-Ray Computed Tomography Images of Laboratory Mice", IEEE TRANSACTIONS ON NUCLEAR SCIENCE, vol. 46, no. 4, 1 August 1999 (1999-08-01), pages 1081 - 1086, XP011088472
MARC KACHELRIEB ET AL: "Empirical Dual Energy Calibration (EDEC) for Cone-Beam Computed Tomography", NUCLEAR SCIENCE SYMPOSIUM CONFERENCE RECORD, 2006. IEEE,, 1 October 2006 (2006-10-01), pages 2456 - 2550, XP031083875
STENNER P. ET AL: "Empirical dual energy calibration (EDEC) for cone-beam computed tomography", MEDICAL PHYSICS, vol. 34, no. 9, 1 September 2007 (2007-09-01), pages 3630 - 2405, XP002462936
STENNER P.: "Empirical dual energy calibration (EDEC) for cone-beam computed tomography", MEDICAL PHYSICS, vol. 34, no. 9, 1 September 2007 (2007-09-01), pages 3630 - 3641, XP002462936

Also Published As

Publication number Publication date
WO2008046498A1 (en) 2008-04-24
DE102006049664A1 (en) 2008-05-08

Similar Documents

Publication Publication Date Title
DE2916486C2 (en)
DE10036142B4 (en) X-ray computed tomography device
DE102008028387B4 (en) A tomographic image reconstruction method for generating an image of an examination object and an imaging device operating according to this method
EP2150179B1 (en) Selection method for two contrast media for use in a dual-energy ct examination, contrast media combination and generation of ct images using a contrast media combination and different energy spectra
DE102006046191B4 (en) Stray radiation correction in radiography and computed tomography with area detectors
DE102012204980B4 (en) Method for the reconstruction of CT images with scatter correction, in particular for dual source CT devices
DE102016221658B4 (en) Scattered radiation compensation for a medical imaging device
DE102008030552A1 (en) A method for generating image data on a virtually predefinable x-ray tube voltage from first and second CT image data
DE102005008767A1 (en) Method for an X-ray device and computer tomograph for reducing beam hardening artifacts from a generated image of an object
DE102004029010A1 (en) Device and method for scattered radiation correction in projection radiography, in particular mammography
DE10356116A1 (en) Method and apparatus for facilitating artifact reduction
DE102005028216A1 (en) Computer tomography medical examination device comprises radiation source, detector, data memory, and evaluation unit connected downstream of detector for correcting radiation hardening of projection image
DE102009015773A1 (en) Method for processing reconstructed CT image data sets and CT system
DE102015207107A1 (en) Method for generating a virtual X-ray projection on the basis of an image data set obtained by means of an X-ray image recording device, computer program, data carrier and X-ray image recording device
EP1415179B1 (en) Device and method for creating a correction characteristic curve for reducing artefacts in tomography
DE102017200282B3 (en) Method for reducing image artifacts
DE102015225395A1 (en) Determination of a spatial distribution of a material property value on the basis of a single energy image acquisition
DE102004022332A1 (en) Method for the post-reconstructive correction of images of a computer tomograph
EP2080015A1 (en) Calibration method for two-spectra or multispectra tomography
DE102006026945B4 (en) Computed tomographic image recording method, method for determining the location-dependent concentration of a number of predetermined substances in an examination subject and associated computed tomography
DE102011005161B4 (en) Method, image data set reconstruction device, X-ray system and computer program for artifact correction
DE102016204709A1 (en) Image reconstruction method
DE102005018660B4 (en) Cupping correction for X-ray computed tomography
DE102009043213A1 (en) Efficient correction of polychromism effects during image reconstruction
DE102020212382B3 (en) Reconstruction of images with filter matrices

Legal Events

Date Code Title Description
PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

17P Request for examination filed

Effective date: 20090428

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IS IT LI LT LU LV MC MT NL PL PT RO SE SI SK TR

RAP1 Party data changed (applicant data changed or rights of an application transferred)

Owner name: CT IMAGING GMBH

TPAC Observations filed by third parties

Free format text: ORIGINAL CODE: EPIDOSNTIPA

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE APPLICATION IS DEEMED TO BE WITHDRAWN

18D Application deemed to be withdrawn

Effective date: 20100401