- Research article
- Open Access
Quantitative analysis of computed tomography images and early detection of cerebral edema for pediatric traumatic brain injury patients: retrospective study
BMC Medicine volume 12, Article number: 186 (2014)
The purpose of this study was to identify whether the distribution of Hounsfield Unit (HU) values across the intracranial area in computed tomography (CT) images can be used as an effective diagnostic tool for determining the severity of cerebral edema in pediatric traumatic brain injury (TBI) patients.
CT images, medical records and radiology reports on 70 pediatric patients were collected. Based on radiology reports and the Marshall classification, the patients were grouped as mild edema patients (n = 37) or severe edema patients (n = 33). Automated quantitative analysis using unenhanced CT images was applied to eliminate artifacts and identify the difference in HU value distribution across the intracranial area between these groups.
The proportion of pixels with HU =17 to 24 was highly correlated with the existence of severe cerebral edema (P <0.01). This proportion was also able to differentiate patients who developed delayed cerebral edema from mild TBI patients. A significant difference between deceased patients and surviving patients in terms of the HU distribution came from the proportion of pixels with HU = 19 to HU = 23 (P <0.01).
The proportion of pixels with an HU value of 17 to 24 in the entire cerebral area of a non-enhanced CT image can be an effective basis for evaluating the severity of cerebral edema. Based on this result, we propose a novel approach for the early detection of severe cerebral edema.
Increased intracranial pressure (ICP) is often observed in traumatic brain injury (TBI). Elevated ICP results in decreased cerebral perfusion pressure (CPP), which can cause cerebral ischemia. One of the most frequent effects of cerebral ischemia is cerebral edema, which is a major cause of brain swelling -. Swollen parenchymal tissue may compress capillaries, aggravating cerebral ischemia by further decreasing the cerebral blood flow (CBF) -. Computed tomography (CT) and magnetic resonance imaging (MRI) are widely used to assess the extent of brain injury in TBI. MRI offers greater resolution and sensitivity for the detection of cerebral edema than CT . By taking the apparent diffusion coefficient (ADC) in diffusion-weighted MRI (DWI), it is possible to achieve even higher sensitivity in detecting cerebral edema -. However, the higher sensitivity of MRI does not lead to significant changes in the management of TBI -. MRI also takes significantly more time than CT and requires patients to lie still during the examination, which can be particularly challenging in pediatric patients. For this reason, sedation or even generalized anesthesia for MRI is often used in pediatric patients . Two major complications of generalized anesthesia in pediatric patients are cardiac and respiratory events . Combined with the damage caused by the primary injury, these complications may further aggravate the secondary insults.
Compared to MRI, CT offers greater accessibility with faster image acquisition time, which can be particularly advantageous in acute settings with unstable patients. Despite these advantages, the utility of CT has been questioned due to the interpretation of its images, particularly in regard to the detection of cerebral edema, as it relies heavily on the radiologist’s subjective impressions of subtle differences in attenuation . To standardize CT interpretation, several CT classifications for brain injury have been devised, mainly to assess the severity of stroke. Some of these systems may be feasible for assessing the severity of TBI because ischemic stroke and TBI share similar mechanisms . More than ten TBI classification systems using CT have been proposed to date . However, no classifications have yet been as heavily evaluated as the Marshall classification . The Marshall classification uses CT findings regarding structural abnormalities, such as the compression of basal cisterns, a shift of the midline and the volume of traumatic lesion. This classification gives simple ordinal scores ranging from 1 to 6 and is known for showing good correlation with the outcome . The Marshall classification provides a significant amount of prognostic information in TBI. Nevertheless, the Marshall classification may need further refinement . The Marshall classification has been criticized for its low sensitivity to neurophysiological changes and significant intra- and inter-rater variability; most of all, it was not designed for assessing the severity of cerebral edema. Another major limit of existing CT classifications is that they may not reflect characteristic aspects of the pediatric brain.
More than 60% of pediatric TBI patients exhibit increased ICP . Immature brains are more vulnerable to increased ICP than adults as a result of their lower intracranial compliance . As a consequence, the mortality rate of pediatric patients is threefold higher than in adults . While an increase in CBF also causes brain swelling, recent studies suggest hyperemia is not a major cause of brain swelling in pediatric patients ,. In fact, decreased CBF is common in pediatric TBI patients . Cerebral edema is thus considered the major cause of increased ICP in pediatric TBI cases. The presence of cerebral edema is assessed by parenchymal hypoattenuation in a CT image. In this study, we hypothesized that the proportion of pixels with specific attenuation levels can be used to detect the presence of cerebral edema. A quantitative standard can compensate for the weakness of the Marshall classification, which relies on the subjective interpretation of images. Furthermore, certain TBI patients may later develop delayed cerebral edema. Early quantitative signs of delayed cerebral edema may exist in the initial CT images. The early detection of cerebral edema from initial CT images may be possible using the quantitative standard. We performed this study to investigate whether the proportion of pixels with specific attenuation levels in brain CT images can be used to detect the presence of cerebral edema and to develop a quantitative standard for rapid assessment of the severity of cerebral edema in pediatric patients.
The difference in density distribution between CT images from patients with severe cerebral edema and patients with mild cerebral edema was investigated. For an accurate comparison, the CT images had to be standardized. A tool for automated artifact elimination and densitometric measurement was devised. After the images were processed and the density distribution had been obtained, statistical analyses were employed to find the pixels that caused the difference between density distributions of the two groups. The proportion of these pixels in a set of CT images was defined as the cerebral edema score (CES). The CES from deceased patients was defined as mortality-related CES (mCES). Additionally, it is well known that the brain water content varies with age ,. The difference caused by age is clearly recognizable after the first two years of life ,. The age-induced difference in the density distribution was also investigated to determine whether the difference between the severe edema group and the mild edema group was affected by the age factor.
The medical records and CT images of 70 pediatric TBI patients were collected for this study (Table 1). The subjects had been admitted to the Seoul National University Hospital (SNUH), Seoul, Korea from 2003 to 2013. The use of the patient material was approved by SNUH (IRB H-1402-020-553). The primary classification of subjects was performed via radiology reports. The existence of cerebral edema and the extent of the resulting damage were confirmed by reports and CT images. Pediatric TBI patients are highly susceptible to the development of cerebral edema. Furthermore, the subtle hypoattenuation caused by ischemic change is often difficult to detect . To minimize methodological complications, detailed severity measures were not taken. Instead, this study dichotomized subjects into two groups by the Marshall classification. Marmarou et al. reported that the brain water content of TBI patients with Marshall scores I and II is close to normal . Based on their findings, we applied the Marshall Classification system to differentiate the edematous subjects into severely edematous (Marshall score III or higher) and mildly edematous (Marshall score I or II). Some subjects with no clear evidence of cerebral edema on initial CT developed severe cerebral edema 12 hours later; they were designated as delayed cerebral edema subjects and considered severely edematous.
The scans were performed on a Brilliance 64 CT scanner (Philips Healthcare, Best, Netherlands). As in many institutions, the brain images were acquired in non-helical mode , thereby effectively avoiding unnecessarily high exposure to radiation and image artifacts. Contiguous, unenhanced 5 mm-slice axial CT images were used for analysis. The images used were obtained within 24 hours after hospital arrival. The image matrix sizes were 512 by 512; X-ray tube voltage, 120 kVp; X-ray tube current, 120 to 250 mA; exposure, 120 to 250 mAs. The median CT dose index volume (CTDIvol) was 21.3 mGy (range = 16.0 to 35.4 mGy), and the median dose length product (DLP) was 372.4 mGy . cm (range = 216.7 to 805.3 mGy . cm). The reconstruction kernel was set as `UB’ (for standard brain imaging). Continuous axial sections from the vertex to the caudal end of the cerebrum were collected for each patient. Artifact elimination and densitometric measurements were performed by in-house software built using MATLAB (MathWorks Inc., Natick, MA, USA).
A tool was designed for the rapid densitometry analysis of CT images. The algorithm regards every pixel outside the cranium as an artifact and can thus quickly dispose of pixels unnecessary for the analysis (Figure 1). To analyze only the intracranial components in CT images, the pixels representing the cranium, soft tissues outside the cranium, and other artifacts, such as foreign bodies or artifacts generated by the CT machine itself, should be eliminated. The cranium pixels and foreign body pixels can be easily eliminated by using their Hounsfield Unit (HU) values. Eliminating the soft-tissue pixels outside the cranium can be more complicated, as the HU value range of soft-tissue pixels is similar to the range of intracranial component pixels. Although they can be differentiated by their relative positions (outside versus inside the cranium), employing such a method for CT images of a fractured cranium can cause an erroneous result. As a preemptive measure to avert such an error, this algorithm seals such fractures by temporarily generating pseudo-cranial pixels. Once a CT image is read (Figure 1A), the algorithm assigns four quadrants to the intracranial area (Figure 1B). These quadrants identify fractures generated by the impact of the primary injury, or craniotomy. Before proceeding to the next step, pixels representing the CT machine itself are eliminated (Figure 1C). Multiple fractures along the cranium can be developed due to the impact of the primary injury. Such fractures are defined as a discontinuation of the cranial pixels and can be sealed off by creating pseudo-cranial pixels (Figure 1D). Sizable fractures, common in CT images after craniotomy, cannot be handled by this method. The algorithm seals them off using the pseudo-axisymmetric property of the human cranium (Figure 1D). It then removes any pixels with a density outside the pre-designated boundary condition (Figure 1E). More details of the process are given in the Appendix. The algorithm allows further post-processing such as densitometric measurements and color mapping. The results are modified CT images composed of pixels depicting only parenchyma, cerebrospinal fluid (CSF) and blood (Figure 1F). The use of this algorithm allows fully automated rapid densitometry analysis for vast numbers of CT images.
The HU values of major intracranial components in a CT image generally do not exceed 80 nor fall below 0 ,. Based on this fact, attenuation threshold limits of 0 to 79 HU were used. Any pixels with an HU value outside this boundary condition were excluded from the analysis. Each pixel in a CT image has a certain HU level, λ, which is between 0 and 79 in this study. The proportion of specific pixels with λ HU in a CT image can be obtained by simple arithmetic. Let λ p be the proportion of pixels with λ HU in a CT image; then,
The proportion (%) of λ in a CT image is 0 to 100, with the sum of every λ p as 100, as per the definition of λ p . In reality, a CT examination on a TBI patient produces a series of images. If a CT examination generates k number of CT images, there will be k number of λ p for each image. Let be the proportion of pixels which have HU =λ in every pixel in the k th CT image in a given set. Then, the proportion of pixels with HU =λ in an entire pixel of the set can be expressed as = the average λ p in a number n of CT images, where = 100, and n is the total number of images generated by the CT examination. The density distribution of a subject can be obtained by plotting , for every λ. If there is a significant difference in density distribution between two sets of CT images, the difference can be expressed as the proportion of specific pixels. The proportion of pixels with HU =λ to β in a set of CT images can be easily denoted as . For convenience, is hereafter denoted as .
Statistical analyses were conducted using commercial software (SPSS version 21; IBM) to obtain the baseline characteristics of the subjects and compare averaged density distributions between the subject groups to find the difference in terms of the proportion of specific pixels. Shapiro-Wilk tests showed that normality cannot be assumed. Non-parametric statistical analyses were used for this reason. The median value, range and percentages were selected as descriptive statistics. The Kruskal-Wallis and Mann-Whitney U tests were conducted to compare the average density distributions of subject groups. With pixels having a specific HU value that caused a significant difference between the density distributions identified, the average proportions of these pixels were obtained. A receiver operator characteristics (ROC) curve determined the cutoff value of the proportions. Pearson’s Chi-square or Fisher’s exact test were applied for categorical variables. A statistically significant result was defined as a P-value of less than 0.05.
A total of 70 TBI patients’ data were used, and 42 (60%) subjects were males. The median age of the subjects was 37 months (range = 2 to 224 months); 33 (47.1%) subjects were classified as severely edematous. Among those subjects, eight (24.20%) were not initially diagnosed as edematous but later developed delayed cerebral edema (mean onset time of edema detection after admission = 36.33 hours). Six patients died during hospitalization. Every deceased patient had CT signs of cerebral edema. Only one patient had had surgery before being subjected to CT examination. The remaining details are listed in Table 1.
Comparison of mild and severe edema subjects
The average density distributions of mild (n = 37) and severe edema subjects (n = 33) were compared. A density histogram depicts the difference between these two groups (Figure 2). The overall CT image density of severe edema subjects was lower in mild edema subjects. A major difference between the two groups was the proportion of pixels with 17 to 24 HU (, P <0.0001). The median values of of the mild edema group and the severe edema group were 14.31 and 22.48, respectively. Subjects with severe edema showed a higher proportion of in their CT image than mild edema subjects. An ROC curve and error bar chart were plotted to verify the determinative power of (Figure 3). The area under the curve was 0.85, with an asymptotic P = 0.005. The cutoff point was determined to have the highest predictive power. greater than or equal to 16.03% (sensitivity = 0.82, specificity = 0.68) was highly correlated with severe cerebral edema. To further verify the credibility of as an indicator for severe edema, additional analyses were applied to confirm whether the age related changes in density distribution affected . There were significant differences in density distribution between patients at or younger than 24 months of age (n = 22, median age = 8.50 months) and patients who were older than 24 months (n = 48, median age = 53.50 months), particularly caused by the proportion of pixels with 7 to 11 HU (, P <0.05), 23 to 28 HU (, P <0.05) and 32 to 47 HU (, P <0.05). Within the severely edematous patient group and the mildly edematous patient group, comparisons regarding age difference resulted in no significant difference in (P >0.5). Tests within each age group showed that the difference in density distribution between severe edema and mild edema was caused by the proportion of pixels having HU 21 to 23 (, P <0.05) in patients at or younger than 24 months of age and HU 13 to 24 (, P <0.05) in patients older than 24 months of age.
Density distribution of subjects with delayed cerebral edema
Among 33 severe edema subjects, 25 exhibited early CT signs of cerebral edema, while the other eight subjects did not. These subjects developed delayed cerebral edema, which was confirmed by later CT examinations due to their neurological changes. Density histograms of the initial CT images from these patients were very similar to the severely edematous group. The Kruskal-Wallis test proved a significant difference between mildly TBI, the severe edema group, and subjects with delayed cerebral edema. To further test the predictive power of , a Mann-Whitney U test was applied to subjects with delayed cerebral edema and the mild edema group. Consequently, was identified as the cause of the difference between these two groups (P <0.03). Furthermore, the average proportions of in subjects with delayed cerebral edema were similar to the proportions in patients who had severe cerebral edema (Figure 4). greater than or equal to 16.29% (area under the curve = 0.75, sensitivity = 0.88, specificity = 0.68) was highly correlated with the development of delayed cerebral edema (Figure 4). The negative and positive predictive values of the threshold were 0.96 and 0.37, reflecting a low prevalence of delayed edema among the study subjects.
Comparison of surviving and deceased subjects
proved to be highly significant in differentiating deceased subjects (n = 6) from surviving subjects. Asymptotic P values were less than 0.01 for . greater than or equal to 15.27% (area under the curve = 0.85, sensitivity = 0.83, specificity = 0.78) was highly correlated with mortality (Figure 5). The peak of the average density histogram of deceased subjects was near HU24, whereas for surviving subjects it was near HU29.
Early studies suggested that immature brains were more prone to develop cerebral edema after TBI than adults ,,. Compared to adult brains, immature brains have a more fragile blood brain barrier and a higher brain water content. These characteristics may contribute to the rapid development of secondary insults from TBI . This implies that the window of successful intervention for cerebral edema can be narrower in pediatric TBI patients than in adult TBI patients. The rapid detection and anticipation of cerebral edema is even more important in pediatric TBI patients.
The difficulty of interpreting the result has always been a problem for CT imaging. Quantitative analysis for CT was devised to compensate for this difficulty. Quantitative CT (QCT) methods have beeen mainly used in other specialties, such as cardiology and pulmonology -. Recent neuroradiological studies, which used densitometry analysis, focused on finding tumors or regional pathological changes -. Only a handful of studies used densitometric methods for whole brain images of abnormal brains -. Although few in number, these studies show promising results. One of the main purposes of this study was to test the feasibility of automated densitometry analysis on CT images.
Conventional problems regarding automated CT densitometry analysis
There are some prerequisite conditions for automating CT densitometry analysis. Diagnostically irrelevant objects, such as artifacts, the cranium, catheter/ICP monitor, air cells, image background, or the CT machine itself, should be eliminated prior to the densitometry analysis. These unwanted data can complicate quantitative analysis and result in misinterpretation. The use of different CT scanners, reconstruction kernels, or adjustments of the tube voltage or currents that can affect the HU value itself in CT imaging should also be addressed . While the effect of scanner variation on quantitative CT analysis is controversial ,, the fluctuation in HU caused by acquisition parameters, such as the tube current and voltage, can be substantial . To minimize HU fluctuation in CT, the tube voltage and current must be properly adjusted . Although the HU values of soft tissues (HU in the range 0 to 100) do not show significant fluctuation, this CT-imaging characteristic can be an obstacle for automated CT image analysis . The designation of ROI is another major concern in quantitative CT analysis. While the HU fluctuation from the use of different CT acquisition parameters can be prevented, the designation of ROI may not. The manual designation of ROI can be susceptible to inter-observer variability. If pre-analytic CT images are to be standardized, the designation of ROI should be automated.
Automation of densitometric measurements through standardization
To address complications regarding automated densitometry analysis, the tool devised for this study has several distinctive functions. In CT imaging, the density of a material is defined as its relative ability to absorb X-ray photons. This ability is expressed as HU. The HU of certain important substances is well known. For instance, the HU of CSF is between 0 and approximately 15, normal cerebral tissue has an HU of less than 40, and the HU of hemorrhagic cells is below 80 ,. For the purpose of this study, any objects other than these three are considered artifacts. For this reason, attenuation threshold limits of 0 to 79 were used in this study. With these limits, artifacts were automatically eliminated. The threshold limits of 0 to 79 enabled the standardization of pre-analytic CT images without compromising the data integrity. Standardization of the CT images allows automated densitometry analysis. The tool devised for this study reads a set of images, eliminates artifacts and leaves only pixels with HU = 0 to 79. The tool then calculates the average proportions of pixels with HU 0, 1, - 79 in the set of images. The result is the density distribution of the entire cerebral area of a subject. Inter-observer variability is minimized with the automation of densitometric measurements. The algorithm employed by the tool enables continuous, automated CT densitometry analysis for groups of images, which allows the density distribution of a group of subjects to be obtained within a few minutes. After the density distributions of groups A and B are calculated, statistical analysis of the distributions can reveal the source of the significant difference between groups A and B. This difference was expressed as , the proportion of pixels with HU =λ to β in a set of CT images. By comparing proportions, the influence of HU variability was minimized.
The exact mechanisms of cerebral edema are not well understood. Understanding cerebral edema from TBI is even more challenging given its heterogeneity. Researchers have focused on the molecular mechanisms of cerebral edema -, which are often difficult to translate into a clinical setting . Currently, the Brain Trauma Foundation recommends ICP-oriented therapy for pediatric TBI patients. Repeat CT scans are considered when ICP therapy does not seem to be effective . Brain CT densitometry analysis can provide a rapid assessment of pathological changes in an abnormal brain.
Quantitative severity assessment for cerebral edema
While there are notable quantitative brain CT analysis studies on cerebral edema or TBI ,, most prior studies on the quantitative assessment of cerebral edema used mean HU values in a pre-designated edematous region as the subject of analysis. However, such approach may not have the desired diagnostic value for several reasons. The designation of ROI is susceptible to inter-observer variability. Inaccurate segmentation of an edematous lesion would result in erroneous interpretation of the CT image. Moreover, the designation of ROI itself can be complicated in diffuse cerebral edema. The `mean HU value of the ROI’ approach, for these reasons, may not be appropriate for quantitative CT analysis of pediatric cerebral edema. Alternatively, other studies have used densitometry analysis, based on predefined values, on pixels of edematous tissue. Rosza et al. separated diffuse brain swelling patients into edema and hyperemia groups according to the number of pixels of specific HU . In this study, pixels having an HU between 11 and 20 were considered edematous pixels. Mangel et al. made a similar assumption, where an edematous pixel was defined as having a HU of 10 to 28 . While intuitive, results from these predefined definitions may not reflect the actual damage caused by cerebral edema. This inherent limitation becomes more problematic when applied to pediatric patients, as these patients have a higher brain water content than adults. To address these problems, standardization of CT images and whole brain analysis were employed in this study. Rather than analyzing pre-designated pixels of interest, we compared the density distributions of the whole cerebrum for severely edematous patients and mildly or moderately edematous patients. This method revealed the actual difference between CT images of severely edematous brains and mildly edematous brains, allowing a standardized global assessment of cerebral edema.
Cerebral edema score
Subjects with severe cerebral edema showed a lower overall density than mild edema subjects. The average density histogram for the two groups clearly illustrates this difference (Figure 2). Compared to the mild edema group, subjects with severe cerebral edema exhibited an increased proportion of pixels having HU 17 to 24 () in an entire set of CT images for the CT examination.
It is well known that CT attenuation and brain water content are inversely correlated . The age related difference in image attenuation is highly correlated with the progression of myelination ,. Thus, the results derived by our method (that is, the quantitative measurement of density distribution) should be affected by an age factor. As prior studies indicated, the change in brain material properties becomes clearest after the first two years of life ,. Indeed, there were significant differences in the density distribution between patients at or younger than 24 months of age and patients older than 24 months; however the two groups could not be differentiated by . Similar results were given within the severely edematous group and mildly edematous group. The results indicate that age does affect the CT density distribution; however, it has little impact on the utility of as an indicator of severe edema, thus allowing the proposed method to be effective. The average HU of human edematous cerebral tissue varies from 18 to 29 . The numerical value of , therefore, implies the severity of insult from cerebral edema: 16% or greater can be considered as a threshold for differentiating mildly edematous subjects from the severely edematous group (Figure 3). Given the definition of the two groups (the existence of severe edema), this threshold can be inferred as a potential quantitative standard for cerebral edema.
If can be used to detect the presence of severe cerebral edema, the remaining question is whether it can be used to predict the development of severe cerebral edema. Delayed cerebral edema can cause delayed intracranial hypertension, which can be fatal. There were eight subjects who had no signs of severe cerebral edema in their initial CT images but later developed severe cerebral edema. If the images from the initial CT examinations of those eight subjects share a similar density distribution with the rest of the severe edema subjects, the early detection of cerebral edema would be possible. The basis for this assumption is that the initial non-enhanced brain CT images possess significant prognostic value regarding TBI . Tests on subjects revealed that was able to differentiate the two groups successfully (Figures 4 and 6). Apparent signs of cerebral edema in a CT image become visible several hours after TBI in delayed edema subjects. However, the abnormal density distribution, expressed in HU, was already observable in their initial CT images (Figure 6). Although cytotoxic edema causes major damage in TBI, the initial edematous insult is vasogenic . Evolving vasogenic edema affects the white matter. According to an earlier study, a white matter density decrease of more than 5 HU is associated with symptoms of cerebral edema . However, such a subtle decrease is nearly impossible to perceive with the human eye . Higher in the initial CT images from delayed cerebral edema subjects may reflect the early, elusive signs of vasogenic insult from TBI. This study concluded that the value of can predict or assess the severity of cerebral edema. Hence, is, hereafter, designated as the CES.
Determinative power of cerebral edema score on mortality
Pediatric mortality is highly correlated with cerebral edema. This study hypothesized, therefore, that the CES of deceased subjects would be significantly higher than the CES of surviving subjects. Five subjects died after arrival at the hospital. The average density distribution in the initial CT images of the deceased subjects was obtained to find a correlation between mortality and CES. The CES was able to differentiate the two groups, but showed an even higher statistical significance. Hence, was considered the mCES. The difference between mCES and CES was found to be small (the proportion of pixels with HU = 18 and 24). This difference can be considered as the effect of the variance caused by a small sample size. The predictive value of in initial CT images for mortality should further be tested. Nevertheless, the method used in this study proved its validity in differentiating the deceased and the surviving subjects.
Cerebral edema score as a complementary CT classification system for pediatric TBI patients
Currently, there is no CT classification system devised for pediatric patients. Instead, the Marshall classification is being employed for pediatric patients  as an alternative. However, this classification system does not consider cerebral edema. This limitation of the Marshall classification, combined with its high inter-observer variability , might result in an erroneous therapeutic approach when applied to pediatric patients. CES and mCES were developed based on CT images of pediatric TBI patients. Most importantly, these scores provide quantitative information about the level of edematous insult. These features can compensate for the previously discussed limitations of the current CT classification systems. The acquisition of CES and mCES is a simple, reproducible process that can be performed within a minute.
Limitations and suggestions
The method devised and used for this study can be a valuable research tool in studying different aspects of TBI, including the progression of cerebral edema after initial insult, correlation between the severity of cerebral edema and clinical presentation and changes in cerebral edema in response to different therapies. Automated densitometry analysis is advisable for large-scale image analysis studies. Nevertheless, several limitations exist in the findings of this study. The patient sample size is relatively small and is limited to pediatric patients. The application of the methods used for this study for adult patients remains untested. Additionally, the method proposed in the study does not solve the problem of HU fluctuations due to the use of different CT acquisition parameters or different CT scanners. Further research is necessary to address these issues. Given these limitations, we advise using this standard to supplement existing classification systems, such as the Marshall classification, for pediatric patients.
This study showed the possibility of utilizing the proportion (not number) of pixels for quantitative CT image analysis. Automated artifact elimination and densitometric measurement were introduced as tools for this analysis. Density distributions of individual patients or groups of patients can be obtained in only seconds or in a few minutes. Statistical analyses can be applied to the average density distributions between patient groups. We derived a quantitative standard (CES) for assessing the severity of cerebral edema by utilizing these methods. The proportion of pixels with HU = 17 to 24 was revealed as the key factor and, hence, was defined as the CES. A CES higher than 16% in a set of images was highly correlated with severe edema. This variable has potential as the standard for damage assessment in TBI and may even be able to serve as an early predictor of severe cerebral edema. A mCES (the proportion of pixels with HU = 19 to 23) higher than 15% can be an effective predictor for mortality. CES and mCES can be obtained by simple densitometric measurements. The acquisition of these variables can be fully automated and be performed rapidly. The application of these standards in clinical and research settings will help to better manage and understand pediatric TBI.
HSK drafted and revised the manuscript. GDK, KWK and YHC carried out the collection of the data and reviewed the manuscript. BCY, BJK, MC and BMO reviewed and participated in the revising of the initial manuscript. DJK conceived and designed the study, established the methodology and oversaw the creation of the final manuscript. All authors read and approved the final manuscript.
This section presents more details of sub-processes in the image artifact eliminating algorithm. Through a series of sub-processes, the algorithm creates an image containing only the major intracranial components (parenchyma, CSF and blood). Assigning quadrants to the intracranial area (Figure 1B), eliminating artifacts developed by the CT machine (Figure 1C), sealing the fractures across the cranium (Figure 1D) and eliminating pixels with an HU higher than 79 (Figure 1E) are the core functions of the algorithm.
Assigning quadrants to the intracranial area
All objects outside the cranium, and even the cranium itself, are considered artifacts. This definition assumes the intactness of the cranium, and hence, fractures and punctures in the cranium needed to be sealed off. This sub-process (Figure 7) exists to locate the fractures.
Elimination of pixels representing the CT machine
Pixels representing the CT machine are irrelevant to the densitometry analysis and are hence considered artifacts. To remove these pixels, the algorithm uses the pixels depicting the empty space between the CT machine and the head of the patient (Figure 8). The HU of air is approximately -1000. A set of pixels representing the CT machine is surrounded by air pixels. By this logic, the algorithm can eliminate the CT machine from a CT image.
Sealing minor fractures of the cranium
Fractures have varying sizes and locations. The algorithm separates minor and major fractures according to their sizes (Figure 9). Minor fractures are defined as less than 25 pixels in longitude and are located between the edges of the fractured cranium. The algorithm seals these fractures by creating pseudo-cranial pixels. Multiple minor fractures can exist in a single CT image, so this sub-process may be executed several times.
Sealing major fractures of the cranium
Major fractures are defined as more than 25 pixels in longitude and stretch across two quadrants. These fractures are mostly caused by craniotomy. The algorithm estimates the coordinates of the lost cranial pixels based on the coordinates of the intact cranial pixels located at the opposite sides of the horizontal or vertical axis. Pseudo-cranial pixels are created along the estimated coordinates (Figure 10).
Elimination of remaining artifacts
The sole purpose of the previous sub-processes is to create a CT image without artifacts outside the intact cranium. The final sub-process can be performed only if the cranium in the CT image is intact. The algorithm eliminates any pixels with an HU outside a predefined range of interest (in this study, 0 to 79) (Figure 11). With this last step of the automated artifact elimination algorithm completed, a CT image is standardized for later analyses.
apparent diffusion coefficient
cerebral blood flow
cerebral edema score
cerebral perfusion pressure
magnetic resonance imaging
receiver operator characteristics
region of interest
traumatic brain injury
Unterberg AW, Stover J, Kress B, Kiening KL: Edema and brain trauma. Neuroscience. 2004, 129: 1021-1029. 10.1016/j.neuroscience.2004.06.046.
Ware AJ, D’Agostino AN, Combes B: Cerebral edema: a major complication of massive hepatic necrosis. Gastroenterology. 1971, 61: 877-884.
Barzo P, Marmarou A, Fatouros P, Hayasaki K, Corwin F: Contribution of vasogenic and cellular edema to traumatic brain swelling measured by diffusion-weighted imaging. J Neurosurg. 1997, 87: 900-907. 10.3171/jns.1997.87.6.0900.
Adelson PD, Kochanek PM: Head injury in children. J Child Neurol. 1998, 13: 2-15. 10.1177/088307389801300102.
Claassen J, Carhuapoma JR, Kreiter KT, Du EY, Connolly ES, Mayer SA: Global cerebral edema after subarachnoid hemorrhage: frequency, predictors, and impact on outcome. Stroke. 2002, 33: 1225-1232. 10.1161/01.STR.0000015624.29071.1F.
Aldrich EF, Eisenberg HM, Saydjari C, Luerssen TG, Foulkes MA, Jane JA, Marshall LF, Marmarou A, Young HF: Diffuse brain swelling in severely head-injured children. A report from the NIH Traumatic Coma Data Bank. J Neurosurg. 1992, 76: 450-454. 10.3171/jns.1992.76.3.0450.
Adelson PD, Clyde B, Kochanek PM, Wisniewski SR, Marion DW, Yonas H: Cerebrovascular response in infants and young children following severe traumatic brain injury: a preliminary report. Pediatr Neurosurg. 1997, 26: 200-207. 10.1159/000121192.
Rosand J, Schwamm LH: Management of brain edema complicating stroke. J Intensive Care Med. 2001, 16: 128-141. 10.1046/j.1525-1489.2001.00128.x.
Doelken M, Lanz S, Rennert J, Alibek S, Richter G, Doerfler A: Differentiation of cytotoxic and vasogenic edema in a patient with reversible posterior leukoencephalopathy syndrome using diffusion-weighted MRI. Diagn Interv Radiol. 2007, 13: 125-128.
Marmarou A, Portella G, Barzo P, Signoretti S, Fatouros P, Beaumont A, Jiang T, Bullock R: Distinguishing between cellular and vasogenic edema in head injured patients with focal lesions using magnetic resonance imaging. Acta Neurochir Suppl. 2000, 76: 349-351.
Schaefer PW, Grant PE, Gonzalez RG: Diffusion-weighted MR imaging of the brain. Radiology. 2000, 217: 331-345. 10.1148/radiology.217.2.r00nv24331.
Hagen T, Ahlhelm F, Reiche W: Apparent diffusion coefficient in vasogenic edema and reactive astrogliosis. Neuroradiology. 2007, 49: 921-926. 10.1007/s00234-007-0272-7.
Provenzale JM: Imaging of traumatic brain injury: a review of the recent medical literature. AJR Am J Roentgenol. 2010, 194: 16-19. 10.2214/AJR.09.3687.
Manolakaki D, Velmahos GC, Spaniolas K, de Moya M, Alam HB: Early magnetic resonance imaging is unnecessary in patients with traumatic brain injury. J Trauma. 2009, 66: 1008-1012. 10.1097/TA.0b013e31819adba2. discussion 1012-1004
Kemp AM, Rajaram S, Mann M, Tempest V, Farewell D, Gawne-Cain ML, Jaspan T, Maguire S: What neuroimaging should be performed in children in whom inflicted brain injury (iBI) is suspected? A systematic review. Clin Radiol. 2009, 64: 473-483. 10.1016/j.crad.2008.11.011.
Serafini G, Zadra N: Anaesthesia for MRI in the paediatric patient. Curr Opin Anaesthesiol. 2008, 21: 499-503. 10.1097/ACO.0b013e328304115b.
Lee C, Mason L: Complications in paediatric anaesthesia. Curr Opin Anaesthesiol. 2006, 19: 262-267. 10.1097/01.aco.0000192787.93386.9c.
Zimmerman RD: Stroke wars: episode IV CT strikes back. AJNR Am J Neuroradiol. 2004, 25: 1304-1309.
Werner C, Engelhard K: Pathophysiology of traumatic brain injury. Br J Anaesth. 2007, 99: 4-9. 10.1093/bja/aem131.
Zhu GW, Wang F, Liu WG: Classification and prediction of outcome in traumatic brain injury based on computed tomographic imaging. J Int Med Res. 2009, 37: 983-995. 10.1177/147323000903700402.
Marshall LF, Marshall SB, Klauber MR, Clark MB, Eisenberg HM, Jane JA, Luerssen TG, Marmarou A, Foulkes MA: A new classification of head injury based on computerized tomography. Spec Suppl. 1991, 75: S14-S20.
Saatman KE, Duhaime AC, Bullock R, Maas AI, Valadka A, Manley GT: Classification of traumatic brain injury for targeted therapies. J Neurotrauma. 2008, 25: 719-738. 10.1089/neu.2008.0586.
Maas AI, Hukkelhoven CW, Marshall LF, Steyerberg EW: Prediction of outcome in traumatic brain injury with computed tomographic characteristics: a comparison between the computed tomographic classification and combinations of computed tomographic predictors. Neurosurgery. 2005, 57: 1173-1182. 10.1227/01.NEU.0000186013.63046.6B. discussion 1173-1182
Ward JD: Pediatric head injury: a further experience. Pediatr Neurosurg. 1994, 20: 183-185. 10.1159/000120784.
Bauer R, Walter B, Fritz H, Zwiener U: Ontogenetic aspects of traumatic brain edema-facts and suggestions. Exp Toxicol Pathol. 2009, 51: 143-150. 10.1016/S0940-2993(99)80088-8.
Zwienenberg M, Muizelaar JP: Severe pediatric head injury: the role of hyperemia revisited. J Neurotrauma. 1999, 16: 937-943. 10.1089/neu.1999.16.937.
Udomphorn Y, Armstead WM, Vavilala MS: Cerebral blood flow and autoregulation after pediatric traumatic brain injury. Pediatr Neurol. 2008, 38: 225-234. 10.1016/j.pediatrneurol.2007.09.012.
Huh JW, Raghupathi R: New concepts in treatment of pediatric traumatic brain injury. Anesthesiol Clin. 2009, 27: 213-240. 10.1016/j.anclin.2009.05.006.
Dzialowski I, Klotz E, Goericke S, Doerfler A, Forsting M, von Kummer R: Ischemic brain tissue water content: CT monitoring during middle cerebral artery occlusion and reperfusion in rats. Radiology. 2007, 243: 720-726. 10.1148/radiol.2432060137.
Marmarou A, Signoretti S, Aygok G, Fatouros P, Portella G: Traumatic brain edema in diffuse and focal injury: cellular or vasogenic?. Acta Neurochir Suppl. 2006, 96: 24-29. 10.1007/3-211-30714-1_6.
Jaffe TA, Hoang JK, Yoshizumi TT, Toncheva G, Lowry C, Ravin C: Radiation dose for routine clinical adult brain CT: Variability on different scanners at one institution. AJR Am J Roentgenol. 2010, 195: 433-438. 10.2214/AJR.09.3957.
Schneider W, Bortfeld T, Schlegel W: Correlation between CT numbers and tissue parameters needed for Monte Carlo simulations of clinical dose distributions. Phys Med Biol. 2000, 45: 459-478. 10.1088/0031-9155/45/2/314.
Creasy JL: Dating Neurological Injury; a Forensic Guide for Radiologists, Other Expert Medical Witnesses, and Attorneys.Springer New York; 2011. doi:10.1007/978-1-60761-250-6..
Bruce DA, Alavi A, Bilaniuk L, Dolinskas C, Obrist W, Uzzell B: Diffuse cerebral swelling following head injuries in children: the syndrome of ‘malignant brain edema’. J Neurosurg. 1981, 54: 170-178. 10.3171/jns.1981.54.2.0170.
Reithmeier T, Speder B, Pakos P, Brinker G, Lohr M, Klug N, Ernestus RI: Delayed bilateral craniectomy for treatment of traumatic brain swelling in children: case report and review of the literature. Childs Nerv Syst. 2005, 21: 249-253. 10.1007/s00381-004-1044-x. discussion 254
Kawata Y, Niki N, Ohmatsu H, Kusumoto M, Tsuchida T, Eguchi K, Kaneko M, Moriyama N: Quantitative classification based on CT histogram analysis of non-small cell lung cancer: correlation with histopathological characteristics and recurrence-free survival. Med Phys. 2012, 39: 988-1000. 10.1118/1.3679017.
Kalef-Ezra J, Karantanas A, Tsekeris P: CT measurement of lung density. Acta Radiol. 1999, 40: 333-337. 10.3109/02841859909175564.
Biernacki W, Redpath AT, Best JJ, MacNee W: Measurement of CT lung density in patients with chronic asthma. Eur Respir J. 1997, 10: 2455-2459. 10.1183/09031936.97.10112455.
Heremans A, Verschakelen JA, Van Fraeyenhoven L, Demedts M: Measurement of lung density by means of quantitative CT scanning. A study of correlations with pulmonary function tests. Chest. 1992, 102: 805-811. 10.1378/chest.102.3.805.
Yamashiro T, Matsuoka S, Estepar RS, Bartholmai BJ, Diaz A, Ross JC, Murayama S, Silverman EK, Hatabu H, Washko GR: Kurtosis and skewness of density histograms on inspiratory and expiratory CT scans in smokers. COPD. 2011, 8: 13-20. 10.3109/15412555.2010.541537.
Nomori H, Ohtsuka T, Naruke T, Suemasu K: Histogram analysis of computed tomography numbers of clinical T1 N0 M0 lung adenocarcinoma, with special reference to lymph node metastasis and tumor invasiveness. J Thorac Cardiovasc Surg. 2003, 126: 1584-1589. 10.1016/S0022-5223(03)00885-7.
Abe S, Tani S, Nishikawa M, Mizuno N, Yabunaka K, Yamamoto K, Katsuda T, Sanada S: Histogram analysis of ventricular enlargement on CT images of pediatric hydrocephaly patients. 5th European Conference of the International Federation for Medical and Biological Engineering. Edited by: Jobbágy A. 2012, Springer, Berlin Heidelberg, 555-558. 10.1007/978-3-642-23508-5_144. doi: 10.1007/978-3-642-23508-5_144
Kato S, Nakamoto T, Iizuka M: Early diagnosis and estimation of pulmonary congestion and edema in patients with left-sided heart diseases from histogram of pulmonary CT number. Chest. 1996, 109: 1439-1445. 10.1378/chest.109.6.1439.
Gardeur D, Sablayrolles JL, Klausz R, Metzger J: Histographic studies in computed tomography of contrast-enhanced cerebral and orbital tumors. J Comput Assist Tomogr. 1977, 1: 231-240. 10.1097/00004728-197704000-00014.
Kondo T, Ikeda T, Go J: Image processing and quantitative assessment in computer tomography for non-surgical treatment of brain tumours. Acta Neurochir. 1982, 64: 19-37. 10.1007/BF01405616.
Imanaga H, Yamamoto M, Jimbo M, Kobayashi N, Saito Y, Kitamura K: Histographic analysis of computed tomography of cerebral hemispheric tumors. No Shinkei Geka. 1983, 11: 481-488. [in Japanese]
Yuh EL, Cooper SR, Ferguson AR, Manley GT: Quantitative CT improves outcome prediction in acute traumatic brain injury. J Neurotrauma. 2012, 29: 735-746. 10.1089/neu.2011.2008.
Mangel L, Vonoczky K, Hanzely Z, Kiss T, Agoston P, Somogyi A, Nemeth G: CT densitometry of the brain: a novel method for early detection and assessment of irradiation induced brain edema. Neoplasma. 2002, 49: 237-242.
Rozsa L, Grote EH, Egan P: Traumatic brain swelling studied by computerized tomography and densitometry. Neurosurg Rev. 1989, 12: 133-140. 10.1007/BF01741486.
Bendszus M, Urbach H, Meyer B, Schultheiss R, Solymosi L: Improved CT diagnosis of acute middle cerebral artery territory infarcts with density-difference analysis. Neuroradiology. 1997, 39: 127-131. 10.1007/s002340050379.
Barras CD, Tress BM, Christensen S, Collins M, Desmond PM, Skolnick BE, Mayer SA, Davis SM: Quantitative CT densitometry for predicting intracerebral hemorrhage growth. AJNR Am J Neuroradiol. 2013, 34: 1139-1144. 10.3174/ajnr.A3375.
Degos V, Lescot T, Puybasset L: Quantitative CT scan and CT-estimated brain specific gravity in TBI. Yearbook of Intensive Care and Emergency Medicine 2010. Edited by: Vincent JL. 2010, Springer, Berlin Heidelberg, 399-410. 10.1007/978-3-642-10286-8_38. doi:10.1007/978-3-642-10286-8_38
Zach JA, Newell JD, Schroeder J, Murphy JR, Curran-Everett D, Hoffman EA, Westgate PM, Han MK, Silverman EK, Crapo JD, Lynch DA: Quantitative computed tomography of the lungs and airways in healthy nonsmoking adults. Invest Radiol. 2012, 47: 596-602. 10.1097/RLI.0b013e318262292e.
Yuan R, Mayo JR, Hogg JC, Paré PD, McWilliams AM, Lam S, Coxson HO: The effects of radiation dose and CT manufacturer on measurements of lung densitometry. Chest. 2007, 132: 617-623. 10.1378/chest.06-2325.
Nakayama Y, Awai K, Funama Y, Hatemura M, Imuta M, Nakaura T, Ryu D, Morishita S, Sultana S, Sato N: Abdominal CT with low tube voltage: preliminary observations about radiation dose, contrast enhancement, image quality, and noise 1. Radiology. 2005, 237: 945-951. 10.1148/radiol.2373041655.
Kahle KT, Simard JM, Staley KJ, Nahed BV, Jones PS, Sun D: Molecular mechanisms of ischemic cerebral edema: role of electroneutral ion transport. Physiology. 2009, 24: 257-265. 10.1152/physiol.00015.2009.
Liang D, Bhatta S, Gerzanich V, Simard JM: Cytotoxic edema: mechanisms of pathological cell swelling. Neurosurg Focus. 2009, 22: E2.
Nag S, Manias JL, Stewart DJ: Pathology and new players in the pathogenesis of brain edema. Acta Neuropathol. 2009, 118: 197-217. 10.1007/s00401-009-0541-0.
Kochanek PM, Carney N, Adelson PD, Ashwal S, Bell MJ, Bratton S, Carson S, Chesnut RM, Ghajar J, Goldstein B, Grant GA, Kissoon N, Peterson K, Selden NR, Tasker RC, Tong KA, Vavilala MS, Wainwright MS, Warden CR: Guidelines for the acute medical management of severe traumatic brain injury in infants, children, and adolescents-second edition. Pediatr Crit Care Med. 2012, 13: S1-S82. 10.1097/PCC.0b013e31823f435c.
Metter RB, Rittenberger JC, Guyette FX, Callaway CW: Association between a quantitative CT scan measure of brain edema and outcome after cardiac arrest. Resuscitation. 2011, 82: 1180-1185. 10.1016/j.resuscitation.2011.04.001.
Dzialowski I, Weber J, Doerfler A, Forsting M, von Kummer R: Brain tissue water uptake after middle cerebral artery occlusion assessed with CT. J Neuroimaging. 2004, 14: 42-48. 10.1177/1051228403258135.
Dietrich RB, Bradley WG, Zaragoza EJ, Otto RJ, Taira RK, Wilson GH, Kangarloo H: MR evaluation of early myelination patterns in normal and developmentally delayed infants. AJR Am J Roentgenol. 1988, 150: 889-896. 10.2214/ajr.150.4.889.
Staudt M, Schropp C, Staudt F, Obletter N, Bise K, Breit A, Weinmann HM: MRI assessment of myelination: an age standardization. Pediatr Radiol. 1994, 24: 122-127. 10.1007/BF02020169.
Camilla Lindan EG, Jerome B: Pediatric Neuroimaging. Fundamentals of Diagnostic Radiology. Edited by: Brant WE, Clyde A. 2012, Wolters Kluwer/Lippincott Williams & Wilkins, Philadelphia, 1420-4
Murakami JW, Weinberger E, Shaw DW: Normal myelination of the pediatric brain imaged with fluid-attenuated inversion-recovery (FLAIR) MR imaging. AJNR Am J Neuroradiol. 1999, 20: 1406-1411.
Clasen RA, Huckman MS, Von Roenn KA, Pandolfi S, Laing I, Lobick JJ: A correlative study of computed tomography and histology in human and experimental vasogenic cerebral edema. J Comput Assist Tomogr. 1981, 5: 313-327. 10.1097/00004728-198106000-00001.
Donkin JJ, Vink R: Mechanisms of cerebral edema in traumatic brain injury: therapeutic developments. Curr Opin Neurol. 2010, 23: 293-299. 10.1097/WCO.0b013e328337f451.
Thalacker U, Liszka G, Somogyi A, Nemeth G: The observation of edema in the substantia alba during postoperative brain irradiation. The role of computed tomographic studies. Strahlenther Onkol. 1998, 174: 14-18. 10.1007/BF03038222. [in German]
Hiler M, Czosnyka M, Hutchinson P, Balestreri M, Smielewski P, Matta B, Pickard JD: Predictive value of initial computerized tomography scan, intracranial pressure, and state of autoregulation in patients with traumatic brain injury. J Neurosurg. 2006, 104: 731-737. 10.3171/jns.2006.104.5.731.
Moppett IK: Traumatic brain injury: assessment, resuscitation and early management. Br J Anaesth. 2007, 99: 18-31. 10.1093/bja/aem128.
This research was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT & Future Planning (2013R1A1A1004827).
The authors declare that they have no competing interests.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Kim, H., Kim, G., Yoon, B.C. et al. Quantitative analysis of computed tomography images and early detection of cerebral edema for pediatric traumatic brain injury patients: retrospective study. BMC Med 12, 186 (2014). https://doi.org/10.1186/s12916-014-0186-2
- Cerebral edema
- Computed tomography
- Traumatic brain injury