Next Article in Journal
Performance Assessment of a New Variable Stiffness Probing System for Micro-CMMs
Next Article in Special Issue
Acute Sleep Deprivation Induces a Local Brain Transfer Information Increase in the Frontal Cortex in a Widespread Decrease Context
Previous Article in Journal
Geiger-Mode Avalanche Photodiode Arrays Integrated to All-Digital CMOS Circuits
Previous Article in Special Issue
RheoStim: Development of an Adaptive Multi-Sensor to Prevent Venous Stasis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Vascular Structure Identification in Intraoperative 3D Contrast-Enhanced Ultrasound Data

by
Elisee Ilunga-Mbuyamba
1,
Juan Gabriel Avina-Cervantes
1,*,
Dirk Lindner
2,
Ivan Cruz-Aceves
3,
Felix Arlt
2 and
Claire Chalopin
4
1
Telematics (CA), Engineering Division (DICIS), University of Guanajuato, Campus Irapuato-Salamanca, Carr. Salamanca-Valle km 3.5 + 1.8, Com. Palo Blanco, Salamanca, Gto. 36885, Mexico
2
Department of Neurosurgery, University Hospital Leipzig, Leipzig 04103, Germany
3
CONACYT Research-Fellow, Center for Research in Mathematics (CIMAT), A.C., Jalisco S/N, Col. Valenciana, Guanajuato, Gto. 36000, Mexico
4
Innovation Center Computer Assisted Surgery (ICCAS), University of Leipzig, Leipzig 04103, Germany
*
Author to whom correspondence should be addressed.
Sensors 2016, 16(4), 497; https://doi.org/10.3390/s16040497
Submission received: 22 February 2016 / Revised: 19 March 2016 / Accepted: 31 March 2016 / Published: 8 April 2016
(This article belongs to the Special Issue Noninvasive Biomedical Sensors)

Abstract

:
In this paper, a method of vascular structure identification in intraoperative 3D Contrast-Enhanced Ultrasound (CEUS) data is presented. Ultrasound imaging is commonly used in brain tumor surgery to investigate in real time the current status of cerebral structures. The use of an ultrasound contrast agent enables to highlight tumor tissue, but also surrounding blood vessels. However, these structures can be used as landmarks to estimate and correct the brain shift. This work proposes an alternative method for extracting small vascular segments close to the tumor as landmark. The patient image dataset involved in brain tumor operations includes preoperative contrast T1MR (cT1MR) data and 3D intraoperative contrast enhanced ultrasound data acquired before (3D-iCEUS s t a r t ) and after (3D-iCEUS e n d ) tumor resection. Based on rigid registration techniques, a preselected vascular segment in cT1MR is searched in 3D-iCEUS s t a r t and 3D-iCEUS e n d data. The method was validated by using three similarity measures (Normalized Gradient Field, Normalized Mutual Information and Normalized Cross Correlation). Tests were performed on data obtained from ten patients overcoming a brain tumor operation and it succeeded in nine cases. Despite the small size of the vascular structures, the artifacts in the ultrasound images and the brain tissue deformations, blood vessels were successfully identified.

Graphical Abstract

1. Introduction

Intraoperative ultrasound imaging is nowadays commonly used in neurosurgery during brain tumor operations [1]. At the beginning of the intervention, the ultrasound images show the surgeon the intraoperative state of the tumor (Figure 1b) [2]. The tumor size or position can be possibly slightly different at the operation time point from the preoperative state depicted in the preoperative MR data (Figure 1a). During the operation ultrasound imaging is a valuable tool to detect the residuals of tumor with the goal to optimize the tumor removal (Figure 1c,d) [3]. However, the interpretation of the ultrasound can be complex [4,5]. The acquisition of the images through the skull opening, called craniotomy, requires skill and experience. Firstly, the sweep of the ultrasound probe is limited by the small opening. The presence of air between the probe and the brain surface blocks the propagation of ultrasound burst. Secondly, the scanned brain area in the images is limited. The orientation of the image and the interpretation of the information it shows can be complex, mostly for images acquired at the end of the operation.
A common method to assist the neurosurgeon in the visual analysis of the intraoperative ultrasound images is to visualize them overlapped to the preoperative MR data in the navigation system because the preoperative image data offers an overview of the entire head [2,6,7]. This technique has limitations for ultrasound images acquired at the end of the operation [5]. The brain structures are so largely deformed, mainly because of tissue resection, so that the information included in the visualization planes does not correspond at all anymore [8].
The evaluation and correction of the brain shift has been extensively studied. Methods are based on similarity measures based on image intensities [9,10,11,12], anatomical landmarks [13,14,15] or biomedical deformation models [16,17]. An interesting previous work proposed to use the vascular structures surrounding the tumor as landmarks in order to estimate the brain deformations [14,18]. In this work, the extraction of the blood vessels was performed in preoperative MR angiographic (MRA) data and in intraoperative Doppler ultrasound images. MRA examination is in routine not acquired in case of brain tumor. Doppler ultrasound can be performed intraoperatively but has limitations to depict the diameters with accuracy and to provide 3D visualization. Moreover, information about brain tissue itself is missing. Several medical studies were conducted to evaluate contrast-enhanced ultrasound imaging for the visualization of brain tumors and the detection of residuals of tumor [19,20,21]. This imaging modality enables to enhance vascularized structures, like lesions and tumors, but also blood vessels [22,23].
In this paper, an alternative method of vascular structures identification in contrast-enhanced ultrasound data is presented with a possible application in the brain shift estimation. Instead of considering the preoperative MRA and intraoperative Doppler ultrasound as previously mentioned, this work introduces the possibility of using cT1MR and 3D-iCEUS data which are routinely involved in brain tumor operations. The method has three main steps i.e., the selection of a vascular segment in cT1MR, the identification of this pattern firstly in 3D-iCEUS s t a r t , then in the 3D-iCEUS e n d . Vesselness-based segmentation methods and rigid registrations techniques were employed, ensuring a computing time compatible with an intraoperative use. Another contribution of this work is the evaluation of the method on CEUS images acquired at the end of the operation, when tissue deformations are large.

