Importance of Veins for Neurosurgery as Landmarks Against Brain Shifting Phenomenon: An Anatomical and 3D-MPRAGE MR Reconstruction of Superficial Cortical Veins

Modern neurosurgery uses preoperative imaging daily. Three-dimensional reconstruction of the cortical anatomy and of the superficial veins helps the surgeons plan and perform neurosurgical procedures much more safely. The target is always to give the patient maximum benefit in terms of outcome and minimize intraoperative and postoperative complications. This study aims to develop a method for the combined representation of the cerebral cortex anatomy and the superficial cerebral veins, whose integration is beneficial in daily practice. Only those patients who underwent surgical procedures with craniotomy and a large opening of the dura mater were included in this study, for a total of 23 patients, 13 females (56.5%) and 10 males (43.5%). The average age was 50.1 years. We used a magnetic resonance tomograph Magnetom Vision® 1.5T (Siemens AG). Two sequences were applied: a strongly T1-weighted magnetization-prepared rapid acquisition with gradient echo (MPRAGE) sequence to visualize cerebral anatomical structures, and a FLASH-2D-TOF angiography sequence to visualize the venous vessels on the cortical surface after the administration of a paramagnetic contrast agent. The two data sets were superimposed manually, co-registered in an interactive process, and merged to create a combined data set, segmented and visualized as a three-dimensional reconstruction. Furthermore, we present our method for visualizing superficial veins, which helps manage brain shift (BS). We also performed anatomical observations on the reconstructions. The reconstructions of the cortical and venous anatomy proved to be a valuable tool in surgical planning and positively influenced the surgical procedure. Due to the good correlation with the existing surgical site, this method should be validated on a larger cohort or in a multicentric study.

Introduction

Modern neurosurgery uses preoperative imaging daily. Three-dimensional reconstruction of the cortical anatomy and of the superficial veins helps the surgeons plan and perform neurosurgical procedures and minimize complications. This study aims to develop a method for the combined representation of the cerebral cortex anatomy and the superficial cerebral veins, whose integration is beneficial in daily practice. Two sequences were applied: a strongly T1-weighted magnetization-prepared rapid acquisition with gradient echo (MPRAGE) sequence to visualize cerebral anatomical structures, and a FLASH-2D-TOF angiography sequence to visualize the venous vessels on the cortical surface after the administration of a paramagnetic contrast agent. The two data sets were superimposed manually, co-registered in an interactive process, and merged to create a combined data set, then segmented and visualized as a three-dimensional reconstruction. We provide photo documentation of eachpatient's surgical site. We present our method for visualizing superficial veins, which helps manage brain shift (BS), as the veins on the cortical surface shift together with the underlying brain tissue and therefore do not—or only minimally—change their relative location to the cortical landmarks. Our study includes x-rays, computer tomography, magnetic resonance tomography, and nuclear medicine examinations, as well as ultrasound images and photographs of the surgical site in the case of a previous surgical intervention. We also performed anatomical observations on the reconstructions. The reconstructions of the cortical and venous anatomy proved to be a valuable tool in surgical planning, and it influenced the surgical procedure positively.

Materials and Methods

Patient Population

Our study is an anatomic, observational study, conducted exclusively in the clinic and laboratory of the first (S.O.T.) and senior (P.A.W.) authors, who also performed the anatomical and surgical procedures. The study includes only those patients who underwent surgical procedures with craniotomy and a large opening of the dura mater. Twenty-three patients were enrolled, 13 females (56.5%) and 10 males (43.5%). The average age was 50.1 years. We divided the patient population into two groups (Figure 1A). The first group consisted of 16 patients affected by brain tumors, 11 females and five males. The average age in this first group was 54.5 years, with a minimum age of 30 and a maximum of 66. The second group, composed of three female and four male patients, had pharmacoresistant focal epilepsy. Their age ranged between 9 and 53 years, with a mean of 40.0 years (Figure 1B). Knowledge of the location of important blood vessels—especially veins—is of great significance, both in planning the surgery itself and implanting subdural electrodes, as cortical stimulation entails the risk of venous complications.

FIGURE 1

Figure 1. (A) Distribution of patients for sex and diagnosis. (B) Age distribution in years in both groups (tumor patients and epilepsy patients). (C) Frequency of tumor types in the tumor patients group and localization of epileptogenic focus in the epilepsy patients group.

Brain Shift

We compared the preoperative reconstruction of the superficial cortical vessels with intraoperative photographs of the cortical surface. For the comparative evaluation of the two representations, it is crucial to consider brain shift, which is the deformation of brain tissue during, and due to, neurosurgical intervention (Holmes and Sargent, 1915; Merwarth, 1942; Huber, 1982; Seeger, 1984; Rhoton, 1989, 2002; Koperna et al., 1992; Apuzzo, 1993; Sakaki et al., 1995; Roberts et al., 1998; Choi et al., 2004). Besides obvious causes, such as the resection of tissues or the drainage of CSF from the ventricles, cerebral blood volume (CBV)—depending on carbon dioxide (CO₂) concentration and/or administration of osmotic agents (e.g., Mannitol), anesthetics or steroids—also plays an important role in the brain shift process. In some cases, an initial swelling of the brain tissue occurs shortly after the opening of the dura mater, without signs of increased intracranial pressure.

Moreover, after removing a brain tumor, there is always a distortion of the brain tissue. Therefore, a brain shift of some centimeters can occur. This situation usually shows a discrepancy between the preoperative navigation imagines and the “real-time” situation in the operating room. Multiple intraoperative MRIs can control this dynamic process. Patient population and related pathologies are summarized in Figure 1C.

Magnetic Resonance Imaging

We used a magnetic resonance tomograph Magnetom Vision^® 1.5T (Siemens AG) (Andeweg, 1996; Siemens, 2000). We examined the patients with two different sequences: MPRAGE sequences, which are TurboFLASH sequences extended to 3D and allow high-contrast imaging of the brain, were acquired natively, and FLASH-2D-TOF angiography, to visualize the venous vessels in as high a contrast as possible (Perese, 1958; Stephens and Stilwell, 1969; Strother et al., 1994; Stevenson et al., 1995; Imai et al., 1996; Hunerbein et al., 1997; Liang et al., 2001; Immonen et al., 2004). We oriented the layers to be as perpendicular to the brain surface as possible in the relevant area.

