DE69931750T2 - Verfahren und gerät zur verkalkungsmessung - Google Patents
Verfahren und gerät zur verkalkungsmessung Download PDFInfo
- 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
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: 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 in1 veranschaulichten Systems. -
3 zeigt eine anschauliche Ansicht von Spalten zwischen Bilddaten, gemäß einem Ausführungsbeispiel der vorliegenden Erfindung. - DETAILLIERTE BESCHREIBUNG DER ERFINDUNG
- In
1 und2 ist ein Computertomographie-(CT)-Bildgebungssystem10 gezeigt, das einen Gantryrahmen12 aufweist, wie er für CT-Scanner der "dritten Generation" typisch ist. Die Gantry12 weist eine Röntgenstrahlenquelle14 auf, die ein Bündel von Röntgenstrahlen16 in Richtung eines auf der entgegengesetzten Seite der Gantry12 angeordneten Detektorarrays18 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 Detektormatrix18 wird durch Detektorelemente20 gebildet, die die projizierten Röntgenstrahlen, die einen Patienten22 durchqueren, gemeinsam erfassen. Die Detektormatrix20 kann ein Einzelschicht-Detektor oder ein Multischicht-Detektor sein. Jedes Detektorelement20 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 Patienten22 erfährt. Während des Scannens zum Gewinnen von Röntgenstrahlprojektionsdaten kreisen die Gantry12 und die daran befestigten Komponenten um eine Rotationsachse24 . - Die Rotation des Gantryrahmens
12 und der Betrieb der Röntgenstrahlenquelle14 werden durch eine Steuervorrichtung26 des CT-Systems10 gesteuert. Der Steuermechanismus26 enthält einen Röntgencontroller28 , der Energie und Zeittaktsignale an die Röntgenstrahlenquelle14 liefert, und einen Gantryantriebscontroller30 , der die Rotationsgeschwindigkeit und die Position der Gantry12 steuert. Ein Datenakquirierungssystem (DAS)32 in dem Steuermechanismus26 tastet von Detektorelementen20 stammende analoge Daten ab und wandelt die Daten in digitale Signale um, um diese anschließend zu verarbeiten. Ein Bildrekonstruktor34 nimmt die von dem DAS32 ankommenden abgetasteten und digitalisierten Röntgenstrahldaten entgegen und führt eine Hochgeschwindigkeits-Bildrekonstruktion durch. Das rekonstruierte Bild wird einem Computer36 als Eingabe zugeführt, der das Bild in einem Massenspeichergerät38 speichert. - Der Rechner
36 nimmt ferner über eine Konsole40 , die eine Tastatur aufweist, von einem Anwender Steuerbefehle und Scanparameter entgegen. Ein zugehöriges Datensichtgerät42 ermöglicht es dem Anwender, das rekonstruierte Bild und sonstige von dem Rechner36 ausgegebene Daten zu beobachten. Die durch den Anwender eingegebenen Steuerbefehle und Parameter werden von dem Rechner36 verwendet, um Steuerungssignale und Daten an das DAS32 , den Röntgencontroller28 und den Gantryantriebscontroller30 auszugeben. Darüber hinaus betätigt der Rechner36 einen Liegenantriebscontroller44 , der eine motorisch angetriebene Liege46 steuert, um den Patienten22 innerhalb des Gantryrahmens12 zu positionieren. Insbesondere bewegt die Liege46 den Patienten22 abschnittsweise durch den Gantrytunnel48 . - 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 Bildrekonstruktor34 implementiert sein und dem Computer36 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 Patienten22 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: 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.
-
- 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: 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)
- 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.
- 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.
- Verfahren nach Anspruch 1, bei dem das Ermitteln des Schichtbildkalziumwerts für jeden interessierenden Bereich durchgeführt wird anhand der Gleichung: 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.
- Verfahren nach Anspruch 1, bei dem das Ermitteln des gewichtetes Gesamtkalziumwerts TCS durchgeführt wird nach der Formel: 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.
- 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. - 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.
- System nach Anspruch 6, bei dem die Bestimmung des Schichtbildkalziumwerts für jeden interessierenden Bereich nach der folgenden Formel durchgeführt wird: 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.
- System nach Anspruch 6, bei dem der gewichtete Gesamtkalziumwert bestimmt wird nach der Formel: 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.
- System nach Anspruch 6, bei dem die Detektormatrix ein Multischicht-Detektorarray ist.
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)
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)
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 |
-
1998
- 1998-11-25 US US09/199,513 patent/US6233304B1/en not_active Expired - Fee Related
-
1999
- 1999-11-23 IL IL13741399A patent/IL137413A0/xx not_active IP Right Cessation
- 1999-11-23 JP JP2000585825A patent/JP4367884B2/ja not_active Expired - Fee Related
- 1999-11-23 WO PCT/US1999/028072 patent/WO2000033252A1/en active IP Right Grant
- 1999-11-23 DE DE69931750T patent/DE69931750T2/de not_active Expired - Fee Related
- 1999-11-23 EP EP99962882A patent/EP1048008B1/de not_active Expired - Lifetime
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 |