2. Materials and Methods

2.1. Patient Image Dataset

Tumor operations were guided by using a neuronavigation system (SonoNavigator, Localite, Sankt Augustin, Germany), including an AplioXG ultrasound device (Toshiba Medical Systems Europe, Zoetermeer, Netherlands) with 2D ultrasound transducers. At the beginning of the intervention, preoperative 3D contrast T1 Magnetic Resonance (cT1MR) data were registered with the patient’s head, based on anatomical landmarks. The procedure was then improved using a surface-based registration technique. During the operation and in addition to the acquisition of intraoperative B-mode ultrasound volumes (noted 3D-iUS), two 3D intraoperative contrast enhanced ultrasound (3D-iCEUS) data were acquired after the injection of 4.8 milliliter of an ultrasound contrast agent (SonoVue, Bracco s.p.a., Milan, Italy) at a rate of 3.0 mL/min using a syringe pump (ACIST VueJect, Bracco s.p.a, Milano, Italy) and using the contrast harmonic imaging (CHI) method. The size and density of the micro-bubbles in the SonoVue contrast agent was 1.9 ± 0.1 micromillimeter and 3.4 ± 0.5 micro-bubbles/mL as recently reported in [24]. The first volume, denoted 3D-iCEUS s t a r t , was obtained transdurally immediately after the craniotomy when the tumor was still entire, and the second one, denoted 3D-iCEUS e n d , was acquired after removing the tumor at the end of the surgery. The surgeon scanned the cerebral region of interest (ROI) with the ultrasound probe (large linear array transducer with a field of view of 38 mm, a range of frequencies of 4.8 MHz to 11.0 MHz and an average frequency of 8 MHz was selected for the acquisitions). The transducer was tracked. The neuronavigation system reconstructs a 3D volume from the 2D-iCEUS slices and overlaid it onto the cT1MR data. The pixel size in the original 2D ultrasound images is 0.422 mm × 0.422 mm and the voxel size of the reconstructed 3D-iCEUS data is 1 × 1 × 1 mm 3 . The data involved in the surgery process are shown in Figure 1.

2.2. Vascular Structure Segmentation

A common solution to distinguish the blood vessels in the 3D-iCEUS data is their extraction by using segmentation methods. However, the segmentation of vascular structures in ultrasound image data is generally a complex task because of the speckle and the blood vessel diameters less than 2 millimeters (Figure 1b,d). Model-based techniques, which include a priori knowledge of the object shape or its data representation, are required to guide the segmentation process and improve the success of algorithms [25,26]. In this work, the extraction of blood vessels is performed in two steps: the computation of a Hessian based vesselness measure for enhancing the vascular structures, followed by the Otsu thresholding algorithm [27]. Moreover, the segmentation is performed in limited target regions of interest in the 3D-iCEUS data, as it will be explained in Section 2.3. The vesselness measure, which is used, was introduced by Sato [28,29] to describe tubular structures within an image using the Hessian matrix. The Hessian matrix H of an image I is defined by:
H = 2 I x 2 2 I x y 2 I x z 2 I x y 2 I y 2 2 I y z 2 I x z 2 I y z 2 I z 2
The second order derivative of the image I is computed from convolution with a Gaussian kernel of standard deviation σ.
I σ = I * G σ
where G σ is a Gaussian function with a standard deviation σ.
G σ = 1 ( 2 π σ 2 ) 3 exp ( - x 2 + y 2 + z 2 2 σ 2 )
Then, the multiscale Hessian matrix will be computed as:
H σ = 2 I σ
Based on experimental analysis of an ideal tube model, Sato et al. [28,30,31] proposed a vessel enhancement filter function given by
v ( σ ) = 0 i f   λ c = 0 exp ( - λ 1 2 2 ( α 1 λ c ) 2 ) . λ c i f   λ 1 0 , λ c 0 exp ( - λ 1 2 2 ( α 2 λ c ) 2 ) . λ c i f   λ 1 > 0 , λ c 0
where α 1 < α 2 and λ c = min ( - λ 2 , - λ 3 ) . The values of α 1 and α 2 are 0.5 and 2, respectively; and λ 1 λ 2 λ 3 (eingenvalues of the Hessian matrix).
Moreover, it has been shown that the variation of sigma values σ enables to describe vascular structures with different radii sizes. Therefore, v ( σ ) is computed at each voxel position in the image for different values of the standard deviations σ m i n σ σ m a x and the maximum response kept v I is described by:
v I = max σ m i n σ σ m a x v ( σ )
Finally the obtained image of vesselness responses is thresholded using the Otsu method [27] in order to extract the vascular structures. Figure 2 presents an illustration of blood vessel segmentation process from cT1MR, 3D-iCEUS s t a r t and 3D-iCEUS e n d . An additional step of 3D representation is added here.

2.3. Vascular Structure Identification

The proposed method consists in identifying the vascular structures in the image data involved in the patient treatment. During the planning stage, some vascular structures close to the tumor are selected by the neurosurgeon. These patterns are automatically recognized in the 3D-iCEUS data by using image registration techniques. The identification method of vascular structures is depicted in Figure 3 and is described as follows:
  • Step 1: Selecting a vascular segment pattern in cT1MR
    Interaction with the application in the operating room has to be limited because of sterilization constraints and restricted time. Due to this, during the operation planning, the user delineates interactively a region of interest including a blood vessel near to the tumor in the cT1MR data. The vascular structure is segmented using the method described in Section 2.2. It performs well because the blood vessels are enhanced in the cT1MR data due to the contrast agent, and any other anatomical structure represented with similar intensities is included in the region of interest. The segmented blood vessel represents the pattern (white frame, Step 1, Figure 3) that is searched for in the 3D-iCEUS s t a r t and 3D-iCEUS e n d data.
  • Step 2: Blood vessel identification in 3D-iCEUS s t a r t
    The pattern is firstly searched in the 3D-iCEUS s t a r t data acquired before resection. In order to reduce the computing time, the search space (large white frame, Step 2, Figure 3) is smaller than the entire 3D-iCEUS s t a r t , but large enough to take the tissue deformations into account. It has twice the volume of the region of interest defined in the cT1MR data, and is centered on the same image position. The enhanced structures in this region are segmented and then a rigid registration method is used to find the sample in the 3D-iCEUS s t a r t data (yellow frame, Step 2, Figure 3) which corresponds best to the pattern. A rigid transformation is sufficient here, since the goal is the identification of the position of the blood vessel in the 3D-iCEUS s t a r t data, which looks like the vascular segment pattern selected in the cT1MR data. The blood vessel detected in the 3D-iCEUS s t a r t then becomes the new pattern, which has to be identified in the 3D-iCEUS e n d data after resection.
  • Step 3: Blood vessel identification in the 3D-iCEUS e n d
    A similar method as described in the step 2 enables to find the position of the target blood vessel (red frame, Step 3, Figure 4) within the search space (yellow frame, Step 3, Figure 3) defined in the 3D-iCEUS e n d data acquired after tumor resection.
