DE69931750T2 - Verfahren und gerät zur verkalkungsmessung - Google Patents

Verfahren und gerät zur verkalkungsmessung Download PDF

Info

Publication number
DE69931750T2
DE69931750T2 DE69931750T DE69931750T DE69931750T2 DE 69931750 T2 DE69931750 T2 DE 69931750T2 DE 69931750 T DE69931750 T DE 69931750T DE 69931750 T DE69931750 T DE 69931750T DE 69931750 T2 DE69931750 T2 DE 69931750T2
Authority
DE
Germany
Prior art keywords
interest
region
value
slice
image data
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
DE69931750T
Other languages
English (en)
Other versions
DE69931750D1 (de
Inventor
Hui Waukesha Hu
Alan Robert Dousman BUCHANAN
Theophano Danvers MITSA
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.)
General Electric Co
Original Assignee
General Electric Co
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 General Electric Co filed Critical General Electric Co
Publication of DE69931750D1 publication Critical patent/DE69931750D1/de
Application granted granted Critical
Publication of DE69931750T2 publication Critical patent/DE69931750T2/de
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Quality & Reliability (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Health & Medical Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Description

  • HINTERGRUND ZU DER ERFINDUNG
  • Diese Erfindung betrifft allgemein die Computertomographie-(CT)-Bildgebung und insbesondere die Erzeugung eines CT-Bild-Verkalkungswerts.
  • In wenigstens einer bekannten CT-System-Konfiguration projiziert eine Röntgenstrahlenquelle einen fächerförmigen Strahl, der so gebündelt ist, dass er innerhalb einer im Allgemeinen als die "Bildgebungsebene" bezeichneten xy-Ebene eines Kartesischen Koordinatensystems verläuft. Der Röntgenstrahl durchquert das abzubildende Objekt, beispielsweise einen Patienten. Der Strahl trifft, nachdem er durch das Objekt geschwächt wurde, auf eine Matrix von Strahlungsdetektoren. Die Intensität des an der Detektormatrix aufgefangenen geschwächten Strahls hängt von der Schwächung ab, die der Röntgenstrahl durch das Objekt erfährt. Jedes Detektorelement in der Matrix erzeugt ein eigenes elektrisches Signal, das die Schwächung des Strahles an der Detektorposition kennzeichnet. Die Schwächungsmesswerte sämtlicher Sensoren werden getrennt akquiriert, um ein Übertragungsprofil zu erzeugen.
  • In bekannten CT-Bildgebungssystemen der dritten Generation werden die Röntgenstrahlenquelle und die Detektormatrix mittels eines Gantryrahmens innerhalb der Bildgebungsebene und um das abzubildende Objekt herum gedreht, so dass sich der Winkel, unter dem der Röntgenstrahl auf das Objekt fällt, ständig ändert. Eine Gruppe von Röntgenstrahlschwächungsmesswerten, d. h. Projektionsdaten des De tektorarrays, unter ein und demselben Gantrywinkel, wird als eine "Ansicht" bezeichnet. Ein "Scandurchgang" des Objekts umfasst einen Satz von Ansichten, der während einer Umdrehung der Röntgenstrahlenquelle und des Detektors unter verschiedenen Gantrywinkeln erzeugt wird. Im Falle eines axialen Scannens werden die Projektionsdaten verarbeitet, um ein Bild zu konstruieren, das einem zweidimensionalen, durch das Objekt geführten Schnitt entspricht.
  • Ein Verfahren zum Rekonstruieren eines Bildes aus einem Satz von Projektionsdaten ist in der Fachwelt als die gefilterte Rückprojektionstechnik bekannt. Dieses Verfahren wandelt die Schwächungsmesswerte eines Scandurchgangs in als "CT-Zahlen" oder "Hounsfield-Einheiten" bezeichnete Integerzahlen um, die verwendet werden, um die Helligkeit eines entsprechenden Pixels auf einem Kathodenstrahlröhrenbildschirm zu steuern.
  • Um die für mehrere Schichten erforderliche Gesamtscanzeit zu reduzieren, kann ein "Spiral"-Scannen durchgeführt werden. Um ein "Spiral"-Scannen durchzuführen, wird der Patient bewegt, während die Daten für die vorgegebene Anzahl von Schichtbildern akquiriert werden. Ein derartiges System erzeugt anhand eines Spiralscandurchgangs mit nur einem Fächerstrahl eine einzelne Helix. Die durch den Fächerstrahl abgebildete Spirale liefert Projektionsdaten, aus denen sich Bilder in jeder vorgegebenen Schicht rekonstruieren lassen. Zusätzlich zu einer reduzierten Scanzeit weist das Spiralscannen weitere Vorteile auf, beispielsweise verbesserte Bildqualität und besser Kontraststeuerung.
  • Beim Spiralscannen wird, wie auch oben erläutert, an jedem Schichtort lediglich eine Ansicht von Daten gesam melt. Um ein Bild einer Schicht zu rekonstruieren, werden die sonstigen Ansichtsdaten für die Schicht basierend auf den für andere Ansichten gesammelten Daten erzeugt. Spiral-Rekonstruktionsalgorithmen sind bekannt und beispielsweise in C. Crawford und K. King, "Computed Tomography Scanning with Simultaneous Patient Translation", Med. Phys. 17(6), Nov./Dez. 1990 beschrieben.
  • Zumindest ein als Elektronenstrahl-Bildgebungssystem bekanntes Bildgebungssystem wird verwendet, um den Nachweis koronarer Atherosklerose zu führen, indem eine Koronararterienverkalkung (CAC = Coronary Artery Calcification) erfasst wird. Bei der Identifizierung von CAC in den Bilddaten wird ein Kalkablagerungsspiegel bestimmt. Allerdings sind Elektronenstrahlbildgebungssysteme sehr kostspielig und stehen nur an einer beschränkten Anzahl von geographischen Orten zur Verfügung. Bis heute lassen sich mittels Universal-CT-Bildgebungssystemen keine stabilen und konsistenten Kalkablagerungsspiegel erzeugen. Ein Faktor der inkonsistenten Ergebnisse ist der ungleichmäßige und/oder nicht zusammenhängende Abstand der Bilddaten.
  • Es wäre erwünscht, einen stabilen und konsistenten Verkalkungswert unter Verwendung von Bilddaten zu erzeugen, die von einem Universal-CT-Bildgebungssystem stammen. Es wäre ebenfalls erwünscht, den Verkalkungswert mit Blick auf ungleichmäßig und/oder nicht zusammenhängend beabstandete Schichten von Bilddaten zu korrigieren.
  • Agatston, A.S. et al. beschreiben in "Quantification of coronary artery calcification using ultrafast computed tomography", Journal Am. Coll. Cardiol, Bd. 15, Nr. 1990, Seiten 827–832 (XP000884744) ein Verfahren, bei dem eine Bewertung jedes interessierenden Bereichs durch Multiplikation eines Dichtewerts mit einer Fläche berechnet wird. Ein koronarer Gesamtkalziumwert wird berechnet, indem jeder dieser Werte für sämtliche Schichten aufsummiert wird.
  • Die US-A-4 837 686 beschreibt ein Verfahren, bei dem eine Kalziumgesamtmasse eines Knotens berechnet wird anhand der Gleichung:
    Figure 00040001
    mit P gleich einem Knochen-Bildpixelwert, K gleich einer Konstanten, die für ein spezielles verwendetes Bildgebungssystem charakteristisch ist, AP gleich einer Fläche eines Pixels, und M gleich einer Knotenvergrößerung.
  • Ukai, Y. et al. beschreiben in "A coronary calcification diagnosis system based on helical CT images", IEEE Nuclear Science Symposium Conference Record, Alberquerque, NM, USA, 9. November 1997, Seiten 1208 bis 1212 (XP002131567), ein Verfahren, bei dem eine CT-Wert-Gleichförmigkeit als ein mittlerer CT-Wert, dividiert durch einen maximalen CT-Wert in einer Herzregion berechnet wird.
  • Schließlich beschreiben Wilson, M. et al. in "Automated Detection of Microcalcifications in Mammograms Through Application of Image Pixel Remapping and Statistical Filter", Proc. 11th IEEE Symp. On Computer-based medical systems, Lubbock, TX, USA, 12.–14. Juni 1998, Seiten 270 bis 274 (XP002131568), einige Bildverfeinerungstechniken.
  • Die Techniken betreffen Bearbeitungen von Histogrammen, Bildkonvolution und statistisches Filter.
  • Erfindungsgemäß sind ein Verfahren zum Ermitteln einer Koronararterienverkalkung anhand von Computertomographiebilddaten, wie es in dem unabhängigen Anspruch 1 definiert ist, und ein entsprechendes System geschaffen, wie es in dem unabhängigen Anspruch 6 definiert ist.
  • KURZDARSTELLUNG DER ERFINDUNG
  • Diese und andere Aufgaben können durch einen Kalkablagerungsbewertungsalgorithmus erreicht werden, der einen Kalkablagerungsspiegel mittels eines Universal-CT-Bildgebungssystems erzeugt. In Übereinstimmung mit einem Ausführungsbeispiel der vorliegenden Erfindung wird der Patient gescannt, um Projektionsdaten zu erzeugen. Die Projektionsdaten werden verarbeitet, um Bilddaten zu erzeugen. Die Bilddaten werden anschließend verarbeitet, um einen Verkalkungswert zu ermitteln, der eine Koronararterienverkalkung repräsentiert. Der Verkalkungswert wird ermittelt, indem eine bewertbare Region in den Bilddaten identifiziert wird, die wenigstens einen interessierenden Bereich in der bewertbaren Region definiert, und indem ein Dichtewert ermittelt wird. Nach der Ermittlung eines Kalziumwerts für jeden interessierenden Bereich, wird ein Gesamtkalziumwert bestimmt, indem der von jedem interessierenden Bereich stammende Kalziumwert aufsummiert wird. Der Gesamtkalziumspiegel repräsentiert den Kalkablagerungsspiegel der Bilddaten.
  • In einem Ausführungsbeispiel wird der Gesamtkalziumwert mit Blick auf ungleichmäßig und/oder nicht zusammen hängend beabstandete Schichten von Bilddaten korrigiert. Insbesondere wichtet der Algorithmus dort, wo die Bilddaten ungleichmäßig und/oder nicht zusammenhängend beabstandete Datenschichten repräsentieren, den aus jedem interessierenden Bereich stammenden Kalziumwert, um das Schichtabstandsproblem zu kompensieren.
  • Der oben beschriebene Algorithmus erzeugt einen Verkalkungswert von Bilddaten, die mittels eines Universal-CT-Bildgebungssystems gewonnen wurden. Darüber hinaus führt der Algorithmus eine Korrektur mit Blick auf nicht gleichmäßig und/oder nicht zusammenhängend beabstandete Schichten von Bilddaten durch.
  • KURZBESCHREIBUNG DER ZEICHNUNGEN
  • 1 zeigt eine anschauliche Ansicht eines CT-Bildgebungssystems.
  • 2 zeigt ein Blockschaltbild des in 1 veranschaulichten Systems.
  • 3 zeigt eine anschauliche Ansicht von Spalten zwischen Bilddaten, gemäß einem Ausführungsbeispiel der vorliegenden Erfindung.
  • DETAILLIERTE BESCHREIBUNG DER ERFINDUNG
  • In 1 und 2 ist ein Computertomographie-(CT)-Bildgebungssystem 10 gezeigt, das einen Gantryrahmen 12 aufweist, wie er für CT-Scanner der "dritten Generation" typisch ist. Die Gantry 12 weist eine Röntgenstrahlenquelle 14 auf, die ein Bündel von Röntgenstrahlen 16 in Richtung eines auf der entgegengesetzten Seite der Gantry 12 angeordneten Detektorarrays 18 projiziert. Ein Röntgenstrahl wird durch einen (nicht gezeigten) Kollimator gebündelt, um innerhalb einer im Allgemeinen als "Bildgebungsebene" bezeichneten xy-Ebene eines kartesischen Koordinatensystems zu verlaufen. Die Detektormatrix 18 wird durch Detektorelemente 20 gebildet, die die projizierten Röntgenstrahlen, die einen Patienten 22 durchqueren, gemeinsam erfassen. Die Detektormatrix 20 kann ein Einzelschicht-Detektor oder ein Multischicht-Detektor sein. Jedes Detektorelement 20 erzeugt ein elektrisches Signal, das die Intensität eines einfallenden Röntgenstrahls, und dementsprechend die Schwächung des Strahles kennzeichnet, die dieser bei einer Durchquerung des Patienten 22 erfährt. Während des Scannens zum Gewinnen von Röntgenstrahlprojektionsdaten kreisen die Gantry 12 und die daran befestigten Komponenten um eine Rotationsachse 24.
  • Die Rotation des Gantryrahmens 12 und der Betrieb der Röntgenstrahlenquelle 14 werden durch eine Steuervorrichtung 26 des CT-Systems 10 gesteuert. Der Steuermechanismus 26 enthält einen Röntgencontroller 28, der Energie und Zeittaktsignale an die Röntgenstrahlenquelle 14 liefert, und einen Gantryantriebscontroller 30, der die Rotationsgeschwindigkeit und die Position der Gantry 12 steuert. Ein Datenakquirierungssystem (DAS) 32 in dem Steuermechanismus 26 tastet von Detektorelementen 20 stammende analoge Daten ab und wandelt die Daten in digitale Signale um, um diese anschließend zu verarbeiten. Ein Bildrekonstruktor 34 nimmt die von dem DAS 32 ankommenden abgetasteten und digitalisierten Röntgenstrahldaten entgegen und führt eine Hochgeschwindigkeits-Bildrekonstruktion durch. Das rekonstruierte Bild wird einem Computer 36 als Eingabe zugeführt, der das Bild in einem Massenspeichergerät 38 speichert.
  • Der Rechner 36 nimmt ferner über eine Konsole 40, die eine Tastatur aufweist, von einem Anwender Steuerbefehle und Scanparameter entgegen. Ein zugehöriges Datensichtgerät 42 ermöglicht es dem Anwender, das rekonstruierte Bild und sonstige von dem Rechner 36 ausgegebene Daten zu beobachten. Die durch den Anwender eingegebenen Steuerbefehle und Parameter werden von dem Rechner 36 verwendet, um Steuerungssignale und Daten an das DAS 32, den Röntgencontroller 28 und den Gantryantriebscontroller 30 auszugeben. Darüber hinaus betätigt der Rechner 36 einen Liegenantriebscontroller 44, der eine motorisch angetriebene Liege 46 steuert, um den Patienten 22 innerhalb des Gantryrahmens 12 zu positionieren. Insbesondere bewegt die Liege 46 den Patienten 22 abschnittsweise durch den Gantrytunnel 48.
  • Die folgende Erörterung, die einen Algorithmus zum Erzeugen eines Verkalkungswerts beschreibt, bezieht sich gelegentlich speziell auf einen Spiralscan. Der Kalziumbewertungsalgorithmus ist jedoch nicht ausschließlich auf die Verwendung in Verbindung mit Spiralscannen beschränkt und kann im Zusammenhang mit sonstigen Scanverfahren, z.B. Axialscannen, verwendet werden. Weiter ist zu beachten, dass der unten beschriebene Algorithmus in dem Computer 36 imp lementiert sein kann und beispielsweise rekonstruierte Bilddaten verarbeiten könnte. In einer Abwandlung könnte der Algorithmus in dem Bildrekonstruktor 34 implementiert sein und dem Computer 36 korrigierte Bilddaten liefern.
  • Wie oben beschrieben werden bei der Durchführung eines CT-Scans, von Detektorelementen 20 ausgegebene Daten gewonnen. Solche Daten werden in der Fachwelt im Allgemeinen als Projektionsdaten bezeichnet. In einem Ausführungsbeispiel beinhaltet das CT Projektionsdaten, die wenigstens eine Schicht des Patienten 22 kennzeichnen. In diesem Fall wird zur Erzeugung der Bilddaten eine Hochgeschwindigkeitsbildrekonstruktion durchgeführt. Hinsichtlich der Bildrekonstruktion sind viele Bildrekonstruktionsalgorithmen gegenwärtig in im Handel erhältlichen CT-Geräten implementiert, und der erfindungsgemäße Kalkablagerungsbewertungsalgorithmus könnte in Verbindung mit vielen derartiger Rekonstruktionsalgorithmen durchgeführt werden.
  • In Übereinstimmung mit einem Ausführungsbeispiel der vorliegenden Erfindung, wird der Verkalkungswert ermittelt, indem eine bewertbare Region in den Bilddaten identifiziert wird, die wenigstens einen interessierenden Bereich in der bewertbaren Region definiert, und indem ein Dichtewert ermittelt wird. Nach der Ermittlung eines Kalziumwerts für jeden interessierenden Bereich, wird ein Gesamtkalziumwert bestimmt.
  • In einem Ausführungsbeispiel wird wenigstens eine bewertbare Region identifiziert, indem an den Bilddaten eine Schwellwertoperation durchgeführt wird. Insbesondere wird jede CT-Zahl in den Bilddaten mit einem ausgewählten Schwellwert oder gewünschten Schwellwertbereich verglichen.
  • Im Allgemeinen weisen unterschiedliche Materialien verschiedene CT-Zahlen auf. Beispielsweise weist Knochengewebe eine CT-Zahl von über 200 auf, Wasser hat eine CT-Zahl von 0, Weichgewebe (im Gehirn) weist eine CT-Zahl von etwa 20–50 auf, und Luft hat eine CT-Zahl von –1000. Da die CT-Zahlen für vielfältige Materialien verschieden sind, kann der Schwellwert so ausgewählt werden, dass gewisse Arten oder Klassen von Objekten von einer weiteren Analyse ausgeschlossen werden. In einem Ausführungsbeispiel dient eines der vielen aus dem Stand der Technik bekannten Schwellwertoperationsverfahren zur Identifizierung der bewertbaren Regionen. Beispielsweise wird bei Verwendung eines ausgewählten Schwellwerts von 90 jeder Teil der Bilddaten, der eine CT-Zahl aufweist, die größer oder gleich 90 ist, als eine bewertbare Region identifiziert.
  • Nach der Ermittlung eines interessierenden Bereichs (ROI) für jede bewertbare Region wird, wie aus dem Stand der Technik bekannt, ein Dichtewert bestimmt. Der Dichtewert wird unter Verwendung einer Gewichtungsreferenztabelle bestimmt. In einem Ausführungsbeispiel ist die Gewichtungsreferenztabelle in dem Computer 36 gespeichert, und zwar in der Form:
    Figure 00100001
    mit:
    • PS
    = CT-Zahl eines Pixels P,
    LT
    = eine untere Schwellwertzahl,
    UT
    = eine obere Schwellwertzahl,
    S
    = ein Skalierungsfaktor und
    C
    = eine Konstante.
  • Der Dichtewert DS jedes interessierenden Bereichs wird anschließend unter Verwendung der Gewichtungsreferenztabelle ermittelt. Insbesondere wird für jedes Pixel jedes interessierenden Bereichs jeder bewertbaren Region ein Dichtewert DSP bestimmt.
  • Beispielsweise wird in einem Ausführungsbeispiel der Dichtewert DSP jedes Pixels mittels eines Skalierungsfaktor S von 0,01, einer Konstanten C von 0,5, eines unteren Schwellwerts LT von 90 und eines oberen Schwellwerts UT von 400 Hounsfield-Einheiten bestimmt. Insbesondere wird im Falle von Bilddaten, die (genau) eine bewertbare Region mit einem interessierenden Bereich von 100 Pixeln aufweisen, der Dichtewert jedes der 100 Pixel bestimmt. Selbstverständlich können andere Werte von S, C, LT und UT verwendet werden, um den Dichtewert jedes Pixels zu ermitteln. Beispielsweise kann der untere Schwellwert 130 betragen.
  • Als Nächstes wird der Kalziumwert jedes interessierenden Bereichs bestimmt. In einem Ausführungsbeispiel wird der Kalziumwert jedes ROI bestimmt als:
    Figure 00110001
    mit:
    • CSR
    = Kalziumwert eines interessierenden Bereichs R,
    AP
    = Fläche eines Pixels P,
    PR
    = eine Anzahl von bewertbaren Pixeln in dem interessierenden Bereich und
    DSP
    = Dichtewert eines Pixels P.
  • Insbesondere ist in einem Ausführungsbeispiel die Fläche AP eines Pixels P durch den Rekonstruktionsalgorithmus definiert. Beispielsweise wird im Falle der Verwendung einer 512 × 512 Matrix zur Abdeckung eines Kreises mit einem Durchmesser von 20 cm die Pixelabmessung bestimmt.
  • Nach Ermittlung des Kalziumwerts für jeden interessierenden Bereich wird dann der Gesamtkalziumwert ermittelt. Insbesondere wird für sämtliche interessierende Regionen ein Gesamtkalziumwert bestimmt, indem der aus jedem interessierenden Bereich stammende Kalziumwert aufaddiert wird.
  • In einem Ausführungsbeispiel, in dem Projektionsdaten keine gleichmäßig und/oder zusammenhängend beabstandete Schichten repräsentieren, wird ein gewichteter Gesamtkalziumwert bestimmt. Der gewichtete Gesamtkalziumwert kompensiert oder berücksichtigt die ungleichmäßig und/oder nicht zusammenhängend beabstandeten Schichten, d.h. ungleichmäßige Abstände zwischen Schichten und überlappende Schichten. Insbesondere errechnet sich der Gesamtkalziumwert durch:
    Figure 00130001
    mit:
    • N
    = ein Normierungsfaktor,
    F
    = ein Skalierungsfaktor,
    G
    = G1 + G2,
    G1
    = ein Abstand zwischen einer vorausgehenden Schicht und einer aktuellen Schicht,
    G2
    = ein Abstand zwischen der aktuellen Schicht und einer nachfolgenden Schicht.
  • Falls entweder eine vorausgehende Schicht oder eine nachfolgende Schicht nicht vorhanden ist, wird ersatzweise für G1 oder G2 die Dicke der aktuellen Schicht verwendet. Beispielsweise wird im Falle einer einzelnen Schicht, die eine Dicke von 5 mm aufweist, der Gesamtkalziumwert unter Verwendung von G = 5 mm bestimmt.
  • In einer ähnlichen Weise lässt sich der Gesamtkalziumwert für überlappende Schichten von Bilddaten ermitteln. Selbstverständlich können andere Werte für N bestimmt werden, beispielsweise N = 0,33.
  • Aus der vorhergehenden Beschreibung vielfältiger Ausführungsbeispiele der vorliegenden Erfindung geht offensichtlich hervor, dass die Aufgaben der Erfindung erreicht sind. Obwohl die Erfindung in Einzelheiten beschrieben und veranschaulicht wurde, ist klar hervorzuheben, dass die Beschreibung lediglich der Veranschaulichung und als Beispiel dient und nicht als beschränkend zu verstehen ist. Dementsprechend ist der Schutzumfang der Erfindung lediglich durch die Begriffe der beigefügten Patentansprüche zu beschränken.

Claims (11)

  1. Verfahren zum Ermitteln einer Koronararterienverkalkung anhand von Computertomographiebilddaten, wobei zu dem Verfahren die Schritte gehören: Sammeln von mindestens zwei Schichten von Bilddaten; Identifizieren von mindestens zwei interessierenden Regionen anhand von mindestens zwei Schichten von Bilddaten; Ermitteln eines Dichtewerts für jeden interessierenden Bereich; Ermitteln eines Verkalkungswerts mittels des Dichtewerts durch Ermitteln eines Schichtbildkalziumwerts jedes interessierenden Bereichs; und Ermitteln eines gewichteten Gesamtkalziumwerts unter Verwendung jedes Schichtbildkalziumwerts und basierend auf einem Abstand zwischen mindestens zwei benachbarten Schichtbildern der mindestens zwei Schichtbilder.
  2. Verfahren nach Anspruch 1, bei dem für jeden interessierenden Bereich der Dichtewert DS bestimmt wird anhand der Formel:
    Figure 00160001
    mit: HP = eine CT-Zahl von Pixel P, LT = eine untere Schwellwertzahl, UT = eine obere Schwellwertzahl, S = ein Skalierungsfaktor, und C = eine Konstante.
  3. Verfahren nach Anspruch 2, ferner mit dem Schritt, vor einem Identifizieren wenigstens eines interessierenden Bereichs basierend auf einer Computertomographie-Zahl des bewertbaren Bereich eine bewertbare Region als einen Schwellwert zu identifizieren.
  4. Verfahren nach Anspruch 1, bei dem das Ermitteln des Schichtbildkalziumwerts für jeden interessierenden Bereich durchgeführt wird anhand der Gleichung:
    Figure 00160002
    mit: R = ein interessierender Bereich, CSR = ein Kalziumwert jedes interessierenden Bereichs, AP = eine Fläche eines Pixels P in dem interessierenden Bereich R, PR = eine Anzahl von bewertbaren Pixeln in dem interessierenden Bereich R, und DSP = ein Dichtewert von Pixel P.
  5. Verfahren nach Anspruch 1, bei dem das Ermitteln des gewichtetes Gesamtkalziumwerts TCS durchgeführt wird nach der Formel:
    Figure 00170001
    mit: N = ein Normierungsfaktor, F = ein Skalierungsfaktor, G = G1 + G2 G1 = ein Abstand zwischen einem vorhergehenden Schichtbild und einem aktuellen Schichtbild, G2 = ein Abstand zwischen dem aktuellen Schichtbild und einem nachfolgenden Schichtbild, R = eine Anzahl von interessierenden Bereichen.
  6. System zum Ermitteln des Grades einer Verkalkung anhand von Bilddaten, nämlich Bilddaten von einem tomographischen Scandurchgang, wobei das System eine Röntgenstrahlenquelle 14 und eine Detektormatrix (20) enthält, wobei das System dazu eingerichtet ist, mindestens zwei Schichten von Bilddaten zu sammeln; anhand von mindestens zwei Schichten von Bilddaten mindestens zwei interessierende Regionen zu identifizieren; für jeden interessierenden Bereich einen Dichtewert zu bestimmen; unter Verwendung des Dichtewert durch Ermittlung eines Schichtbildkalziumwerts jedes interessierenden Bereichs einen Verkalkungswert zu bestimmen; und unter Verwendung jedes Schichtbildkalziumwerts und basierend auf einem Abstand zwischen mindestens zwei benachbarten Schichtbildern der mindestens zwei Schichtbilder einen gewichteten Gesamtkalziumwert zu bestimmen.
  7. System nach Anspruch 6, bei dem der Dichtewert für jeden interessierenden Bereich nach folgender Formel bestimmt wird:
    Figure 00190001
    mit: Hp = eine CT-Zahl von Pixel P, LT = eine untere Schwellwertzahl, UT = eine obere Schwellwertzahl, S = ein Skalierungsfaktor, und C = eine Konstante.
  8. System nach Anspruch 6, das ferner dazu eingerichtet ist, vor dem Identifizieren wenigstens eines interessierenden Bereichs basierend auf einer Computertomographie-Zahl der bewertbaren Region als einen Schwellwert eine bewertbare Region zu identifizieren.
  9. System nach Anspruch 6, bei dem die Bestimmung des Schichtbildkalziumwerts für jeden interessierenden Bereich nach der folgenden Formel durchgeführt wird:
    Figure 00190002
    mit: R = ein interessierender Bereich, CSR = ein Kalziumwert von jeden interessierenden Bereich, AP = eine Fläche eines Pixels P in dem interessierenden Bereich R, PR = eine Anzahl von bewertbaren Pixeln in dem interessierenden Bereich R und DSP = ein Dichtewert von Pixel P.
  10. System nach Anspruch 6, bei dem der gewichtete Gesamtkalziumwert bestimmt wird nach der Formel:
    Figure 00200001
    mit: N = eine Normierungsfaktor, F = eine Skalierungsfaktor, G = G1 + G2, G1 = ein Abstand zwischen einem vorhergehenden Schichtbild und einem aktuellen Schichtbild, G = ein Abstand zwischen einem vorhergehenden Schichtbild und einem nachfolgenden Schichtbild, und R = eine Anzahl von interessierenden Regionen.
  11. System nach Anspruch 6, bei dem die Detektormatrix ein Multischicht-Detektorarray ist.
DE69931750T 1998-11-25 1999-11-23 Verfahren und gerät zur verkalkungsmessung Expired - Fee Related DE69931750T2 (de)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US199513 1988-05-27
US09/199,513 US6233304B1 (en) 1998-11-25 1998-11-25 Methods and apparatus for calcification scoring
PCT/US1999/028072 WO2000033252A1 (en) 1998-11-25 1999-11-23 Methods and apparatus for calcification scoring

Publications (2)

Publication Number Publication Date
DE69931750D1 DE69931750D1 (de) 2006-07-20
DE69931750T2 true DE69931750T2 (de) 2007-05-16

Family

ID=22737848

Family Applications (1)

Application Number Title Priority Date Filing Date
DE69931750T Expired - Fee Related DE69931750T2 (de) 1998-11-25 1999-11-23 Verfahren und gerät zur verkalkungsmessung

Country Status (6)

Country Link
US (1) US6233304B1 (de)
EP (1) EP1048008B1 (de)
JP (1) JP4367884B2 (de)
DE (1) DE69931750T2 (de)
IL (1) IL137413A0 (de)
WO (1) WO2000033252A1 (de)

Families Citing this family (65)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6421552B1 (en) 1999-12-27 2002-07-16 Ge Medical Systems Global Technology Company, Llc Methods and apparatus for estimating cardiac motion using projection data
US6674834B1 (en) 2000-03-31 2004-01-06 Ge Medical Systems Global Technology Company, Llc Phantom and method for evaluating calcium scoring
US8565860B2 (en) 2000-08-21 2013-10-22 Biosensors International Group, Ltd. Radioactive emission detector equipped with a position tracking system
US8909325B2 (en) 2000-08-21 2014-12-09 Biosensors International Group, Ltd. Radioactive emission detector equipped with a position tracking system and utilization thereof with medical systems and in medical procedures
WO2005119025A2 (en) 2004-06-01 2005-12-15 Spectrum Dynamics Llc Radioactive-emission-measurement optimization to specific body structures
US8489176B1 (en) 2000-08-21 2013-07-16 Spectrum Dynamics Llc Radioactive emission detector equipped with a position tracking system and utilization thereof with medical systems and in medical procedures
US6697451B2 (en) * 2001-09-05 2004-02-24 Ge Medical Systems Global Technology Company, Llc Dynamic phantom and method for evaluating calcium scoring
EP1293925A1 (de) * 2001-09-18 2003-03-19 Agfa-Gevaert Verfahren zur Bestimmung von RÖntgenbildern
US7286866B2 (en) * 2001-11-05 2007-10-23 Ge Medical Systems Global Technology Company, Llc Method, system and computer product for cardiac interventional procedure planning
US7127096B2 (en) * 2001-11-20 2006-10-24 Accuimage Diagnostics Corp. Method and software for improving coronary calcium scoring consistency
US6990222B2 (en) * 2001-11-21 2006-01-24 Arnold Ben A Calibration of tissue densities in computerized tomography
US8649843B2 (en) * 2001-11-24 2014-02-11 Ben A. Arnold Automated calcium scoring of the aorta
EP1451753A2 (de) * 2001-11-24 2004-09-01 Image Analysis, Inc. Automatische detektion und quantifizierung des koronaren und aortischen kalzium
US7311705B2 (en) 2002-02-05 2007-12-25 Medtronic, Inc. Catheter apparatus for treatment of heart arrhythmia
US20030190063A1 (en) * 2002-03-08 2003-10-09 Acharya Kishore C. Method and system for performing coronary artery calcification scoring
US7499743B2 (en) * 2002-03-15 2009-03-03 General Electric Company Method and system for registration of 3D images within an interventional system
US7346381B2 (en) * 2002-11-01 2008-03-18 Ge Medical Systems Global Technology Company Llc Method and apparatus for medical intervention procedure planning
US7778686B2 (en) * 2002-06-04 2010-08-17 General Electric Company Method and apparatus for medical intervention procedure planning and location and navigation of an intervention tool
US6904118B2 (en) * 2002-07-23 2005-06-07 General Electric Company Method and apparatus for generating a density map using dual-energy CT
JP4619781B2 (ja) 2002-07-23 2011-01-26 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー プラークの成分を検出するためのシステム
US6922462B2 (en) * 2002-07-31 2005-07-26 Ge Medical Systems Global Technology Company, Llc Method, system and computer product for plaque characterization
US20040133100A1 (en) * 2002-08-23 2004-07-08 Morteza Naghavi Novel risk assessment method based upon coronary calcification distribution pattern imaged by computed tomography
US7149331B1 (en) 2002-09-03 2006-12-12 Cedara Software Corp. Methods and software for improving thresholding of coronary calcium scoring
US7747047B2 (en) 2003-05-07 2010-06-29 Ge Medical Systems Global Technology Company, Llc Cardiac CT system and method for planning left atrial appendage isolation
US7343196B2 (en) * 2003-05-09 2008-03-11 Ge Medical Systems Global Technology Company Llc Cardiac CT system and method for planning and treatment of biventricular pacing using epicardial lead
US7565190B2 (en) * 2003-05-09 2009-07-21 Ge Medical Systems Global Technology Company, Llc Cardiac CT system and method for planning atrial fibrillation intervention
US7813785B2 (en) 2003-07-01 2010-10-12 General Electric Company Cardiac imaging system and method for planning minimally invasive direct coronary artery bypass surgery
US7344543B2 (en) 2003-07-01 2008-03-18 Medtronic, Inc. Method and apparatus for epicardial left atrial appendage isolation in patients with atrial fibrillation
US7308299B2 (en) 2003-10-22 2007-12-11 General Electric Company Method, apparatus and product for acquiring cardiac images
US7308297B2 (en) * 2003-11-05 2007-12-11 Ge Medical Systems Global Technology Company, Llc Cardiac imaging system and method for quantification of desynchrony of ventricles for biventricular pacing
EP1694212B1 (de) * 2003-11-25 2013-08-07 Koninklijke Philips Electronics N.V. Automatische dosisanpassung als funktion des body mass index eines patienten bei der ct-calcium-bewertung
US7176466B2 (en) 2004-01-13 2007-02-13 Spectrum Dynamics Llc Multi-dimensional image reconstruction
US9470801B2 (en) 2004-01-13 2016-10-18 Spectrum Dynamics Llc Gating with anatomically varying durations
US7968851B2 (en) 2004-01-13 2011-06-28 Spectrum Dynamics Llc Dynamic spect camera
US8571881B2 (en) 2004-11-09 2013-10-29 Spectrum Dynamics, Llc Radiopharmaceutical dispensing, administration, and imaging
US8586932B2 (en) 2004-11-09 2013-11-19 Spectrum Dynamics Llc System and method for radioactive emission measurement
WO2007010534A2 (en) 2005-07-19 2007-01-25 Spectrum Dynamics Llc Imaging protocols
WO2005118659A2 (en) * 2004-06-01 2005-12-15 Spectrum Dynamics Llc Methods of view selection for radioactive emission measurements
US9040016B2 (en) 2004-01-13 2015-05-26 Biosensors International Group, Ltd. Diagnostic kit and methods for radioimaging myocardial perfusion
JP4206044B2 (ja) * 2004-01-20 2009-01-07 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー カルシウム・スコア測定方法および装置
US7454248B2 (en) * 2004-01-30 2008-11-18 Ge Medical Systems Global Technology, Llc Method, apparatus and product for acquiring cardiac images
GB2416223A (en) * 2004-07-15 2006-01-18 Medicsight Plc Quantification of coronary artery calcification
US8515527B2 (en) * 2004-10-13 2013-08-20 General Electric Company Method and apparatus for registering 3D models of anatomical regions of a heart and a tracking system with projection images of an interventional fluoroscopic system
US7327872B2 (en) * 2004-10-13 2008-02-05 General Electric Company Method and system for registering 3D models of anatomical regions with projection images of the same
US8615405B2 (en) 2004-11-09 2013-12-24 Biosensors International Group, Ltd. Imaging system customization using data from radiopharmaceutical-associated data carrier
US9316743B2 (en) 2004-11-09 2016-04-19 Biosensors International Group, Ltd. System and method for radioactive emission measurement
US9943274B2 (en) 2004-11-09 2018-04-17 Spectrum Dynamics Medical Limited Radioimaging using low dose isotope
US8000773B2 (en) * 2004-11-09 2011-08-16 Spectrum Dynamics Llc Radioimaging
EP1827505A4 (de) 2004-11-09 2017-07-12 Biosensors International Group, Ltd. Radiodarstellung
WO2008059489A2 (en) 2006-11-13 2008-05-22 Spectrum Dynamics Llc Radioimaging applications of and novel formulations of teboroxime
US8111886B2 (en) * 2005-07-19 2012-02-07 Spectrum Dynamics Llc Reconstruction stabilizer and active vision
US8837793B2 (en) 2005-07-19 2014-09-16 Biosensors International Group, Ltd. Reconstruction stabilizer and active vision
US8894974B2 (en) 2006-05-11 2014-11-25 Spectrum Dynamics Llc Radiopharmaceuticals for diagnosis and therapy
US9275451B2 (en) 2006-12-20 2016-03-01 Biosensors International Group, Ltd. Method, a system, and an apparatus for using and processing multidimensional data
US8521253B2 (en) 2007-10-29 2013-08-27 Spectrum Dynamics Llc Prostate imaging
US8186880B1 (en) 2008-11-27 2012-05-29 Arnold Ben A Extended and fixed INTable simultaneously imaged calibration and correction methods and references for 3-D imaging devices
US8649577B1 (en) 2008-11-30 2014-02-11 Image Analysis, Inc. Automatic method and system for measurements of bone density and structure of the hip from 3-D X-ray imaging devices
US8338788B2 (en) * 2009-07-29 2012-12-25 Spectrum Dynamics Llc Method and system of optimized volumetric imaging
JP4981938B2 (ja) * 2010-03-08 2012-07-25 富士フイルム株式会社 診断支援装置、冠動脈解析プログラムおよび冠動脈解析方法
US8517608B1 (en) 2011-08-03 2013-08-27 Ben A. Arnold System and method for calibration of CT scanners and display of images in density units without the use of water phantoms
US8867822B2 (en) 2011-10-14 2014-10-21 Fujifilm Corporation Model-based coronary artery calcium scoring
KR101984247B1 (ko) 2012-03-15 2019-05-30 삼성전자 주식회사 관상동맥석회화 수준 변화 예측장치 및 예측방법
CN113362404B (zh) * 2020-03-05 2024-03-22 上海西门子医疗器械有限公司 计算机断层扫描的散射校正方法、装置和存储介质
CN112927196B (zh) * 2021-02-09 2024-03-22 深圳睿心智能医疗科技有限公司 钙化积分方法及装置
WO2022185231A1 (en) * 2021-03-04 2022-09-09 Multict Imaging Ltd X-ray computed tomography (ct) scanner

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4837686A (en) 1985-11-15 1989-06-06 Picker International Substance quantification in animal bodies
US5864146A (en) * 1996-11-13 1999-01-26 University Of Massachusetts Medical Center System for quantitative radiographic imaging
IL134557A0 (en) 1997-08-28 2001-04-30 Qualia Computing Inc Method and system for automated detection of clustered microcalcifications from digital mammograms

Also Published As

Publication number Publication date
JP2002531199A (ja) 2002-09-24
DE69931750D1 (de) 2006-07-20
IL137413A0 (en) 2001-07-24
WO2000033252A1 (en) 2000-06-08
WO2000033252A9 (en) 2002-08-22
EP1048008B1 (de) 2006-06-07
JP4367884B2 (ja) 2009-11-18
US6233304B1 (en) 2001-05-15
EP1048008A1 (de) 2000-11-02

Similar Documents

Publication Publication Date Title
DE69931750T2 (de) Verfahren und gerät zur verkalkungsmessung
DE69934475T2 (de) Verfahren und Gerät zur Bewertung des Herzens mit Hilfe eines Mehrstrahlabtasters
DE60034748T2 (de) Verfahren und Vorrichtung zur bewegungsfreien kardiologischen Computertomographie
DE60024907T2 (de) Verfahren und Vorrichtung zur Kompensation von Artefakten mittels veränderlicher Winkelabtastung
DE60224770T2 (de) Verfahren und Vorrichtung zur Rauschverminderung in Computertomographen
DE60014001T2 (de) Verfahren und Vorrichtung zur auf Scout basierten Verkalkungsmessung
DE69935413T2 (de) Herzsignal-triggerung eines computertomographen
DE60026681T2 (de) Verfahren und Gerät zur Mehrschichttomografie mit partieller Abtastung
DE10356116A1 (de) Verfahren und Vorrichtung zum Erleichtern eines Verringerns von Artefakten
DE60219097T2 (de) Verfahren und Vorrichtung zur bewegungsfreien Computertomographie
DE10355383A1 (de) Verfahren und Vorrichtung zum Erfassen von Perfusionsdaten
DE19921116A1 (de) Artefaktkorrektur für stark dämpfende Objekte
DE10043725A1 (de) Verfahren und Vorrichtung zur Zwei-Durchlauf-Kegelstrahlbildrekonstruktion
DE10064785A1 (de) Verfahren und Vorrichtung zum Schätzen einer Herzbewegung unter Verwendung von Projektionsdaten
DE10211581A1 (de) Verfahren und Vorrichtung unter Verwendung eines verallgemeinerten spiralenförmigen Interpolationsalgorithmus
DE4421845A1 (de) Adaptive Filter zur Verminderung von Streifen-Artefakten in Röntgen-Tomographiebildern
DE19743220A1 (de) Verfahren und Vorrichtung zur dreidimensionalen Bildrekonstruktion mit maximaler Intensitätsprojektion bei einem Computer-Tomographie-System
DE3689166T2 (de) Substanzquantifizierung in Tierkörpern.
DE102005005839A1 (de) Verfahren und Vorrichtung zur Verminderung von Artefakten bei Kegelstrahl-CT-Bildrekonstruktionen
DE19743217A1 (de) Verfahren und Gerät zur Verringerung von Teilvolumen-Bildartefakten
DE10393159T5 (de) Verfahren und Anordnung zum medizinischen Röntgen
DE102004048402A1 (de) Verfahren und Vorrichtung für Bildgebung mittels dynamischen Spiralscannens
DE19635017A1 (de) Verfahren und Vorrichtung zur Erhöhung der Bildschärfe
DE19650529A1 (de) Verfahren und Vorrichtungen zur Mehrfach-Schnitt-Wendel-Bildrekonstruktion in einem Computer-Tomographie-System
DE102004029474A1 (de) System und Verfahren zum Scannen eines Objekts in Tomosynthese-Anwendungen

Legal Events

Date Code Title Description
8364 No opposition during term of opposition
8339 Ceased/non-payment of the annual fee