DE102010014211A1 - Auswertungsverfahren für eine Folge von Bilddatensätzen mit zweifacher elastischer Verzerrung der Bilddatensätze - Google Patents

Auswertungsverfahren für eine Folge von Bilddatensätzen mit zweifacher elastischer Verzerrung der Bilddatensätze Download PDF

Info

Publication number
DE102010014211A1
DE102010014211A1 DE102010014211A DE102010014211A DE102010014211A1 DE 102010014211 A1 DE102010014211 A1 DE 102010014211A1 DE 102010014211 A DE102010014211 A DE 102010014211A DE 102010014211 A DE102010014211 A DE 102010014211A DE 102010014211 A1 DE102010014211 A1 DE 102010014211A1
Authority
DE
Germany
Prior art keywords
image data
data sets
distorted
evaluation method
computing device
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Withdrawn
Application number
DE102010014211A
Other languages
English (en)
Inventor
Dr. Deimling Michael
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Siemens Healthcare GmbH
Original Assignee
Siemens AG
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 Siemens AG filed Critical Siemens AG
Priority to DE102010014211A priority Critical patent/DE102010014211A1/de
Priority to US13/080,096 priority patent/US8463009B2/en
Publication of DE102010014211A1 publication Critical patent/DE102010014211A1/de
Withdrawn legal-status Critical Current

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/503—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of the heart
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5217—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
    • A61B6/5264—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to motion
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0012—Biomedical image inspection
    • G06T7/0014—Biomedical image inspection using an image reference approach
    • G06T7/0016—Biomedical image inspection using an image reference approach involving temporal comparison
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/20—Analysis of motion
    • G06T7/246—Analysis of motion using feature-based methods, e.g. the tracking of corners or segments
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5229—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
    • A61B6/5235—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/54—Control of apparatus or devices for radiation diagnosis
    • A61B6/541—Control of apparatus or devices for radiation diagnosis involving acquisition triggered by a physiological signal
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/54—Control of the diagnostic device
    • A61B8/543—Control of the diagnostic device involving acquisition triggered by a physiological signal
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G06T2207/10088—Magnetic resonance imaging [MRI]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G06T2207/10108—Single photon emission computed tomography [SPECT]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20048—Transform domain processing
    • G06T2207/20056—Discrete and fast Fourier transform, [DFT, FFT]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20076—Probabilistic image processing
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30048—Heart; Cardiac

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Public Health (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Biomedical Technology (AREA)
  • Pathology (AREA)
  • Surgery (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Optics & Photonics (AREA)
  • Veterinary Medicine (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Quality & Reliability (AREA)
  • Multimedia (AREA)
  • Physiology (AREA)
  • Cardiology (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Databases & Information Systems (AREA)
  • Data Mining & Analysis (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

Es wird eine Folge von Bilddatensätzen (B) eines Bereichs (2) eines Lebewesens (3) erfasst, während sich der erfasste Bereich (2) des Lebewesens (3) atmungs- und herzschlagbedingt bewegt. In den erfassten Bilddatensätzen (B) wird eine erste Kontur (9) bestimmt, die sich atmungsbedingt bewegt. Die erfassten Bilddatensätze (B) werden derart elastisch zu einfach verzerrten Bilddatensätzen (B') verzerrt, dass die ersten Konturen (9) der einfach verzerrten Bilddatensätze (B') örtlich miteinander korrespondieren. In den erfassten Bilddatensätzen (B) oder in den einfach verzerrten Bilddatensätzen (B') wird eine zweite Kontur (11) bestimmt, die sich herzschlagbedingt bewegt. Die einfach verzerrten Bilddatensätze (B') werden derart elastisch zu zweifach verzerrten Bilddatensätzen (B'') verzerrt, dass die zweiten Konturen (11) der zweifach verzerrten Bilddatensätze (B'') örtlich miteinander korrespondieren. Örtlich miteinander korrespondierende Bereiche (12) der zweifach verzerrten Bilddatensätze (B'') werden bilddatensatzübergreifend ausgewertet.

Description

  • Die vorliegende Erfindung betrifft ein Auswertungsverfahren für eine Folge von Bilddatensätzen eines Bereichs eines Lebewesens,
    • – wobei die Folge von Bilddatensätzen erfasst wurde, während sich der erfasste Bereich des Lebewesens atmungs- und herzschlagbedingt bewegte.
  • Die vorliegende Erfindung betrifft weiterhin ein Computerprogramm, das Maschinencode umfasst, der von einer Recheneinrichtung unmittelbar ausführbar ist und dessen Ausführung durch die Recheneinrichtung bewirkt, dass die Recheneinrichtung ein derartiges Auswertungsverfahren ausführt.
  • Die vorliegende Erfindung betrifft weiterhin eine Recheneinrichtung, die derart ausgebildet ist, dass sie im Betrieb ein derartiges Auswertungsverfahren ausführt.
  • Eine wichtige diagnostische Anwendung der Herzbildgebung ist es, Areale im Myocard, die eine verminderte Funktionsfähigkeit aufweisen, sichtbar zu machen. Im Stand der Technik hat sich hierfür insbesondere die nuklearmedizinische Bildgebung mit der so genannten SPECT Methode (Single Photon Emission Computed Tomography) etabliert. Es ist weiterhin zukünftig beabsichtigt, mehr und mehr die Kernspintomographie für derartige diagnostische Anwendungen einzusetzen. Die Kernspintomographie bietet insbesondere die Vorteile einer erheblich besseren Ortsauflösung sowie die Lieferung zusätzlicher funktioneller Informationen, beispielsweise eine Prognose über die Erholung des regionalen Kontraktionsvermögens des Herzens, der Beweglichkeit oder des Entzündungszustands der Herzwand usw.. Ein weiterer Vorteil der Kernspintomographie liegt darin, nicht mit ionisierender Strahlung arbeiten zu müssen.
  • Für die Darstellung der Funktionsfähigkeit wird im Stand der Technik ein Kontrastmittel (beispielsweise ein Gadoliniummittel) appliziert und nach einer gewissen, im Minutenbereich liegenden Zeit mit einer entsprechenden Messsequenz das Gewebe untersucht. An den Stellen des Herzmuskels, an denen die Durchblutung des Herzmuskels stark eingeschränkt ist, wird das Kontrastmittel nur stark verzögert ausgeschwemmt. Das erkrankte Areal wird daher mit höherem Signal dargestellt als gesundes Gewebe. Die entsprechende Vorgehensweise ist im Stand der Technik als „delayed enhancement” bekannt.
  • Ein anderes, auf Kernspinverfahren basierendes Erkennungsverfahren wird SPAMM genannt (SPAtial Modulation of Magnetization). Bei diesem Verfahren wird dem Herzbereich ein geordnetes Gitter- oder Streifenmuster aufgeprägt. Eine verminderte Beweglichkeit des Myocards zeigt sich in diesem Fall in einer Deformation der Streifen- bzw. Gitterstruktur innerhalb des zeitaufgelösten Intervalls zwischen zwei Herzschlägen.
  • Die Auswertungen des Standes der Technik sind problematisch, weil der zu untersuchende Bereich des Lebewesens bewegt ist. Insbesondere sind in der Regel zwei überlagerte Bewegungen vorhanden, nämlich zum einen die durch die Bewegung des Zwerchfells induzierte Lungenbewegung, zum anderen der Herzschlag und die dadurch bedingten Bewegungen. Fast alle Aufnahmeverfahren des Standes der Technik basieren daher auf einer bewegungssynchronisierten Datenerfassung wie beispielsweise prospektives oder retrospektives Triggern und/oder Gaten. Eine Messzeitverkürzung durch eine Aufnahme unter Atemstillstand ist in manchen Fällen möglich, aber nicht in allen Fällen.
  • Die Aufgabe der vorliegenden Erfindung besteht darin, ein Auswertungsverfahren für eine Folge von Bilddatensätzen eines Bereichs eines Lebewesens zu schaffen, bei dem weder eine Synchronisierung mit der atmungsbedingten Bewegung noch eine Synchronisierung mit der herzschlagbedingten Bewegung erforderlich ist.
  • Die Aufgabe wird durch ein Auswertungsverfahren mit den Merkmalen des Anspruchs 1 gelöst. Vorteilhafte Ausgestaltungen des erfindungsgemäßen Auswertungsverfahrens sind Gegenstand der abhängigen Ansprüche 2 bis 5.
  • Erfindungsgemäß ist vorgesehen, ein Auswertungsverfahren der eingangs genannten Art dadurch auszugestalten,
    • – dass in den erfassten Bilddatensätzen eine erste Kontur bestimmt wird, die sich entsprechend der atmungsbedingten Bewegung des erfassten Bereichs des Lebewesens bewegt,
    • – dass die erfassten Bilddatensätze derart elastisch zu einfach verzerrten Bilddatensätzen verzerrt werden, dass die ersten Konturen der einfach verzerrten Bilddatensätze örtlich miteinander korrespondieren,
    • – dass in den erfassten Bilddatensätzen oder in den einfach verzerrten Bilddatensätzen eine zweite Kontur bestimmt wird, die sich entsprechend der herzschlagbedingten Bewegung des erfassten Bereichs des Lebewesens bewegt,
    • – dass die einfach verzerrten Bilddatensätze derart elastisch zu zweifach verzerrten Bilddatensätzen verzerrt werden, dass die zweiten Konturen der zweifach verzerrten Bilddatensätze örtlich miteinander korrespondieren, und
    • – dass örtlich miteinander korrespondierende Bereiche der zweifach verzerrten Bilddatensätze bilddatensatzübergreifend ausgewertet werden.
  • Die erste Kontur kann insbesondere die Lungenkontur des Lebewesens oder eine Rippenkontur des Lebewesens sein. Die zweite Kontur kann insbesondere die Herzkontur des Lebewesens sein.
  • Die Art der Bilddatensätze kann nach Bedarf bestimmt sein. Insbesondere können die Bilddatensätze alternativ zweidimensionale Bilddatensätze oder dreidimensionale Bilddatensätze sein. Dreidimensionale Bilddatensätze sind bevorzugt. Auch die Art der physikalischen Erfassung der Bilddatensätze kann nach Bedarf bestimmt sein. In Frage kommen insbesondere Röntgentechniken, Ultraschalltechniken und – bevorzugt – Magnetresonanztechniken.
  • Die einzelnen Bilddatensätze weisen in der Regel einen konstanten zeitlichen Versatz auf. Der zeitliche Versatz kann nach Bedarf bestimmt sein. Er sollte aber ungleich der Herzperiode und auch ungleich der Atmungsperiode sein. Vorzugsweise werden die Bilddatensätze mit einem Zeitversatz erfasst, der kleiner als der zeitliche Abstand zwischen zwei Herzschlägen des Lebewesens ist. Beispielsweise können die Bilddatensätze mittels einer multislice 2D true FISP Messsequenz erfasst werden. Eine derartige Messsequenz kann mit einer temporären Auflösung von 300 ms Bilddatensätze liefern, die ein FoV (FoV = Field of View) von ca. 300 mm aufweisen. Wird die Messsequenz über längere Zeit wiederholt, fällt die Folge von Bilddatensätzen mit einer Rate von ca. 3 Bilddatensätzen pro Sekunde an. Wird zusätzlich eine Vielkanalarrayspule oder eine andere parallele Aufnahmetechnik verwendet, verwendet, ist eine Beschleunigung um (beispielsweise) den Faktor 4 möglich, so dass in einem Herzzyklus etwa 10 bis 15 Bilddatensätze anfallen. Eine derartige Bildrate ist hoch genug, um die Bewegung quasi kontinuierlich darstellen zu können. Falls erforderlich, kann durch zeitliche Interpolation die Zeitauflösung noch weiter verbessert werden. Pro Atemzyklus werden bei dieser Vorgehensweise 40 Bilder und mehr erfasst.
  • Die Anzahl an insgesamt erfassten Bilddatensätzen kann nach Bedarf sein, beispielsweise über 100.
  • Die Art der Auswertung kann ebenfalls nach Bedarf bestimmt sein. Beispielsweise kann die Auswertung die Ermittlung eines statistischen Schwankungsmaßes von Bilddatenwerten der örtlich miteinander korrespondierenden Bereiche umfassen. Ein Beispiel eines geeigneten statistischen Schwankungsmaßes ist die Standardabweichung. Alternativ oder zusätzlich kann die Auswertung eine Auswertung des Fourierspekturms von Bilddatenwerten der örtlich miteinander korrespondierenden Bereiche umfassen. Insbesondere deutet ein starker niederfrequenter Anteil des Fourierspektrums auf eine verringerte Beweglichkeit hin.
  • Das erfindungsgemäße Auswertungsverfahren arbeitet besonders zuverlässig, wenn die ausgewerteten örtlich miteinander korrespondierenden Bereiche innerhalb der zweiten Kontur liegen.
  • Die vorliegende Erfindung wird weiterhin durch ein Computerprogramm der eingangs genannten Art gelöst, wobei die Ausführung des Computerprogramms bewirkt, dass die Recheneinrichtung ein erfindungsgemäßes Auswertungsverfahren ausführt. Das Computerprogramm kann beispielsweise auf einem Datenträger in maschinenlesbarer Form gespeichert sein. Der Datenträger kann insbesondere Bestandteil der Recheneinrichtung sein.
  • Die Aufgabe wird weiterhin durch eine Recheneinrichtung gelöst, die derart ausgebildet ist, dass sie im Betrieb ein erfindungsgemäßes Auswertungsverfahren ausführt.
  • Weitere Vorteile und Einzelheiten ergeben sich aus der nachfolgenden Beschreibung von Ausführungsbeispielen in Verbindung mit den Zeichnungen. Es zeigen in Prinzipdarstellung:
  • 1 eine Auswertungsanordnung,
  • 2 ein Zeitdiagramm,
  • 3 ein Ablaufdiagramm,
  • 4 schematisch einen Bilddatensatz und
  • 5 ein Ablaufdiagramm.
  • Gemäß 1 wird mittels einer geeigneten Erfassungseinrichtung 1 eine Folge von Bilddatensätzen B eines Bereichs 2 eines Lebewesens 3 erfasst. Der erfasste Bereich 2 ist derart gewählt, dass sich der erfasste Bereich 2 des Lebewesens 3 sowohl atmungsbedingt als auch herzschlagbedingt bewegt. Beispielsweise kann der erfasste Bereich 2 des Lebewesens 3 der Herz-Lungen-Region des Lebewesens 3 entsprechen. Die beiden genannten Bewegungen des Bereichs 2 des Lebewesens 3 erfolgen auch während der Erfassung der Bilddatensätze B.
  • Die Art der Erfassungseinrichtung 1 und die Art der Bilddatensätze B kann nach Bedarf gewählt sein. Beispielsweise können die Bilddatensätze B dreidimensionale Bilddatensätze sein, die mittels einer Magnetresonanzanlage erfasst werden. Alternativ können die Bilddatensätze B zweidimensionaler Natur sein. Auch kann die Erfassungseinrichtung 1 eine andere Art der Erfassung realisieren, die beispielsweise auf Röntgentechniken oder auf Ultraschalltechniken basiert. Die erfassten Bilddatensätze B können zweidimensional oder – bevorzugt – dreidimensional sein.
  • Die Bilddatensätze B werden gemäß 2 in der Regel mit einem festen zeitlichen Versatz T voneinander erfasst. Es ist in Einzelfällen möglich, dass der zeitliche Versatz T größer als die Dauer T' einer Atmungsperiode ist. In aller Regel ist der zeitliche Versatz T jedoch kleiner als die Dauer T' der Atmungsperiode, oftmals sogar erheblich kleiner. In der Regel ist der zeitliche Versatz T sogar kleiner als die Dauer T'' einer Herzperiode. Beispielsweise können pro Herzschlag 10 bis 15 Bilddatensätze B erfasst werden. Die hierfür erforderlichen Techniken sind Fachleuten bekannt. Beispielsweise kann bei Magnetresonanzanlagen die bereits erwähnte multislice 2D true FISP Technik in Verbindung mit einer parallelen Aufnahmetechnik (Vielkanalarrayspule oder dergleichen) verwendet werden.
  • Die Anzahl an Bilddatensätzen B kann nach Bedarf bestimmt werden. In vielen Fällen wird die Anzahl größer als 100 sein.
  • Die erfassten Bilddatensätze B werden einer Recheneinrichtung 4 zugeführt. Die Recheneinrichtung 4 wertet die Bilddatensätze B aus. Zum Auswerten der Bilddatensätze B führt die Recheneinrichtung 4 ein Computerprogramm 5 aus, mit dem die Recheneinrichtung 4 programmiert ist und das die Wirkungsweise der Recheneinrichtung 4 bestimmt.
  • Das Computerprogramm 5 ist innerhalb der Recheneinrichtung 4 in einem Datenträger 6 gespeichert, beispielsweise auf der Festplatte der Recheneinrichtung 4. Der Datenträger 6 ist Bestandteil der Recheneinrichtung 4. Die Speicherung erfolgt in maschinenlesbarer Form, meist sogar in ausschließlich maschinenlesbarer Form.
  • Das Computerprogramm 5 kann der Recheneinrichtung 4 über einen mobilen Datenträger 7 zugeführt werden, beispielsweise einem USB-Memorystick (dargestellt) oder eine SD-Speicherkarte (nicht dargestellt) oder einen anderen geeigneten mobilen Datenträger 7. Auch auf dem mobilen Datenträger 7 ist das Computerprogramm 5 in maschinenlesbarer Form (meist sogar in ausschließlich maschinenlesbarer Form) gespeichert. Alternativ ist es möglich, das Computerprogramm 5 der Recheneinrichtung 4 über eine Rechnernetzanbindung zuzuführen.
  • Das Computerprogramm 5 umfasst Maschinencode 8, der von der Recheneinrichtung 4 unmittelbar ausführbar ist. Das Ausführen des Maschinencodes 8 durch die Recheneinrichtung 4 bewirkt, dass die Recheneinrichtung 4 ein erfindungsgemäßes Auswertungsverfahren ausführt. Dieses Auswertungsverfahren wird nachfolgend in Verbindung mit 3 näher erläutert.
  • Gemäß 3 nimmt die Recheneinrichtung 4 in einem Schritt S1 zunächst die erfassten Bilddatensätze B entgegen.
  • In einem Schritt S2 wird in den erfassten Bilddatensätzen B eine erste Kontur 9 (siehe 4) bestimmt. Die erste Kontur 9 ist dadurch charakterisiert, dass sie sich entsprechend der atmungsbedingten Bewegung des erfassten Bereichs 2 des Lebewesens 3 bewegt.
  • Die erste Kontur 9 kann nach Bedarf bestimmt sein. In vielen Fällen wird es sich entsprechend der Darstellung von 4 um die Kontur der Lunge und/oder der Rippen des Lebewesens 3 handeln. Der Schritt S2 kann von der Recheneinrichtung 4 alternativ vollautomatisch oder interaktiv unter Mitwirkung eines Anwenders 10 (siehe 1) ausgeführt werden.
  • In einem Schritt S3 werden die erfassten Bilddatensätze B elastisch verzerrt. Die elastisch verzerrten Bilddatensätze werden nachstehend als einfach verzerrte Bilddatensätze bezeichnet und mit dem Bezugszeichen B' versehen, um sie von den erfassten, unverzerrten Bilddatensätzen B unterscheiden zu können. Das Verzerren der erfassten Bilddatensätze B erfolgt im Schritt S3 derart, dass die ersten Konturen 9 der einfach verzerrten Bilddatensätze B' örtlich miteinander korrespondieren. Elastische Verzerrungsverfahren als solche sind Fachleuten bekannt. Auf Verzerrungsverfahren muss im Rahmen der vorliegenden Erfindung daher nicht näher eingegangen werden.
  • Zur Implementierung des Schrittes S3 wird in der Regel einer der erfassten Bilddatensätze B zu einem Referenzbilddatensatz bestimmt. Dieser eine erfasste Bilddatensatz B (das heißt der Referenzbilddatensatz) bleibt unverzerrt. Die anderen Bilddatensätze B werden entsprechend der oben genannten Vorschrift verzerrt.
  • Der Referenzbilddatensatz kann prinzipiell beliebig bestimmt sein. In der Regel wird der Referenzbilddatensatz jedoch derart bestimmt, dass es sich in Bezug auf die Bewegung der ersten Kontur 9 um einen mittleren Bilddatensatz handelt. Beispielsweise kann insbesondere einer derjenigen Bilddatensätze B als Referenzbilddatensatz herangezogen werden, bei dem sich die atmungsbedingte Bewegung in etwa in der Mitte zwischen eingeatmetem und ausgeatmetem Zustand befindet. Die Auswahl des entsprechenden erfassten Bilddatensatzes B kann alternativ automatisch oder interaktiv unter Mitwirkung des Anwenders 10 erfolgen. Das Verzerren der übrigen Bilddatensätze B wird von der Recheneinrichtung 4 in aller Regel vollautomatisch ausgeführt.
  • Nach Ausführung des Schrittes S3 korrespondieren die ersten Konturen 9 der einfach verzerrten Bilddatensätze B' örtlich miteinander.
  • Nunmehr wird in einem Schritt S4 eine zweite Kontur 11 (siehe 4) bestimmt. Die Bestimmung des Schrittes S4 kann alternativ in den erfassten Bilddatensätzen B oder in den einfach verzerrten Bilddatensätzen B' vorgenommen werden. Der Schritt S4 kann – analog zum Schritt S2 – alternativ automatisch oder interaktiv unter Mitwirkung des Anwenders 10 ausgeführt werden. Er kann, falls die zweite Kontur 11 in den erfassten Bilddatensätzen B bestimmt wird, vor dem Schritt S3 und sogar vor dem Schritt S2 ausgeführt werden.
  • Die zweite Kontur 11 ist dadurch bestimmt, dass sie sich entsprechend der herzschlagbedingten Bewegung des erfassten Bereichs 2 des Lebewesens 3 bewegt. Insbesondere kann es sich bei der zweiten Kontur 11 um die Herzkontur selbst handeln. Alternativ kann jedoch auch eine andere Kontur herangezogen werden.
  • In einem Schritt S5 werden die einfach verzerrten Bilddatensätze B' erneut verzerrt und so zweifach verzerrte Bilddatensätze B'' erzeugt. Die zweifach verzerrten Bilddatensätze werden nachfolgend zur Unterscheidung von den erfassten Bilddatensätzen B und den einfach verzerrten Bilddatensätzen B' mit dem Bezugszeichen B'' versehen.
  • Das Verzerren des Schrittes S5 erfolgt derart, dass die zweiten Konturen 11 der Bilddatensätze B örtlich miteinander korrespondieren. Seine Implementierung kann analog zum Schritt S3 sein.
  • Im Rahmen des Schrittes S5 wird oftmals analog zum Schritt S3 – einer der einfach verzerrten Bilddatensätze B' zu einem Referenzdatensatz bestimmt, bezüglich dessen die anderen einfach verzerrten Bilddatensätze B' verzerrt werden. Analog zur Vorgehensweise im Rahmen des Schrittes S3 wird in der Regel ein einfach verzerrter Bilddatensatz B' zum Referenzdatensatz bestimmt, bei dem ein mittlerer Bewegungszustand der zweiten Kontur 11 vorliegt, also zwischen einer maximalen Kontraktion des Herzens und einer maximalen Relaxation des Herzens. Der Referenzbilddatensatz für die Ausführung des Schrittes S5 kann mit dem Referenzbilddatensatz für das Ausführen des Schrittes S3 identisch sein. Zwingend ist dies jedoch nicht erforderlich. Es kann sich alternativ um einen anderen Bilddatensatz handeln, der im Rahmen des Schrittes S3 bereits verzerrt wurde.
  • Der Schritt S5 wird in aller Regel von der Recheneinrichtung 4 vollautomatisch ausgeführt. Lediglich das Bestimmen des Referenzbilddatensatzes kann gegebenenfalls unter Mitwirkung des Anwenders 10 erfolgen.
  • In einem Schritt S6 werden örtlich miteinander korrespondierende Bereiche 12 der zweifach verzerrten Bilddatensätze B'' (siehe 4) bilddatensatzübergreifend ausgewertet.
  • Die Art der Auswertung des Schrittes S6 kann nach Bedarf bestimmt sein. Beispielsweise kann zur Implementierung des Schrittes S6 entsprechend 5 in einem Schritt S11 die Ermittlung eines statistischen Schwankungsmaßes – beispielsweise der Standardabweichung – von Bilddatenwerten der örtlich miteinander korrespondierenden Bereiche 12 vorgenommen werden. Die Bereiche 12 können nach Bedarf bestimmt sein. Es kann sich um einzelne Pixel bzw. Voxel handeln oder um Gruppen benachbarter Pixel bzw. Voxel. Liegt das statistische Schwankungsmaß oberhalb einer bestimmten Schwelle, kann der entsprechende Bereich 12 beispielsweise als „gesund” klassifiziert werden. Liegt das statistische Schwankungsmaß unterhalb der Schranke, kann der entsprechende Bereich 12 als krankhaft verändert klassifiziert werden.
  • Der Schritt S11 kann, soweit es die Ermittlung des statistischen Schwankungsmaßes und den Vergleich mit der bestimmten Schwelle betrifft, von der Recheneinrichtung 4 automatisch ausgeführt werden. Die diagnostische Einstufung selbst muss selbstverständlich der Anwender 10 vornehmen. Die Schwelle kann fest vorgegeben sein oder vom Anwender 10 einstellbar sein.
  • Alternativ oder zusätzlich kann zur Implementierung des Schrittes S6 gemäß 5 in einem Schritt S12 eine Auswertung des Fourierspektrums der Bilddatenwerte der entsprechenden Bereiche 12 vorgenommen werden. In diesem Fall kann insbesondere ermittelt werden, welches Gewicht einer Nullfrequenz zukommt, also einem Gleichanteil entspricht. Je größer dieses Gewicht ist, desto größer ist die Wahrscheinlichkeit, dass eine krankhafte Veränderung vorliegt.
  • Der Schritt S12 kann analog zum Schritt S11 bezüglich der Ermittlung des Fourierspektrums und des Gewichtes einzelner Frequenzbereiche von der Recheneinrichtung 4 automatisch ausgeführt werden. Die diagnostische Einstufung selbst muss jedoch der Anwender 10 vornehmen.
  • Die Bereiche 12, die in den zweifach verzerrten Bilddatensätzen B'' ausgewertet werden, können prinzipiell beliebig angeordnet sein. In der Regel werden sie jedoch entsprechend der Darstellung von 4 innerhalb der zweiten Kontur 11 liegen.
  • Die vorliegende Erfindung weist viele Vorteile auf. Insbesondere ist es im Rahmen des erfindungsgemäßen Auswertungsverfahrens nicht erforderlich, ein Atem-Gating vorzunehmen. Auch ist es nicht erforderlich, die Bilddatensätze B mit dem Herzschlag zu synchronisieren. Je nach Art der Datenerfassung kann es sogar möglich sein, ohne Kontrastmittel zu arbeiten. Weiterhin ist es nicht erforderlich, ionisierende Strahlung zu verwenden.
  • Die obige Beschreibung dient ausschließlich der Erläuterung der vorliegenden Erfindung. Der Schutzumfang der vorliegenden Erfindung soll hingegen ausschließlich durch die beigefügten Ansprüche bestimmt sein.
  • Bezugszeichenliste
    • 1
    Erfassungseinrichtung
    2
    Bereich des Lebewesens
    3
    Lebewesen
    4
    Recheneinrichtung
    5
    Computerprogramm
    6, 7
    Datenträger
    8
    Maschinencode
    9, 11
    Konturen
    10
    Anwender
    12
    Bereich der zweifach verzerrten Bilddatensätze
    B, B', B''
    Bilddatensätze
    S1–S12
    Schritte
    T
    zeitlicher Versatz
    T', T''
    Dauern

Claims (9)

  1. Auswertungsverfahren für eine Folge von Bilddatensätzen (B) eines Bereichs (2) eines Lebewesens (3), – wobei die Folge von Bilddatensätzen (B) erfasst wurde, während sich der erfasste Bereich (2) des Lebewesens (3) atmungs- und herzschlagbedingt bewegte, – wobei in den erfassten Bilddatensätzen (B) eine erste Kontur (9) bestimmt wird, die sich entsprechend der atmungsbedingten Bewegung des erfassten Bereichs (2) des Lebewesens (3) bewegt, – wobei die erfassten Bilddatensätze (B) derart elastisch zu einfach verzerrten Bilddatensätzen (B') verzerrt werden, dass die ersten Konturen (9) der einfach verzerrten Bilddatensätze (B') örtlich miteinander korrespondieren, – wobei in den erfassten Bilddatensätzen (B) oder in den einfach verzerrten Bilddatensätzen (B') eine zweite Kontur (11) bestimmt wird, die sich entsprechend der herzschlagbedingten Bewegung des erfassten Bereichs (2) des Lebewesens (3) bewegt, – wobei die einfach verzerrten Bilddatensätze (B') derart elastisch zu zweifach verzerrten Bilddatensätzen (B'') verzerrt werden, dass die zweiten Konturen (11) der zweifach verzerrten Bilddatensätze (B'') örtlich miteinander korrespondieren, und – wobei örtlich miteinander korrespondierende Bereiche (12) der zweifach verzerrten Bilddatensätze (B'') bilddatensatzübergreifend ausgewertet werden.
  2. Auswertungsverfahren nach Anspruch 1, dadurch gekennzeichnet, dass die Bilddatensätze (B) dreidimensionale Bilddatensätze sind.
  3. Auswertungsverfahren nach Anspruch 1 oder 2, dadurch gekennzeichnet, dass die Auswertung die Ermittlung eines statistischen Schwankungsmaßes, beispielsweise der Standardabweichung, von Bilddatenwerten der örtlich miteinander korrespondierenden Bereiche (12) umfasst.
  4. Auswertungsverfahren nach Anspruch 1, 2 oder 3, dadurch gekennzeichnet, dass die Auswertung eine Auswertung des Fourierspektrums von Bilddatenwerten der örtlich miteinander korrespondierenden Bereiche (12) umfasst.
  5. Auswertungsverfahren nach einem der obigen Ansprüche, dadurch gekennzeichnet, dass die ausgewerteten örtlich miteinander korrespondierenden Bereiche (12) innerhalb der zweiten Kontur (11) liegen.
  6. Computerprogramm, das Maschinencode (8) umfasst, der von einer Recheneinrichtung (4) unmittelbar ausführbar ist und dessen Ausführung durch die Recheneinrichtung (4) bewirkt, dass die Recheneinrichtung (4) ein Auswertungsverfahren nach einem der obigen Ansprüche ausführt.
  7. Computerprogramm nach Anspruch 6, dadurch gekennzeichnet, dass es auf einem Datenträger (6, 7) in maschinenlesbarer Form gespeichert ist.
  8. Computerprogramm nach Anspruch 7, dadurch gekennzeichnet, dass der Datenträger (6) Bestandteil der Recheneinrichtung (4) ist.
  9. Recheneinrichtung, dadurch gekennzeichnet, dass die Recheneinrichtung derart ausgebildet ist, dass sie im Betrieb ein Auswertungsverfahren nach einem der Ansprüche 1 bis 5 ausführt.
DE102010014211A 2010-04-08 2010-04-08 Auswertungsverfahren für eine Folge von Bilddatensätzen mit zweifacher elastischer Verzerrung der Bilddatensätze Withdrawn DE102010014211A1 (de)

Priority Applications (2)

Application Number Priority Date Filing Date Title
DE102010014211A DE102010014211A1 (de) 2010-04-08 2010-04-08 Auswertungsverfahren für eine Folge von Bilddatensätzen mit zweifacher elastischer Verzerrung der Bilddatensätze
US13/080,096 US8463009B2 (en) 2010-04-08 2011-04-05 Evaluation method for a series of image data sets with two-fold elastic distortion of the image data sets

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
DE102010014211A DE102010014211A1 (de) 2010-04-08 2010-04-08 Auswertungsverfahren für eine Folge von Bilddatensätzen mit zweifacher elastischer Verzerrung der Bilddatensätze

Publications (1)

Publication Number Publication Date
DE102010014211A1 true DE102010014211A1 (de) 2011-10-13

Family

ID=44658026

Family Applications (1)

Application Number Title Priority Date Filing Date
DE102010014211A Withdrawn DE102010014211A1 (de) 2010-04-08 2010-04-08 Auswertungsverfahren für eine Folge von Bilddatensätzen mit zweifacher elastischer Verzerrung der Bilddatensätze

Country Status (2)

Country Link
US (1) US8463009B2 (de)
DE (1) DE102010014211A1 (de)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110507359B (zh) * 2014-08-28 2022-06-07 深圳迈瑞生物医疗电子股份有限公司 剪切波成像方法及系统
US9526468B2 (en) * 2014-09-09 2016-12-27 General Electric Company Multiple frame acquisition for exposure control in X-ray medical imagers

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102005010093A1 (de) * 2005-03-04 2006-09-14 Siemens Ag Verfahren zur Darstellung von Respirationsmustern
DE102005037426A1 (de) * 2005-08-08 2007-02-15 Siemens Ag Vorrichtung und Verfahren zum Fusionieren von Bilddatensätzen
US20080226149A1 (en) * 2005-08-04 2008-09-18 Hans-Aloys Wischmann Motion Compensation in Functional Imaging
DE102008008601A1 (de) * 2008-02-12 2009-08-20 Siemens Aktiengesellschaft Verfahren zur Verarbeitung von medizinischen Bilddaten und Magnetresonanzanlage zur Aufnahme und Verarbeitung medizinischer Bilddaten
US20100061611A1 (en) * 2008-09-11 2010-03-11 Siemens Corporate Research, Inc. Co-registration of coronary artery computed tomography and fluoroscopic sequence

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102006011242B4 (de) * 2006-03-10 2012-03-29 Siemens Ag Verfahren zur Rekonstruktion einer 3D-Darstellung
US7539533B2 (en) * 2006-05-16 2009-05-26 Bao Tran Mesh network monitoring appliance

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102005010093A1 (de) * 2005-03-04 2006-09-14 Siemens Ag Verfahren zur Darstellung von Respirationsmustern
US20080226149A1 (en) * 2005-08-04 2008-09-18 Hans-Aloys Wischmann Motion Compensation in Functional Imaging
DE102005037426A1 (de) * 2005-08-08 2007-02-15 Siemens Ag Vorrichtung und Verfahren zum Fusionieren von Bilddatensätzen
DE102008008601A1 (de) * 2008-02-12 2009-08-20 Siemens Aktiengesellschaft Verfahren zur Verarbeitung von medizinischen Bilddaten und Magnetresonanzanlage zur Aufnahme und Verarbeitung medizinischer Bilddaten
US20100061611A1 (en) * 2008-09-11 2010-03-11 Siemens Corporate Research, Inc. Co-registration of coronary artery computed tomography and fluoroscopic sequence

Also Published As

Publication number Publication date
US20110249874A1 (en) 2011-10-13
US8463009B2 (en) 2013-06-11

Similar Documents

Publication Publication Date Title
DE102014219915B3 (de) Verfahren zur Kompensation von durch eine zyklische Bewegung eines Patienten auftretenden Ortszuordnungsfehlern von PET-Daten und Bildaufnahmesystem
DE102013205830B4 (de) Verfahren und Bilddatenerzeugungseinrichtung zur Erzeugung von Bilddaten eines sich bewegenden Objekts, Magnetresonanzanlage und Computerprogrammprodukt
DE102010036538A1 (de) System und Verfahren zum Kompensieren einer respiratorischen Bewegung in akquirierten radiographischen Bildern
DE102005017492A1 (de) Verfahren zum rechnerischen Kompensieren einer periodischen Bewegung eines Organs
DE112010001224T5 (de) Auf Bewegungskompensation basierende CT-Vorrichtung und CT-Verfahren
DE102005023907A1 (de) Verfahren zur Ermittlung von Positronen-Emissions-Messinformationen im Rahmen der Positronen-Emissions-Tomographie
DE202017106016U1 (de) Medizinische Informationsverarbeitungsvorrichtung, Röntgen-CT-Vorrichtung und computerlesbares Speichermedium mit einem Programm für ein medizinisches Informationsverarbeitungsverfahren
DE102014200303B4 (de) Verfahren zur Durchführung einer Positronenemissionstomographie in einer Hybridanlage und entsprechende Hybridanlage
DE102016219887A1 (de) Verfahren und System zur Nutzung von Messdaten
DE112015002791T5 (de) Intra-Rekonstruktions-Bewegungskorrektur
DE102009017439A1 (de) Verfahren und Vorrichtung zur Bildgebung eines vorbestimmten Volumenabschnitts mittels PET-Daten
DE202020005686U1 (de) MRT-abgeleitete Strain-basierte Messungen und damit zusammenhängende Bilddatenerfassung, Bilddatenverarbeitung, Patientenbewertungen und Überwachungssysteme
DE102009030110A1 (de) Verfahren zur Ermittlung der Ventilation einer Lunge
DE102014210218B4 (de) Verfahren zum Generieren von medizinischen Bilddatensätzen, sowie medizinisches Bildgebungsgerät
DE102012222073B4 (de) Verfahren zur Auswertung von Bilddatensätzen und Kombinationsbildaufnahmeeinrichtung
DE102015215584B4 (de) Verfahren und System zur Rekonstruktion von Planungsbildern
DE102008009190A1 (de) Verfahren zur automatischen Bestimmung der Ruhephase des Herzens
DE102013209939B4 (de) Verfahren zur Bewegungskorrektur von Emissionscomputertomographiedaten mittels Magnetresonanztomographiedaten
DE112014005888T5 (de) Verfahren und Vorrichtung zur bereichsspezifischen MR-basierten Abschwächungskorrektur von zeitbereichsgefilterten PET-Studien
DE102015200850B4 (de) Verfahren zur Auswertung von medizinischen Bilddaten
DE102014201499A1 (de) Verfahren und medizinische Bildgebungseinrichtung zur Bestimmung der Perfusion
DE102010014211A1 (de) Auswertungsverfahren für eine Folge von Bilddatensätzen mit zweifacher elastischer Verzerrung der Bilddatensätze
DE102007009182A1 (de) Verfahren und Vorrichtung zur Bilddarstellung von sich zyklisch bewegenden Objekten
DE102013218047B3 (de) Verfahren zur automatischen Anzeige und/oder Messung von Knochenveränderungen in medizinischen Bilddaten, sowie medizinisches Bildgebungsgerät und elektronisch lesbarer Datenträger
DE102015205270A1 (de) Verfahren zur Ermittlung eines Röntgenbilddatensatzes eines bewegten Zielgebiets und Röntgeneinrichtung

Legal Events

Date Code Title Description
R016 Response to examination communication
R079 Amendment of ipc main class

Free format text: PREVIOUS MAIN CLASS: A61B0019000000

Ipc: A61B0034000000

R079 Amendment of ipc main class

Free format text: PREVIOUS MAIN CLASS: A61B0034000000

Ipc: A61B0090000000

R081 Change of applicant/patentee

Owner name: SIEMENS HEALTHCARE GMBH, DE

Free format text: FORMER OWNER: SIEMENS AKTIENGESELLSCHAFT, 80333 MUENCHEN, DE

R119 Application deemed withdrawn, or ip right lapsed, due to non-payment of renewal fee