Geautomatiseerde Midline Shift en intracraniale druk Schatting op basis van Brain CT afbeeldingen

Summary

Een geautomatiseerde middellijn shift schatting en intracraniële druk (ICP) pre-screening op basis van computertomografie (CT) beelden voor patiënten met traumatisch hersenletsel (TBI) wordt voorgesteld met behulp van beeldverwerking en machine learning technieken.

Abstract

In dit artikel presenteren we een geautomatiseerd systeem vooral gebaseerd op de computertomografie (CT) beelden, bestaande uit twee hoofdonderdelen: de middellijn shift schatting en intracraniële druk (ICP) pre-screening systeem. De middellijn shift eerste raming een schatting van het ideale middellijn uitgevoerd op basis van de symmetrie van de schedel en anatomische kenmerken in de hersenen CT scan. Vervolgens wordt segmentatie van de ventrikels van de CT-scan en gebruikt als richtlijn voor de identificatie van de werkelijke middellijn door matching vorm. Deze processen bootsen de meting door artsen en hebben veelbelovende resultaten bij de evaluatie. In de tweede component, meer mogelijkheden geëxtraheerd verband met ICP, zoals de textuurinformatie, bloed hoeveelheid van CT scans en andere opgenomen kenmerken, zoals leeftijd, letselernst score van de ICP schatten zijn ook opgenomen. Machine learning technieken, waaronder feature selectie en classificatie, zoals ondersteuningVector Machines (SVMs), worden gebruikt om het voorspellingsmodel met RapidMiner te bouwen. De evaluatie van de voorspelling toont potentiële nut van het model. De geschatte ideale middellijn shift en voorspelde ICP niveaus kan worden gebruikt als een snelle pre-screening stap voor artsen om beslissingen te nemen, om zo aan te bevelen voor of tegen invasieve ICP monitoring.

Introduction

Elk jaar zijn er ongeveer 1,4 miljoen traumatisch hersenletsel (TBI) gerelateerde afdeling spoedeisende gevallen in de Verenigde Staten, waarvan meer dan 50.000 de dood tot gevolg ^1. Ernstige TBI gaat meestal gepaard met een toename van de intracraniale druk (ICP) met symptomen zoals hematomen en zwelling hersenweefsel. Deze resulteren in een verminderde cerebrale perfusie druk en cerebrale doorbloeding, het plaatsen van de gewonde hersenen in extra risico. Ernstige ICP stijging fataal kan zijn, zodat de controle ICP voor patiënten met TBI is cruciaal. Dit vereist typisch plaatsing van katheters direct in de hersenen voor de monitoring van druk, een riskante procedure voor patiënten die alleen kan worden uitgevoerd in gespecialiseerde medische centra. Deze procedure houdt ook risico zoals een infectie. Echter, sommige tekenen van verhoogde ICP merkbaar zijn in de medische beeldvorming. In het bijzonder wordt middellijn shift vaak geassocieerd met een verhoging van de ICP en kunnen worden overgebracht naar de hersenen berekende tomography (CT) beelden. Als zodanig zijn deze beelden een gelegenheid voor niet-invasieve detectie van verhoogde ICP die kunnen worden gebruikt als pre-screening stap voor craniale trepanatie. CT is nog steeds de gouden standaard voor TBI eerste beoordeling van alle andere beeldvormingsmodaliteiten zoals MRI, vanwege de hoge snelheid en relatief lage kosten ^2. Bovendien is een CT-onderzoek niet strikte patiënt immobiliteit, en heeft voordeel bij het onthullen van ernstige afwijkingen, zoals botbreuken en hematomen. Terwijl CT wordt gebruikt voor detectie van letsels in de hersenen, op basis van de huidige technologie is middellijn verschuiving niet automatisch gemeten en derhalve artsen moeten deze belangrijke factor bepalen door visuele inspectie. Onnauwkeurig of niet CT interpretatie wordt vaak geassocieerd met de aard van het menselijke visuele systeem en de complexe structuur van de hersenen. Terwijl kleine middellijn verschuivingen zijn ongrijpbaar, ze zijn vaak van onschatbare waarde voor de beoordeling van hersenletsel, in particulaire vroeg stadium van schade voordat toestand van de patiënt ernstiger. Aan de andere kant van het spectrum, grote middellijn shift suggereert sterk verhoogde ICP en ernstige TBI. Het is echter een zeer uitdagende taak voor mensen visueel inspecteren CT beelden en kwantitatief te voorspellen het niveau van ICP. Door de vooruitgang in geautomatiseerde computationele technieken kunnen features uit CT beelden, zoals middellijn verschuiving hematoom volume en textuur van hersenen CT beelden, nauwkeurig worden gemeten en automatisch met behulp van geavanceerde beeldverwerkingsmethoden. De relatie tussen ICP en middellijn shift en andere functies zoals de mate van bloeding, is de structuur van CT beelden niet onderzocht. In dit artikel is een computationeel raamwerk voorgesteld aan de middellijn shift meting en andere fysiologische / anatomische kenmerken op de hersenen CT-beelden vervolgens te meten en de mate van ICP niet-opdringerig met behulp van machine learning technieken voorspellen.