Overlay, Co-registration, and Processing of Records

The MRI data were transferred to a VoxelQ workstation (Picker International, Cleveland, OH, USA) for further processing. The MPRAGE sequence was defined as the primary data set and the FLASH-2D-TOF angiography as the secondary data set (active data set).

Interactive Overlay of the Data Sets

The active data set was color-coded and transparently superimposed onto the primary data set. We assigned a color to each gray value of the active data set. We were interested in larger venous and other non-vascular—structures, especially the SSS, transverse sinus, rectus sinus, larger cortical and bridging veins, the vessels of the interhemispheric gap, and the tentorium cerebelli (Figure 2).

FIGURE 2

Figure 2. Interactive overlay of the data sets. Representation of the primary and secondary data sets in the interactive module of VoxelQ. The active dataset can be moved and rotated in the selected window (yellow border). The views are updated automatically and immediately in all windows (left). The secondary data set has been aligned with the primary data set by displacements in different planes and rotations about different axes (successful co-registration) (right).

Co-registration of the Data Sets and Adjustment of the Threshold Value

In the Combiner module of the VoxelQ, we merged the two data sets into a single, combined data set. Before we created the new data set, however, we increased the contrast of the secondary data set to a maximum so that the pixels above a specific brightness value would be displayed in monochrome blue, while the pixels with a brightness value below this threshold would not be displayed at all.

Segmentation of the Data Set

To obtain an unobstructed view of the cortical surface and the vessels running on it, the overlaying structures—i.e., cranium, and, if possible, dura mater—must be removed from the data volume. In most cases, this was done using the Volume Sculpting module of the VoxelQ workstation. We used a commercial PC with Amira (Indeed—Visual Concepts GmbH, Berlin, Germany). The combined data set is displayed on the screen; the parts to be removed are marked and deleted from the data volume. Since the parts of the 3D data set to be deleted are marked on a two-dimensional projection, the data set must be rotated successively about one or more axes in a process similar to peeling an apple (Figure 3).

FIGURE 3

Figure 3. Procedure for volume sculpting. Marked extracerebral structures (1) are cut away perpendicularly to the viewing plane (2). Rotation of the data set around one or more axes again results in well-exposed extracerebral structures, which are marked (3) and removed. This process is repeated until undisturbed supervision on the cerebrum is possible (4). In most cases, this was done using the Volume Sculpting module of the VoxelQ workstation (Amira, Indeed-Visual Concepts GmbH, Berlin, Germany).

The veins are difficult to separate from the adjacent structures that have been removed, especially in the parasagittal and temporal regions. In particularly difficult cases, we used a different segmentation method, with which we divided the undesired structures from those desired in the 2D layer images before creating the 3D reconstruction (Figure 4). This manual segmentation method—which, however, requires more time—was particularly useful in those cases where veins ran near contrast-absorbing tumors, or where segmentation by volume sculpting did not lead to a satisfactory result. Depending on the goal, segmentation took place either in the primary data set, in the secondary data set, or in the combined data set. In our study, the tumor was often also segmented to hide it later in the combined data set or display it in a color-coded manner.

FIGURE 4

Figure 4. Manual segmentation in the two-dimensional slice images of the data sets. Left: Segmentation of the skull base in the primary data set (MPRAGE sequence) by marking and removing the unwanted region before combining the records. Middle: Segmentation of the secondary data set (venous MR angiography) in a patient with a convexity meningioma. The veins are marked with a stamping tool. Only the marked voxels are later combined with the primary data set. A safe demarcation of the veins from the contrast medium-absorbing tumor is therefore possible. Right: Segmentation of the combined data set in the temporal area, to create an undisturbed view of the cortical surface in the 3D reconstruction of this region.

Determination of the Geometric Error of the Interactive Co-registration

Since there is no binding gold standard, and external reference systems are themselves subject to errors, we examined the reproducibility of co-registration as a measure of the accuracy achieved. To determine the error in the interactive co-registration, we assumed that the secondary data set's position in relation to the primary data set would statistically fluctuate around the correct position after repeated geometric correlation. To obtain a gold standard for the secondary data set's correct position, we registered the morphological and angiographic MR data sets of three patients in 10 consecutive trials. We extracted the parameters of the absolute spatial position and orientation of the active data set from all the individual results. The geometric center of the 10 individual positions was determined, and it served as the gold standard against which we compared the results of the individual experiments.

To calculate the deviation from the individual overlay results' gold standard, we started by limiting the image volume to an average adult's cerebral cortex volume relevant to our study. We calculated this volume by applying a three-dimensional binary mask, which we obtained by increasing the contrast from an FDG-PET examination of a healthy adult's cortex. The image volume thus limited was then determined with the help of a Monte Carlo simulation, using 5,000 random pixels (Duvernoy et al., 1981). Then, we calculated the expected value of the normalized volume integral between the source and the image points of the geometric transformation applied by us, and the average value of the spatial distance between each of the 5,000 corresponding pixels of the gold standard and a single overlay experiment. Here we considered the specific properties of the VoxelQ workstation we used, such as the sequence of the rotations, axis orientation, and direction of rotation. The selection of 5,000 points proved to be sufficiently accurate (variation < 0.01 mm) to make representative statements for the entire image volume. The mean value of the spatial error results, obtained from the 10 registration experiments performed on one patient, indicates the average error—for this patient—of the interactive registration method used by us. To avoid influencing this result by selecting a patient with conditions favorable to superimposition (e.g., good quality MR imaging), we calculated this mean value for three randomly selected patients. The three values were averaged, and the result would correspond to the average spatial deviation of our procedure for one operator (intraobserver variability).

We also investigated the extent to which the registration method's spatial error depends on the user and their previous knowledge. For this purpose, we had five test operators each perform five overlays of the two data sets of a randomly selected patient. Two of the operators were well acquainted with the methodology of interactive overlapping used by us, one had little experience with it, and the other two had no previous training or experience in this field. The test subjects with little or no prior knowledge received a brief introduction to the software and the anatomical landmarks useful for the superposition, but they received no assistance in the actual geometric alignment of the data sets during the co-registration process. With the same method described above, we determined the average spatial error for each subject. The mean value over all five test subjects would correspond to our method's expected geometric error when performed by any one operator, i.e., experts and laypersons alike (interobserver variability).

