WO2006094493A2 - Korrektur von nichtlinearitäten eines abbildungssytems durch ein a-priori wissen bei durchstrahlungsaufnahmen - Google Patents

Korrektur von nichtlinearitäten eines abbildungssytems durch ein a-priori wissen bei durchstrahlungsaufnahmen Download PDF

Info

Publication number
WO2006094493A2
WO2006094493A2 PCT/DE2006/000420 DE2006000420W WO2006094493A2 WO 2006094493 A2 WO2006094493 A2 WO 2006094493A2 DE 2006000420 W DE2006000420 W DE 2006000420W WO 2006094493 A2 WO2006094493 A2 WO 2006094493A2
Authority
WO
WIPO (PCT)
Prior art keywords
reconstruction
correction
data
registration
specimen
Prior art date
Application number
PCT/DE2006/000420
Other languages
English (en)
French (fr)
Other versions
WO2006094493A3 (de
WO2006094493A8 (de
Inventor
Stefan Kasperl
Ingo Bauscher
Matthias Franz
Stefan Schroepfer
Original Assignee
Fraunhofer-Gesellschaft Zur Foerderung Der Angewandten Forschung E.V.
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 Fraunhofer-Gesellschaft Zur Foerderung Der Angewandten Forschung E.V. filed Critical Fraunhofer-Gesellschaft Zur Foerderung Der Angewandten Forschung E.V.
Priority to US11/908,267 priority Critical patent/US20080212734A1/en
Priority to CA002600648A priority patent/CA2600648A1/en
Priority to EP06722577A priority patent/EP1861734A2/de
Priority to DE112006001354T priority patent/DE112006001354A5/de
Publication of WO2006094493A2 publication Critical patent/WO2006094493A2/de
Publication of WO2006094493A3 publication Critical patent/WO2006094493A3/de
Publication of WO2006094493A8 publication Critical patent/WO2006094493A8/de

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/29—Measurement performed on radiation beams, e.g. position or section of the beam; Measurement of spatial distribution of radiation
    • G01T1/2914—Measurement of spatial distribution of radiation
    • G01T1/2985—In depth localisation, e.g. using positron emitters; Tomographic imaging (longitudinal and transverse section imaging; apparatus for radiation diagnosis sequentially in different planes, steroscopic radiation diagnosis)

Definitions

  • the technical field of application of the invention is an industrial quality control of specimens with regard to quantitative statements, e.g. Surveying tasks.
  • CT computed tomography
  • various physical effects cause artifacts in the reconstructed tomograms, which reduce image quality.
  • CT reconstructions must be as free of artifacts as possible, cf. WO-A 2003/062856 (Fraunhofer).
  • Previous correction methods for example for beam hardening or scattered beam correction according to the above WO-A, considerably reduce the artifacts and the image quality thus achieved already enables meaningful dimensional stability analyzes.
  • these methods work iteratively and require the presence of the complete projection data.
  • An initial CT reconstruction initially provides artifact-laden 3D voxel data from the candidate.
  • Postprocessing image processing steps determine therefrom correction parameters for an improved second CT reconstruction. If necessary, it will go through another iteration.
  • the input data required for the correction methods can no longer be correctly determined from the test sample itself.
  • the object of the invention is to provide a method for the on-line correction of non-linearities of the imaging system during data acquisition in industrial computed tomography (CT).
  • CT computed tomography
  • Claim 1 or claim 10 alternatively claim 20 solve this task, with the aid of target data of the specimen.
  • the X-ray tubes used in CT emit polychromatic radiation.
  • the interaction of X-radiation when passing through matter is e.g. energy-dependent.
  • Real system characteristics therefore have a nonlinear course caused by effects such as beam hardening, stray radiation and detector nonlinearities. This causes artifacts in the reconstructed layers, such as stripes, blurred edges, barrel distortions and cupping effects, which degrade image quality and make it difficult or even prevent surveying tasks.
  • the method claimed here corrects nonlinearities of the imaging system in computed tomography already during the data acquisition (claim 1) or calculates parameters used at least before the end of the data acquisition (acquisition process or "data acquisition 11" ).
  • the advantage over the prior art is that the claimed method manages with a single CT reconstruction. Time-consuming iterative post-processing steps (IAR) are eliminated.
  • the correction methods can resort to better input data, resulting in better quality CT reconstructions.
  • the procedure uses the target data of the DUT and provides input data for CT reconstruction correction methods.
  • correction The correction parameters are determined during data acquisition. A correction is made; either here or later.
  • initialization is meant a rough screened registration of the device under test.
  • a rough screening is therefore a registry whose accuracy
  • singular pairs of points are searched, with a singular point is one that is measurably highlighted by its environment.
  • Singular points can be ones that have a maximum or a minimum, both two-dimensional and one-dimensional. Measurable is the singular point emphasizing its environment.
  • Other possibilities of singular points to be understood are those that are edge points of the object shadow or intersections of edges.
  • a point of a digital model of a test specimen (usually a CAD model) is projected on the detector.
  • the singular point of the model and the singular point of the figure form a pair of points called "feature point”.
  • feature points can also be extracted. This extraction of the singular points mentioned (in terms of the most unique feature points possible) is carried out with search algorithms from the measurements.
  • the search algorithms are adapted to the simulated projection of the digital model.
  • the position can be registered at the beginning of the CT scan. This registration is made from a projection. Possible algorithms that can be used to perform this registration are the SoftPOSIT process, cf. DeMenthon et al., SoftPOSIT, Simultaneous Pose and Correspondence Determination, International Journal of Computer Vision, 59 (3), 2004, pages 259 to 284. This ability to register the start position is relatively insensitive to mismatched feature points when using the well-known SoftPOSIT process (couples), if they do not get too much.
  • the procedure for intensity-based registrations is to determine the similarity between the reference image and the template image. Similarities are obtained with statistical methods, all pixel information serves as the basis, cf. Penney et al., "A Comparison of Similarity Measurements for Use in 2-D-3-D Medical Image Registration," IEEE Transactions on Medical Imaging, 17 (4), 1998, pages 586 to 595. Intensity-based 2D or 3D registration algorithms optimize, starting from a sufficiently good starting value, the similarity of reference and transformed template, cf. Pluim, IEEE Transactions on Medical Imaging, 22 (8), 2003, pages 986 to 1004.
  • the CT model as target data of the test object and the a-priori knowledge used thereby can be used on several projections in different positions of the test object. Each layer is characterized by a different angle of rotation which the specimen assumes relative to a rotation axis.
  • the registration as a 2D registration or 3D registration takes place alternatively and prompted by the application.
  • From a 2D fan-beam CT can be easily generalized to a 3D cone-beam CT.
  • the manner of the detector which is designed either as a line detector in a 2D-CT or as an area detector in a 3D-CT, is tracked in each case. In both cases, attenuated intensities are imaged on the detector by the object and by the radiation with the measuring radiation from the punctiform source, in each case as a projection in each case at a rotation angle of the test object.
  • the ideal case is a perfectly adjusted CT imaging system.
  • the position of the axis of rotation needs to be known, around which the test object is rotated in angular increments.
  • the registration on some projections makes it possible to use the CT on remaining projections in such a way that the position of the test object can be calculated for further projections.
  • a simulation in the form of a virtual CT can be done on the basis of the above knowledge. It provides the required input data for correction procedures during the reconstruction. A correction, at least a provision of correction parameters, takes place during the data acquisition.
  • radiated lengths are associated with arbitrary detector locations (pixels) for each assumed incremental rotational position of the test object. Each irradiated length and an associated measured intensity at the detector are combined into pairs of values. In order to determine the correction data during the data acquisition, data from all projections is not necessary.
  • Some projections are sufficient, for example a representative selection that covers an angle range of less than 360 °, in particular clearly below it.
  • the determination of the correction data can already be started when this representative selection of projections has been recorded. So at least part of the calculation of the correction parameters and the further recording process run in parallel.
  • the calculation of the correction parameters can be completed or completed essentially with the end of the acquisition process, that is also those projections that are not required for the representative selection.
  • the reconstruction can take place in the temporal range after / at the conclusion of the recording, thus allowing a lesser delay until the result is available.
  • FIG. 1 is a schematic side view of an imaging system with a symbolization of a radiographic image, caused by a radiation source Q, measuring beams q, a test object 10 and a detector 31.
  • FIG. 1a is a plan view of the arrangement of FIG. 1, showing the turntable with its axis 100; the two edge points of the test object form the boundary rays of the fan of the measuring beams q, representing an intensity distribution at the detector 31; Layer represents, but can also image in three-dimensional CT a volume of the specimen in the form of a planar x / y extent, wherein the detector 31 is formed correspondingly flat.
  • FIG. 2 illustrates the stepwise change of the angular position of the
  • Figure 3 illustrates not to scale, but symbolically and for clarity greatly increases the registration of a specimen 11, which is shown in solid in its actual position 11 and is shown in its inaccurate coarse determination position with 11 'by dashed lines.
  • the difference angle as a registration error is denoted by ⁇ .
  • the beam source Q is much farther away from the device under test than represented by the symbolic distance z1, and the device under test 11 is further away from the detector than symbolically shows the distance z2.
  • FIG. 3 a is the intensity profile or the associated intensity profile in the x-direction (from top to bottom in FIG. 3) with respect to a point-shaped radiation source having a fan-shaped beam as measuring beams.
  • Significant feature points whose position is designated by xa, xe and xf, and which belong to the edge points 11a, 11e and 11f of the test object 11 from FIG. 3, become apparent therefrom.
  • Figure 4 illustrates a flow chart for the implementation of
  • the side view of Figure 1 shows a test piece 10 in an L-shape (in side view) and a radiation source Q, which can emit X-rays or neutron beams. These rays are denoted by q, either cone-shaped or fan-shaped for a 2D or a 3D tomography.
  • the axis 100 is the axis of rotation of a plate 20, which drives a shaft 21 via a drive 22 with a gear, which is rotationally rigidly coupled to the turntable 20.
  • the rotation is characterized by ⁇ (omega), the shaft 21 is rotatably mounted on a base surface 25.
  • the axis 100 is perpendicular to the radiation axis emanating from the source Q, penetrates the device under test 10 and extends to a screen 31, which is used as a detector.
  • an intensity distribution I is shown which has a two-dimensional shape in the case of a three-dimensional tomography with a weakened intensity course corresponding to the shape, shape and the material of the test object 10, as l (x, y).
  • l (x, y) In the case of a layer transmission and a fan-shaped steel q, only one height direction is to be measured, for example, which has an intensity distribution l (y).
  • FIG. 1a (without the test object 10) with a turntable 20 which is rotatable about the axis 100.
  • the marginal rays of the radiation source Q are shown, which just touch the turntable, the central beam axis as well as the intensity distribution l (x) in the horizontal direction on the detector 31st
  • a driving beam q1 is shown, which would irradiate the test piece 10 when placed on the turntable 20 and is located within the two object shadow lines (boundary beams).
  • the turntable 20 is gradually rotated by the drive 22 by angular increments ⁇ , as Figure 2 shows clearly.
  • Each time interval T1, T2 or T3, an angle increment is valid for a transmission of radiation from the radiation source Q.
  • the angle increments are each symbolized as equalized as 20a, 20b and 20c.
  • Figure 3 illustrates symbolically, but not to scale, the specimen here in a similar shape and designated as specimen 11 in the coarse registration.
  • An orientation of the test piece 11 is roughly determined in a first, rapid registration.
  • the specimen is in the position shown in solid lines with corner points 11a, 11e and 11f and is irradiated by the radiation source Q with the example. Fan beam q.
  • the beam axis is perpendicular to the detector plane 31, in a surface detector, in a line detector is only a dependence of x.
  • the position of the test object 11 is determined precisely via feature point pairs.
  • Other possibilities which will be described separately, are statistical methods which also achieve a positioning of the test piece that is more accurate than the first rough (fast) registration which identifies the coarse layer 11 'of the test piece.
  • an angular error ⁇ shown between the actual ply 11 and the registered ply 11 '.
  • the angle error ⁇ is above one degree.
  • a translational error may also occur which lies in the range above 1 mm to 2 mm (or is dimensioned on the test object as at least 1% of its largest, in particular typical length).
  • the distances z1, z2 are not to scale, but to understand symbolically.
  • FIG. 3a An intensity distribution is shown in Figure 3a, which results in a fan-shaped beam q.
  • the course of the fan-shaped beams from top to bottom, starting from the corner point 11a to the corner point 11f (respective boundary beam) shows the profile of Figure 3a corresponding to the increasing or decreasing thickness of the radiation-absorbing specimen 11.
  • the curve of the intensity l (x ) shows some singular points at the positions xa, xe and xf which correspond to the vertices 11a, 11e and 11f of the position of the specimen.
  • the function profile of FIG. 3a shifts in the x direction by a small amount.
  • Each singular point forms a pair of points with a corresponding model point in a digital model, usually on a CAD model of the test object.
  • Several such pairs of points can each achieve a more accurate registration of the candidate in a projection.
  • the measurements of the singular points on the detector can be understood as an extraction, in any case they lead to the positioning of the specimen beyond the coarse registration more accurate.
  • statistical Methods are used, as explained above.
  • the similarity between the reference image and the template image plays a role here, cf. Penney in IEEE Transactions as mentioned above. These statistical methods are intensity-based and allow more accurate registration.
  • the position of the test piece 11 relative to the axis of rotation possibly can be calculated with translational error for at least one further projection.
  • Influence of the target data of the test specimen from the digital model allows the said improvement of the coarse positioning of the specimen.
  • at least one further projection of the test object can be calculated. This can be done relative to the axis of rotation and / or with a translational displacement.
  • a virtual CT can be made with the help of the acquired knowledge.
  • This is a simulated CT that provides input data for a reconstruction reconstruction procedure. This is only possible when the rough registration was made.
  • Use of the correction data resulting from the simulation can either begin during the data acquisition process or only after completion of this data acquisition in the temporal area of the end of the acquisition process.
  • the necessary correction data which has already been determined during the data acquisition, are available at the end of the acquisition process, so that a rapid correction can be carried out, which does not first have to calculate the parameters for the correction from the data acquisition, but rather for a reconstruction at the conclusion of the acquisition Data collection already available.
  • a large time saving of the calculation method is the result.
  • the correction and thus the reconstruction at the end of the acquisition process can provide improved CT reconstruction.
  • the first reconstruction can thus work with correction data, which are available immediately at the end of the recording process, after they were previously determined during the data acquisition.
  • a correction can also be made during the recording process (data acquisition), even if only partially with regard to some of the artifacts that occur during the recording.
  • a reconstruction of the surveyed test object is therefore carried out with corrected measurement data and is not only available faster in time, but also qualitatively better.
  • FIG. 4 illustrates a symbolic signal flow diagram or flowchart of a data acquisition 70 that can be viewed in time, starting with its beginning on the left and its end on the right.
  • A-priori knowledge 69 is first given and allows a registration 71 which is coarse and can be made more accurate by using, for example, feature point pairs, each of which is a measurable singular point (e) on the detector 31 and associated with each one singular point (s) in the digital model.
  • the successful registration then allows a simulation 72 which is a virtual CT. It is used to provide input data for CT reconstruction correction procedures.
  • correction data are determined which can lead to a correction of the data of the data acquisition 70, which is symbolized by the arrows 73a.
  • a correction 73b can take place only when the data acquisition has been completed and the projection or data recording is transferred to the calculations "correction of the measured data" 74. From this correction, which can take place very quickly in terms of time, there results a reconstruction 75, which likewise can proceed very fast in order to obtain the corrected volume 11 * which forms the reconstruction.
  • the right edge of the block "recording" 70 symbolizes both the area before the immediate end by influencing the correction parameters via the influences 73a on the data acquisition, and / or the area 74,73b, which is downstream and the correction and the reconstruction is concerned.
  • Industrial quality control is a preferred field of application, in particular in the field of automotive engineering and based on castings as test items 10,11. X-rays were mentioned as preferred measuring beams.
  • the projections required for the parameter determination are less than all the images provided for a rotation angle of 360 °. Which are recorded in the increments ⁇ .

Abstract

Es witd ein Verfahren zur Online-Korrektur von Nichtlinearitäten des Abb ildungs Systems während der Datenaufnahme in der industriellen Computertomographie (CT) angegeben. Vorgeschlagen wird ein Verfahren zur Bereitstellung von korrigierten Projektionsdaten als eine verbesserte CT-Rekonstruktion, bei der Messstrahlen (g) ausgesendet werden, die den Prüfling (10,11) durchstrahlen und auf einem Detektor (31) in ihrer Intensität erfasst werden. Folgende Schritte sind vorgesehen. Eine erste Initialisierung, wobei eine Orientierung des Prüflings (10) mit einer schnellen Registrierung nur grobbestimmt wird. Eine Registrierung, wobei die Positionierung des Prüflings (10) genauer bestimmt wird, insbesondere über Merkmalspunktepaare. Eine Bewegung, wobei nach einer erfolgreich durchgeführten Registrierung an zumindest einigen Projektionen die Lage des Prüflings (10,11) für zumindest eine weitere Projektion berechnet wird. Eine Simulation, wobei mit Hilfe der Ergebnisse der vorigen Schritte eine virtuelle CT durchgeführt wird, welche Eingangsdaten für ein zu erfolgendes Korrekturverfahren der CT-Rekonstruktion liefert. Die Vornehmung einer Korrektur, wobei während der Datenaufnahme (70) vom Detektor (31) aus den Korrekturdaten Parameter bestimmt werden und eine Korrektur erfolgt (73a, 73b). Die Rekonstruktion, wobei am Endes des Aufnahmeprozesses als Datenaufnahme (70) korrigierte Projektionsdaten (11*) als verbesserte CT-Rekonstruktion vorlie en 74 75.

Description

Korrektur von Nichtlinearitaeten eines Abbildungssystems durch ein a-priori Wissen bei Durchstrahlungsaufnahmen
Das technische Anwendungsgebiet der Erfindung ist eine industrielle Qualitätskontrolle an Prüflingen im Hinblick auf quantitative Aussagen, z.B. Vermessungsaufgaben.
Bei der Computertomographie (CT) verursachen verschiedene physikalische Effekte Artefakte in den rekonstruierten Tomogrammen, welche die Bildqualität vermindern. Um Vermessungsaufgaben mit der gewünschten Genauigkeit und automatisiert durchführen zu können, müssen CT Rekonstruktionen aber möglichst artefaktfrei sein, vgl. WO-A 2003/062856 (Fraunhofer).
Bisherige Korrekturverfahren, beispielsweise zur Strahlaufhärtungs- bzw. Streustrahlkorrektur nach obiger WO-A reduzieren die Artefakte erheblich und die damit erzielte Bildqualität ermöglicht bereits sinnvolle Maßhaltigkeitsanalysen. Allerdings arbeiten diese Verfahren iterativ und setzen das Vorliegen der vollständigen Projektionsdaten voraus. Eine erste CT-Rekonstruktion liefert zunächst artefaktbehaftete 3D Voxeldaten des Prüflings. Nachverarbeitende Bildverarbeitungsschritte ermitteln daraus Korrekturparameter für eine verbesserte zweite CT-Rekonstruktion. Falls nötig, wird eine weitere Iteration durchlaufen. Bei stark artefaktbehafteten CT-Rekonstruktionen können die für die Korrekturverfahren benötigen Eingangsdaten nicht mehr korrekt aus dem Prüfling selbst ermittelt werden.
Aufgabe der Erfindung ist es, ein Verfahren zur Online-Korrektur von Nichtlinearitäten des Abbildungssystems während der Datenaufnahme in der industriellen Computertomographie (CT) anzugeben.
Anspruch 1 oder Anspruch 10 alternativ Anspruch 20 lösen diese Aufgabe, unter Zuhilfenahme von Solldaten des Prüflings.
Einer der wichtigsten Anwendungsfälle ist eine Gussteilefertigung (Anspruch 5 oder 15) in der Automobilindustrie. Qualitätskontrolle bei Gussteilen beinhaltet vor allem das Auffinden von Fehlstellen und das Prüfen von Abmessungen. Hauptaufgaben in der Vorserienentwicklung ist die schnelle Überprüfung der Maßhaltigkeit von Gussteilen mit komplexer Geometrie, sowie die Analyse von Abweichungen von den Solldaten. Unter dem Aspekt industrieller Anwendbarkeit im Vergleich mit anderen Quellen (Synchroton oder Gammastrahler) werden Röntgenröhren als Strahlungsquelle bevorzugt (Anspruch 4 oder Anspruch 14).
Allerdings emittieren die in der CT verwendeten Röntgenröhren polychromatische Strahlung. Die Wechselwirkung der Röntgenstrahlung beim Durchgang durch Materie ist z.B. energieabhängig. Reale Systemkennlinien haben daher einen nichtlinearen Verlauf, verursacht durch Effekte, wie Strahlaufhärtung, Streustrahlung und Nichtlinearitäten des Detektors. Dies verursacht in den rekonstruierten Schichten Artefakte, wie Streifen, unscharfe Kanten, tonnenförmige Verzeichnungen und Cupping- Effekte, welche die Bildqualität verschlechtern und Vermessungsaufgaben erschweren oder sogar verhindern.
Das hier beanspruchte Verfahren korrigiert Nichtlinearitäten des Abbildungssystems in der Computertomographie bereits während der Datenaufnahme (Anspruch 1) oder berechnet zumindest dazu verwendete Parameter vor dem Ende der Datenerfassung (Aufnahmeprozess oder kurz "Datenaufnahme11).
Dadurch wird die Bildqualität der Rekonstruktion erhöht und quantitative Aussagen ermöglicht, z.B. Vermessungsaufgaben erreicht, beinhaltend die Überprüfung der Maßhaltigkeit, oder Soll-Ist-Vergleiche des Prüfkörpers mit Solldaten zum Beispiel aus einem CAD-System.
Der Vorteil gegenüber dem Stand der Technik ist, dass das beanspruchte Verfahren mit einer einzigen CT-Rekonstruktion auskommt. Zeitaufwendige iterative Nachverarbeitungsschritte (IAR) entfallen. Durch die Verwendung der Solldaten (des Prüflings) als a-priori Wissen können die Korrekturverfahren auf bessere Eingangsdaten zurückgreifen, was sich in qualitativ besseren CT-Rekonstruktionen äußert.
Das Verfahren verwendet die Solldaten des Prüflings und liefert Eingangsdaten für Korrekturverfahren der CT-Rekonstruktion.
Es ist ein mehrstufiges Verfahren. Einzelne Stufen sind:
• Initialisierung: Die Orientierung des Prüflings wird mit einer ersten schnellen Registrierung grob bestimmt. • Registrierung: Ausgehend von der groben Positionierung wird eine auf Merkmale und/oder auf Intensitäten gestützte Registrierung durchgeführt. Dies ist eine genauere Registrierung.
• Bewegung: Nach erfolgreich durchgeführter Registrierung an einigen Projektionen kann die Lage des Prüflings, z.B. relativ zur Drehachse für weitere Projektionen berechnet werden.
• Simulation: Mit Hilfe dieses Wissens kann eine virtuelle CT simuliert werden, welche die benötigten Eingangsdaten für die Korrekturverfahren der CT- Rekonstruktion liefert.
• Korrektur: Die Korrekturparameter werden während der Datenaufnahme bestimmt. Eine Korrektur erfolgt; entweder hier oder später.
• Rekonstruktion: Am Ende des Aufnahmeprozesses liegen korrigierte Projektionsdaten für eine verbesserte CT-Rekonstruktion des Prüflings vor.
Mit der Initialisierung ist eine grob gerasterte Registrierung des Prüflings gemeint. Eine grobe Rasterung ist also eine Registrierung, deren Genauigkeit
• bezüglich der Rotation wenige Grad, insbesondere oberhalb eines Winkelfehlers von etwa einem Grad; und/oder
• bezüglich der Translation von etwa 1 mm bis 2 mm oder in einem Bereich von 1 % einer typischen Prüflingsabmessung (Anspruch 16,19)
beträgt. Damit bildet sich ein Startwert, der für eine im folgenden ablaufende genauere Registrierung benutzt wird. Dazu werden bspw. Merkmalspunkte verwendet. Diese sind bestimmte Paare.
Die genaue Registrierung erfolgt in Abhängigkeit vom Prüfling merkmals- und/oder intensitäts"basiert" im Sinne einer Auswertung dieser Messwerte.
Merkmalsbasierte Registrierung:
Nach erfolgter grober Registrierung, beispielsweise ein Festlegen eines grob gerasterten Winkelwerts eines drehbar gelagerten Prüflings, ggf. auch mit einer zugehörigen Translation, werden singuläre Punktepaare gesucht, wobei ein singulärer Punkt ein solcher ist, der sich messbar von seiner Umgebung hervorhebt. Diese singulären Punkte können einmal solche sein, die ein Maximum oder ein Minimum besitzen, sowohl zweidimensional wie auch eindimensional. Messbar ist der sich von seiner Umgebung hervorhebende singuläre Punkt. Andere Möglichkeiten von zu verstehenden singulären Punkten sind solche, die Randpunkte des Objektschattens sind oder Schnittpunkte von Kanten darstellen.
Ein Punkt eines digitalen Modells eines Prüflings (meist eines CAD-Modells) wird bei Projektion auf dem Detektor abgebildet. Der singuläre Punkt des Modells und der singuläre Punkt der Abbildung bilden ein Punktepaar, das als "Merkmalspunkt" bezeichnet wird.
Ist das Modell grob registriert, können Projektionen simuliert werden. Durch diese Simulation ergeben sich - entsprechend der groben Registrierung - ungefähre Positionen von Projektionen von Modell-Merkmalspunkten. Diese sind der Berechnung bekannt. Eine solche Kenntnis kann aber erst durch die grobe Registrierung des CAD- Modells erfolgen, welche in Folge die Simulation zur ungefähren Position der Projektion bringt.
Aus den Messungen lassen sich ebenfalls Merkmalspunkte extrahieren. Diese Extraktion von den genannten singulären Punkten (im Sinne möglichst eindeutiger Merkmalspunkte) erfolgt mit Suchalgorithmen aus den Messungen. Die Suchalgorithmen sind an die simulierte Projektion des digitalen Modells angepasst.
Bestehen nunmehr Merkmalspunkte (als Punktepaare), kann die Position zu Beginn der CT-Aufnahme registriert werden. Diese Registrierung erfolgt aus einer Projektion. Mögliche verwendbare Algorithmen, um diese Registrierung vorzunehmen, sind der Prozess SoftPOSIT, vgl. DeMenthon et al., SoftPOSIT, Simultaneous Pose and Correspondence Determination, International Journal of Computer Vision, 59 (3), 2004, Seiten 259 bis 284. Diese Möglichkeit der Registrierung der Startposition ist bei Anwendung des bekannten Prozesses SoftPOSIT relativ unempfindlich gegen falsch zugeordnete Merkmalspunkte(paare), wenn diese nicht zu sehr überhand nehmen.
Intensitätsbasierte Registrierung:
Das Vorgehen bei intensitätsbasierter Registrierungen ist es, die Ähnlichkeit zwischen Referenz- und Schablonenbild zu bestimmen. Hierbei werden Ähnlichkeiten mit statistischen Methoden gewonnen, als Grundlage dienen alle Pixelinformationen, vgl. Penney et al., "A Comparison of Similarity Measures for Use in 2-D-3-D Medical Image Registration", IEEE Transactions on Medical Imaging, 17(4), 1998, Seiten 586 bis 595. Intensitätsbasierte 2D oder 3D-Registrierungsalgorithmen optimieren, ausgehend von einem hinreichend guten Startwert, die Ähnlichkeit von Referenz und transformierter Schablone, vgl. Pluim, IEEE Transactions on Medical Imaging, 22(8), 2003, Seiten 986 bis 1004.
So kann ein A-priori Wissen ausgenutzt werden, eine Registrierung erfolgreich durchzuführen. Das CT-Modell als Soll-Daten des Prüflings und das dadurch eingesetzte A-priori Wissen können an mehreren Projektionen in unterschiedlichen Lagen des Prüflings eingesetzt werden. Jede Lage ist durch einen anderen Drehwinkel charakterisiert, den der Prüfling relativ zu einer Drehachse einnimmt.
Die Registrierung als 2D-Registrierung oder 3D-Registrierung erfolgt alternativ und veranlasst vom Anwendungsfall. Aus einer 2D-Fächerstrahl-CT kann problemlos auf eine 3D-Konusstrahl-CT generalisiert werden. Jeweils nachgeführt wird die Art und Weise des Detektors, der entweder als Zeilendetektor bei einer 2D-CT ausgebildet ist, oder als Flächendetektor bei einer 3D-CT. In beiden Annahmefällen werden durch das Objekt und durch die Durchstrahlung mit der Messstrahlung von der punktförmigen Quelle aus abgeschwächte Intensitäten auf dem Detektor abgebildet, jeweils als eine Projektion bei jeweils einem Drehwinkel des Prüflings.
Der Idealfall ist ein perfekt justiertes CT-Abbildungssystem. Hier braucht lediglich noch die Lage der Drehachse bekannt zu sein, um den der Prüfling in Winkelinkrementen verdreht wird.
Diese Winkelinkremente zwischen den Aufnahmepositionen des Prüflings sind hinreichend genau bekannt. Die Registrierung ebenfalls. Mit einem registrierten digitalen Modell des Prüflings ist es jetzt möglich, eine CT-Simulation vorzunehmen. Diese kann für beliebige Detektorpixel auf dem Detektor bei einer beliebigen Drehlage des Prüflings die zugehörig eine von einem von der punktförmigen Quelle jeweils gedacht ausgehenden Messstrahl durchstrahlte Länge liefert.
Die Registrierung an einigen Projektionen erlaubt es, an verbliebenen Projektionen die CT so einzusetzen, dass die Lage des Prüflings für weitere Projektionen berechenbar wird.
Eine Simulation in Form einer virtuellen CT kann auf Grund des vorstehenden Wissens erfolgen. Sie liefert die benötigten Eingangsdaten für Korrekturverfahren bei der Rekonstruktion. Eine Korrektur, zumindest eine Bereitstellung von Korrekturparametern, erfolgt während der Datenaufnahme. Bei der virtuellen CT entstehen für beliebige Detektorstelien (Pixel) bei jeder angenommenen inkrementellen Drehlage des Prüflings zugehörig durchstrahlte Längen. Eine jeweils durchstrahlte Länge und eine zugehörig gemessene Intensität am Detektor werden zu Wertepaaren kombiniert. Um die Korrekturdaten bei der Datenaufnahme zu bestimmen, sind nicht Daten von allen Projektionen notwendig.
Einige Projektionen reichen, beispielsweise eine repräsentative Auswahl, die einen Winkelbereich von unter 360° abdeckt, insbesondere deutlich darunter. Nachdem die Korrekturdaten schon während der Datenaufnahme bestimmt werden und nicht alle Projektionen als Eingangsgröße notwendig sind, die Korrekturparameter zu bestimmen, kann mit der Bestimmung der Korrekturdaten bereits begonnen werden, wenn diese repräsentative Auswahl von Projektionen aufgenommen wurde. So läuft zumindest ein Teil der Berechnung der Korrekturparameter und der weitere Aufnahmeprozess parallel. Bevorzugt kann die Berechnung der Korrekturparameter abgeschlossen sein, oder abgeschlossen werden im Wesentlichen mit dem Ende des Aufnahmeprozesses, also auch derjenigen Projektionen, die für die repräsentative Auswahl nicht benötigt werden. Die Rekonstruktion kann im zeitlichen Bereich nach/bei Abschluss der Aufnahme erfolgen, erlaubt also eine geringere Verzögerung bis zum Vorliegen des Ergebnisses.
Als Korrekturverfahren können solche Verfahren angewendet werden, die in "Quality Improvements for Cone-beam CT using Beam Hardening and Scattering Correction", Third World Congress on Industrial Process Tomography, Banff, Canada, 2002, Seiten 90 bis 95 herangezogen werden. Für die Rekonstruktion liegen bereits korrigierte Projektionsdaten vor, so dass die erste Rekonstruktion bereits mit Korrekturdaten arbeiten kann. Eine Rekonstruktion kann bereits auf korrigierte Messdaten zurückgreifen. Bereits die erste Rekonstruktion liefert ein vollständig korrigiertes Volumen des rekonstruierten Prüflings. Es ergibt sich eine verbesserte CT-Rekonstruktion.
Die für die Korrektur verwendeten Eingangsdaten sind besser, was sich in einer qualitativ besseren CT-Rekonstruktion niederschlägt. Die Erfindung wird nachfolgend an Hand von mehreren Ausführungsbeispielen erläutert und ergänzt und dabei näher beschrieben.
Figur 1 ist eine schematische Seitenansicht eines Abbildungssystems mit einer Symbolisierung einer Durchstrahlungsaufnahme, veranlasst von einer Strahlenquelle Q, Messstrahlen q, einem Prüfling 10 und einem Detektor 31.
Figur 1 a ist eine Aufsicht auf die Anordnung von Figur 1 , aus der der Drehteller mit seiner Achse 100 hervorgeht, die beiden Randpunkte des Prüflings bilden die Grenzstrahlen des Fächers der Messstrahlen q, zur Abbildung einer Intensitätsverteilung am Detektor 31 , der für eine Ebene eine Schicht darstellt, der aber auch bei dreidimensionaler CT ein Volumen des Prüflings in Form einer flächigen x/y-Erstreckung abbilden kann, wobei der Detektor 31 entsprechend flächig ausgebildet ist.
Figur 2 veranschaulicht die stufenweise Veränderung der Winkellage des
Prüflings jeweils um einen Differenzwinkel Δα.
Figur 3 veranschaulicht nicht maßstabgerecht, sondern symbolisch und zur Verdeutlichung stark vergrößert die Registrierung eines Prüflings 11 , der in seiner tatsächlichen Lage 11 durchgezogen dargestellt ist und in seiner ungenau liegenden Grobbestimmungslage mit 11' strichliniert dargestellt ist. Der Differenzwinkel als Registrierfehler ist mit γ bezeichnet. Die Strahlquelle Q ist wesentlich weiter von dem Prüfling entfernt, als mit dem symbolischen Abstand z1 dargestellt, auch ist der Prüfling 11 weiter vom Detektor entfernt, als dass symbolisch der Abstand z2 zeigt.
Figur 3a ist das Intensitätsprofil oder der zugehörige Intensitätsverlauf in x-Richtung (in Figur 3 von oben nach unten), bezogen auf eine punktförmige Strahlenquelle mit einem fächerförmigen Strahl als Messstrahlen. Es werden daraus wesentliche Merkmalspunkte deutlich, deren Position mit xa, xe und xf bezeichnet ist, und die zu den Randpunkten 11a, 11e und 11f des Prüflings 11 aus Figur 3 gehören. Figur 4 veranschaulicht eine Ablaufskizze für die Durchführung der
Rekonstruktion mit teilweise parallel laufender Bestimmung von Korrekturparametern, sodass die korrigierten Messdaten der ersten Rekonstruktion bereits ein vollständig korrigiertes Volumen 11* rekonstruieren kann.
Die Seitenansicht der Figur 1 zeigt einen Prüfling 10 in einer L-förmigen Gestalt (in Seitenansicht) und eine Strahlenquelle Q, die Röntgenstrahlen oder Neutronenstrahlen abgeben kann. Diese Strahlen sind mit q bezeichnet, entweder konusförmig oder fächerförmig für eine 2D- oder eine 3D-Tomografie. Die Achse 100 ist die Drehachse eines Tellers 20, der über einen Antrieb 22 mit einem Getriebe einen Schaft 21 antreibt, der mit dem Drehteller 20 drehsteif gekoppelt ist. Die Drehung ist mit ω (omega) gekennzeichnet, der Schaft 21 ist auf einer Sockelfläche 25 drehbar gelagert.
Die Achse 100 steht senkrecht auf der Strahlungsachse, welche von der Quelle Q ausgeht, den Prüfling 10 durchdringt und zu einem Schirm 31 reicht, der als Detektor eingesetzt wird. In Höhenrichtung der Darstellung ist eine Intensitätsverteilung I gezeigt, die bei einer dreidimensionalen Tomografie mit einem entsprechend der Form, Gestalt und dem Werkstoff des Prüflings 10 geschwächten Intensitätsverlauf eine zweidimensionale Form hat, als l(x,y). Bei einer Schichtdurchstrahlung und einem fächerförmigen Stahl q ist nur eine Höhenrichtung bspw. zu messen, welche eine Intensitätsverteilung l(y) aufweist. Dieses sind die zu erfassenden Daten einer Durchstrahlungsaufnahme.
In einer Ansicht von Oben zeigt sich diese Anordnung als Figur 1a (ohne den Prüfling 10) mit einem Drehteller 20 der um die Achse 100 drehbar ist. Die Randstrahlen der Strahlenquelle Q sind eingezeichnet, welche den Drehteller gerade noch berühren, die zentrale Strahlenachse ebenso, wie auch die Intensitätsverteilung l(x) in horizontaler Richtung auf dem Detektor 31.
Ein Fahrstrahl q1 ist eingezeichnet, welcher den Prüfling 10 bei Aufsetzen auf den Drehteller 20 durchstrahlen würde und innerhalb der beiden Objektschattenlinien (Grenzstrahlen) gelegen ist.
Der Drehteller 20 ist mit dem Antrieb 22 stufenweise um Winkelinkremente Δα verdrehbar, wie Figur 2 anschaulich zeigt. Jeweils eine bestimmte Zeitspanne T1 , T2 oder T3 ist ein Winkelinkrement gültig für eine Durchstrahlungsaufnahme von der Strahlenquelle Q. Die Winkelinkremente sind als 20a, 20b und 20c jeweils gleich gestuft symbolisiert. Figur 3 veranschaulicht symbolisch, aber nicht maßstabsgerecht den Prüfling hier in einer ähnlichen Gestalt und als Prüfling 11 bezeichnet bei der groben Registrierung.
Eine Orientierung des Prüflings 11 wird in einer ersten, schnellen Registrierung grob bestimmt. Dabei liegt der Prüfling in der durchgezogen dargestellten Lage mit Eckpunkten 11a, 11e und 11f und wird von der Strahlenquelle Q mit dem bspw. Fächerstrahl q durchstrahlt. Die Strahlenachse ist senkrecht auf der Detektorebene 31 , bei einem Flächendetektor, bei einem Zeilendetektor handelt es sich nur um eine Abhängigkeit von x. Die Position des Prüflings 11 wird über Merkmalspunktepaare genau bestimmt. Andere Möglichkeiten, die gesondert beschrieben werden, sind statistische Methoden, welche ebenfalls eine Positionierung des Prüflings erreichen, die genauer ist als die erste grobe (schnelle) Registrierung, welche die Groblage 11' des Prüflings identifiziert. Hier ist ein Winkelfehler γ anzunehmen, der zwischen der tatsächlichen Lage 11 und der registrierten Lage 11' gezeigt ist. Der Winkelfehler γ liegt oberhalb von einem Grad. Es kann zusätzlich auch ein translatorischer Fehler hinzutreten, der im Bereich oberhalb von 1mm bis 2mm liegt (oder am Prüfling bemessen ist als zumindest 1 % seiner größten, insbesondere typischen Länge).
Die Abstände z1 , z2 sind nicht maßstabsgetreu, sondern symbolisch zu verstehen.
Eine Intensitätsverteilung ist in Figur 3a gezeigt, welche sich bei einem fächerförmigen Strahl q ergibt. Der Verlauf der fächerförmigen Strahlen von oben nach unten, beginnend von dem Eckpunkt 11a bis zum Eckpunkt 11f (jeweiliger Grenzstrahl) zeigt den Verlauf von Figur 3a entsprechend der stärker zunehmenden oder abnehmenden Dicke des die Strahlung absorbierenden Prüflings 11. Die Verlaufskurve der Intensität l(x) zeigt einige singuläre Punkte an den Positionen xa, xe und xf, welche den Eckpunkten 11a, 11e und 11f der Position des Prüflings entsprechen. Bei entsprechend ungenauer Registrierung 11' verschiebt sich der Funktionsverlauf von Figur 3a in x- Richtung um einen geringen Betrag.
Jeder singuläre Punkt bildet ein Punktepaar mit einem entsprechenden Modellpunkt in einem digitalen Modell, meist an einem CAD-Modell des Prüflings. Mehrere solche Punktepaare können jeweils eine genauere Registrierung des Prüflings in einer Projektion erreichen.
Die Messungen der singulären Punkte auf dem Detektor können als Extraktion verstanden werden, jedenfalls führen sie dazu, dass die Positionierung des Prüflings über die grobe Registrierung hinaus genauer wird. Alternativ können statistische Methoden angewendet werden, wie eingangs erläutert. Die Ähnlichkeit zwischen Referenzbild und Schablonenbild spielt hierbei eine Rolle, vgl. Penney in IEEE Transactions wie eingangs angegeben. Diese statistischen Methoden arbeiten intensitätsbasiert und erlauben die genauere Registrierung.
Wenn insoweit von einer "hinreichend genauen Bestimmung" die Rede ist, ist diese jedenfalls genauer als die schnelle grobe Registrierung und die Grobbestimmung der Lage des Prüflings, die Ausgangspunkt war.
Nach einer erfolgreich durchgeführten Registrierung an zumindest einigen Projektionen kann die Lage des Prüflings 11 relativ zur Drehachse, ggf. auf mit translatorischem Fehler für zumindest eine weitere Projektion berechnet werden.
Einfluss der Solldaten des Prüflings aus dem digitalen Modell erlaubt die genannte Verbesserung der Grobpositionierung des Prüflings. Nach einer solchen durchgeführten Registrierung kann zumindest eine weitere Projektion des Prüflings berechnet werden. Dies kann relativ zur Drehachse und/oder mit einer translatorischen Verschiebung geschehen.
Nach einer Registrierung des Prüflings kann mit Hilfe des erworbenen Wissens eine virtuelle CT erfolgen. Diese ist eine simulierte CT, mit der Eingangsdaten für ein Korrekturverfahren für die Rekonstruktion geliefert werden. Dies ist erst dann möglich, wenn die grobe Registrierung erfolgte. Eine Verwendung der Korrekturdaten, die aus der Simulation entstehen, kann entweder schon während der erfolgenden Datenerfassung einsetzen, oder erst nach Abschluss dieser Datenerfassung im zeitlichen Bereich des Endes des Aufnahmeprozesses.
Die notwendigen Korrekturdaten, die schon während der Datenaufnahme bestimmt worden sind, stehen beim Ende des Aufnahmeprozesse zur Verfügung, so dass eine schnelle Korrektur erfolgen kann, die nicht erst die Parameter zur Korrektur aus der Datenerfassung errechnen muss, sondern diese für eine Rekonstruktion beim Abschluss der Datenerfassung schon verfügbar hat. Eine große zeitliche Einsparung des Berechnungsverfahrens ist die Folge.
Aus den Korrekturdaten, die bereits während der Datenaufnahme bestimmt wurden, kann die Korrektur und damit die Rekonstruktion zum Ende des Aufnahmeprozesses eine verbesserte CT-Rekonstruktion liefern. Bereits die erste Rekonstruktion kann so mit Korrekturdaten arbeiten, die unmittelbar zum Ende des Aufnahmeprozesse zur Verfügung stehen, nachdem sie zuvor während der Datenaufnahme bestimmt wurden. Eine Korrektur kann aber auch bereits während des Aufnahmeprozesses (der Datenaufnahme) erfolgen, und sei es auch nur teilweise hinsichtlich einiger der Artefakte, die bei der Aufzeichnung entstehen. Eine Rekonstruktion des vermessenen Prüflings erfolgt also mit korrigierten Messdaten und ist nicht nur zeitlich schneller verfügbar, sondern auch qualitativ besser.
Die Figur 4 veranschaulicht einen symbolischen Signalflussplan oder Ablaufplan einer Datenaufnahme 70, die zeitlich betrachtet werden kann, beginnend links mit ihrem Anfang und rechts mit ihrem Ende. A-priori Wissen 69 wird zunächst vorgegeben und erlaubt eine Registrierung 71 die grob ist und genauer ausgestaltet werden kann, durch Einsatz von bspw. Merkmalspunkte-Paaren, die jeweils ein messbarer singulärer Punkt(e) auf dem Detektor 31 sind und mit einem jeweils zugehörigen singulären Punkt(en) im digitalen Modell gepaart werden. Die erfolgreiche Registrierung erlaubt dann eine Simulation 72 die eine virtuelle CT ist. Mit ihr werden Eingangsdaten für Korrekturverfahren der CT-Rekonstruktion geliefert.
Bei der Bestimmung der Korrekturdaten 73, die bereits während der Datenaufnahme 70 abläuft, werden Korrekturdaten bestimmt, die zu einer Korrektur der Daten der Datenerfassung 70 führen kann, was durch die Pfeile 73a symbolisiert wird. Ebenso kann in einem alternativen Ausführungsbeispiel, oder auch kumulativ, eine Korrektur 73b erst anschließend erfolgen, wenn die Datenaufnahme abgeschlossen ist und die Projektion oder Datenaufnahme den Berechnungen "Korrektur der Messdaten" 74 übergeben wird. Aus dieser Korrektur, die zeitlich sehr schnell ablaufen kann, ergibt sich eine Rekonstruktion 75, die ebenfalls sehr schnell ablaufen kann, um das korrigierte Volumen 11* zu erhalten, welche die Rekonstruktion bildet.
Zum Ende des Aufnahmeprozesses symbolisiert den rechten Rand des Blocks "Aufnahme" 70, sowohl den Bereich vor dem unmittelbaren Ende durch Einflussnahme der Korrekturparameter über die Einflüsse 73a auf die Datenerfassung, und/oder den Bereich 74,73b, der nachgelagert ist und die Korrektur sowie die Rekonstruktion betrifft.
Die industrielle Qualitätskontrolle ist ein bevorzugtes Anwendungsgebiet, insbesondere im Bereich des Automobilbaus und bezogen auf Gussteile als Prüflinge 10,11. Röntgenstrahlen waren als bevorzugte Messstrahlen erwähnt. Durch die Art und Weise des Aufbaus nach Figur 4 können Artefakte auch ohne Iteration reduziert werden, und dies mit großem zeitlichen Vorteil. Die für die Parameterbestimmung erforderlichen Projektionen sind weniger als alle für einen Drehwinkel von 360° zur Verfügung gestellten Aufnahmen. Die in den Inkrementen Δα erfasst werden.
^P ^P ^^ ^^

Claims

Ansprüche:
1. Mehrstufiges Verfahren zur Bereitstellung von korrigierten Projektionsdaten als eine verbesserte CT-Rekonstruktion, mit folgenden Schritten
• Initialisierung, wobei eine Orientierung des Prüflings (10,11) mit einer ersten schnellen Registrierung grob bestimmt wird;
• Registrierung, wobei die Position des Prüflings (10) über Merkmalspunktepaare und/oder statistische Methoden hinreichend genau bestimmt wird;
• Bewegung, wobei nach einer erfolgreich durchgeführten Registrierung an zumindest einigen Projektionen die Lage des Prüflings insbesondere relativ zur Drehachse für zumindest eine weitere Projektion berechnet wird;
• Merkmalspunkte, wobei die grobe Registrierung zur Extraktion von bevorzugt möglichst eindeutigen Merkmalspunkten genutzt wird; o Startposition, wobei eine 2D oder 3D-Registrierung an Solldaten des Prüflings mit den Messdaten anhand der extrahierten Merkmalspunkte erfolgt; o Bewegung, wobei nach einer erfolgreich durchgeführten Registrierung an zumindest einigen Projektionen die Lage des Prüflings relativ zur Drehachse für zumindest eine weitere Projektion berechnet wird;
• Simulation, wobei mit Hilfe dieses Wissens eine virtuelle CT simuliert wird, um Eingangsdaten für ein Korrekturverfahren der CT-Rekonstruktion zu liefern;
• Korrektur, wobei aus den Korrekturdaten während der Datenaufnahme (70) Parameter bestimmt werden (73) und eine Korrektur erfolgt (73a;73b);
• Rekonstruktion, wobei zum Ende des Aufnahmeprozesses korrigierte Projektionsdaten für eine verbesserte CT-Rekonstruktion vorliegen (74,75).
2. Verfahren nach Anspruch 1 , wobei die Bereitstellung bei der industriellen Qualitätskontrolle erfolgt.
3. Verfahren nach Anspruch 2, wobei zumindest eine Vermessung am Prüfling erfolgt.
4. Verfahren nach Anspruch 1 , wobei Röntgenstrahlen bei der CT-Erstellung verwendet werden.
5. Verfahren nach Anspruch 3 oder 1 , wobei der Prüfling ein Gussteil ist, insbesondere im Automobilbau.
6. Verfahren nach Anspruch 1 , bei dem keine Iteration zur Reduzierung der Artefakte in den Schichtrekonstruktionen eingesetzt wird
7. Verfahren nach Anspruch 1 , wobei die Eingangsdaten für das Korrekturverfahren Wertepaare sind, die aus jeweils durchstrahlter Länge und zugehörige gemessener Intensität auf dem Detektor (31) aufgebaut sind.
8. Verfahren nach Anspruch 7, wobei zur Bestimmung der Wertepaare Projektionen aus einem Winkelbereich von weniger als 360°, insbesondere deutlich weniger als 360° verwendet werden.
9. Verfahren nach Anspruch 1 , wobei die Merkmalspunkte jeweils ein singuläres Punktepaar aus Modellpunkt und zugehörigem Punkt der Projektion sind.
10. Mehrstufiges Verfahren zur Bereitstellung von korrigierten Projektionsdaten als eine verbesserte CT-Rekonstruktion, bei der von einer Strahlenquelle (Q) insbesondere fächerförmige Messstrahlen (q) ausgesendet werden, die den Prüfling (10,11) durchstrahlen und auf einem Detektor (31) in ihrer Intensität erfasst werden, mit folgenden Schritten des Verfahrens,
(a) einer ersten Initialisierung, wobei eine erste Orientierung des Prüflings (10) mit einer ersten schnellen Registrierung nicht mehr als grobbestimmt wird;
(b) Registrierung, wobei die Positionierung des Prüflings (10) genauer bestimmt wird, insbesondere über Merkmalspunktepaare;
(c) Bewegung, wobei nach einer erfolgreich durchgeführten Registrierung an zumindest einigen Projektionen die Lage des Prüflings (10, 11 ) für zumindest eine weitere Projektion berechnet wird;
(d) Simulation, wobei mit Hilfe der Ergebnisse der vorigen Schritte eine virtuelle CT durchgeführt wird, welche Eingangsdaten für ein zu erfolgendes Korrekturverfahren der CT-Rekonstruktion liefert; (e) Vornehmen einer Korrektur, wobei während der Datenaufnahme (70) vom Detektor (31) aus den Korrekturdaten Parameter bestimmt werden und eine Korrektur erfolgt (73a,73b);
(f) Rekonstruktion, wobei im zeitlichen Bereich des Endes des Aufnahmeprozesses als Datenaufnahme (70) korrigierte Projektionsdaten (11*) als eine verbesserte CT-Rekonstruktion vorliegen (74,75).
11. Verfahren nach Anspruch 10, wobei Nichtlinearitäten des Abbildungssystems aus Quelle und Detektor (Q;31) mit dazwischen gelegenem Prüfling korrigiert werden.
12. Verfahren nach Anspruch 10, wobei die Bereitstellung bei der industriellen Qualitätskontrolle erfolgt.
13. Verfahren nach Anspruch 12, wobei zumindest eine Vermessung am Prüfling erfolgt, im Wege einer Durchstrahlungsaufnahme.
14. Verfahren nach Anspruch 10, wobei Röntgenstrahlen als Messstrahlen bei der Tomogramm Erstellung als Rekonstruktion des Prüflings verwendet werden.
15. Verfahren nach Anspruch 13 oder 10, wobei der Prüfling ein Gussteil ist, insbesondere im Automobilbau.
16. Verfahren nach Anspruch 10, wobei die Initialisierung mit einem Winkelfehler von wenigen Grad, insbesondere oberhalb von 1 ° oder/und mit einem translatorischen Fehler oberhalb von im wesentlichen 1 mm erfolgt.
17. Verfahren nach Anspruch 16, wobei die Drehachse (100) des Prüflings gegeben ist um welche der Prüfling in einzelnen indexierten Schritten eines gegebenen Winkelinkrements (Δα) bei der Durchstrahlung gedreht wird.
18. Verfahren nach Anspruch 10, wobei die Merkmalspunkte extrahiert werden und dabei ein jeweiligen singulärer Punkt aus einen digitalen Modell, insbesondere CAD-Modell, auf dem Detektor (31) als ein jeweils abgebildeter Punkt erscheint, und beide korrespondierenden Punkte ein Merkmalspunktepaar bilden.
19. Verfahren nach Anspruch 10 oder 16, wobei die Initialisierung mit einem translatorischen Fehler im Bereich von im Wesentlichen 1 % einer typischen Abmessung des Prüflings 10 erfolgt.
20. Mehrstufiges Verfahren zur Bereitstellung von korrigierten Projektionsdaten als eine verbesserte CT-Rekonstruktion, mit folgenden Schritten
• Initialisierung, wobei eine Orientierung des Prüflings (10,11) mit einer ersten schnellen Registrierung grob bestimmt wird;
• Registrierung, wobei die Position des Prüflings (10) genauer insbesondere hinreichend genau bestimmt wird;
• Bewegung, wobei nach einer erfolgreich durchgeführten Registrierung an zumindest einigen Projektionen die Lage des Prüflings für zumindest eine weitere Projektion berechnet wird;
• Simulation, wobei mit Hilfe dieses Wissens eine virtuelle CT simuliert wird, um Eingangsdaten für ein Korrekturverfahren der CT-Rekonstruktion zu liefern;
• Korrektur, wobei aus den Daten während der Datenaufnahme Parameter zur Korrektur bestimmt werden und eine Korrektur erfolgt;
• Rekonstruktion, wobei im Bereich des zeitlichen Endes des Aufnahmeprozesses korrigierte Projektionsdaten für eine verbesserte CT-Rekonstruktion vorliegen.
21. Verfahren nach Anspruch 20, wobei die Position des Prüflings über Merkmalspunktepaare erfolgt.
22. Verfahren nach Anspruch 20 oder 21 , wobei die Position des Prüflings über eine intensitätsgestützte statistische Methode erfolgt.
* A ¥ ♦
PCT/DE2006/000420 2005-03-09 2006-03-09 Korrektur von nichtlinearitäten eines abbildungssytems durch ein a-priori wissen bei durchstrahlungsaufnahmen WO2006094493A2 (de)

Priority Applications (4)

Application Number Priority Date Filing Date Title
US11/908,267 US20080212734A1 (en) 2005-03-09 2006-03-09 Correction of Non-Linearities in an Imaging System by Means of a Priori Knowledge in Radiography
CA002600648A CA2600648A1 (en) 2005-03-09 2006-03-09 Correction of non-linearities in an imaging system by means of a priori knowledge in radiography
EP06722577A EP1861734A2 (de) 2005-03-09 2006-03-09 Korrektur von nichtlinearitäten eines abbildungssytems durch ein a-priori wissen bei durchstrahlungsaufnahmen
DE112006001354T DE112006001354A5 (de) 2005-03-09 2006-03-09 Korrektur von Nichtlinearitäten eines Abbildungssystems durch ein a-priori Wissen bei Durchstrahlungsaufnahmen

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
DE102005011161 2005-03-09
DE102005011161.0 2005-03-09

Publications (3)

Publication Number Publication Date
WO2006094493A2 true WO2006094493A2 (de) 2006-09-14
WO2006094493A3 WO2006094493A3 (de) 2006-11-16
WO2006094493A8 WO2006094493A8 (de) 2010-08-05

Family

ID=36649618

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/DE2006/000420 WO2006094493A2 (de) 2005-03-09 2006-03-09 Korrektur von nichtlinearitäten eines abbildungssytems durch ein a-priori wissen bei durchstrahlungsaufnahmen

Country Status (5)

Country Link
US (1) US20080212734A1 (de)
EP (1) EP1861734A2 (de)
CA (1) CA2600648A1 (de)
DE (1) DE112006001354A5 (de)
WO (1) WO2006094493A2 (de)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010094774A2 (de) 2009-02-20 2010-08-26 Werth Messtechnik Gmbh Verfahren zum messen eines objekts
DE102009038505A1 (de) 2009-08-21 2011-03-03 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Konzept zum Ermitteln von Korrekturwerten für eine Reduktion von Artefakten bei einer Computertomographie eines zu untersuchenden Teilbereichs eines Prüfkörpers
AT509965A1 (de) * 2010-05-25 2011-12-15 Fh Ooe Forschungs & Entwicklungs Gmbh Verfahren zur bestimmung der optimalen lage eines objektes für eine 3d-ct-durchstrahlung
DE102013104720A1 (de) 2012-05-07 2013-11-07 Werth Messtechnik Gmbh Computertomografie-Verfahren und Anordnung zur Bestimmung von Merkmalen an einem Messobjekt
EP3798684A1 (de) * 2019-09-26 2021-03-31 Siemens Healthcare GmbH Datenkorrektur in der röntgenbildgebung
CN112763517A (zh) * 2020-12-24 2021-05-07 中国原子能科学研究院 一种堆内辐照试验装置

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102007016370A1 (de) * 2007-04-03 2008-10-09 Carl Zeiss Industrielle Messtechnik Gmbh Verfahren und eine Messanordnung zum Erzeugen von dreidimensionalen Bildern von Messobjekten mittels invasiver Strahlung
JP6131606B2 (ja) * 2013-01-21 2017-05-24 株式会社島津製作所 放射線撮影装置およびそれにおける画像処理方法
CN107533018B (zh) * 2015-04-24 2021-06-04 株式会社尼康 X射线检查装置、x射线检查方法及构造物的制造方法
EP3537137A1 (de) * 2016-11-01 2019-09-11 Shimadzu Corporation Verfahren zur kalibrierung einer bildgebungsvergrösserung für eine strahlungstomografievorrichtung
DE102017205113A1 (de) * 2017-03-27 2018-09-27 Siemens Aktiengesellschaft Ermitteln der Pose einer Röntgeneinheit relativ zu einem Objekt anhand eines digitalen Modells des Objekts
CN112964738B (zh) * 2021-01-29 2022-11-22 山东大学 一种工业ct快速扫描系统及方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4302675A (en) * 1980-01-21 1981-11-24 Technicare Corporation Method of multiplanar emission tomography and apparatus therefor
US4969110A (en) * 1988-08-01 1990-11-06 General Electric Company Method of using a priori information in computerized tomography
US6047041A (en) * 1997-09-08 2000-04-04 Scientific Measurement System Apparatus and method for comparison
US20040066908A1 (en) * 2000-10-11 2004-04-08 Randolf Hanke Method and device for the representing an object by means of an irradiation and for reconstructing said object

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4920491A (en) * 1988-05-16 1990-04-24 General Electric Company Enhancement of image quality by utilization of a priori information
US5390112A (en) * 1993-10-04 1995-02-14 General Electric Company Three-dimensional computerized tomography scanning method and system for imaging large objects with smaller area detectors
FR2820822B1 (fr) * 2001-02-14 2003-09-05 Peugeot Citroen Automobiles Sa Dispositif et procede de manipulation d'un produit et de traitement d'images radiocospiques du produit pour obtenir des coupes tomographiques et utilisations
DE10202732A1 (de) * 2002-01-24 2003-08-07 Fraunhofer Ges Forschung Vorrichtung und Verfahren zum Erstellen einer Korrekturkennlinie für eine Reduktion von Artefakten bei einer Tomographie

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4302675A (en) * 1980-01-21 1981-11-24 Technicare Corporation Method of multiplanar emission tomography and apparatus therefor
US4969110A (en) * 1988-08-01 1990-11-06 General Electric Company Method of using a priori information in computerized tomography
US6047041A (en) * 1997-09-08 2000-04-04 Scientific Measurement System Apparatus and method for comparison
US20040066908A1 (en) * 2000-10-11 2004-04-08 Randolf Hanke Method and device for the representing an object by means of an irradiation and for reconstructing said object

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
GRAEME P PENNEY ET AL: "A Comparison of Similarity Measures for Use in 2-D-3-D Medical Image Registration" IEEE TRANSACTIONS ON MEDICAL IMAGING, IEEE SERVICE CENTER, PISCATAWAY, NJ, US, Bd. 17, Nr. 4, August 1998 (1998-08), XP011035760 ISSN: 0278-0062 *
KOENIG A ET AL: "OBJECT POSE ESTIMATION USING A SET OF LOCAL RADIOGRAPHS OF A PART AND ITS CAD MODEL" REVIEW OF PROGRESS IN QUANTITATIVE NONDESTRUCTIVE EVALUATION, Nr. 16, 28. Juli 1996 (1996-07-28), Seiten 829-836, XP001069375 *
See also references of EP1861734A2 *

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2665035A2 (de) 2009-02-20 2013-11-20 Werth Messtechnik GmbH Verfahren zum Messen eines Objekts
DE102010000473A1 (de) 2009-02-20 2010-08-26 Werth Messtechnik Gmbh Verfahren zum Messen eines Objektes
WO2010094774A2 (de) 2009-02-20 2010-08-26 Werth Messtechnik Gmbh Verfahren zum messen eines objekts
US9025855B1 (en) * 2009-02-20 2015-05-05 Werth Messtechnik Gmbh Method for measuring an object
EP2654016A2 (de) 2009-02-20 2013-10-23 Werth Messtechnik GmbH Verfahren zum Messen eines Objekts
EP2654017A2 (de) 2009-02-20 2013-10-23 Werth Messtechnik GmbH Verfahren zum Messen eines Objekts
EP2665034A2 (de) 2009-02-20 2013-11-20 Werth Messtechnik GmbH Verfahren zum Messen eines Objekts
DE102009038505A1 (de) 2009-08-21 2011-03-03 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Konzept zum Ermitteln von Korrekturwerten für eine Reduktion von Artefakten bei einer Computertomographie eines zu untersuchenden Teilbereichs eines Prüfkörpers
DE102009038505A8 (de) * 2009-08-21 2011-06-01 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Konzept zum Ermitteln von Korrekturwerten für eine Reduktion von Artefakten bei einer Computertomographie eines zu untersuchenden Teilbereichs eines Prüfkörpers
AT509965A1 (de) * 2010-05-25 2011-12-15 Fh Ooe Forschungs & Entwicklungs Gmbh Verfahren zur bestimmung der optimalen lage eines objektes für eine 3d-ct-durchstrahlung
AT509965B1 (de) * 2010-05-25 2012-06-15 Fh Ooe Forschungs & Entwicklungs Gmbh Verfahren zur bestimmung der optimalen lage eines objektes für eine 3d-ct-durchstrahlung
WO2013167616A2 (de) 2012-05-07 2013-11-14 Werth Messtechnik Gmbh Computertomografie-verfahren und anordnung zur bestimmung von merkmalen an einem messobjekt
DE102013104720A1 (de) 2012-05-07 2013-11-07 Werth Messtechnik Gmbh Computertomografie-Verfahren und Anordnung zur Bestimmung von Merkmalen an einem Messobjekt
EP3798684A1 (de) * 2019-09-26 2021-03-31 Siemens Healthcare GmbH Datenkorrektur in der röntgenbildgebung
US11666301B2 (en) 2019-09-26 2023-06-06 Siemens Healthcare Gmbh Data correction in X-ray imaging
CN112763517A (zh) * 2020-12-24 2021-05-07 中国原子能科学研究院 一种堆内辐照试验装置

Also Published As

Publication number Publication date
CA2600648A1 (en) 2006-09-14
DE112006001354A5 (de) 2008-03-06
US20080212734A1 (en) 2008-09-04
WO2006094493A3 (de) 2006-11-16
EP1861734A2 (de) 2007-12-05
WO2006094493A8 (de) 2010-08-05

Similar Documents

Publication Publication Date Title
WO2006094493A2 (de) Korrektur von nichtlinearitäten eines abbildungssytems durch ein a-priori wissen bei durchstrahlungsaufnahmen
DE102011004598B4 (de) Verfahren und Computersystem zur Streustrahlkorrektur in einem Multi-Source-CT
EP2847620B1 (de) Computertomografie-verfahren und anordnung zur bestimmung von merkmalen an einem messobjekt
EP3548876B1 (de) Dunkelfeldtensortomographieverfahren
EP3111417B1 (de) Rauschreduktion in tomogrammen
EP1415179B1 (de) Vorrichtung und verfahren zum erstellen einer korrekturkennlinie für eine reduktion von artefakten bei einer tomographie
DE102016221658A1 (de) Streustrahlungskompensation für ein bildgebendes medizinisches Gerät
DE102007016370A1 (de) Verfahren und eine Messanordnung zum Erzeugen von dreidimensionalen Bildern von Messobjekten mittels invasiver Strahlung
EP1861822B1 (de) Verfahren und vorrichtung zur konturfeinermittlung eines objekts bei bildgebenden untersuchungsverfahren
DE102006022104B4 (de) Vorrichtung zur dreidimensionalen Vermessung eines Festkörpers
WO2006122756A1 (de) Verfahren und vorrichtung zum bestimmen von materialgrenzen eines prüfobjekts
EP1899714B1 (de) Verfahren und anordnung zum untersuchen eines messobjekts mittels invasiver strahlung
DE102006022103B4 (de) Verfahren zum Vermessen eines Festkörpers
DE102006011008A1 (de) Korrektur von Nichtlinearitäten eines Abbildungssystems durch ein a-priori-Wissen bei Durchstrahlungsaufnahmen
DE102015007934A1 (de) Verfahren und Computerprogrammprodukt zum Erzeugen eines artefaktreduzierten Voxeldatensatzes
DE102015218596A1 (de) Verbesserung der Bildqualität bei der Computertomographie unter Ausnutzung redundanter Information in Projektionsdatensätzen
DE102021204628B3 (de) Verfahren zum Betreiben eines Computertomographen beim Vermessen einer Interessensregion eines Objekts und Computertomograph
WO2011045351A1 (de) Verfahren zur zerstörungsfreien 3d-untersuchung eines prüfkörpers mittels computertomographie
EP3992620A1 (de) Computerimplementiertes verfahren zur ermittlung mindestens eines für eine auswertung von messdaten benötigten geometrischen parameters
DE112020003891T5 (de) Röntgen-bildgebungsvorrichtung und röntgenbehandlungsvorrichtung
DE69932037T2 (de) Grössenbestimmung eines objektdetails
DE102010022285A1 (de) Verfahren zur Bestimmung der optimalen Lage eines Messobjektes bei der 3D-Computertomografie
DE102022103888A1 (de) Verfahren und Vorrichtung zur Computertomografiemessung
DE102009038505A1 (de) Konzept zum Ermitteln von Korrekturwerten für eine Reduktion von Artefakten bei einer Computertomographie eines zu untersuchenden Teilbereichs eines Prüfkörpers
DE102018132011A1 (de) Verfahren zur Korrektur von mit einem Detektor eines Computertomographen aufgenommenen Projektionsdaten zur Abbildung eines Objekts

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application
DPE1 Request for preliminary examination filed after expiration of 19th month from priority date (pct application filed from 20040101)
WWE Wipo information: entry into national phase

Ref document number: 2600648

Country of ref document: CA

WWE Wipo information: entry into national phase

Ref document number: 11908267

Country of ref document: US

WWE Wipo information: entry into national phase

Ref document number: 2006722577

Country of ref document: EP

NENP Non-entry into the national phase

Ref country code: RU

WWW Wipo information: withdrawn in national office

Country of ref document: RU

WWE Wipo information: entry into national phase

Ref document number: 1120060013548

Country of ref document: DE

WWP Wipo information: published in national office

Ref document number: 2006722577

Country of ref document: EP

REF Corresponds to

Ref document number: 112006001354

Country of ref document: DE

Date of ref document: 20080306

Kind code of ref document: P

DPE1 Request for preliminary examination filed after expiration of 19th month from priority date (pct application filed from 20040101)