Finally, the preselected vessel in the cT1MR is overlapped on the 3D-iCEUS s t a r t and 3D-iCEUS e n d images for visualization purpose. Figure 4 illustrates the workflow of the entire method and illustrations of the blood vessel identification on patient data is shown in Section 3.
It should be noted that the method is so far not able to recognize an incorrect identification of the selected blood vessel. In this case, the user has to reposition interactively the area including the vascular segment corresponding to the pattern. However, the approach is to minimize the interaction between surgeon and computer in operating room because of the sterilization requirement.

2.4. Validation

One key point in image registration is the similarity measure. The Normalized Cross Correlation (NCC) and the Normalized Mutual Information (NMI) are typical similarity measures commonly used [32,33]. The NCC computes pixel-wise normalized cross-correlation between two images and usually aims at registering images of the same modality. The NMI describes how well one image can predict the other one and is, therefore, suitable for multi-modality registration. Robust and easier to interpret, the Normalized Gradient Field (NGF) measures a normalized distance between the gradients of the images to be registered [34,35]. Some works proposed calculating the similarity measures from binary data, arguing that the registration process is, in this way, accelerated [36,37,38,39]. These three similarity measures were used for comparison in the registration algorithm. In order to validate the method, a neurosurgeon visually checked that the vascular segments, automatically identified in the 3D-iCEUS s t a r t and 3D-iCEUS e n d , corresponded to the blood vessel selected in the cT1MR data. Moreover, the results of the algorithm were quantitatively compared with a registration, interactively performed, using the Dice Similarity Index (DSI) and the Hausdorff distance. The DSI represents the percentage of overlapping of two objects. Its value is 1 if they overlap perfectly and 0 if they do not intersect at all. The Hausdorff distance represents the maximal distance between two objects. The computing time of the registration process was also measured for each similarity measure. The NCC, NMI and NCC measures were alternately computed by using the following equations:
N C C = 1 σ x σ y i ( x i - x ¯ ) ( y i - y ¯ )
N M I = H ( X ) + H ( Y ) H ( X , Y )
where x i and y i are respectively the set of points of images X (fixed image) and Y (moving image). In addition, x ¯ and y ¯ are the mean of X and Y; and σ x and σ y are their standard deviation. H ( X ) and H ( Y ) are the entropies of random variables X and Y, and H ( X , Y ) represents their joint entropy.
n ( I , x ) : = I ( x ) | | I ( x ) | | I ( x ) 0 0 o t h e r w i s e
d c ( Y , X ) = | | n ( X , x ) × n ( Y , x ) | | 2
d d ( Y , X ) = < n ( X , x ) , n ( Y , x ) > 2
where n ( I , x ) is a regularized normalized gradient field of a given image I ( x ) . Then, d c ( Y , X ) and d d ( Y , X ) represent the distance between the gradients of images X and Y.
Besides, the Dice Similarity Index (DSI) and the Hausdorff distance are described as follows:
D ( X , Y ) = 2 | X Y | | X | + | Y |
d H ( X , Y ) = max max x X min y Y x - y , max y Y min x X y - x
where X is the fixed image and Y the moving image.

3. Results

The collection of data has been performed at the University Hospital, Department of Neurosurgery, University of Leipzig, Germany, in the context of a previous research project funded by the German Research Society (Deutsche Forschungsgemein-schaft) and accepted by the ethic commission of the University of Leipzig. The implementation was done with an Intel Celeron, 1.5 GHz and 2 GB of memory using MeVisLab tool. The manual registration was validated by the neurosurgeon of this institution. The identification algorithm was tested on ten datasets of patients who overcame a brain tumor operation. A simple vascular segment near to the tumor was selected in the preoperative cT1MR data.
Table 1 includes the DSI, Hausdorff distances and the processing time values computed between the segmented blood vessels in the cT1MR and 3D-iCEUS s t a r t data and between the segmented blood vessels in the 3D-iCEUS s t a r t and 3D-iCEUS e n d data, both normalized with the scores obtained by the expert registration. The values followed by a star indicate that the algorithm failed to find the correct vascular segment. If it occurred before resection, the identification in the 3D-iCEUS e n d data was not performed. It is clearly observed that the proposed approach was able to find correctly the targeted vascular structure with at least one similarity measure used.
Besides, Table 2 includes the mean DSI values and the mean Hausdorff distances, averaged on the successful cases, while Table 3 presents the mean computing time. The outcomes show that the NGF achieves the highest mean rate of registration than the NMI and NCC. Nevertheless, the last ones performs the registration in a shortest time.
In Figure 5, the results obtained step by step for five cases are illustrated. First, the ROI is defined in the cT1MR by the surgeon (column (a)) in the planning stage and it encloses the selected vascular structure so as to reduce the space of segmentation. Second, the vessel is extracted by using the process described in Section 2.2 (column (b)). Third, the vascular structure extracted from the cT1MR data is overlaid on the 3D-iCEUS s t a r t (column (c)). It should be observed that the extracted structure is not spatially aligned with its corresponding in the 3D-iCEUS s t a r t . Fourth, after the segmentation of the vascular structure in 3D-iCEUS s t a r t , column (d) presents the registration results, in which, vascular structures are aligned. At this step, the blood vessel found in 3D-iCEUS s t a r t becomes the new pattern and it is superimposed on the 3D-iCEUS e n d image in column (e). In most of cases, there is no matching between the new pattern with its corresponding vessel in the 3D-iCEUS e n d because of large tissue deformations at this stage of the operation. Finally, the registration of vascular structures from the 3D-iCEUS s t a r to 3D-iCEUS e n d is carried out (column (f)).
An illustration of using vascular structures from cT1MR and 3D-iCEUS s t a r t for correcting the brain shift is depicted in Figure 6. First, the blood vessels are extracted from the both modalities and by overlapping them it is clearly observed that they are not aligned. Second, a registration process is performed to align the segmented structures. By registering the targets, a matrix denoted T which describes the spatial transformation necessary for aligning structures is obtained. Then, this transformation matrix is applied on the cT1MR images so as to achieve its matching with the 3D-iCEUS s t a r t data. Finally, the superimposing of these image modalities is presented for the input before registration and for the output after registration.