Time Required for Interactive Co-registration

To allow for a comparison with other data sets co-registration methods, we measured the time required for the method we used. The timing began at the opening of the VoxelQ workstation's interactive module, after both data sets had been loaded into the workspace, and ended when the geometric correlation was completed to our satisfaction.

The average duration of the interactive overlay of new records was 4.7 min, with a minimum of 3.0 min and a maximum of 8.7 min. The mean overlay duration of the records of two randomly selected patients was, for the first patient, 4.2 min, with a minimum of 2.8 min and a maximum of 6.2 min, and 4.6 min for the second patient, with a minimum of 3.3 min and a maximum of 8.7 min. Globally, in multiple overlays of the same records, the overlaying time shortens slightly the more overlay operations are carried out (Figure 5).

FIGURE 5

Figure 5. Image volume to the cerebral cortex volume. We calculated the image volume to the cerebral cortex volume of an average adult relevant to our study by applying a three-dimensional binary mask, which we obtained by increasing the contrast from an FDG-PET examination of a healthy adult patient's cortex.

We timed the first overlay of 12 patients' data sets and used the mean as a benchmark for the first overlay of new data sets. We also timed 10 consecutive overlay attempts on the data sets of two randomly selected patients, then calculated the two mean values and examined whether there was a trend in them.

Validation of the Vein Representation

To measure the reliability of the reconstructions we made, we took intraoperative photographs of each patient's exposed cerebral cortex and of the vessels running on it. Based on these photographs, we then created views of our reconstructions from the same angle we took intraoperatively. The reconstruction and the photo were then aligned by rotating and scaling them in a commercially available image processing program (Adobe Photoshop 7.0) (Figure 6).

FIGURE 6

Figure 6. Validation of vein display. Transparent superimposition (right) of the reconstruction (left) and photograph of the surgical site (middle) for the validation of the vein display.

The brain shift (BS) phenomenon is the basis for a geometric comparison of the reconstruction and the photograph, down to the millimeter, as a mode for correlation. Therefore, we performed a quantitative comparison focusing on neurosurgically relevant veins. For this purpose, relevant veins were marked both in the reconstruction and in the photo. Since it is difficult to determine which vein is the stem and which the confluence vessel in junctions and forks, we decided to mark each subarea individually to facilitate the evaluation.

To mark all the structures we identified as veins in the reconstruction, we desaturated the photo and traced in yellow, on a separate layer, the veins' course. Similarly, the veins found in the photograph of the surgical site were marked in green on a new layer, in cooperation with the neurosurgeons who operated, as it was essential for us to mark all the vessels relevant for the surgical procedure. Fine lateral or end branches without great relevance to perioperative complications were not marked (Figures 7, 8).

FIGURE 7

Figure 7. Marking the veins shown on different drawing levels and evaluation of the levels by superimposition. Left: In the reconstruction, vein sections identified as such are marked in yellow. Middle: In the surgical site photograph, important vein sections for the neurosurgeon are marked in green. Right: The yellow and green markings are displayed one above the other for evaluation.

FIGURE 8

Figure 8. Distribution of the bridging veins and their junctions with the SSS. Anastomosing veins are highlighted in red. This illustration was made by Prof. Peter A. Winkler, based on the results of our anatomical considerations.

The layers with the green and yellow markings were displayed simultaneously to assess the reconstruction's quality, and the results were evaluated using the scheme shown in Table 1. We classified each marked vein section as one of this scheme's five categories and recorded each category's frequency as a percentage. The vein sections marked in the same way in photograph and reconstruction were classified as Category 0; therefore, this category's frequency is a measure of our vein presentation method's reliability.

TABLE 1

Table 1. Overview of veins categories.

Anatomical Examinations

In order to use the reconstructions created within the framework of this manuscript for the performance of anatomical observations, we created several images of each patient, showing the brain and the cortical vessels running on it from different fixed angles. The views we used are shown in Table 2.

TABLE 2

Table 2. Overview of the anatomical analysis and the standardized views of reconstructions.

Determination of the Number of Bridging Veins at the SSS

We examined 23 patients' reconstructions in the views 1–3 presented in Table 1 for the number of bridging veins entering the SSS. This procedure was performed separately for the left and right hemispheres, as well as for the frontal, middle and occipital thirds of the superior sagittal sinus. We summed the mean values for the individual sinus sections to obtain the average number of bridging veins on the corresponding side of the sinus.

Determination of Drainage Types Frequency

We investigated the occurrence of large anastomotic veins in 19 of our patients and classified them according to the types of venous drainage of the lateral cerebral hemispheres presented by Erös et al. (1999). For this purpose, we used the views 4–6 of the reconstructions, presented in Table 1. In four patients, the arrangement of the layer stack of the MR angiography prevented an assessment of the venous anastomoses.

Determination of the Bridging Veins' Confluence Angle at the SSS

In order to determine the bridging veins' confluence angle at the superior sagittal sinus from our reconstructions' views 1–3, it was necessary to think about the underlying type of imaging and the resulting geometric conditions (see Appendix).

We divided the bridging veins at the SSS into six different groups, based on the localization of their confluence: anterior frontal, medial frontal, posterior frontal, anterior parietal, posterior parietal, and occipital. As far as possible, we determined the confluence angles for two veins per group in each patient. We calculated 267 angles: 46 for the anterior frontal group, 45 for the medial frontal, 45 for the posterior frontal group, 44 for the anterior parietal, 46 for the posterior parietal, 41 for the occipital group. For each group, we calculated the mean and mean absolute deviation values. For the distance and angle measurements, we used Adobe Photoshop 7.0. Whenever possible, we compared and verified the anatomical considerations made on the reconstructions with the corresponding intraoperative photos.

Results

Confluence Angle of the Bridging Veins at the SSS