Protocol

1. Methodologie Overzicht

Het voorgestelde kader verwerkt de hersenen CT beelden van traumatisch hersenletsel (TBI) patiënten automatisch berekenen middellijn verschuiving in pathologische gevallen en evenals andere geëxtraheerde informatie gebruiken om intracraniale druk (ICP) te voorspellen. Figuur 1 toont het schema van de gehele kader. De geautomatiseerde middellijn shift meting kan worden onderverdeeld in drie stappen. Ten eerste, de ideale middellijn van de hersenen, de middellijn voor letsel dat wil zeggen, is te vinden via een hiërarchische zoeken op basis van de schedel symmetrie en weefsel functies ^3. Ten tweede, het ventriculaire systeem is gesegmenteerd voor elk hersenen CT-beeld ^4. Ten derde is de werkelijke middellijn op basis van de gesegmenteerde vervormd ventriculaire systeem met een vorm matching methode ^5. De horizontale verschuiving van het ventriculaire systeem wordt dan geschat op basis van de schatting van de ideale middellijn en de werkelijke middellijn. Na thij middellijn shift achtereenvolgens geschat functies zoals middellijn shift, textuurinformatie van CT beelden, en andere demografische gegevens worden gebruikt om ICP voorspellen. Machine learning algoritmen worden gebruikt om de relatie tussen de ICP en de geëxtraheerde kenmerken ⁶ model.

2. Ideaal Midline Schatting

1. Deze stap detecteert de geschatte ideale middellijn met behulp van symmetrie van de schedel. Eerst met grijstinten Thresholding, de schedel moet gesegmenteerd van de rest van de inhoud van het CT beeld. Waarna voeren een grondig onderzoek naar de rotatiehoeken in het massamiddelpunt van de schedel identificeren. De optimale rotatiehoek wordt gedefinieerd als de hoek die de symmetrie van de verkregen helften van de schedel maximaliseert. De geschatte ideale middellijn is de lijn die door het massamiddelpunt punt en heeft de optimale rotatie angel ten opzichte van de oorspronkelijke verticale richting van het CT beeld.
2. Deze stap detecteert de posterior falx cerebri anterior en falx gehechtheid aan de marge van de sagittale sulcus. Deze stap wordt uitgevoerd om de geïdentificeerde benadering ideale middellijn verfijnen. Eerste definiëren twee zoeken rechthoeken gecentreerd op de twee snijpunten tussen de geschatte ideale middellijn en de calvarium. Kies vervolgens de grootte van de rechthoeken empirisch zodat ze betrekking hebben op de anatomische kenmerken te detecteren zoals hieronder uiteengezet. De voorste falx beslag wordt gedetecteerd als de piek punt van de nok op de calvarium en de falx cerebri wordt gedetecteerd als de grijze lijn in het posterieure gebied ^3.
3. Deze stap maakt gebruik van deze functies boven gedetecteerd om de ideale middellijn positie te verfijnen. Zodra de piek punt van de voorste falx gehechtheid en het verste punt in het achterste falx cerebri van de calvarum zijn gespecificeerd, de verfijnde ideale middellijn is de lijn tussen de twee punten.

3. Ventrikel Segmentatie

1. Breng eerst een laag-level segmentatie met behulp van een Gauss-mengsel model (GGM) voor elke CT slice ^4,7. CT-beelden kunnen worden onderverdeeld in 4 soorten weefsel: bot / bloed, cerebrospinale vloeistof (CSF), grijs en witte stof. Voor het initialiseren van de Gaussian Mixture Model worden de parameters geschat op basis van een iteratieve K-means segmentatie resultaat van het CT beeld. Maak dan gebruik van de verwachting-maximalisatie (EM) methode om iteratief te optimaliseren GGM beter vertegenwoordigen de CT-beeld. De harde segmentatie resultaat kan worden verkregen door afscherming van de CT beeld in gebieden gebaseerd op de maximale waarschijnlijkheid van het lidmaatschap dat elke pixel behoort tot verschillende types regio.
2. Na het lage segmentatie van GMM, bijkomende beperkingen gelden de segmentatie resultaat om ventriculaire regio herkennen. Alleen ventriculaire gebieden met grootte boven een bepaalde drempel blijven behouden. Beperkingen van toepassing op de locatie van ventriculaire regio en met de hersenen kader en de reeks ventricle sjablonen. Pak de set van ventrikel sjablonen van een standaard hersenen MR beelden en dan vergroten met behulp van de morfologische dilatatie van variaties tussen de verschillende onderwerpen en pathologische gevallen tegemoet te komen.

4. Werkelijke Midline Schatting

1. Geef de functie punten op de ventrikel sjablonen uit MR.
2. Stel de Meerdere Regio's Shape Matching ^5,8 tussen gesegmenteerde ventrikels en de MR sjabloon.
3. Raming van de werkelijk middellijn gebaseerd op de geïdentificeerde kenmerkende punten op de ventrikel vormen die vorm matching. Dan het gemiddelde van de gemiddelde linker en rechterzijde gemiddelde van de x-coördinaten van de kenmerkende punten van de x-coördinaat van de middellijn definiëren.