4. Discussion

4.1. Visual Validation

The visual checking for identification a selected blood vessel was performed without problem in the 3D-iCEUS s t a r t . The ultrasound images were easy to interpret because the anatomical structures are, in general, similarly represented in the 3D-iCEUS s t a r t and in the preoperative cT1MR data. Brain tissues were already moved at this stage of the operation (after the craniotomy) by comparing the images. The corresponding planes in the 3D-iCEUS s t a r t and preoperative cT1MR data were coarsely only translated and the visual comparison of anatomical structures in these planes were still possible. After tumor resection, the comparison of the 3D-iCEUS s t a r t and the 3D-iCEUS e n d was much more complex for the expert. Firstly, the loss of tumor in the 3D-iCEUS e n d deprives of a valuable reference structure for the comparison of the images. Secondly, the image quality was by trend lower in the 3D-iCEUS e n d . The acquisition was performed with the resection cavity filled with liquid in order to conduct the ultrasound waves. The liquid outflow outside the cavity caused artifacts in the images. Thirdly, brain tissues deformed largely between the beginning of the operation and the stage after tumor resection. Corresponding anatomical structures were elastically deformed. They were represented in different slices and the mental reconstruction of the information was difficult. Therefore, the postoperative MR image were important for assisting the visual evaluation.

4.2. Quantitative Validation

Additionally, the Dice similarity index and the Hausdorff distance were used for quantitatively assess the results obtained from the experiments. It was observed that the NGF similarity provided the highest mean DSI values and the lowest mean Hausdorff distances (Table 2). In contrast, it was found as the slowest (Table 2). Compared to the NMI and NCC, the particularity of the NGF measure is its robustness in cases of curved vascular segments. A T-test was performed to compare the DSI values obtained with the different similarity measures and it showed that the NGF and NCC are statistically different (p-value: 0.0304). As well as, the p-values obtained by comparing the NGF and NMI, then the NMI and NCC were 0.0611 and 0.1747, respectively. In general, the computing time is lower than 15 s on average. The calculation of the NCC and NMI involves the image intensities and was more quickly performed (computing time lower than 10 s on average, Table 3) while the NGF requires the previous computation of the image gradients (computing time lower than 13 s on average, Table 3). However, these mean computing time values are still acceptable for use in the operating theater. The size of the ROI plays an important role on the value of the computing time. It should be noticed that the results from the Table 1 show that the method will always perform a registration even if the correct vascular structure was not found. Despite the tissues deformation which tends to decrease the DSI value, a higher value than 0.60 was found acceptable compared to the neurosurgeon validation. However, the DSI close to 1 indicate the perfect matching obtained.

4.3. Limitations and Future Improvements

4.3.1. Image Quality

The quality of the ultrasound images limits the performance of the registration algorithm. Mainly, it is affected by factors such as: the size of the craniotomy, the location and depth of the tumor in the brain that may complicate the acquisition process. Also, a precise sweeping of the 2D US probe has to be performed within a restricted time window, when the contrast agent enhances optimally the structures. The time window for an optimal acquisition is about 30 s [21]. Nevertheless, a technical solution to improve the image quality is the reduction of the size of voxel in the 3D ultrasound volume, for example 0.5 mm × 0.5 mm × 0.5 mm. The use a 3D ultrasound transducer will enable to perform the acquisition faster to overcome the washed in and out of contrast agent [40].

4.3.2. Misidentification of Vascular Structures

The algorithm succeeded better if a vessel with a specific shape (e.g., curve) is selected. A single straight vascular segment can be in certain cases, confused with other vascular segments or with the resection cavity edges in the 3D-iCEUS e n d data if they are close. Possible approaches to solve this problem would be to segment no-vascular structures in the images in order remove them in the registration process. For example, the tumor can be extracted using a model-based segmentation technique [41]. The resection cavity can be easier identified in the B-mode ultrasound images. Another interesting possibility to test in the future could be the use of the Scale-Invariant Feature Transform (SIFT) proposed by David Lowe [42,43] for registering images as applied in [44,45]. However, in spite of the difficulties mentioned above, the proposed methodology was able to find the correct vascular structures via at least one similarity measure used and, furthermore it was capable to carry out successfully the registration.

4.3.3. Vascular Structures Segmentation

The traditional vesselness based on Hessian matrix was used in this study for vascular structure segmentation. The future intent could be, for instance, the use of alternative methods based on Gabor filter and multiobjective optimization applied in [46] on automatic segmentation of coronary arteries.

5. Conclusions