The veins draining the frontal lobes show approximately a right angle at the SSS level (Table 3). The values we measured for the bridging veins' mean confluence angle in the various sections of the SSS, and the mean absolute deviation (MAD) from these mean values as a measure of the values' range, are shown in Table 4.

TABLE 3

Table 3. Average confluence angle of superficial cortical veins at the SSS level.

TABLE 4

Table 4. Values-measurements for the bridging veins' mean confluence angle in the various sections of the SSS, and the mean absolute deviation from these mean values (MAD) as a measure of the range of the values.

Results of the Correlation

In our reconstructions, we identified a total of 158 vein sections. The intraoperative photographs showed 161 relevant veins. The total number of veins categorized by us was 179. In 140 cases, we found a match (Category 0 in Figure 1), which corresponds to 78.2%. Twenty-one patients had veins visible in the photograph that were not shown in our reconstructions (Category 1—11.7%). In eight cases, veins identified in the reconstructions were not visible in the photographs (Category 2—4.5%). In four patients, structures identified as veins in the reconstruction were recognized intraoperatively as arteries (Category 3—2.2%). Finally, in six cases (3.4%), veins identified in the reconstructions could not be assessed in the surgical site (Category 4), as they were covered by cotton wool, already coagulated, or resected together with the surrounding tissue before we took the photograph.

Geometric Error of the Interactive Co-registration

In the three patients, the geometric deviation from the 10 overlay attempts' gold standard was 0.90 mm in the first patient, 1.20 mm in the second patient, and 0.58 mm in the third patient. The intraobserver variability of the interactive co-registration procedure used by us, calculated from these three values, is 0.89 mm (Table 5A).

TABLE 5A

Table 5A. Determination of intra-observer variability.

In our interobserver variability studies, the two experienced users (users 1 and 2) achieved average geometric deviations of 0.71 and 1.32 mm from the gold standard. The less experienced user (user 3) achieved a value of 0.61 mm, while the inexperienced users (users 4 and 5) achieved average deviations of 0.57 and 0.87 mm from the gold standard. The average deviation of each user was 0.82 mm, which is the interobserver variability value of our method (Table 5B).

TABLE 5B

Table 5B. Determination of the inter-observer variability.

Results of the Anatomical Considerations

Number of Bridging Veins

We found an average of 10.9 bridging veins on the left and 10.5 bridging veins on the right of the superior sagittal sinus during the examination of the 23 reconstructions in our study. We observed an average of:

• 4.3 bridging veins on the left and 4.1 on the right of the frontal third of the sinus;

• 4.4 bridging veins on the left and 4.5 on the right of the middle third of the sinus;

• 2.3 bridging veins on the left and 2.0 on the right of the occipital third of the sinus.

Frequency of Drainage Types

We examined 19 reconstructions of the cortical venous anatomy for the drainage type using the classification by Erös et al. We found a Type I configuration in nine cases. We assigned three patients each to Type II, III, and IV, and one patient to Type V cortical venous drainage (Figure 9 and Table 6).

FIGURE 9

Figure 9. Frequency of the different drainage types (according to Erös et al.) in our study's patient population.

TABLE 6

Table 6. Distribution of the found drainage types (according to Erös et al.) in the reconstructions created by us.

Results: Adjustment of the Threshold Value

If the threshold value is set too high, the veins are separated clearly from the brain parenchyma in the reconstruction, but they tend to present themselves discontinuous and broken, and smaller veins and branches do not appear at all, thus increasing the number of false negatives. If, on the other hand, the threshold value is set too low, the main stems are continuous, and the veins are more finely branched, but it is often difficult to distinguish the vessels from the blue “background noise” of the brain parenchyma, and the veins from artifacts. In this case, one would experience an increase in false positives.

Discussion

Preoperative knowledge of the cerebral cortex anatomy, and of the veins running on it, helps the surgeons plan and execute a neurosurgical procedure and avoid complications. Several methods have been developed, in recent years, to visualize the course of the cerebral vessels, often limiting the visualization only to the veins by appropriate methods (conventional angiography, CT angiography, MR angiography) or with imaging reconstruction (Di Chiro, 1962; Kaplan and Browder, 1974; Katada, 1990; Lüders, 1992; Imai et al., 1996; Kikinis et al., 1996; Nakajima et al., 1997; Hawkes, 1998; Tsuchiya et al., 1998, 1999a,b, 2002; Kettenbach et al., 1999; Kubota et al., 2001; Nabavi et al., 2001; Kaminogo et al., 2002; Murphy et al., 2004). To reconstruct the cortical anatomy of the cerebral hemispheres, we used a strongly T1-weighted MPRAGE sequence (Padget, 1956; Sugita et al., 1982; Lanzieri et al., 1987; Liang et al., 2001). These sequences have been used for years (Tsutsumi et al., 1991; Weber and Ivanovic, 1994; Guppy et al., 1997).

Digital Subtraction Angiography (DSA) was the gold standard for intracranial vascular imaging until a few years ago, but it is an invasive procedure with a non-negligible rate of complications (von Lanz and Wachsmuth, 1985; Laub, 1994; West et al., 1997). The method of choice for presenting and evaluating cerebral veins has been described in the literature (Huber, 1982; Lee et al., 1997). The bridging veins and the large cortical anastomotic veins examined in this study can be visualized with CT angiography in 76–97% of cases. Before beginning our study and including patients, we decided to use MR angiography to visualize the cortical veins and examined volunteers with different MRA sequences and contrast agent dosages. We first tried to image the cortical veins using contrast-enhanced MPRAGE sequences, which was successful thanks to the independence of layer orientation and flow velocity (Philips Medical Systems, 1996; Liang et al., 2001). However, it turned out that the smaller superficial veins were only recognizable after the administration of higher doses of paramagnetic contrast medium. The resulting uptake of contrast medium by the meninges posed a problem in the differentiation of the cortical veins from the meninges, significantly limiting the cortical surface's virtual segmentation. Therefore, we turned to flux-sensitive MRA sequences. Slow blood flow in the smaller cortical veins limits their representation in the MRA (Cline et al., 1991; Gandhe et al., 1994; Pant et al., 1997; Mamata et al., 1999; Reichenbach et al., 2000; Wilms et al., 2001). Also, changes in flow direction and velocity, caused by various pathological factors (tumor, edema), can make the representation of venous structures more difficult or impossible (Noachtar et al., 2003). For example, the time-of-flight method might be inadequate to distinguish vessels from stationary tissue at low flow velocities (Bartels and van Overbeeke, 1997). The situation is similar in phase-contrast MRA, where the motion-induced phase shifts may not be detected in slow-flowing blood (Bartels and van Overbeeke, 1997). The loss of signal in vessels that run mostly parallel to the layer orientation can be misread as occlusion or absence of blood vessels (Bartels and van Overbeeke, 1997). For this reason, we oriented the layer stack at an angle of at least 10° to the SSS's course to obtain a good representation of the sinus. Furthermore, we tried to orient the layers in the relevant areas as perpendicular to the cortical surface as possible, to represent the veins better and prevent saturation of the signal of smaller vessels in this region.

