Automated Midline Shift and Intracranial Pressure Estimation based on Brain CT Images

Podsumowanie

An automated midline shift estimation and intracranial pressure (ICP) pre-screening system based on computed tomography (CT) images for patients with traumatic brain injury (TBI) is proposed using image processing and machine learning techniques.

Streszczenie

In this paper we present an automated system based mainly on the computed tomography (CT) images consisting of two main components: the midline shift estimation and intracranial pressure (ICP) pre-screening system. To estimate the midline shift, first an estimation of the ideal midline is performed based on the symmetry of the skull and anatomical features in the brain CT scan. Then, segmentation of the ventricles from the CT scan is performed and used as a guide for the identification of the actual midline through shape matching. These processes mimic the measuring process by physicians and have shown promising results in the evaluation. In the second component, more features are extracted related to ICP, such as the texture information, blood amount from CT scans and other recorded features, such as age, injury severity score to estimate the ICP are also incorporated. Machine learning techniques including feature selection and classification, such as Support Vector Machines (SVMs), are employed to build the prediction model using RapidMiner. The evaluation of the prediction shows potential usefulness of the model. The estimated ideal midline shift and predicted ICP levels may be used as a fast pre-screening step for physicians to make decisions, so as to recommend for or against invasive ICP monitoring.

Wprowadzenie

Every year there are about 1.4 million traumatic brain injuries (TBI) related emergency department cases in the United States, of which, over 50,000 result in death¹. Severe TBI is usually accompanied by an increase in intracranial pressure (ICP) with symptoms such as hematomas and swelling brain tissue. These result in reduced cerebral perfusion pressure and cerebral blood flow, placing the injured brain in additional risk. Severe ICP increase can be fatal, so monitoring ICP for patients with TBI is crucial. This typically requires placement of indwelling catheters directly into the brain for monitoring of pressure, a risky procedure for patients that can only be performed at specialized medical centers. The procedure also involves risk such as infection. However, some signs of elevated ICP may be noticeable in medical imaging. In particular, midline shift is often associated with an increase in the ICP and can be captured from the brain computed tomography (CT) images. As such, these images provide an opportunity for non-invasive detection of elevated ICP which can be used as a pre-screening step before cranial trepanation. CT imaging is still the gold standard for initial TBI assessment among all other imaging modalities, e.g. MRI, because of its high speed and relative low cost². Furthermore, a CT examination does not require strict patient immobility, and has advantage in revealing severe abnormalities such as bone fractures and hematomas. While CT is commonly used for detection of injuries in the brain, based on the current technology, midline shift is not automatically measured and therefore physicians must assess this important factor by visual inspection. Inaccurate or inconsistent CT interpretation is often associated with the nature of the human visual system and the complex structure of the brain. While small midline shifts are elusive, they are often invaluable for assessment of brain injury, in particular at early stages of injury before a patient's condition becomes more severe. On the other side of the spectrum, large midline shift suggests highly elevated ICP and more severe TBI. However, it is a very challenging task for humans to visually inspect CT images and predict the level of ICP quantitatively. Due to advances in automated computational techniques, features extracted from CT images, such as midline shift, hematoma volume, and texture of brain CT images, can be measured accurately and automatically using advanced image processing methods. However, the relationship between ICP and midline shift as well as other features such as degree of bleeding, the texture from CT images has not been explored. In this paper, a computational framework has been proposed to measure the midline shift measurement as well as other physiological / anatomical features on brain CT images and then predict the degree of ICP non-intrusively using machine learning techniques.

Protokół

1. Methodology Overview

The proposed framework processes the brain CT images of traumatic brain injury (TBI) patients to automatically calculate midline shift in pathological cases and use it as well as other extracted information to predict intracranial pressure (ICP). Figure 1 shows the schematic diagram of the entire framework. The automated midline shift measurement can be divided into three steps. First, the ideal midline of the brain, i.e. the midline before injury, is found via a hierarchical search based on skull symmetry and tissue features³. Secondly, the ventricular system is segmented for each brain CT image⁴. Thirdly, the actual midline is estimated from the segmented deformed ventricular system using a shape matching method⁵. The horizontal shift of the ventricular system is then estimated based on the estimation of the ideal midline and the actual midline. After the midline shift is successively estimated, features including midline shift, texture information of CT images, as well as other demographic information are used to predict ICP. Machine learning algorithms are used to model the relationship between the ICP and the extracted features⁶.

2. Ideal Midline Estimation

1. This step detects the approximate ideal midline using symmetry of the skull. First, using grayscale thresholding, the skull needs to be segmented from the rest of the content in the CT image. After which, perform an exhaustive search to identify the rotation angles around the mass center of the skull. The optimal rotation angle is defined as the angle that maximizes the symmetry of the resulting halves of the skull. The approximate ideal midline is the line passing through the mass center point and has the optimal rotation angel with respect to the original vertical direction of the CT image.
2. This step detects the posterior falx cerebri and anterior falx attachment to the margins of the sagittal sulcus. This step is performed to refine the identified approximate ideal midline. First define two searching rectangles centered at the two intersection points between the approximate ideal midline and the calvarium. Next, choose the sizes of the rectangles empirically so that they cover the anatomical features to be detected as explained below. The anterior falx attachment is detected as the peak point of the ridge on the calvarium and the falx cerebri is detected as the gray line in the posterior region³.
3. This step uses those features detected above to refine the ideal midline position. Once the peak point of the anterior falx attachment and the furthest point in the posterior falx cerebri from the calvarum are specified, the refined ideal midline is the line connecting the two points.