Despite the small size of the vascular structures surrounding the brain tumors, the low signal to noise ratio in ultrasound image and the brain tissue deformation, it was possible to correctly identify vascular segments in 3D-iCEUS patient’s data. Moreover, we showed that the NGF is more robust, especially in cases of vascular segments with specific shapes. However, the NCC and NMI have a lower computation time than the former. But it is important to notice that the computation time achieved by the algorithm by using these three similarity measures is still compatible with intraoperative use. Besides, an application of this work is the use of the vascular structures as landmarks for the estimation and correction of brain shift. The future works should be the improving of the vascular segmentation in 3D-iCEUS due to the low signal noise ratio in this modality in general. On the other hand, a solution has to be found in cases where the blood vessel is not visible in the intraoperative ultrasound data. For instance, the combination of information included in the B-mode and contrast-enhanced ultrasound data should make the method more robust.

Acknowledgments

This work has been supported by the National Council of Science and Technology of Mexico (CONACYT) under the grant number 493442 and under the project Cátedras-CONACYT No. 3150-3097. Also, the authors would like to thank the POA-DAIP, Engineering Division, University of Guanajuato, Campus Irapuato-Salamanca for providing the necessary funds for covering the costs to publish in open access.

Author Contributions

Claire Chalopin and Dirk Lindner conceived and designed the experiments; Felix Arlt performed the experiments; Elisee Ilunga, Ivan Cruz-Aceves and Juan Gabriel Avina-Cervantes analyzed numerical data and proposed the image vision methodology; Claire Chalopin and Elisee Ilunga wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Unsgaard, G.; Rygh, O.; Selbekk, T.; Muller, T.; Kolstad, F.; Lindseth, F.; Nagelhus Hernes, T. Intra-operative 3D ultrasound in neurosurgery. Acta Neurochir. 2006, 148, 235–253. [Google Scholar] [CrossRef] [PubMed]
  2. Unsgaard, G.; Selbekk, T.; Muller, T.; Ommedal, S.; Torp, S.; Myhr, G.; Bang, J.; Nagelhus Hernes, T. Ability of navigated 3D ultrasound to delineate gliomas and metastases—Comparison of image interpretations with histopathology. Acta Neurochir. 2005, 147, 1259–1269. [Google Scholar] [CrossRef] [PubMed]
  3. Selbekk, T.; Jakola, A.; Solheim, O.E.A. Ultrasound imaging in neurosurgery: Approaches to minimize surgically induced image artefacts for improved resection control. Acta Neurochir. 2013, 155, 973–980. [Google Scholar] [CrossRef] [PubMed]
  4. Solheim, O.; Selbekk, T.; Jakola, A.; Unsgard, G. Ultrasound-guided operations in unselected high-grade gliomas-overall results, impact of image quality and patient selection. Acta Neurochir. 2010, 152, 1873–1886. [Google Scholar] [CrossRef] [PubMed]
  5. Selbekk, T.; Jakola, A.; Solheim, O.; Johansen, T.; Lindseth, F.; Reinertsen, I.; Unsgard, G. Ultrasound imaging in neurosurgery: approaches to minimize surgically induced image artefacts for improved resection control. Acta Neurochir. 2013, 155, 973–980. [Google Scholar] [CrossRef] [PubMed]
  6. Trantakis, C.; Meixensberger, J.; Lindner, D.; Straub, G.; Grunst, G.; Schmidtgen, A.; Arnold, S. Iterative neuronavigation using 3D ultrasound. A feasibilty study. Neurol. Res. 2002, 24, 666–670. [Google Scholar] [CrossRef] [PubMed]
  7. Lindner, D.; Trantakis, C.; Renner, C.; Arnold, S.; Schmitgen, A.; Schneider, J.; Meixensberger, J. Application of Intraoperative 3D Ultrasound During Navigated Tumor Resection. Minim. Invasive Neurosurg. 2006, 49, 197–202. [Google Scholar] [CrossRef] [PubMed]
  8. Maurer, C.R.; Hill, D.L.G.; Maciunas, R.J.; Barwise, J.A.; Fitzpatrick, J.M.; Wang, M.Y. Medical Image Computing and Computer-Assisted Interventation. In Proceedings of the MICCAI’98: First International Conference, Cambridge, MA, USA, 11–13 October 1998; pp. 51–62.
  9. Letteboer, M.; Willems, P.; Viergever, M.; Niessen, W. Brain shift estimation in image-guided neurosurgery using 3-D ultrasound. IEEE Trans. Biomed. Eng. 2005, 52, 268–276. [Google Scholar] [CrossRef] [PubMed]
  10. Ji, S.; Wu, Z.; Hartov, A.; Roberts, D.W.; Paulsen, K.D. Mutual-information-based image to patient re-registration using intraoperative ultrasound in image-guided neurosurgery. Med. Phys. 2008, 35, 4612–4624. [Google Scholar] [CrossRef] [PubMed]
  11. Coupe, P.; Hellier, P.; Morandi, X.; Barillot, C. 3D Rigid Registration of Intraoperative Ultrasound and Preoperative MR Brain Images Based on Hyperechogenic Structures. Int. J. Biomed. Imaging 2012, 2012. [Google Scholar] [CrossRef] [PubMed]
  12. Fuerst, B.; Wein, W.; Muller, M.; Navab, N. Automatic ultrasound-MRI registration for neurosurgery using the 2D and 3D {LC2} Metric. Med. Image Anal. 2014, 18, 1312–1319. [Google Scholar] [CrossRef] [PubMed]
  13. Comeau, R.; Sadikot, A.; Fenster, A.; Peters, T. Intraoperative ultrasound for guidance and tissue shift correction in image-guided neurosurgery. Med. Phys. 2000, 27, 787–800. [Google Scholar] [CrossRef] [PubMed]
  14. Reinertsen, I.; Lindseth, F.; Unsgaard, G.; Collins, D. Clinical validation of vessel-based registration for correction of brain-shift. Med. Image Anal. 2007, 11, 673–684. [Google Scholar] [CrossRef] [PubMed]
  15. Hartov, A.; Roberts, D.; Paulsen, K. A comparative analysis of coregistered ultrasound and magnetic resonance imaging in neurosurgery. Neurosurgery 2008, 62, 99–101. [Google Scholar] [CrossRef] [PubMed]
  16. Ferrant, M.; Nabavi, A.; Macq, B.; Jolesz, F.; Kikinis, R.; Warfield, S. Registration of 3-D intraoperative MR images of the brain using a finite-element biomechanical model. IEEE Trans. Med. Imaging. 2001, 20, 1384–1397. [Google Scholar] [CrossRef] [PubMed]
  17. Hawkes, D.; Barratt, D.; Blackall, J.; Chan, C.; Edwards, P.; Rhode, K.; Penney, G.; McClelland, J.; Hill, D. Tissue deformation and shape models in image-guided interventions: A discussion paper. Med. Image Anal. 2005, 9, 163–175. [Google Scholar] [CrossRef] [PubMed]
  18. Reinertsen, I.; Lindseth, F.; Askeland, C.; Iversen, D.H.; Unsgard, G. Intra-operative correction of brain-shift. Acta Neurochir. 2014, 156, 1301–1310. [Google Scholar] [CrossRef] [PubMed]
  19. Hansen, C.; Wilkening, W.; Ermert, H.; Engelhardt, M.; Schmieder, K.; Krogias, C.; Eyding, J. Intraoperative contrast enhanced perfusion imaging of cerebral tumors. In Proceedings of the 2005 IEEE Ultrasonics Symposium, Rotterdam, The Netherlands, 18–21 September 2005; Volume 2, pp. 743–746.
  20. Kanno, H.; Ozawa, Y.; Sakata, K.; Sato, H.; Tanabe, Y.; Shimizu, N.; Yamamoto, I. Intraoperative power Doppler ultrasonography with a contrast-enhancing agent for intracranial tumors. J. Neurosurg. 2005, 102, 295–301. [Google Scholar] [CrossRef] [PubMed]
  21. Prada, F.; Perin, A.; Martegani, A.; Aiani, L.; Solbiati, L.; Lamperti, M.; Casali, C.; Legnani, F.; Mattei, L.; Saladino, A.; et al. Intraoperative contrast-enhanced ultrasound for brain tumor surgery. Neurosurgery 2014, 74, 542–552. [Google Scholar] [CrossRef] [PubMed]
  22. Holscher, T.; Ozgur, B.; Singel, S.; Wilkening, W.; Mattrey, R.; Sang, H. Intraoperative ultrasound using phase inversion harmonic imaging: first experiences. Neurosurgery 2007, 60, 382–387. [Google Scholar] [CrossRef] [PubMed]
  23. Prada, F.; Bene, M.; Saini, M.; Ferroli, P.; DiMeco, F. Intraoperative cerebral angiosonography with ultrasound contrast agents: How I do it. Acta Neurochir. 2015, 157, 1025–1029. [Google Scholar] [CrossRef] [PubMed]
  24. Hyvelin, J.M.; Greis, C.; Gaud, E.; Costa, M.; Helbert, A.; Bussat, P.; Bettinger, T.; Frinking, P. Characteristics and echogenicity of clinical ultrasound contrast agents: An in vitro and in vivo comparison study. In Proceedings of the 21 European Symposium on Ultrasound Contrast Imaging, An ICUS Conference, Erasmus MC Rotterdam, Rotterdam, The Netherlands, 21–22 January 2016; pp. 5–8.
  25. Gill, J.; Ladak, H.; Steinman, D.; Fenster, A. Accuracy and variability assessment of a semiautomatic technique for segmentation of the carotid arteries from three-dimensional ultrasound images. Med. Phys. 2000, 27, 1333–1342. [Google Scholar] [CrossRef] [PubMed]
  26. Chalopin, C.; Krissian, K.; Meixensberger, J.; Muns, A.; Arlt, F.; Lindner, D. Evaluation of a semi-automatic segmentation algorithm in 3D intraoperative ultrasound brain angiography. Biomed. Tech. 2013, 58, 293–302. [Google Scholar] [CrossRef] [PubMed]
  27. Otsu, N. A Threshold Selection Method from Gray-Level Histograms. IEEE Trans. Syst. Man Cybern. 1979, 9, 62–66. [Google Scholar]
  28. Yoshinobu, S.; Shin, N.; Atsumi, H.; Thomas, K.; Guido, G.; Shigeyuki, Y.; Ron, K. CVRMed-MRCAS’97: First Joint Conference Computer Vision. In Proceedings of the Virtual Reality and Robotics in Medicine and Medical Robotics and Computer-Assisted Surgery, Grenoble, France, 19–22 March 1997; pp. 213–222.
  29. Sato, Y.; Nakajima, S.; Shiraga, N.; Atsumi, H.; Yoshida, S.; Koller, T.; Gerig, G.; Kikinis, R. Three-dimensional multi-scale line filter for segmentation and visualization of curvilinear structures in medical images. Med. Image Anal. 1998, 2, 143–168. [Google Scholar] [CrossRef]
  30. Luu, H.M.; Klink, C.; Moelker, A.; Niessen, W.; van Walsum, T. Quantitative evaluation of noise reduction and vesselness filters for liver vessel segmentation on abdominal CTA images. Phys. Med. Biol. 2015, 60, 3905–3926. [Google Scholar] [CrossRef] [PubMed]
  31. Drechsler, K.; Laura, C.O. Comparison of vesselness functions for multiscale analysis of the liver vasculature. In Proceedings of the 2010 10th IEEE International Conference on Information Technology and Applications in Biomedicine (ITAB), Corfu, Greece, 3–5 November 2010; pp. 1–5.
  32. Andronache, A.; von Siebenthal, M.; Székely, G.; Cattin, P. Non-rigid registration of multi-modal images using both mutual information and cross-correlation. Med. Image Anal. 2008, 12, 3–15. [Google Scholar] [CrossRef] [PubMed]
  33. Pluim, J.; Maintz, J.; Viergever, M. Mutual-information-based registration of medical images: A survey. IEEE Trans. Med. Imaging 2003, 22, 986–1004. [Google Scholar] [CrossRef] [PubMed]
  34. Hodneland, E.; Lundervold, A.; Rorvik, J.; Munthe-Kaas, A.Z. Normalized gradient fields for nonlinear motion correction of DCE-MRI time series. Comput. Med. Imaging Graph. 2014, 38, 202–210. [Google Scholar] [CrossRef] [PubMed]
  35. Haber, E.; Modersitzki, J. Bildverarbeitung für die Medizin 2005: Algorithmen—Systeme—Anwendungen Proceedings des Workshops vom 13.–15. März 2005 in Heidelberg; Springer Berlin Heidelberg: Heidelberg, Germany, 2005; pp. 350–354. [Google Scholar]
  36. Bogush, R.; Maltsev, S.; Ablameyko, S.; Uchida, S.; Kamata, S. An efficient correlation computation method for binary images based on matrix factorisation. In Proceedings of the Sixth International Conference on Document Analysis and Recognition, Seattle, WA, USA, 10–13 September 2001; pp. 312–316.
  37. Chanwimaluang, T.; Fan, G.; Fransen, S. Hybrid retinal image registration. IEEE Trans. Inf. Technol. Biomed. 2006, 10, 129–142. [Google Scholar] [CrossRef] [PubMed]
  38. Sun, M.; Qiao, G.; Zhang, R.; Zong, G. Characteristics of Independence on Image Gray Level in NCCO Applications. In Proceedings of the International Conference on Information Technology and Computer Science (ITCS 2009), Kiev, Ukraine, 25–26 July 2009; Volume 2, pp. 275–278.
  39. Crabb, M.G.; Davidson, J.L.; Little, R.; Wright, P.; Morgan, A.R.; Miller, C.A.; Naish, J.H.; Parker, G.J.M.; Kikinis, R.; McCann, H.; et al. Mutual information as a measure of image quality for 3D dynamic lung imaging with EIT. Physiol. Meas. 2014, 35, 863. [Google Scholar] [CrossRef] [PubMed]
  40. Müns, A.; Meixensberger, J.; Arnold, S.; Schmitgen, A.; Arlt, F.; Chalopin, C.; Lindner, D. Integration of a 3D ultrasound probe into neuronavigation. Acta Neurochir. 2011, 153, 1529–1533. [Google Scholar] [CrossRef] [PubMed]
  41. Chalopin, C.; Lindenberg, R.; Arlt, F.; Muns, A.; Meixensberger, J.; Lindner, D. Brain tumor enhancement revealed by 3D intraoperative ultrasound imaging in a navigation system. Biomed. Eng./Biomed. Tech. 2012, 57, 468–471. [Google Scholar] [CrossRef]
  42. Lowe, D.G. Object recognition from local scale-invariant features. In Proceedings of the Seventh IEEE International Conference on Computer Vision, Kerkyra, Greece, 20–27 September 1999; Volume 2, pp. 1150–1157.
  43. Lowe, D.G. Distinctive Image Features from Scale-Invariant Keypoints. Int. J. Comput. Vis. 2004, 60, 91–110. [Google Scholar] [CrossRef]
  44. Ghassabi, Z.; Shanbehzadeh, J.; Sedaghat, A.; Fatemizadeh, E. An efficient approach for robust multimodal retinal image registration based on UR-SIFT features and PIIFD descriptors. EURASIP J. Image Video Process. 2013, 2013, 1–16. [Google Scholar] [CrossRef]
  45. Chen, J.; Tian, J. Real-time multi-modal rigid registration based on a novel symmetric-SIFT descriptor. Progress Natural Sci. 2009, 19, 643–651. [Google Scholar] [CrossRef]
  46. Cruz-Aceves, I.; Oloumi, F.; Rangayyan, R.M.; Avina-Cervantes, J.G.; Hernandez-Aguirre, A. Automatic segmentation of coronary arteries using Gabor filters and thresholding based on multiobjective optimization. Biomed. Signal Process. Control 2016, 25, 76–85. [Google Scholar] [CrossRef]