MRA may be improved by developing MRIs that work with higher magnetic field strengths–MRAs in the range of 3 Tesla can both significantly improve the spatial resolution and shorten the duration of the examination (Gandhe et al., 1994).

We are aware that there is still potential for optimizing the MR angiography sequence we used. The 2D-TOF sequence does not always detect smaller veins, especially in the temporal region. The maximum of 64 layers is also an obstacle, as the sagittal layers—with a maximum thickness of 2 mm and width of 12.8 cm—do not allow simultaneous imaging of both hemispheres' veins. A more suitable sequence could probably also reduce, if not altogether avoid, the use of paramagnetic contrast agents. A shorter examination time would also be desirable to minimize movement artifacts and stress on the patient.

Registration Method

The registration method we use is interactive and retrospective (Chan and Thompson, 1984; Wenz et al., 1994; Liang et al., 2001). Despite the growing applications of automatic registration, the interactive procedures show a comparable or better accuracy. With interactive image superimposition, the operator can limit the result.

Similarly to a study by Pfluger et al., in our investigations no significant difference was noted between people knowledgeable in both clinical anatomy and three-dimensional image superimposition and laypersons without background knowledge in the evaluation of superimpositions (Wong et al., 1997). The lay group's expected error was also below the value of automatic procedures, such as Wood's algorithm (Woods et al., 1992; Bosmans and Marchal, 1996; Wong et al., 1997). With an average duration of 4.7 min for the primary and secondary data sets' co-registration and a rendering time usually of <30 min, our procedure can be easily integrated into clinical routine Globally, in multiple overlays of the same records, the overlaying time shortens slightly the more overlay operations are carried out.

Validation of the Vein Representation

A reconstruction only has clinical value if it reflects, as reliably as possible, the real anatomical conditions, and intraoperative comparison with the patient's real situs is the best way to compare the two. By marking the veins independently on several layers and subsequently evaluating the marked sections in cooperation with the neurosurgeons performing the procedure, we obtained a concrete quantitative evaluation of the procedure. With a positive correlation of 78.2%, our method for imaging superficial cortical veins and cortical gyral anatomy proved to be a reliable aid in planning and performing neurosurgical interventions. In 11.7% of cases, we could not visualize the relevant veins preoperatively. These were mostly smaller cortical veins, often in the temporal region. This obstacle is an incentive to optimize further the MR angiography sequence we use. The ratio of false-positive veins (4.5%) is associated with the topic of “suitable threshold value.” Identifying cortical arterial vessels as veins, which happened in 2.2% of the marked sections, is unavoidable when using flow-sensitive MRA sequences.

Anatomical Considerations

Every imaging reconstruction provides an incomplete image of reality. Casey et al. reported an average of 6.2 bridging veins to the SSS on each side, while the mean value in our studies was significantly higher, with 10.5 and 10.9 veins per side (Casey et al., 1998). A different sensitivity and specificity of the methods employed, and/or a different interpretation of the reconstructions, could explain this discrepancy. As we were present during the operations and photographed the surgical site, we were able to compare the veins reconstructed with MR imaging to the real ones. We confirmed the results of our anatomical observations in all cases.

The method we have developed for measuring confluence angles has proven to be extremely robust and easy to use in practice. Comparative measurements of the same angles in different images showed only minor deviations in the range of ±5°.

Limits and Strengths of the Study

The technology used in this study could appear dated or old fashioned and—as it requires manual intervention—even difficult to use in the clinical practice, compared with other medical imaging software currently in use in our operating rooms.

However, this limit is not necessarily a weakness. The Brain Shift is particularly important in the navigation systems, which are usually focused on preoperative imaging procedures. Consequently, the intraoperative shift always happens, especially on the brain surface, distorting the superficial venous patterns and causing spatial variations.

In this context, the method for visualizing superficial veins we presented is robust, as the veins on the cortical surface experience the brain shift phenomenon and therefore are not limited to the cortical landmarks or are only minimally displaced.

Furthermore, at present, the study represents the most extensive and scholarly work on the subject, including mathematical and geometrical calculations of the veins flowing into the superior sagittal sinus, which have never been presented before.

Conclusion

Three-dimensional reconstructions of the cortical anatomy and its superficial veins help plan and optimize neurosurgical procedures and minimize complications. Our method for visualizing superficial veins helps manage brain shift (BS), as the veins on the cortical surface shift together with the underlying brain tissue and therefore do not—or only minimally—change their relative location to the cortical landmarks. We developed a method for the combined representation of the cortical anatomy and the superficial cerebral veins that could represent a valuable tool, and a possible adjunct in the neuronavigation armamentarium, in planning and performing neurosurgical interventions. Due to the good correlation with the existing surgical site, it should be validated on a larger cohort or in a multicentric study.

Data Availability Statement

The original contributions presented in the study are included in the article/supplementary materials, further inquiries can be directed to the corresponding author/s.

Ethics Statement

Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements. Written informed consent was not obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author Contributions

All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnana.2020.596167/full#supplementary-material

References

Andeweg, J. (1996). The anatomy of collateral venous flow from the brain and its value in aetiological interpretation of intracranial pathology. Neuroradiology 38, 621–628. doi: 10.1007/s002340050321

PubMed Abstract | CrossRef Full Text | Google Scholar