5. Meer functie Extracties

1. Meet de intracraniële hematoom / bloedingen volume op basis van het Gaussiaans mengsel model (GMM) segmentatie resultaten van de CT-beelden. De gesegmenteerde resultaat kan include kleine regio's, waaronder veneuze sinus bloed en falx cerebri, maar ze kunnen meestal worden verwaarloosd ten opzichte van de regio's van hematoom. Dan tel het aantal pixels geclassificeerd als bloed voor elk segment en vatten ze op. De uiteindelijke som kwantificeert de extravasatie van bloed volume in de CT-onderzoek.
2. Extract textuur eigenschappen. Selecteer eerst zes ramen in elke CT-beeld, dat de grijze en witte stof vangt, maar voorkomen dat het bloed en de ventrikel structuren in het CT-beeld, zie Figuur 6. Pak dan de overeenkomstige structuur functies met behulp van de volgende methoden: Gray Level Run lengte ^9, Histogram analyse, Fourier analyse, Dual Tree Complex Wavelet analyse ^10.
3. Demografische gegevens en letselernst score worden verzameld.
4. Alle gewonnen kenmerken van elke CT-beeld worden samengevoegd om de volledige CT-onderzoek vertegenwoordigen. In het bijzonder, min (f), max (f), mediaan (f), gemiddelde (f), std (f) worden berekend onder all de geselecteerde functies die behoren tot de specifieke CT-onderzoek, bijvoorbeeld een functie ƒ van de middellijn shift of een textuur functie. Voor de intracraniale bloeding bedrag kenmerk, naast de 5 operators hierboven genoemde som (f) wordt ook toegevoegd aan het totale bloedvolume te nemen.

6. ICP Schatting

Het belangrijkste idee van ICP schatting is machine learning technieken toe te passen om een ​​model op basis van een set van training monsters op te bouwen. Is de ingebouwde model wordt geëvalueerd op de resterende monsters. Vanwege de hoge dimensie van geëxtraheerde kenmerken, zoals die van de CT-scans en demografische informatie, feature selectie is van belang voor niet-verbonden functies te verwijderen voor een relatief eenvoudige manier stabiel model. Daarom zijn er twee stappen die moeten worden uitgevoerd voor ICP schatting / voorspelling. Selecteer eerst de relatieve kenmerken die zijn informatief in het voorspellen van ISP's. De tweede stap is het Support Vector Machines (SVM) te gebruiken als deleeralgoritme te ontwikkelen en evaluatie van de opleidings-model. Software zoals RapidMiner ¹¹ is ideaal voor deze taak, want het is een zeer goed ontwikkeld instrument voor de meeste van machine learning algoritmen en biedt een zeer krachtige interfaces om te trainen en te evalueren modellen.

1. Voer de functie selectie in twee fasen. Gebruik eerst de informatie winst verhouding criterium in de eerste etappe naar de top 50 functies te selecteren. In de tweede fase met een genetisch algoritme waarin SVM verder optimaliseren van de functie selectie. Maak dan gebruik van de laatste set van de geselecteerde functies om het model voor ICP voorspelling te bouwen in de volgende stap.
2. Voer ICP classificatie en evaluatie door middel van machine learning technieken. Het hoogste niveau diagram van de opleiding en het testen zijn weergegeven in RapidMiner in figuur 8. Breng een 10-voudige cross-validatie regeling als de buitenste laag voor evaluatie. Om een ​​model te bouwen met een betere generalisatie en voorkomen dat over-fitting naar thij trainingsgegevens, nest een andere laag van cross-validatie in elke training vouw. In deze modules, eerst van toepassing de functie selectieproces hierboven beschreven en vervolgens een SVM te gebruiken voor classificatie met een eigen parameter selectie module ^6.

Results

Het testen CT datasets werden verstrekt door de Carolina Healthcare System (CHS) onder Institutional Review Board goedkeuring. Alle proefpersonen werden gediagnosticeerd met een lichte tot ernstige TBI bij de eerste opname in het ziekenhuis. Voor elke patiënt werd de ICP-waarde per uur met ICP probes in de ventrikel gebied voor en na CT scans werden verkregen. De ICP waarde koppelen aan elke CT scan, gemiddelde de twee grootste afmetingen van ICP de tijd van CT scan, die beide binnen een uur na de CT scan. Wijs de gemid...

Discussion

In deze studie wordt een intuïtieve en flexibel framework voorgesteld om twee uitdagende problemen aan te pakken: de schatting van de middellijn verschuiving in CT-beelden en ICP-niveau voorspelling op basis van geëxtraheerde kenmerken. De evaluatie resultaten tonen de effectiviteit van de voorgestelde werkwijze. Voor zover wij weten is dit de eerste keer dat van een systematische studie in de aanpak van deze twee problemen. We merken dat op basis van het algemene kader, zijn er vele mogelijke verbeteringen die kunnen w...