Research Article  Open Access
MengYin Ho, WeiLung Tseng, Furen Xiao, "Estimation of the Craniectomy Surface Area by Using Postoperative Images", International Journal of Biomedical Imaging, vol. 2018, Article ID 5237693, 8 pages, 2018. https://doi.org/10.1155/2018/5237693
Estimation of the Craniectomy Surface Area by Using Postoperative Images
Abstract
Decompressive craniectomy (DC) is a neurosurgical procedure performed to relieve the intracranial pressure engendered by brain swelling. However, no easy and accurate method exists for determining the craniectomy surface area. In this study, we implemented and compared three methods of estimating the craniectomy surface area for evaluating the decompressive effort. We collected 118 sets of preoperative and postoperative brain computed tomography images from patients who underwent craniectomy procedures between April 2009 and April 2011. The surface area associated with each craniectomy was estimated using the marching cube and quasiMonte Carlo methods. The surface area was also estimated using a simple AC method, in which the area is calculated by multiplying the craniectomy length (A) by its height (C). The estimated surface area ranged from 9.46 to 205.32 cm^{2}, with a median of 134.80 cm^{2}. The rootmeansquare deviation (RMSD) between the marching cube and quasiMonte Carlo methods was 7.53 cm^{2}. Furthermore, the RMSD was 14.45 cm^{2} between the marching cube and AC methods and 12.70 cm^{2} between the quasiMonte Carlo and AC methods. Paired tests indicated no statistically significant difference between these methods. The marching cube and quasiMonte Carlo methods yield similar results. The results calculated using the AC method are also clinically acceptable for estimating the DC surface area. Our results can facilitate additional studies on the association of decompressive effort with the effect of craniectomy.
1. Introduction
Decompressive craniectomy (DC) is a common neurosurgical procedure, and it involves removing part of the cranium, thus relieving an edematous brain and reducing intracranial pressure (ICP) by creating extra space. It is typically performed on patients with head injuries or stroke. Although the effectiveness of DC in controlling ICP has been demonstrated [1, 2], whether this control yields more favorable clinical outcomes remains controversial [2, 3].
A larger craniectomy might result in a more favorable ICP control [4], less delayed intracranial hematoma [5, 6], and more desirable clinical outcomes [6–8]; however, certain studies have reported that the craniectomy size is not correlated with complications [9], ICP control [9], or outcomes [10, 11].
The craniectomy size can be expressed in terms of the diameter [5, 8, 12–14], area [7, 9–13, 15, 16], or volume [17] of the skull flap. Although researchers in many studies have measured the craniectomy size according to its area, most of them have used a simplified formula to estimate the area. However, unlike the established ABC/2 formula for estimating the hematoma volume [18–20], the accuracy of such estimation of craniectomy area was never verified, probably because of the technical difficulties associated with computerassisted area analysis.
In accordance with our previous study on the estimation of the skull defect volume [17], we implemented computer algorithms for surface area estimation and used them to validate a simple manual method of estimating the surface area of the skull flap by using postoperative images. To verify the robustness of these methods, data on small craniectomies, including posterior fossa and nondecompressive craniectomies, were also collected.
2. Methods
2.1. Patients
In this retrospective observational study, we collected 118 sets of preoperative and postoperative brain computed tomography (CT) images from patients who underwent craniectomies between April 2009 and April 2011 at the National Taiwan University Hospital. We included image sets of 75 male and 43 female patients. The indications for craniectomy were trauma in 50 patients, spontaneous cerebral hemorrhage in 29, cerebral infarct in 15, infection in 8, dural sinus thrombosis in 1, and brain swelling after various procedures in 15. Posterior fossa craniectomy was performed on 19 patients, and unilateral frontotemporoparietal craniectomy was conducted on the remaining patients.
2.2. Data Acquisition And Preprocessing
The preoperative and postoperative imaging procedures were performed according to a standard protocol of the hospital. All axial images were downloaded from the picture archiving and communication system (PACS) of the hospital to a personal computer for further processing. To obtain the craniectomy size, we first registered the postoperative image with the preoperative one. The missing skull flap was then generated by subtraction between the registered image pair. The processing is illustrated in Figure 1. Image segmentation and registration tasks were developed using C++ and the Insight Toolkit [21]. After the image analysis, the computed missing skull flap in each case was confirmed by boardcertified neurosurgeons on a computer screen, slice by slice.
(a)
(b)
(c)
(d)
2.3. Surface Area Estimation
Three methods were used to estimate the surface area of the skull flap: marching cube, quasiMonte Carlo, and our simple AC method (Figure 2).
(a)
(b)
The marching cube method is a computer graphics algorithm that creates polygonal models of constantdensity surfaces from threedimensional (3D) data [22]. We used the Visualization Toolkit [23] to create a triangulated representation of the missing skull flap and then summed the areas of small triangles. A skull flap mainly has two sides (inner and outer); therefore, we considered only triangles demonstrating a surface normal pointing outside the image center.
The quasiMonte Carlo method is based on the CauchyCrofton formula from integral geometry. Liu et al. published the detailed methodology for applying this method for the surface area estimation of digitized 3D objects [24]. The method essentially entails estimating the surface area of a volumetric object by counting the number of intersection points between the object’s boundary surface and a set of uniformly distributed lines generated using lowdiscrepancy sequences. Liu et al. further found that a desirable surface area estimate can be obtained using few thousand lines if a favorable lowdiscrepancy sequence is used. We followed this method to estimate the surface area of the missing skull flap. After generating a ball with a known surface area () enclosing the skull flap, we used a fourdimensional Sobol sequence [25] implemented in the GNU Scientific Library [26] to generate 20000 lines in the ball. Furthermore, we counted the total intersection points between these lines and the skull flap (). The skull flap surface could then be estimated by , where is the number of intersection points between these lines and the enclosing ball, and its value was in our implementation. The total surface area was then halved as the estimate of the outer surface area.
The proposed manual method can be termed the AC method, according to the ABC convention used to estimate intracerebral hemorrhage volumes [20]. The greatest dimension of the skull defect was identified after reviewing the axial brain CT images, and the linear distance between corners of the outer table of the skull defect was used to determine the craniectomy length (A) (Figure 2(b)). Furthermore, the craniectomy height (C) was determined by adding the interslice distance at which the fullthickness skull defect was visible on the postoperative CT images. The mathematical analysis involved in this method is highly similar to that reported in our previous study [17] and is included in the appendix because the mathematical part is subordinate to the theme of this study. In our study, one neurosurgeon estimated the surface areas of all postoperative images obtained using this method.
2.4. Data Analysis
After surface areas estimated through three methods were recorded, data were analyzed using R [27]. The Pearson correlation coefficient was calculated using linear regression. We used the rootmeansquare deviation (RMSD) to measure the pairwise differences in the surface areas obtained using two methods. BlandAltman analysis [28] was also performed to assess the agreement between these methods. To conduct comparisons with previous estimation studies calculating the hematoma volume [18, 19], we also calculated the mean deviation (in percentage) between the various methods.
3. Results
Table 1 presents a summary of the surface areas estimated using the marching cube (Sm), quasiMonte Carlo (Sq), and AC (Sac) methods. The estimated median surface areas were 134.52 cm^{2} for Sm, 131.73 cm^{2} for Sq, and 133.28 cm^{2} for Sac. Considering their mean as metaestimates, we determined that the craniectomy size ranged from 9.46 to 205.32 cm^{2}.