Apuzzo, M. L. J. (1993). Brain Surgery: Complication Avoidance and Managment. New York, NY: Churchill Livingstone Inc.

Google Scholar

Bartels, R. H., and van Overbeeke, J. J. (1997). Charles Labbe (1851-1889). J. Neurosurg. 87, 477–480. doi: 10.3171/jns.1997.87.3.0477

PubMed Abstract | CrossRef Full Text | Google Scholar

Bosmans, H., and Marchal, G. (1996). Contrast-enhanced MR angiography. Radiologe 36, 115–123. doi: 10.1007/s001170050048

PubMed Abstract | CrossRef Full Text | Google Scholar

Casey, S. O., Rubinstein, D., Lillehei, K. O., Alberico, R. A., Ozsvath, R. R., Cajade-Law, A. G., et al. (1998). Integral and shell-MIP display algorithms in MR and CT three-dimensional models of the brain surface. Am. J. Neuroradiol. 19, 1513–1521.

PubMed Abstract | Google Scholar

Chan, R. C., and Thompson, G. B. (1984). Morbidity, mortality, and quality of life following surgery for intracranial meningiomas. A retrospective study in 257 cases. J. Neurosurg. 60, 52–60. doi: 10.3171/jns.1984.60.1.0052

PubMed Abstract | CrossRef Full Text | Google Scholar

Choi, J. Y., Chang, J. W., Park, Y. G., Kim, T. S., Lee, B. I., and Chung, S. S. (2004). A retrospective study of the clinical outcomes and significant variables in the surgical treatment of temporal lobe tumor associated with intractable seizures. Stereotact. Funct. Neurosurg. 82, 35–42. doi: 10.1159/000076659

PubMed Abstract | CrossRef Full Text | Google Scholar

Cline, H. E., Lorensen, W. E., Souza, S. P., Jolesz, F. A., Kikinis, R., Gerig, G., et al. (1991). 3D surface rendered MR images of the brain and its vasculature. J. Comput. Assist. Tomogr. 15, 344–351. doi: 10.1097/00004728-199103000-00035

PubMed Abstract | CrossRef Full Text | Google Scholar

Di Chiro, G. (1962). Angiographic patterns of cerebral convexity veins and superficial dural sinuses. Am. J. Roentgenol. 87:308–321.

PubMed Abstract | Google Scholar

Duvernoy, H. M., Delon, S., and Vannson, J. L. (1981). Cortical blood vessels of the human brain. Brain Res. Bull. 7, 519–579. doi: 10.1016/0361-9230(81)90007-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Erös, C. A., Herzog, C., and Winkler, P. A. (1999). “Topography of the superficial anastomotic veins of the brain - a microanatomical study,” in Jahrestagung der Deutschen Gesellschaft für Neurochirurgie (München), 71.

Google Scholar

Gandhe, A. J., Hill, D. L., Studholme, C., Hawkes, D. J., Ruff, C. F., Cox, T. C. S., et al. (1994). Combined and three-dimensional rendered multimodal data for planning cranial base surgery: a prospective evaluation. Neurosurgery 35, 463–470; discussion 471. doi: 10.1097/00006123-199409000-00015

CrossRef Full Text | Google Scholar

Guppy, K. H., Origitano, T. C., Reichman, O. H., and Segal, S. (1997). Venous drainage of the inferolateral temporal lobe in relationship to transtemporal/transtentorial approaches to the cranial base. Neurosurgery 41, 615–619; discussion 619–620. doi: 10.1227/00006123-199709000-00021

PubMed Abstract | CrossRef Full Text | Google Scholar

Hawkes, D. J. (1998). Algorithms for radiological image registration and their clinical application. J. Anat. 193(Pt 3), 347–361. doi: 10.1046/j.1469-7580.1998.19330347.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Holmes, G., and Sargent, P. (1915). Injuries of the superior longitudinal sinus. Br. Med. J. 2, 493–498. doi: 10.1136/bmj.2.2857.493

PubMed Abstract | CrossRef Full Text | Google Scholar

Huber, P. (1982). Cerebral Angiography. Stuttgart; New York, NY: Thieme.

Google Scholar

Hunerbein, R., Reuter, P., Meyer, W., and Kuhn, F. P. (1997). [CT angiography of cerebral venous circulation: anatomical visualization and diagnostic pitfalls in interpretation]. Rofo 167, 612–618.

PubMed Abstract | Google Scholar

Imai, F., Ogura, Y., Kiya, N., Zhou, J., Ninomiya, T., Katada, K., et al. (1996). Synthesized surface anatomy scanning (SSAS) for surgical planning of brain metastasis at the sensorimotor region: initial experience with 5 patients. Acta Neurochir. 138, 290–293. doi: 10.1007/BF01411739

PubMed Abstract | CrossRef Full Text | Google Scholar

Immonen, A., Jutila, L., Kalviainen, R., Mervaala, E., Partanen, K., Partanen, J., et al. (2004). Preoperative clinical evaluation, outline of surgical technique and outcome in temporal lobe epilepsy. Adv. Tech. Stand. Neurosurg. 29, 87–132. doi: 10.1007/978-3-7091-0558-0_3

PubMed Abstract | CrossRef Full Text | Google Scholar

Kaminogo, M., Hayashi, H., Ishimaru, H., Morikawa, M., Kitagawa, N., Matsuo, Y., et al. (2002). Depicting cerebral veins by three-dimensional CT angiography before surgical clipping of aneurysms. Am. J. Neuroradiol. 23, 85–91.

PubMed Abstract | Google Scholar

Kaplan, H. A., and Browder, J. (1974). Importance of veins in partial cerebral lobectomy. J. Neurosurg. 41, 360–366. doi: 10.3171/jns.1974.41.3.0360

PubMed Abstract | CrossRef Full Text | Google Scholar

Katada, K. (1990). MR imaging of brain surface structures: surface anatomy scanning (SAS). Neuroradiology 32, 439–448. doi: 10.1007/BF00588477

PubMed Abstract | CrossRef Full Text | Google Scholar

