EP1563460A1 - Verfahren zur bildregistrierung - Google Patents
Verfahren zur bildregistrierungInfo
- Publication number
- EP1563460A1 EP1563460A1 EP03778261A EP03778261A EP1563460A1 EP 1563460 A1 EP1563460 A1 EP 1563460A1 EP 03778261 A EP03778261 A EP 03778261A EP 03778261 A EP03778261 A EP 03778261A EP 1563460 A1 EP1563460 A1 EP 1563460A1
- Authority
- EP
- European Patent Office
- Prior art keywords
- criterion
- displacement
- distance
- control points
- smoothness
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Ceased
Links
- 238000000034 method Methods 0.000 title claims abstract description 57
- 238000006073 displacement reaction Methods 0.000 claims abstract description 22
- 238000011156 evaluation Methods 0.000 claims abstract description 4
- 230000009466 transformation Effects 0.000 claims abstract description 3
- 238000013459 approach Methods 0.000 claims description 4
- 238000005452 bending Methods 0.000 claims description 2
- 230000015572 biosynthetic process Effects 0.000 claims description 2
- 238000011161 development Methods 0.000 claims description 2
- 230000000737 periodic effect Effects 0.000 claims description 2
- 230000002123 temporal effect Effects 0.000 claims description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000004613 tight binding model Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 206010028980 Neoplasm Diseases 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000011478 gradient descent method Methods 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/14—Transformations for image registration, e.g. adjusting or mapping for alignment of images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/33—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
Definitions
- the invention relates to a method for image registration, ie for correcting geometric differences in different representations of an object.
- These processes play e.g. play an important role in medical technology and especially in the analysis of tissue changes in the context of early cancer detection.
- This object is achieved according to the invention by iteratively determining a transformation which is optimal with regard to a predetermined distance and smoothness criterion, in which corresponding control points are guaranteed to be mapped to one another by (1.) initializing an iteration counter and the initial displacement field, (2.) Determination of the numerical solutions of the non-linear partial differential equation (PDE) with the differential operator derived from a given smoothness criterion and the point evaluation functionals located at given control points, (3.) summarizing the interpolation conditions, (4.) calculating a special numerical solution of the PDE with the force determined on the basis of the distance criterion and the current displacement field and the differential operator derived from the smoothness criterion, (5.) evaluating the special solution at the control points, (6.) determining the coefficients to calculate an updated displacement, (7.) update the displacement field and increase the iteration counter, (8.) check the displacement for convergence and (9.) if the convergence criterion is not met, repeat steps (4.) to (8.
- reference image reference image
- template T template
- the reference and template can be in discrete form.
- T u T u (x)
- T u T (x -u (x)
- a minimizer of the abovementioned distance criterion can be determined in an iterative manner by means of a gradient descent method.
- any distance criterion can be selected.
- the forces associated with the common distance criteria can be found in the literature (Modersitzki 2002). However, the specific way in which these forces are calculated is not essential for the registration process.
- any functional known from the literature can be used as a smoothness criterion.
- a partial differential operator A can be derived from the smoothness criterion. These operators are known for the criteria used in the literature (Modersitzki 2002).
- the desired shift u can then be characterized as a solution of a non-linear partial differential equation (PDE).
- PDE non-linear partial differential equation
- This procedure coincides with the procedure for the method based solely on the distance criterion and the smoothness criterion.
- the specific numerical method for the solution of the PDE is irrelevant for the registration procedure.
- the v J , j ⁇ ⁇ , ..., m are Green's functions of the differential operator A, which represent a solution of the PDE for a given single point shift.
- a suitable linear combination of these Green functions therefore ensures that all control points are mapped onto one another as required in the overall process.
- the function v ° is determined using an iterative process so that the distance criterion is minimized while maintaining the required smoothness.
- the weighting factors / L * are adjusted so that the control points are mapped in the required manner.
- the initialization according to the invention is followed by a common iteration procedure, in the course of which a gradient descent is carried out taking into account the control points. Human intervention is not necessary.
- the described method thus combines the advantages of the methods based on distance criteria (in particular the ability to be automated and an optimal registration on average) with those of the checkpoint procedure (guaranteed registration of excellent points) and, when specifying an initial set of checkpoints, provides reproducible, optimal results regardless of the user or computer program.
- the details of the computer code do not play a significant role in the final result of the image registration and only influence the required computing time and the storage requirements.
- the images to be registered can be digital images, pixels, JPEG, wavelet-based objects or acoustic signals.
- linear equation systems occurring in the method can be solved directly, indirectly, iteratively or by means of multigrid and a reference coordinate system can be used for the method, which is represented by Euler or Lagrange coordinates.
- the invention proposes to register one-, two- or three-dimensional objects as well as sequences of one-, two- and three-dimensional objects and to use control points which are anatomical landmarks, fiducial markers or other characteristic parameters.
- One proposed distance criterion is based on intensity, edge, corner, surface normal or level set or on the sum of squared differences, I 2 distance, correlation, variants of the correlation, Mutual information or variants of the mutual information is based.
- the force terms associated with the distance measure should be calculated using finite difference methods or gradient formation and the smoothness criterion used should be physically motivated via an elastic potential or a fluid approach based on temporal or spatial derivations of the displacement via diffusive or curvature approaches.
- boundary conditions of the differential operator should advantageously be given via explicit or implicit, Neumann, Dirichlet, sliding, bending or periodic boundary conditions.
- the type of discretization of the differential operator should be based on finite differences, finite volumes, finite elements, Fourier methods, series development, filter techniques, collocations or multigrid, and the interpolation should be carried out one-dimensionally using splines or wavelets.
- the move can be explicitly updated using the increment of the move or its time derivative.
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Image Processing (AREA)
Abstract
Verfahren zur Bildregistrierung durch iterative Bestimmung einer hinsichtlich eines vorgegebenen Distanz- und Glattheitskriteriums optimalen Transformation, bei dem in den Bildern korrespondierende Kontrollpunkte garantierbar aufeinander abgebildet werden, durch (1) Initialisieren eines Iterationszählers und des initialen Verrückungsfeldes, (2) Bestimmen der numerischen Lösungen der nicht-linearen partiellen Differentialgleichung (PDE) mit dem aus einem vorgegebenen Glattheitskriterium ableitbarem Differentialope rator und den an vorgegebenen Kontrollpunkten lokalisierten Punktauswertefunktionalen, (3) Zusammenfassen der Interpolationsbedingungen, (4) Berechnen einer speziellen numerischen Lösung der PDE mit der auf der Basis des Distanzkriteriums und des aktuellen Verrückungsfeldes bestimmten Kraft und dem aus dem Glattheitskriterium abgeleiteten Differentialoperator, (5) Auswerten der speziellen Lösung an den Kontrollpunkten, (6) Bestimmen der Koeffizienten zur Berechnung einer aktualisierten Verrückung, (7) Aufdatieren des Verrückungsfelds und Erhöhen des Iterationszählers, (8) Überprüfen der Verrückung auf Konvergenz und (9) bei Nichterfüllen des Konvergenzkriteriums erneutes Durchlauf der Schritte (4) bis (8).
Description
Verfahren zur Bildregistrierung
Die Erfindung betrifft ein Verfahren zur Bildregistrierung, also zur Korrektur von geometrischen Differenzen in unterschiedlichen Darstellungen eines Objektes. Diese Verfahren spielen z.B. in der Medizintechnik und insbesondere bei der Analyse von Gewebeveränderungen im Rahmen der Krebsfrüherkennung eine wichtige Rolle.
Es sind bereits Verfahren bekannt, die eine Bildregistrierung aufgrund eines Distanzkriteriums durchführen (Lisa Gottesfeld Brown: A survey of Image registration techniques, ACM Computing Surveys, 24(4): 325-376, 1992, Jan Modersitzki: Numerical Methods or Image Registration, Habilitation, Institute of Mathematics; University of Lübeck, Germany, 2002). Die allgemeine Methodik basiert auf der Optimierung einer anwendungskonform zu wählenden Zielfunktion, die typischerweise auf Bildintensitäten basiert. Bei derartigen Verfahren werden außer der Bildinformation keine weiteren Kenntnisse zur Registrierung herangezogen. Das Ergebnis der Registrierung ist nur im Sinne einer globalen Mittelung optimal. Kommen in einer Anwendung speziellen, charakteristischen Punkten (wie z.B. in medizintechnischen Anwendungen den sogenannten anatomischen Landmarken) eine besondere Bedeutung zu, sind diese Verfahren also nicht empfehlenswert.
Neben der Bildregistrierung auf der Basis eines Distanzkriteriums sind auch Verfahren bekannt, die die Bildregistrierung ausschließlich auf der Basis von Kontrollpunkten durchführen (Karl Rohr: Landmark-based Image Analysis. Computational Imaging and Vision. Kluwer Academic Publishers, Dordrecht, 2001). Bei diesen Verfahren werden den zu registrierenden Ansichten prospektiv oder retrospektiv korrespondierende Kontrollpunkte assoziiert, die dann mittels Registrierung zum Abgleich gebracht werden. Der Nachteil dieser Verfahren besteht darin, dass die Registrierung ausschließlich Kontrollpunkte berücksichtigt. Weitere Bildinformationen wie z.B. Bildintensitäten können bei diesen Verfahren nicht berücksichtigt werden. Bei unbefriedigenden Registrierungsergebnissen kann ein Anwendet nur versuchen, das Ergebnis durch geschicktes Einfügen weiterer Kontrollpunkte zu verbessern. Das Einfügen weiterer Kontrollpunkte basiert auf subjektivem Ausprobieren. Hierzu gibt es keine Richtlinien und insbesondere kein automatisiertes Vorgehen.
Der Erfindung liegt die Aufgabe zugrunde, ein Verfahren zur Bildregistrierung zu entwickeln, das sowohl zu einer perfekten, garantierbaren Korrespondenz einer Anzahl vorgegebener Kontrollpunkte führt, als auch zu einem im Sinne des Distanzkriterium optimalen Ergebnis.
Diese Aufgabe wird erfindungsgemäß gelöst durch iterative Bestimmung einer hinsichtlich eines vorgegebenen Distanz- und Glattheitskriteriums optimalen Transformation, bei dem in den Bildern korrespondierende Kontrollpunkte garantierbar aufeinander abgebildet werden, durch (1.) Initialisieren eines Iterationszählers und des initialen Verrük- kungsfeldes, (2.) Bestimmen der numerischen Lösungen der nicht-linearen partiellen Differentialgleichung (PDE) mit dem aus einem vorgegebenen Glattheitskriterium ableitbarem Differentialoperator und den an vorgegebenen Kontrollpunkten lokalisierten Punkt- auswertefunktionalen, (3.) Zusammenfassen der Interpolationsbedingungen, (4.) Berechnen einer speziellen numerischen Lösung der PDE mit der auf der Basis des Distanzkriteriums und des aktuellen Verruckungsfeldes bestimmten Kraft und dem aus dem Glattheitskriterium abgeleiteten Differentialoperator, (5.) Auswerten der speziellen Lösung an den Kontrollpunkten, (6.) Bestimmen der Koeffizienten zur Berechnung einer aktualisierten Verrückung, (7.) Aufdatieren des Verruckungsfelds und Erhöhen des Iterationszählers, (8.) Überprüfen der Verrückung auf Konvergenz und (9.) bei Nichterfüllen des Konvergenzkriteriums erneutes Durchlauf der Schritte (4.) bis (8.).
Der Verfahrensablauf wird durch das Flußdiagramm der Figur 1 beschreibend abgebildet.
Zur Vereinfachung bezeichnen wir eine Ansicht als Referenzbild (Referenz R) und eine weitere, zu korrigierende Ansicht als Template (Template T). Formal handelt es sich hierbei um Funktionen eines ύf-dimensionalen reellen Raumes bzw. einer Teilmenge Ω a Rrf in die Menge der reellen Zahlen. Jedem -dimensionalen Punkt x e Ω wird also durch i?(x)und T(x) ein Wert zugewiesen, der z.B. als Färb- oder als Grauwert interpretiert werden kann.
In praktischen Anwendungen - insbesondere bei jeder Programmierung des hier erläuterten Verfahrens - können Referenz und Template in diskreter Form vorliegen. Die Bilder sind dann Funktionen auf einem Gitter (z.B. Ω = {l,...,«, }x {l,...,«2 }} ^a die Dimension d = 2) in eine diskrete Menge (z.B. in die Menge {θ,...,255}) und können als aus Pixeln
aufgebaut interpretiert werden. Für das Registrierungsverfahren sind diese Einschränkungen und insbesondere die konkrete Art der Diskretisierung unerheblich und unwesentlich. Die Einschränkungen werden lediglich zum Zwecke einer vereinfachten Beschreibung vorgenommen. Das Verfahren lässt sich analog auf beliebige -dimensionale Datensätze anwenden.
Die Aufgabe der Bildregistrierung besteht in der Bestimmung einer Verrückungsfunktion u, so dass die Forderung R(x) = Tu (x) mit der Abkürzung Tu (x) := T(x -u(x)) für alle x e Ω möglichst gut erfüllt wird. Zur Berechnung des durch u deformierten Templates Tu wird bei diskret vorgegebenen Bildern wie in der Bildverarbeitung allgemein üblich eine Interpolation (z.B. -linear) durcrizuführen sein, da die verrückten Koordinaten x - u(x) nicht notwendig auch auf dem diskreten Gitter liegen müssen. Wie eine solche Interpolation erfolgt, ist für das Registrierungsverfahren nicht wesentlich.
Über die oben formulierte Ähnlichkeit hinaus müssen sowohl Forderungen an die Glattheit der Verrückung gestellt werden sowie an Abbildungseigenschaften in Bezug auf eine Anzahl von vorab gewählten Kontrollpunkten. Im einfachsten Fall sollen die Koordinaten jedes der m Kontrollpunkte KT des Templates auf den jeweils korrespondierenden Kontrollpunkt KR der Referenz abgebildet werden, j = \,...,m . Stimmen die Koordinaten der Kontrollpunkte bereits überein, was ggf. durch eine Vorregistrierung gewährleistet werden kann, so gilt also u = 0 in diesen Punkten.
Wie bei Optimierungsproblemen üblich kann die Bestimmung eines Minimierers des oben genannten Distanzkriteriums mittels eines Gradientenabstiegsverfahren in iterativer Weise erfolgen. Im Prinzip kann hierbei jedes beliebige Distanzkriterium gewählt werden. Die den gängigen Distanzkriterien assoziierten Kräfte finden sich in der Literatur (Modersitzki 2002). Die konkrete Art der Berechnung dieser Kräfte ist für das Registrierungsverfahren jedoch nicht wesentlich.
Als Glattheitskriterium kann wiederum prinzipiell jedes aus der Literatur bekannte Funktional genutzt werden. Aus dem Glattheitskriterium lässt sich ein partieller Differentialoperator A ableiten. Für die in der Literatur verwendeten Kriterien sind diese Operatoren bekannt (Modersitzki 2002). Die gesuchte Verrückung u lässt sich dann als Lösung einer nicht-linearen partiellen Differentialgleichung (PDE) charakterisieren.
Zur Bestimmung einer numerischen Lösung dieser PDE verwenden wir eine finite Differenzen Approximation des Differentialoperators, die dann auf ein Gleichungssystem für die Gitterwerte der Verrückung führt. Die konkrete Diskretisierung der Differentialgleichung ist jedoch für das Registrierungsverfahren ohne Bedeutung.
Diese Vorgehensweise deckt sich mit dem Verfahren für die ausschließlich auf dem Distanzkriterium und dem Glattheitskriterium basierende Methode. Der neue Aspekt besteht in einer geeigneten Einbindung der vorgegebenen Kontrollpunkte in die Berechnung der Verrückung, bei der eine Korrespondenz der Kontrollpunkte garantiert werden kann. Da bereits Verfahren zur Bestimmung der Verrückung auf der Basis des Distanz- und des Glattheitskriteriums bekannt sind, geben wir ein Verfahren an, welches Teillösungen in geeigneter Weise zu einer Gesamtlösung kombiniert, z.B. als m ue(x) = v°(x) + ∑λJ evj (x) , e Ω, £ = l,...,d .
Bezeichnet A den zum Glattheitsterm gehörigen Differentialoperator und / das zum Distanzkriterium gehörige Kraftfeld, dann ist v° eine numerische Lösung von Av° = -f , die Funktionen vJ sind numerische Lösungen der distributioneilen PDE AvJ = δJ , j = \,...,.m wobei δJ das Punktauswertefunktional (Dirac-Stoß) lokalisiert am Kontrollpunkt KT'J bezeichnet. Die konkrete numerische Methode für die Lösung der PDE ist für das Registrierungsverfahren unerheblich.
Aus mathematischer Sicht sind die vJ ,j ~ \,...,m , Greensche Funktionen des Differentialoperators A, die eine Lösung der PDE bei vorgegebener Einzelpunktverrückung darstellen. Eine geeignete Linearkombination dieser Greenschen Funktionen stellt daher sicher, dass bei dem Gesamtverfahren alle Kontrollpunkte wie gefordert aufeinander abgebildet werden.
Die Funktion v° wird über ein iteratives Verfahren so bestimmt, dass das Distanzkriterium unter Einhaltung der geforderten Glattheit minimiert wird. Die Gewichtsfaktoren /L* werden dabei so angepasst, dass die Kontrollpunkte in der geforderten Weise abgebildet werden.
Die Initialisierung des Programms erfordert die Wahl eines Distanz- und eines Glattheitskriteriums bzw. der aus diesen Kriterien ableitbaren Kraft und des Differentialoperators. Auf der Basis der an den Kontrollpunkten lokalisierten Punktauswertefunktio-
nale können dann die Greenschen Funktionen vJ , j = 1,..., m mit einem numerischen Verfahren bestimmt werden. Diese werden im weiteren Verlauf des Verfahrens nicht mehr geändert.
Der erfindungsgemäßen Initialisierung folgt eine gängige Iterationsprozedur, in deren Verlauf ein Gradientenabstieg unter Berücksichtigung der Kontrollpunkte durchgeführt wird. Menschliches Eingreifen ist nicht von Nöten. Das beschriebene Verfahren kombiniert also die Vorteile der auf Distanzkriterien basierenden Methoden (also insbesondere die Automatisierbarkeit und eine im Mittel optimale Registrierung) mit denen des Kontrollpunktverfahrens (garantierte Registrierung ausgezeichneter Punkte) und liefert bei Vorgabe eines initialen Satzes von Kontrollpunkten reproduzierbare, optimale Ergebnisse unabhängig vom Anwender bzw. Computerprogramm. Die Details des Computercodes spielen für das Endergebnis der Bildregistrierung keine wesentliche Rolle und beeinflussen nur die benötigte Rechenzeit und die Speicheranforderungen.
Bei den zu registrierenden Bildern kann es sich um digitale Bilder, Pixel, JPEG, wavelet basierte Objekte oder akustische Signale handeln.
Die in dem Verfahren auftretenden linearen Gleichungssysteme können direkt, indirekt, iterativ oder mittels Multigrid gelöst werden und es kann für das Verfahren ein Referenzkoordinatensystem verwendet werden, das durch Euler- bzw. Lagrange-Koordinaten abgebildet wird.
Im übrigen schlägt die Erfindung vor ein-, zwei- oder drei-dimensionale sowie Sequenzen von ein-, zwei- und drei-dimensionalen Objekten zu registrieren, und Kontrollpunkte zu verwenden, die anatomische Landmarken, fiducial marker oder andere charakteristische Kenngrößen sind.
Als Distanzkriterium wird ein solches vorgeschlagen, daß auf Intensitäts-, Kanten-, Ek- ken-, Oberflächennormalen- oder Level-Set basiert oder auf der ,Sum of Squared Diffe- rences', I2-Distanz, Korrelation, Varianten der Korrelation, Mutual-Information oder Varianten der Mutual-Information beruht.
Die dem Distanzmaß assoziierten Kraftterme sollten über Finite-Differenzen- Verfahren oder Gradientenbildung berechnet werden und das verwendete Glattheitskriterium über ein elastisches Potential oder einen fluidalen Ansatz physikalisch motiviert sein oder
über Diffusive oder Curvature-Ansätze auf zeitlichen oder räumlichen Ableitungen der Verrückung basieren.
Die Randbedingungen des Differentialoperators sollten vorteilhafterweise über explizite oder implizite, Neumann-, Dirichlet-, sliding-, bending- oder periodische Randbedingungen gegeben sein.
Die Art der Diskretisierung des Differentialoperators sollte auf finiten Differenzen, fini- ten Volumen, finiten Elemente, Fourier-Methoden, Reihenentwicklung, Filtertechniken, Kollokationen oder Multigrid beruhen und die Interpolation -dimensional über Splines oder Wavelets durchgeführt werden.
Schließlich kann die Verrückung explizit über das Inkrement der Verrückung oder deren Zeitableitung aufdatiert werden.
Claims
1. Verfahren zur Registrierung von Bildern durch iterative Bestimmung einer hinsichtlich eines vorgegebenen Distanz- und Glattheitskriteriums optimalen Transformation, dadurch gekennzeichnet, dass
in den Bildern korrespondierende Kontrollpunkte garantierbar aufeinander abgebildet werden, durch
(1) initialisieren eines Iterationszählers und des initialen Verruckungsfeldes,
(2) Bestimmen der numerischen Lösungen der nicht-linearen partiellen Differentialgleichung (PDE) mit dem aus einem vorgegebenen Glattheitskriterium ableitbarem Differentialoperator und den an vorgegebenen Kontrollpunkten lokalisierten Punktaus- wertefunktionalen,
(3) Zusammenfassen der Interpolationsbedingungen,
(4) Berechnen einer speziellen numerischen Lösung der PDE mit der auf der Basis des Distanzkriteriums und des aktuellen Verruckungsfeldes bestimmten Kraft und dem aus dem Glattheitskriterium abgeleiteten Differentialoperator,
(5) Auswerten der speziellen Lösung an den Kontrollpunkten,
(6) Bestimmen der Koeffizienten zur Berechnung einer aktualisierten Verrückung,
(7) Aufdatieren des Verruckungsfelds und Erhöhen des Iterationszählers,
(8) Überprüfen der Verrückung auf Konvergenz und
(9) bei Nichterfüllen des Konvergenzkriteriums erneuter Durchlauf der Schritte (4) bis (8).
2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass ein- zwei- oder dreidimensionale sowie Sequenzen von ein-, zwei- und drei-dimensionalen Objekten registriert werden.
3. Verfahren nach einem der vorhergehenden Ansprüche, dadurch gekennzeichnet, dass die Kontrollpunkte anatomische Landmarken, fiducial marker oder andere charakteristische Kenngrößen sind.
4. Verfahren nach einem der vorhergehenden Ansprüche, dadurch gekennzeichnet, dass das Distanzkriterium Intensitäts-, Kanten-, Ecken-, Oberflächennormalen- oder Level-Set basiert oder auf der Sum of Squared Differences, Z2-Distanz, Korrelation, Varianten der Korrelation, Mutual-Information oder Varianten der Mutual-Information beruht.
5. Verfahren nach einem der vorhergehenden Ansprüche, dadurch gekennzeichnet, dass die dem Distanzmaß assoziierten Kraftterme über Finite-Differenzen- Verfahren oder Gradientenbildung berechnet werden.
6. Verfahren nach einem der vorhergehenden Ansprüche, dadurch gekennzeichnet, dass das verwendete Glattheitskriterium über ein elastisches Potential oder einen fluidalen Ansatz physikalisch motiviert oder über Diffusive oder Curvature-Ansätze auf zeitlichen oder räumlichen Ableitungen der Verrückung basiert.
7. Verfahren nach einem der vorhergehenden Ansprüche, dadurch gekennzeichnet, dass die Randbedingungen des Differentialoperators explizite oder implizite, Neumann-, Di- richlet-, sliding-, bending- oder periodische Randbedingungen sind.
8. Verfahren nach einem der vorhergehenden Ansprüche, dadurch gekennzeichnet, dass die Art der Diskretisierung des Differentialoperators auf finiten Differenzen, finiten Volumen, finiten Elemente, Fourier-Methoden, Reihenentwicklung, Filtertechniken, Kollokationen oder Multigrid beruht.
9. Verfahren nach einem der vorhergehenden Ansprüche, dadurch gekennzeichnet, dass die Interpolation -dimensional über Splines oder Wavelets durchgeführt wird.
10. Verfahren nach einem der vorhergehenden Ansprüche, dadurch gekennzeichnet, dass die Verrückung explizit über das Inkrement der Verrückung oder deren Zeitableitung aufdatiert wird.
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE10253784 | 2002-11-19 | ||
DE10253784A DE10253784A1 (de) | 2002-11-19 | 2002-11-19 | Verfahren zur Bildregistrierung |
PCT/DE2003/003805 WO2004047024A1 (de) | 2002-11-19 | 2003-11-18 | Verfahren zur bildregistrierung |
Publications (1)
Publication Number | Publication Date |
---|---|
EP1563460A1 true EP1563460A1 (de) | 2005-08-17 |
Family
ID=32318525
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
EP03778261A Ceased EP1563460A1 (de) | 2002-11-19 | 2003-11-18 | Verfahren zur bildregistrierung |
Country Status (5)
Country | Link |
---|---|
US (1) | US20060008179A1 (de) |
EP (1) | EP1563460A1 (de) |
AU (1) | AU2003285277A1 (de) |
DE (1) | DE10253784A1 (de) |
WO (1) | WO2004047024A1 (de) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7639896B2 (en) | 2004-08-09 | 2009-12-29 | Carestream Health, Inc. | Multimodal image registration using compound mutual information |
WO2006095221A2 (de) * | 2005-03-11 | 2006-09-14 | Philips Intellectual Property & Standards Gmbh | Bildgebendes verfahren |
JP4398919B2 (ja) * | 2005-08-22 | 2010-01-13 | 株式会社東芝 | 画像マッチング装置、画像マッチング方法および画像マッチングプログラム |
US8064664B2 (en) * | 2006-10-18 | 2011-11-22 | Eigen, Inc. | Alignment method for registering medical images |
CN103533973B (zh) * | 2011-05-16 | 2016-12-28 | 泰尔茂株式会社 | 夹子及血液袋系统 |
US9939509B2 (en) * | 2014-01-28 | 2018-04-10 | Ohio State Innovation Foundation | Variable density incoherent spatiotemporal acquisition (VISTA) for highly accelerated magnetic resonance imaging |
WO2016091534A1 (en) | 2014-12-09 | 2016-06-16 | Asml Netherlands B.V. | Method and apparatus for image analysis |
US10607334B2 (en) | 2014-12-09 | 2020-03-31 | Asml Netherlands B.V. | Method and apparatus for image analysis |
CN107437253B (zh) * | 2017-08-07 | 2020-10-23 | 江西农业大学 | 一种基于灭点的条播作物行提取的预处理方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5850486A (en) * | 1996-04-29 | 1998-12-15 | The Mclean Hospital Corporation | Registration of image data |
US6226418B1 (en) * | 1997-11-07 | 2001-05-01 | Washington University | Rapid convolution based large deformation image matching via landmark and volume imagery |
US6633686B1 (en) * | 1998-11-05 | 2003-10-14 | Washington University | Method and apparatus for image registration using large deformation diffeomorphisms on a sphere |
WO2002056241A1 (en) * | 2001-01-12 | 2002-07-18 | University Of Florida | Computational algorithms for image registration |
-
2002
- 2002-11-19 DE DE10253784A patent/DE10253784A1/de not_active Withdrawn
-
2003
- 2003-11-18 WO PCT/DE2003/003805 patent/WO2004047024A1/de not_active Application Discontinuation
- 2003-11-18 EP EP03778261A patent/EP1563460A1/de not_active Ceased
- 2003-11-18 US US10/535,682 patent/US20060008179A1/en not_active Abandoned
- 2003-11-18 AU AU2003285277A patent/AU2003285277A1/en not_active Abandoned
Non-Patent Citations (1)
Title |
---|
See references of WO2004047024A1 * |
Also Published As
Publication number | Publication date |
---|---|
AU2003285277A1 (en) | 2004-06-15 |
WO2004047024A1 (de) | 2004-06-03 |
DE10253784A1 (de) | 2005-06-02 |
US20060008179A1 (en) | 2006-01-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP1191476B1 (de) | Vorrichtung und Verfahren zur Berechnung einer zwei Abbildungen verbindenden Transformation | |
CN104008538B (zh) | 基于单张图像超分辨率方法 | |
DE102008032006B4 (de) | Verfahren zur Steurung der Bildaufnahme bei einer Bildaufnahmeeinrichtung, sowie eine Bildaufnahmeeinrichtung | |
EP1958154B1 (de) | Verfahren zur topographischen darstellung von veränderungen in einem untersuchten gehirn | |
EP2284795A2 (de) | Quantitative Analyse, Visualisierung und Bewegungskorrektur in dynamischen Prozessen | |
DE10392259T5 (de) | Objektkorrespondenzidentifikation ohne volle Volumenregistration | |
WO2005059831A1 (en) | Elastic image registration | |
DE102007019514A1 (de) | Allgemeines Framework für die Bildsegmentierung unter Verwendung geordneter räumlicher Abhängigkeit | |
DE102005003605A1 (de) | Vorwissen, Niveaumengenrepräsentationen und visuelle Gruppierung | |
DE102010016034A1 (de) | Verfahren und Systeme zum Anzeigen quantitativer segmentaler Daten im 4D-Rendering | |
WO2004047024A1 (de) | Verfahren zur bildregistrierung | |
DE102006050364A1 (de) | Verfahren zum Detektieren und Verfolgen deformierbarer Objekte unter Verwendung eines adaptiven zeitvariierenden autoregressiven Modells | |
CN108701360B (zh) | 图像处理系统和方法 | |
CN110945564A (zh) | 基于混合上下文cnn模型的医学图像分割 | |
Attia et al. | Spatially aware melanoma segmentation using hybrid deep learning techniques | |
Rogelj et al. | Validation of a nonrigid registration algorithm for multimodal data | |
Boone et al. | ROOD-MRI: Benchmarking the robustness of deep learning segmentation models to out-of-distribution and corrupted data in MRI | |
CN100550056C (zh) | 在图像序列中进行运动校正的方法及系统 | |
DE19829170C2 (de) | Berechnung und Darstellung 2- oder 3-dimensionaler Deformationen räumlicher Strukturen | |
CN112435284A (zh) | 图像配准方法、电子设备及存储介质 | |
Boukellouz et al. | Evaluation of several similarity measures for deformable image registration using T1-weighted MR images of the brain | |
Escalante-Ramírez et al. | Motion estimation and segmentation in ct cardiac images using the hermite transform and active shape models | |
DE102009015116B4 (de) | Verfahren und Vorrichtung zur Registrierung von Bilddatensätzen und zur Reduktion von lagebedingten Grauwertschwankungen nebst zugehörigen Gegenständen | |
Safran et al. | Curvelet-based sampling for accurate and efficient multimodal image registration | |
CN114255265A (zh) | 单模态医学图像的配准方法、系统及计算机可读存储介质 |
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: 20050510 |
|
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 IT LI LU MC NL PT RO SE SI SK TR |
|
AX | Request for extension of the european patent |
Extension state: AL LT LV MK |
|
DAX | Request for extension of the european patent (deleted) | ||
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: THE APPLICATION HAS BEEN REFUSED |
|
18R | Application refused |
Effective date: 20060623 |