The data plots obtained using the line of equality and BlandAltman analyses were our main results and are shown in Figure 3. Sm and Sq were similar, and their correlation coefficient was 0.990 and RMSD was 7.53 cm^{2}. The mean deviation was 3.65%. Using BlandAltman analysis, we determined that the 95% confidence interval (CI) for the bias was −13.66 to 15.69 cm^{2}.
The differences between Sm and Sac were significantly higher than those between Sm and Sq. The correlation coefficient and RMSD were 0.968 and 14.45 cm^{2}, respectively. In addition, the mean deviation was 12.12%, and the 95% CI for the bias was −25.57 to 29.20 cm^{2}.
The differences between Sq and Sac were slightly smaller than those between Sm and Sac. The correlation coefficient and RMSD were 0.976 and 12.70 cm^{2}, respectively. Moreover, the mean deviation was 11.41%, and the 95% CI for the bias was −25.18 to 24.79 cm^{2}.
The paired test revealed that Sm and Sq demonstrated no statistically significant difference ( = 0.144). No statistically significant difference was observed between Sm and Sac ( = 0.543) or Sq and Sac (p = 0.865). In brief, the surface area measured using these three methods showed no statistically significant differences.
4. Discussion
DC is a neurosurgical procedure to relieve brain swelling. The craniectomy size might be correlated with ICP control [4], complications [5, 6], and clinical outcomes [6–8, 11]; therefore, the measurement or estimation of the craniectomy size is clinically relevant and should be routinely performed for every patient undergoing DC.
The simplest measurement of the craniectomy size is the skull flap diameter, estimated either during the operation or on postoperative images. Because of bone loss engendered by craniotomy and possible additional bone removal induced by the use of a rongeur, size estimation on postoperative images might be a more accurate approach, considering the advancement of PACSs. However, the simple anteroposterior diameter or largest diameter represents only one dimension of the craniectomy size. Therefore, measuring both the anteroposterior (AP) and superoinferior (SI) diameters is also common in clinical practice, although notations, such as 15 × 12 cm, complicate comparisons and statistics.
The most common measurement of the craniectomy size is its area, either the base or surface area. If the craniectomy base is modeled as an ellipse (Figure 4), then its base area can be estimated by , where and represent the AP and SI diameters, respectively [7]. Although this estimation method is straightforward, validating it is difficult because the craniectomy base cannot be precisely defined to facilitate computerassisted planimetry.
Studies have proposed estimating the craniectomy surface area by considering it as a spherical cap [12, 13] (Figure 4). The surface area of such a cap is , where is the AP diameter and h is the longest distance from to the dural flap. The formula ignores the SI dimension and estimates the protruding brain surface area instead of the craniectomy surface area; therefore, it is prone to change over the entire clinical course.
As per our review of the relevant literature, the present study is the first to validate the estimation formula through computerized methods. Our results revealed that the surface areas estimated using the marching cube and quasiMonte Carlo methods were highly close because these two computerized methods are adequately established in the surface area estimation of digitized 3D objects. However, outliers exist, of which the difference may increase to 40 cm^{2}. These images were typically degraded by metallic artifacts, which produced greater errors during the surface area estimation. We did not discard these data because we aimed to examine the robustness of these methods. Our results confirm the high level of agreement between these digital methods.
Even when the computerized methods were used, the estimated surface area of the digitized object is slightly different from that of the original object [24]. Our implementation of these two digital methods was also not perfect. In the marching cube method, we summed the area of small triangles exhibiting outwardpointing surface normals; this would inevitably include triangles on the edge in addition to the outer surface of craniectomy. Similarly, in the quasiMonte Carlo method, dividing the total surface area by two may overestimate the surface area. This should not be a problem in large craniectomies but would contribute to inaccuracies in small craniectomies.
As expected, the agreement between the results of the AC method and those of the other two methods was not as high. Between Sm and Sac, the Pearson value was 0.968 and mean deviation was 12.12%, and between Sq and Sac, was 0.976 and the mean deviation was 11.41%. For comparison, a study reported that the correlation coefficient () between the ABC/2 method and volumetric analysis was 0.842 for subdural hematoma and 0.929 for intraparenchymal hematoma [18]. The mean deviation between the ABC/2 method and volumetric analysis was reported to be 14.54% [19]. Therefore, the effectiveness of the AC method in estimating the craniectomy surface area is similar to that of the ABC/2 method in estimating the hematoma volume, if not more favorable.
Certain factors render the AC method more useful than the digital methods. First, the AC method requires only the PACS for estimating the craniectomy surface area. Even if the CT image is placed in the traditional view box, we can still obtain an estimate from the film by using only axial slices. More accurate digital methods require extra processes, including segmentation and registration, which are typically not available on the commercial PACS.
The AC method requires only postoperative images to estimate the craniectomy surface area. By contrast, if preoperative images are unavailable for digital methods, we can only provide the mirror image of postoperative CT to estimate the surface area. However, the estimation may be error prone if the patient has an asymmetric skull or if pathology or artifact exists on the other side of the skull.
The AC method is the product of two length measurements; therefore, it should not be affected by most motion and metalrelated artifacts in images. Although digital methods are quite accurate, they might fail if severe image artifacts exist.
Our previous simple ABC method for estimating the skull flap volume is the product of the estimated surface area and bone thickness. However, the bone flap thickness might be erroneously measured because the thickness is usually the smallest dimension and is not even around the skull flap. For instance, we can expect a significant difference in the bone thickness measured between the bone and brain windows. We consider the AC method as a more robust method for estimating the craniectomy surface area; however, we also encourage investigators to record the skull flap volume estimated using the ABC method. When available, digital methods can be used to achieve a more accurate estimation of the craniectomy volume and surface area. If they are not readily available, we consider that the AC and ABC methods should be used to report the craniectomy size for associated studies. If only one measurement of the craniectomy size can be obtained, we recommend using the surface area estimated through the AC method, which is a simple, robust, and widely applicable method.
The association between the decompressive effort and craniectomy effect has been described in numerous previous studies. Studies have considered that DC may reduce the medical refractory ICP [29–32]. However, other studies conclude that the adverse effect of DC results in poor patient outcomes [13, 31]. Little evidence describing the necessity of craniectomy, appropriate timing, and clinical outcomes of different craniectomy sizes is available. Future studies regarding the immediate postcraniectomy effect, such as the change in ICP, reduction of the midline shift, or neurological improvement, can determine the correlation between the craniectomy size and treatment outcomes according to our approach. Studies may also evaluate the complication rate through various craniectomy sizes. They can also evaluate the ideal craniectomy size after further clinical investigation, and individualized preoperative plans can be devised in diverse circumstances.
The adequacy of craniectomy does not solely depend on its size. Other factors, including the extent of brain injury, spaceoccupying lesions, and systemic disorders, may also contribute to the effects of this procedure. However, these factors do not mean that a quantification tool for estimating the craniectomy size is not necessary or is not valuable. The methods we presented for estimating the craniectomy surface area provide tools for further investigation on the effects of size.
Our study had limitations. First, we examined the possible inaccuracy in small craniectomy procedures by using digital methods. Only one neurosurgeon evaluated all the craniectomy surface area by using the AC method; therefore, interrater agreement could not be established. Finally, this study focused on the craniectomy size; therefore, the effect of craniectomy and its association with size were not evaluated.
5. Conclusion
In this study, we compared three methods for estimating the craniectomy surface area. We confirmed that the marching cube and quasiMonte Carlo methods were consistent. The accuracy of the simple AC method was also evaluated. These methods provide a quantitative evaluation for postoperative assessment in patients who have undergone craniectomy.
Appendix
Mathematical Analysis of the AC Method
The mathematical analysis of this method is highly similar to that reported in our previous study [17]. To simplify the derivation, we first model the convexity of the skull as a spherical dome with radius and diameter . We also assume that the craniectomy procedure is performed on a round piece of skull bone with diameter . Therefore, the skull defect can be considered a spherical cap of the skull (Figure 5(a)). The apex angle, denoted as 2θ, is defined as follows:An estimate of the surface area of the skull defect (S) is considered as that of the spherical cap with a base diameter of A [33]:where is the cap height. When (A.1) and (A.2) are used, the estimated surface area of the skull defect can be expressed as follows:The formula is similar to that of the area of a disk (), where in our case.
(a)
(b)
The skull is not a true sphere. If we model it as an ellipsoid, we can also obtain a similar estimate of . However, the exact formula for the surface area of the assumed ellipsoid or spheroid cap is extremely complicated [34, 35] and beyond our need for a simple approximate formula. If we can simply obtain two orthogonal base axes of and of the skull flap, as an extremely intuitive generalization of (A.4), we can estimate the surface area of such a defect by using the following formula:where .
Most craniotomy procedures are performed close to the skull base; therefore, we assume that the skull defect reaches the base of an imaginary dome and derive another axis from the coronal section (Figure 5(b)). The height (C) that we obtained from the axial slices of CT is shorter than the second axis () of the ellipsoid cap that we aimed to investigate (Figure 5(b)). From (A.5), we know thatwhere .
The value of β is dependent on θ and φ. It can also be expressed in terms of A/D and C/R. If both ratios are close to 0.8, as in most large craniotomy procedures, the factor β is close to 1 [17]. Even if the skull defect is extremely small or large, β rarely falls out of the range 0.8–1.2. Hence, we can assume that and estimate the surface area by using .
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This work was supported by Taiwan Ministry of Science and Technology (Grant nos. 1012314B002039 and 1062314B002082).
References
 E. BorSengShu, E. G. Figueiredo, R. L. O. Amorim et al., “Decompressive craniectomy: A metaanalysis of influences on intracranial pressure and cerebral perfusion pressure in the treatment of traumatic brain injury,” Journal of Neurosurgery, vol. 117, no. 3, pp. 589–596, 2012. View at: Publisher Site  Google Scholar
 D. James Cooper, J. V. Rosenfeld, L. Murray et al., “Decompressive craniectomy in diffuse traumatic brain injury,” The New England Journal of Medicine, vol. 364, no. 16, pp. 1493–1502, 2011. View at: Publisher Site  Google Scholar
 J. Sahuquillo, F. MartínezRicarte, and M.A. Poca, “Decompressive craniectomy in traumatic brain injury after the DECRA trial. Where do we stand?” Current Opinion in Critical Care, vol. 19, no. 2, pp. 101–106, 2013. View at: Publisher Site  Google Scholar
 C. P. Gao and B. T. Ang, “Biomechanical modeling of decompressive craniectomy in traumatic brain injury,” Acta Neurochirurgica, Supplementum, vol. 102, pp. 279–282, 2008. View at: Publisher Site  Google Scholar
 S. Wagner, H. Schnippering, A. Aschoff, J. A. Koziol, S. Schwab, and T. Steiner, “Suboptimum hemicraniectomy as a cause of additional cerebral lesions in patients with malignant infarction of the middle cerebral artery,” Journal of Neurosurgery, vol. 94, no. 5, pp. 693–696, 2001. View at: Publisher Site  Google Scholar
 J.Y. Jiang, W. Xu, W.P. Li et al., “Efficacy of standard trauma craniectomy for refractory intracranial hypertension with severe traumatic brain injury: A multicenter, prospective, randomized controlled study,” Journal of Neurotrauma, vol. 22, no. 6, pp. 623–628, 2005. View at: Publisher Site  Google Scholar
 J. Chung, O. Y. Bang, Y. C. Lim, S. K. Park, and Y. S. Shin, “Newly suggested surgical method of decompressive craniectomy for patients with middle cerebral artery infarction,” The Neurologist, vol. 17, no. 1, pp. 11–15, 2011. View at: Publisher Site  Google Scholar
 C. L. Sedney, T. Julien, J. Manon, and A. Wilson, “The effect of craniectomy size on mortality, outcome, and complications after decompressive craniectomy at a rural trauma center,” Journal of Neurosciences in Rural Practice, vol. 5, no. 3, pp. 212–217, 2014. View at: Publisher Site  Google Scholar
 L. Tanrikulu, A. OezTanrikulu, C. Weiss et al., “The bigger, the better? About the size of decompressive hemicraniectomies,” Clinical Neurology and Neurosurgery, vol. 135, article no. 4052, pp. 15–21, 2015. View at: Publisher Site  Google Scholar
 W. T. Curry Jr., M. K. Sethi, C. S. Ogilvy, and B. S. Carter, “Factors associated with outcome after hemicraniectomy for large middle cerebral artery territory infarction,” Neurosurgery, vol. 56, no. 4, pp. 681–692, 2005. View at: Publisher Site  Google Scholar
 C. R. Wirtz, T. Steiner, A. Aschoff et al., “Hemicraniectomy with dural augmentation in medically uncontrollable hemispheric infarction,” Neurosurgical Focus, vol. 2, no. 5, p. E7, 1997. View at: Publisher Site  Google Scholar
 P. De Bonis, A. Pompucci, A. Mangiola, L. Rigante, and C. Anile, “Posttraumatic hydrocephalus after decompressive craniectomy: An underestimated risk factor,” Journal of Neurotrauma, vol. 27, no. 11, pp. 1965–1970, 2010. View at: Publisher Site  Google Scholar
 E. Münch, P. Horn, L. Schürer, A. Piepgras, T. Paul, and P. Schmiedek, “Management of severe traumatic brain injury by decompressive craniectomy,” Neurosurgery, vol. 47, no. 2, pp. 315–323, 2000. View at: Publisher Site  Google Scholar
 T. J. Kenning, R. H. Gandhi, and J. W. German, “A comparison of hinge craniotomy and decompressive craniectomy for the treatment of malignant intracranial hypertension: early clinical and radiographic analysis.,” Neurosurgical Focus, vol. 26, no. 6, p. E6, 2009. View at: Publisher Site  Google Scholar
 M. Olivecrona, M. RodlingWahlstrom, S. Naredi, and L.O. D. Koskinen, “Effective ICP reduction by decompressive craniectomy in patients with severe traumatic brain injury treated by an ICPtargeted therapy,” Journal of Neurotrauma, vol. 24, no. 6, pp. 927–935, 2007. View at: Publisher Site  Google Scholar
 C. M. Schirmer, D. A. Hoit, and A. M. Malek, “Decompressive hemicraniectomy for the treatment of intractable intracranial hypertension after aneurysmal subarachnoid hemorrhage,” Stroke, vol. 38, no. 3, pp. 987–992, 2007. View at: Publisher Site  Google Scholar
 F. Xiao, I.J. Chiang, T. M.H. Hsieh et al., “Estimating postoperative skull defect volume from CT images using the ABC method,” Clinical Neurology and Neurosurgery, vol. 114, no. 3, pp. 205–210, 2012. View at: Publisher Site  Google Scholar
 J. M. Gebel, C. A. Sila, M. A. Sloan et al., “Comparison of the ABC/2 estimation technique to computerassisted volumetric analysis of intraparenchymal and subdural hematomas complicating the GUSTO1 trial,” Stroke, vol. 29, no. 9, pp. 1799–1801, 1998. View at: Publisher Site  Google Scholar
 H. B. Huttner, T. Steiner, M. Hartmann et al., “Comparison of ABC/2 estimation technique to computerassisted planimetric analysis in warfarinrelated intracerebral parenchymal hemorrhage,” Stroke, vol. 37, no. 2, pp. 404–408, 2006. View at: Publisher Site  Google Scholar
 R. U. Kothari, T. Brott, and J. P. Broderick, “The ABCs of measuring intracerebral hemorrhage volumes,” Stroke, vol. 27, no. 8, pp. 13041305, 1996. View at: Publisher Site  Google Scholar
 T. S. Yoo, M. J. Ackerman, W. E. Lorensen et al., “Engineering and algorithm design for an image processing API: a technical report on ITKthe insight toolkit,” Studies in health technology and informatics, pp. 586–592, 2002. View at: Publisher Site  Google Scholar
 W. E. Lorensen and H. E. Cline, “Marching cubes: a high resolution 3D surface construction algorithm,” ACM SIGGRAPH Computer Graphics, vol. 21, no. 4, pp. 163–169, 1987. View at: Publisher Site  Google Scholar
 W. Schroeder, K. Martin, and W. Lorensen, The Visualization Toolkit, Kitware, New York, NY, USA, 2004. View at: Publisher Site
 Y.S. Liu, J. Yi, H. Zhang, G.Q. Zheng, and J.C. Paul, “Surface area estimation of digitized 3D objects using quasiMonte Carlo methods,” Pattern Recognition, vol. 43, no. 11, pp. 3900–3909, 2010. View at: Publisher Site  Google Scholar
 I. A. Antonov and V. M. Saleev, “An economic method of computing LPτsequences,” USSR Computational Mathematics and Mathematical Physics, vol. 19, no. 1, pp. 252–256, 1979. View at: Publisher Site  Google Scholar
 B. Gough, GNU scientific library reference manual: Network Theory Ltd, GNU scientific library reference manual, Network Theory Ltd, 2009.
 Team RC, R: A Language And Environment for Statistical Computing, R Foundation for Statistical Computing, Vienna, Austria, 2012.
 J. Martin Bland and D. Altman, “Statistical methods for assessing agreement between two methods of clinical measurement,” The Lancet, vol. 327, no. 8476, pp. 307–310, 1986. View at: Publisher Site  Google Scholar
 T. S. Skoglund, C. ErikssonRitzén, C. Jensen, and B. Rydenhag, “Aspects on decompressive craniectomy in patients with traumatic head injuries,” Journal of Neurotrauma, vol. 23, no. 10, pp. 1502–1509, 2006. View at: Publisher Site  Google Scholar
 J. S. Huh, H. S. Shin, J. J. Shin, T. H. Kim, Y. S. Hwang, and S. K. Park, “Surgical Management of Massive Cerebral Infarction,” Journal of Korean Neurosurgical Society, vol. 42, no. 4, p. 331, 2007. View at: Publisher Site  Google Scholar
 S.H. Chen, Y. Chen, W.K. Fang, D.W. Huang, K.C. Huang, and S.H. Tseng, “Comparison of craniotomy and decompressive craniectomy in severely headinjured patients with acute subdural hematoma,” Journal of Trauma  Injury Infection and Critical Care, vol. 71, no. 6, pp. 1632–1636, 2011. View at: Publisher Site  Google Scholar
 B. M. Eberle, B. Schnüriger, K. Inaba, J. Peter Gruen, D. Demetriades, and H. Belzberg, “Decompressive craniectomy: Surgical control of traumatic intracranial hypertension may improve outcome,” Injury, vol. 41, no. 9, pp. 894–898, 2010. View at: Publisher Site  Google Scholar
 J. W. Harris and H. Stocker, Handbook of Mathematics and Computational Science, Springer, New York, NY, USA, 1998. View at: Publisher Site  MathSciNet
 H. Y. Erbil and R. A. Meric, “Evaporation of sessile drops on polymer surfaces: Ellipsoidal cap geometry,” The Journal of Physical Chemistry B, vol. 101, no. 35, pp. 6867–6873, 1997. View at: Publisher Site  Google Scholar
 S. T. Wong, A. K. Ball, and J. G. Sivak, “Model of retinal surface area and neuron distribution in the avian eye,” Journal of Neuroscience Methods, vol. 123, no. 1, pp. 1–9, 2003. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 MengYin Ho et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.