Kettenbach, J., Wong, T., Kacher, D., Hata, N., Schwartz, R. B., McL Black, P., et al. (1999). Computer-based imaging and interventional MRI: applications for neurosurgery. Comput. Med. Imag. Graph 23, 245–258. doi: 10.1016/S0895-6111(99)00022-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Kikinis, R., Gleason, P. L., Moriarty, T. M., Moore, M. R., Alexander, E. III., Stieg, P. E., et al. (1996). Computer-assisted interactive three-dimensional planning for neurosurgical procedures. Neurosurgery 38, 640–649; discussion 649–651. doi: 10.1227/00006123-199604000-00003

PubMed Abstract | CrossRef Full Text | Google Scholar

Koperna, T., Tschabitscher, M., and Knosp, E. (1992). The termination of the vein of “Labbe” and its microsurgical significance. Acta Neurochir. 118, 172–175. doi: 10.1007/BF01401304

PubMed Abstract | CrossRef Full Text | Google Scholar

Kubota, M., Saeki, N., Yamaura, A., Ono, J., and Ozawa, Y. (2001). Influences of venous involvement on postoperative brain damage following the anterior interhemispheric approach. Acta Neurochir. 143, 321–325; discussion 325–326. doi: 10.1007/s007010170085

PubMed Abstract | CrossRef Full Text | Google Scholar

Lanzieri, C. F., Sacher, M., Duchesneau, P. M., Rosenbloom, S. A., and Weinstein, M. A. (1987). The preoperative venogram in planning extended craniectomies. Neuroradiology 29, 360–365. doi: 10.1007/BF00348915

PubMed Abstract | CrossRef Full Text | Google Scholar

Laub, G. (1994). [Basic principles of MR angiography. An introduction]. Radiologe 34, 416–422.

PubMed Abstract | Google Scholar

Lee, B. C., Hatfield, G., Park, T. S., and Kaufman, B. A. (1997). MR imaging surface display of the cerebral cortex in children. Pediatr. Radiol. 27, 199–206. doi: 10.1007/s002470050102

PubMed Abstract | CrossRef Full Text | Google Scholar

Liang, L., Korogi, Y., Sugahara, T., Onomichi, M., Shigematsu, Y., Yanget, D., et al. (2001). Evaluation of the intracranial dural sinuses with a 3D contrast-enhanced MPRAGE sequence: prospective comparison with 2D-TOF MR venography and digital subtraction angiography. Am. J. Neuroradiol. 22, 481–492.

PubMed Abstract | Google Scholar

Lüders, H. O. (1992). Epilepsy Surgery. New York, NY: Raven Press Ltd.

Google Scholar

Mamata, H., Komiya, T., Muro, I., and Matsuyama, S. (1999). Application and validation of three-dimensional data sets from a phase contrast MR angiography for preoperative computer simulation of brain tumors. J. Magn. Reson. Imag. 10, 102–106. doi: 10.1002/(SICI)1522-2586(199907)10:1<102::AID-JMRI15>3.0.CO;2-C

PubMed Abstract | CrossRef Full Text | Google Scholar

Merwarth, H. R. (1942). The syndrome of the rolandic vein (hemiplegia of venous origin). Am. J. Surg. 56, 526–544. doi: 10.1016/S0002-9610(42)90290-3

CrossRef Full Text | Google Scholar

Murphy, M. A., O'Brien, T. J., Morris, K., and Cook, M. J. (2004). Multimodality image-guided surgery for the treatment of medically refractory epilepsy. J. Neurosurg. 100, 452–462. doi: 10.3171/jns.2004.100.3.0452

PubMed Abstract | CrossRef Full Text | Google Scholar

Nabavi, A., Black, P. M., Gering, D. T., Westin, C. F., Mehta, V., Pergolizzi, J., et al. (2001). Serial intraoperative magnetic resonance imaging of brain shift. Neurosurgery 48, 787–797; discussion 797–798. doi: 10.1227/00006123-200104000-00019

PubMed Abstract | CrossRef Full Text | Google Scholar

Nakajima, S., Atsumi, H., Kikinis, R., Moriarty, T. M., Metcalf, D. C., Jolesz, F. A., et al. (1997). Use of cortical surface vessel registration for image-guided neurosurgery. Neurosurgery 40, 1201–1208; discussion 1208–1210. doi: 10.1097/00006123-199706000-00018

PubMed Abstract | CrossRef Full Text | Google Scholar

Noachtar, S., Winkler, P. A., and Lüders, H. O. (2003). “Surgical therapy of epilepsy,” in Neurological Disorders: Course and Treatment, 2nd Edn, eds T. Brandt, C. Caplan, J. Dichgans, J. Diener, C. Kennard (San Diego: Academic Press), 235–44.

Google Scholar

Padget, D. H. (1956). The cranial venous system in man in reference to development, adult configuration, and relation to the arteries. Am. J. Anat. 98, 307–355. doi: 10.1002/aja.1000980302

PubMed Abstract | CrossRef Full Text | Google Scholar

Pant, B., Sumida, M., Kurisu, K., Arita, K., Ikawa, F., Migita, K., et al. (1997). Usefulness of two-dimensional time-of-flight MR angiography combined with surface anatomy scanning for convexity lesions. Neurosurg. Rev. 20, 108–113. doi: 10.1007/BF01138193

PubMed Abstract | CrossRef Full Text | Google Scholar

Perese, D. M. (1958). Superficial veins of the brain from a surgical point of view. J. Neurosurg. 17, 402–412. doi: 10.3171/jns.1960.17.3.0402

PubMed Abstract | CrossRef Full Text | Google Scholar

Philips Medical Systems (1996). Basic Principles of MR and Imaging. Philips Medical Systems.

Google Scholar

Reichenbach, J. R., Barth, M., Haacke, E. M., Klarhofer, M., Kaiser, W. A., and Moser, E. (2000). High-resolution MR venography at 3.0 Tesla. J. Comput. Assist. Tomogr. 24, 949–957. doi: 10.1097/00004728-200011000-00023

PubMed Abstract | CrossRef Full Text | Google Scholar

Rhoton, A. L. Jr. (1989). Comment to the Article: microsurgical anatomy of the venous drainage into the superior sagittal sinus. Neurosurgery 24, 519–20. doi: 10.1227/00006123-198904000-00005

PubMed Abstract | CrossRef Full Text | Google Scholar

Rhoton, A. L. Jr. (2002). The cerebral veins. Neurosurgery 51(4 Suppl.), S159–S205. doi: 10.1097/00006123-200210001-00005

CrossRef Full Text | Google Scholar

Roberts, D. W., Hartov, A., Kennedy, F. E., Miga, M. I., and Paulsen, K. D. (1998). Intraoperative brain shift and deformation: a quantitative analysis of cortical displacement in 28 cases. Neurosurgery 43, 749–758; discussion 758–760. doi: 10.1097/00006123-199810000-00010

PubMed Abstract | CrossRef Full Text | Google Scholar

Sakaki, T., Kakizaki, T., Takeshima, T., Miyamoto, K., and Tsujimoto, S. (1995). Importance of prevention of intravenous thrombosis and preservation of the venous collateral flow in bridging vein injury during surgery: an experimental study. Surg. Neurol. 44, 158–162. doi: 10.1016/0090-3019(95)00160-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Seeger, W. (1984). Microsurgery of Cerebral Veins. Wien; New York, NY: Springer.

Google Scholar

Siemens, A. (2000). Applikationshandbuch zum MAGNETOM Vision/Vision Plus VB33E. Erlangen: Siemens AG, Medizinische Technik.

Google Scholar

Stephens, R. B., and Stilwell, D. L. (1969). Arteries and Veins of the Human Brain. Springfield, IL: Thomas, C.C.

Google Scholar

Stevenson, J., Knopp, E. A., and Litt, A. W. (1995). MPRAGE subtraction venography: a new technique. J. Magn. Reson. Imag. 5, 239–241. doi: 10.1002/jmri.1880050220

PubMed Abstract | CrossRef Full Text | Google Scholar

Strother, S. C., Anderson, J. R., Xu, X. L., Liow, J. S., Bonar, D. C., and Rottenberg, D. A. (1994). Quantitative comparisons of image registration techniques based on high- resolution MRI of the brain. J. Comput. Assist. Tomogr. 18, 954–962. doi: 10.1097/00004728-199411000-00021

PubMed Abstract | CrossRef Full Text | Google Scholar

Sugita, K., Kobayashi, S., and Yokoo, A. (1982). Preservation of large bridging veins during brain retraction. Technical note. J. Neurosurg. 57, 856–858. doi: 10.3171/jns.1982.57.6.0856

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsuchiya, K., Hachiya, J., Hiyama, T., and Maehara, T. (1999a). A new MRI technique for demonstrating the surface of the brain together with the cortical veins. Neuroradiology 41, 425–427. doi: 10.1007/s002340050776

CrossRef Full Text | Google Scholar

Tsuchiya, K., Hachiya, J., Hiyama, T., Maehara, T., and Kassai, Y. (1998). Combination of surface anatomy MRI and MR venography to demonstrate cerebral cortex and cortical veins on one image. J. Comput. Assist. Tomogr. 22, 972–975. doi: 10.1097/00004728-199811000-00023

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsuchiya, K., Katase, S., Hachiya, J., Hiyama, T., and Shiokawa, Y. (1999b). A new technique of surface anatomy MR scanning of the brain: its application to scalp incision planning. Am. J. Neuroradiol. 20, 515–518.

PubMed Abstract | Google Scholar

Tsuchiya, K., Katase, S., Yoshino, A., Hachiya, J., and Shiokawa, Y. (2002). MR-angiogram-added surface anatomy scanning of superficial cerebral arteriovenous malformations. Eur. Radiol. 12, 2330–2334. doi: 10.1007/s00330-001-1299-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsutsumi, K., Shiokawa, Y., Sakai, T., Aoki, N., Kubota, M., and Saito, I. (1991). Venous infarction following the interhemispheric approach in patients with acute subarachnoid hemorrhage. J. Neurosurg. 74, 715–719. doi: 10.3171/jns.1991.74.5.0715

PubMed Abstract | CrossRef Full Text | Google Scholar

von Lanz, T., and Wachsmuth, W. (1985). Praktische Anatomie - Ein Lehr- und Hilfsbuch der Anatomischen Grundlagen Ärztlichen Handelns. Berlin Heidelberg New York, NY: Springer Verlag.

PubMed Abstract | Google Scholar

Weber, D. A., and Ivanovic, M. (1994). Correlative image registration. Semin. Nucl. Med. 24, 311–323. doi: 10.1016/S0001-2998(05)80021-2

CrossRef Full Text | Google Scholar

Wenz, F., Hess, T., Knopp, M. V., Weisser, G., Blüml, S., Lothar, R., et al. (1994). 3D MPRAGE evaluation of lesions in the posterior cranial fossa. Magn. Reson. Imag. 12, 553–558. doi: 10.1016/0730-725X(94)92449-X

PubMed Abstract | CrossRef Full Text | Google Scholar

West, J., Fitzpatrick, J. M., Wang, M. Y., Dawant, B. M., Maurer, C. R. Jr., Kessler, R. M., et al. (1997). Comparison and evaluation of retrospective intermodality brain image registration techniques. J. Comput. Assist. Tomogr. 21, 554–566.

PubMed Abstract | Google Scholar

Wilms, G., Bosmans, H., Demaerel, P., and Marchal, G. (2001). Magnetic resonance angiography of the intracranial vessels. Eur. J. Radiol. 38, 10–18. doi: 10.1016/S0720-048X(01)00285-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Wong, J. C. H., Studholme, C., Hawkes, D. J., and Maisey, M. N. (1997). Evaluation of the limits of visual detection of image misregistration in a brain fluorine-18 fluorodeoxyglucose PETMRI study. Eur. J. Nucl. Med. 24, 642–650. doi: 10.1007/BF00841402

PubMed Abstract | CrossRef Full Text | Google Scholar

Woods, R. P., Cherry, S. R., and Mazziotta, J. C. (1992). Rapid automated algorithm for aligning and reslicing PET images. J. Comput. Assist. Tomogr. 16, 620–633. doi: 10.1097/00004728-199207000-00024

PubMed Abstract | CrossRef Full Text | Google Scholar