3. Ventricle Segmentation

1. First apply a low-level segmentation using a Gaussian mixture model (GMM) for each CT slice^4,7. CT images can be divided into 4 types of tissue: bone/blood, cerebrospinal fluid (CSF), gray matter and white matter. For the initialization of the Gaussian Mixture Model, the parameters are estimated based on an iterative K-means segmentation result of the CT image. Then use the expectation-maximization (EM) method to optimize GMM iteratively to better represent the CT image. The hard segmentation result can be obtained by partitioning the CT image into regions based on the maximum probability of the membership that each pixel belongs to different region types.
2. After the low-level segmentation of GMM, apply further constraints on the segmentation result in order to recognize ventricular regions. Only ventricular regions with size above a certain threshold are retained. Apply constraints on the location of ventricular regions as well using the brain bounding box and the set of ventricle templates. Extract the set of ventricle templates from a standard brain MR images and then enlarge it using the morphological dilation to accommodate variations among different subjects and pathological cases.

4. Actual Midline Estimation

1. Specify the feature points on the ventricle templates extracted from MR.
2. Do the Multiple Regions Shape Matching^5,8 between segmented ventricles and the MR template.
3. Estimate the actual midline based on the identified feature points on the ventricle shapes using shape matching. Then use the average of left side mean and right side mean of the x-coordinates of those feature points to define the x-coordinate of the midline.

5. More Feature Extractions

1. Measure the intracranial hematoma/bleeding volume based on the Gaussian mixture model (GMM) segmentation results obtained from the CT images. The segmented result may include small regions including venous sinus blood and falx cerebri, but they can usually be neglected in comparison to the regions of hematoma. Then count the number of pixels classified as blood for each slice and sum them up. The final sum quantifies the extravasated blood volume in the CT examination.
2. Extract texture features. First select six windows in each CT image which captures the gray and white matter but avoid the blood and ventricle structures in the CT image, see Figure 6. Then extract the corresponding texture features by using the following methods: Gray Level Run length⁹, Histogram analysis, Fourier analysis, Dual Tree Complex Wavelet analysis¹⁰.
3. Demographic information and injury severity score are also collected.
4. All extracted features of each CT image are aggregated to represent the entire CT examination. Specifically, min(f), max(f), median(f), mean(f), std(f) are calculated among all the selected features belonging to the particular CT examination, e.g. a feature ƒ of the midline shift or a texture feature. For the intracranial hemorrhage amount feature, besides the 5 operators listed above, sum(f) is also added to record the total blood volume.

6. ICP Estimation

The main idea of ICP estimation is to apply machine learning techniques to build a model based on a set of training samples. Then the built model is evaluated on the remaining test samples. Because of the high dimension of extracted features including those from the CT scans and demographic information, feature selection is important to remove unrelated features for a relatively simple thus stable model. Therefore there are two steps to be performed for ICP estimation/prediction. First, select the relative features which are informative in predicting ICPs. The second step is to use Support Vector Machines (SVM) as the learning algorithm to develop and evaluate the training model. Software such as RapidMiner¹¹ is ideal for this task because it is a very well developed tool for most of machine learning algorithms and provides very powerful interfaces to train and evaluate models.

1. Perform the feature selection in two stages. First use the information gain ratio criterion in the first stage to select the top 50 features. In the second stage, use a genetic algorithm incorporating SVM to further optimize the feature selection. Then use the final set of selected features to build the model for ICP prediction in the following step.
2. Perform ICP classification and evaluation through machine learning techniques. The top level diagram of the training and testing are shown in RapidMiner in Figure 8. Apply a 10 fold cross-validation scheme as the outmost layer for evaluation. In order to build a model with better generalization and avoid over-fitting to the training data, nest another layer of cross-validation inside each training fold. In these modules, first apply the feature selection process described above and then use a SVM for classification with its own parameter selection module⁶.

Wyniki

The testing CT datasets were provided by the Carolinas Healthcare System (CHS) under Institutional Review Board approval. All subjects were diagnosed with mild to severe TBI when first admitted to hospital. For each patient, the ICP value was recorded every hour using ICP probes inside the ventricle region both before and after CT scans were obtained. To associate the ICP value with each CT scan, average the two closest measurements of ICP to the time of CT scan, both of which are within an hour of the CT scan. Then assi...

Dyskusje

In this study, an intuitive and flexible framework is proposed to address two challenging problems: the estimation of the midline shift in CT images and ICP level prediction based on extracted features. The evaluation results show the effectiveness of the proposed method. As far as we know, this is the first time of a systematic study in addressing these two problems. We notice that based on the general framework, there are many potential improvements that can be achieved. For example, in the proposed segmentation, the l...