Figure 1. Patient image data acquired during tumor resection. The brain tumor (white arrow) is represented in the preoperative cT1MR data (a) and in the 3D-iCEUS s t a r t (b) acquired at the beginning of the operation. After tumor removing, the resection cavity (indicated by two white arrows in (c) and (d) is well visible in the B-mode ultrasound image (c). In the 3D-iCEUS e n d acquired at the end of the operation (d); the borders of the cavity can be easily interpreted as a blood vessel.
Figure 1. Patient image data acquired during tumor resection. The brain tumor (white arrow) is represented in the preoperative cT1MR data (a) and in the 3D-iCEUS s t a r t (b) acquired at the beginning of the operation. After tumor removing, the resection cavity (indicated by two white arrows in (c) and (d) is well visible in the B-mode ultrasound image (c). In the 3D-iCEUS e n d acquired at the end of the operation (d); the borders of the cavity can be easily interpreted as a blood vessel.
Sensors 16 00497 g001
Figure 2. Illustration of vascular structure segmentation in cT1MR, 3D-iCEUS s t a r t and 3D-iCEUS e n d .
Figure 2. Illustration of vascular structure segmentation in cT1MR, 3D-iCEUS s t a r t and 3D-iCEUS e n d .
Sensors 16 00497 g002
Figure 3. Process of the proposed vascular structure identification.
Figure 3. Process of the proposed vascular structure identification.
Sensors 16 00497 g003
Figure 4. Workflow of the proposed method.
Figure 4. Workflow of the proposed method.
Sensors 16 00497 g004
Figure 5. Illustration of vascular structures identification for five cases: (a) preoperative cT1MR; (b) vascular structure selection in cT1MR; (c) overlapping of vascular structure segmented in cT1MR on CEUS s t a r t before registration; (d) overlapping of vascular structure segmented in cT1MR on CEUS s t a r t after registration; (e) overlapping of vascular structure segmented in CEUS s t a r t on CEUS e n d before registration; (f) overlapping of vascular structure segmented in CEUS s t a r t on CEUS e n d after registration.
Figure 5. Illustration of vascular structures identification for five cases: (a) preoperative cT1MR; (b) vascular structure selection in cT1MR; (c) overlapping of vascular structure segmented in cT1MR on CEUS s t a r t before registration; (d) overlapping of vascular structure segmented in cT1MR on CEUS s t a r t after registration; (e) overlapping of vascular structure segmented in CEUS s t a r t on CEUS e n d before registration; (f) overlapping of vascular structure segmented in CEUS s t a r t on CEUS e n d after registration.
Sensors 16 00497 g005
Figure 6. Application of vascular structure identification in brain shift correction between cT1MR and 3D-iCEUS s t a r t .
Figure 6. Application of vascular structure identification in brain shift correction between cT1MR and 3D-iCEUS s t a r t .
Sensors 16 00497 g006
Table 1. Quantitative evaluation of the performance of three different similarity measures in the blood vessel registration algorithm using the DSI and the Hausdorff distance. The values in bold are the best obtained among the three similarity measures. * Indicate that the algorithm failed to find the correct vascular segment.
Table 1. Quantitative evaluation of the performance of three different similarity measures in the blood vessel registration algorithm using the DSI and the Hausdorff distance. The values in bold are the best obtained among the three similarity measures. * Indicate that the algorithm failed to find the correct vascular segment.
PatientSimilarity MeasureROI in cT1MR (voxels)Blood Vessel Identification in the 3D-iCEUS s t a r t Using the Pattern from the cT1MRBlood Vessel Identification in the 3D-iCEUS e n d Using the Pattern from the 3D-iCEUS s t a r t
Processing Time (s)DSIHausdorff Distance (mm)Processing Time (s)DSIHausdorff Distance (mm)
NGF2.80.898.0833.50.82412.207
1NMI50 × 29 × 820.9037.8742.80.93711.180
NCC1.80.9827.4831.70.99410.770
NGF11.10.64521.772140.81311.180
2NMI19 × 37 × 263.50.94617.9172.30.93610.677
NCC5.90.67618.8152.10.0 *14.071 *
NGF3.10.9765.3851.30.98612.042
3NMI35 × 16 × 132.90.9975.3850.30.263 *21.237 *
NCC0.90.9795.3850.50.0 *22.159 *
NGF5.70.76612.3691.50.8299.487
4NMI15 × 17 × 220.70.9378.1240.70.9258.602
NCC0.90.8399.2747.40.8259.487
NGF2.10.74813.403130.76329.172
5NMI61 × 55 × 82.20.422 *28.792 *4.90.3834.015 *
NCC1.20.408 *29.682 *6.70.67432.016
NGF0.40.0 *22.023 *---
6NMI8 × 39 × 261.80.8329.1104.90.9810.863
NCC2.40.8329.11016.20.91310.863
NGF17.60.95328.54756.6117.234
7NMI49 × 56 × 2911.30.82729.65330.20.051 *19.519 *
NCC9.10.98328.5308.80.0 *19.519 *
NGF3.20.9884.6908.30.84910.050
8MI32 × 14 × 192.90.031 *17.720 *---
NCC1.30.097 *15.524 *---
NGF3.10.98814.0004.10.87914.000
9NMI31 × 22 × 222.50.475 *24.083 *---
NCC3.50.458 *25.729 *---
NGF2.70.9522.4491.50.9613.162
10NMI12 × 18 × 120.70.8182.5000.70.0 *8.464 *
NCC2.40.9522.4491.80.0 *8.718 *
Table 2. Mean DSI values and Hausdorff distances computed forthe three similarity measures and by comparison with the expert registration.
Table 2. Mean DSI values and Hausdorff distances computed forthe three similarity measures and by comparison with the expert registration.
Comparative StudiesMean DSI (Algorithm, Expert) Computed after Registration in the:Mean Hausdorff Distances in mm (Algorithm, Expert) Computed after Registration in the:
3D-iCEUS s t a r t 3D-iCEUS e n d 33D-iCEUS s t a r t 3D-iCEUS e n d
NGF vs. Expert0.8700.86813.27213.170
NMI vs. Expert0.8380.83215.11615.570
NCC vs. Expert0.8150.85215.19815.950
Table 3. Mean processing time in s during the registration computed for the successful cases (patient 10 excluded).
Table 3. Mean processing time in s during the registration computed for the successful cases (patient 10 excluded).
Similarity MeasurecT1MR − 3D-iCEUS s t a r t 3D-iCEUS s t a r t − 3D-iCEUS e n d
NGF6.112.8
NMI3.53.1
NCC3.28.0

Share and Cite

MDPI and ACS Style

Ilunga-Mbuyamba, E.; Avina-Cervantes, J.G.; Lindner, D.; Cruz-Aceves, I.; Arlt, F.; Chalopin, C. Vascular Structure Identification in Intraoperative 3D Contrast-Enhanced Ultrasound Data. Sensors 2016, 16, 497. https://doi.org/10.3390/s16040497

AMA Style

Ilunga-Mbuyamba E, Avina-Cervantes JG, Lindner D, Cruz-Aceves I, Arlt F, Chalopin C. Vascular Structure Identification in Intraoperative 3D Contrast-Enhanced Ultrasound Data. Sensors. 2016; 16(4):497. https://doi.org/10.3390/s16040497

Chicago/Turabian Style

Ilunga-Mbuyamba, Elisee, Juan Gabriel Avina-Cervantes, Dirk Lindner, Ivan Cruz-Aceves, Felix Arlt, and Claire Chalopin. 2016. "Vascular Structure Identification in Intraoperative 3D Contrast-Enhanced Ultrasound Data" Sensors 16, no. 4: 497. https://doi.org/10.3390/s16040497

APA Style

Ilunga-Mbuyamba, E., Avina-Cervantes, J. G., Lindner, D., Cruz-Aceves, I., Arlt, F., & Chalopin, C. (2016). Vascular Structure Identification in Intraoperative 3D Contrast-Enhanced Ultrasound Data. Sensors, 16(4), 497. https://doi.org/10.3390/s16040497

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop