2D myocardial deformation imaging based on RF-based non-rigid image registration.

Myocardial deformation imaging is a well-established echocardiographic technique for the assessment of myocardial function. Although some solutions make use of speckle tracking of the reconstructed B-mode images, others apply block matching on the underlying radio-frequency (RF) data in order to increase sensitivity to small inter-frame motion and deformation. However, for both approaches, lateral motion estimation remains a challenge due to the relatively poor lateral resolution of the ultrasound image in combination with the lack of phase information in this direction. Hereto, non-rigid image registration (NRIR) of B-mode images has previously been proposed as an attractive solution. However, hereby, the advantages of RF-based tracking were lost. The aim of this study was therefore to develop an NRIR motion estimator adopted to RF data sets. The accuracy of this estimator was quantified using synthetic data and was contrasted against a state of the art block matching solution. The results show that RF-based NRIR outperforms BM in terms of tracking accuracy particularly, as hypothesized, in the lateral direction. Finally, this RF-based NRIR algorithm was applied clinically, illustrating its ability to estimate both in-plane velocity components in-vivo.

2-D Myocardial Deformation Imaging Based on RF-Based Nonrigid Image Registration.

Myocardial deformation imaging is a well-established echocardiographic technique for the assessment of myocardial function. Although some solutions make use of speckle tracking of the reconstructed B-mode images, others apply block matching (BM) on the underlying radio frequency (RF) data in order to increase sensitivity to small interframe motion and deformation. However, for both approaches, lateral motion estimation remains a challenge due to the relatively poor lateral resolution of the ultrasound image in combination with the lack of phase information in this direction. Hereto, nonrigid image registration (NRIR) of B-mode images has previously been proposed as an attractive solution. However, hereby, the advantages of RF-based tracking were lost. The aim of this paper was, therefore, to develop an NRIR motion estimator adapted to RF data sets. The accuracy of this estimator was quantified using synthetic data and was contrasted against a state-of-the-art BM solution. The results show that RF-based NRIR outperforms BM in terms of tracking accuracy, particularly, as hypothesized, in the lateral direction. Finally, this RF-based NRIR algorithm was applied clinically, illustrating its ability to estimate both in-plane velocity components in vivo.

Realistic Vendor-Specific Synthetic Ultrasound Data for Quality Assurance of 2-D Speckle Tracking Echocardiography: Simulation Pipeline and Open Access Database.

Two-dimensional (2-D) echocardiography is the modality of choice in the clinic for the diagnosis of cardiac disease. Hereto, speckle tracking (ST) packages complement visual assessment by the cardiologist by providing quantitative diagnostic markers of global and regional cardiac function (e.g., displacement, strain, and strain-rate). Yet, the reported high vendor-dependence between the outputs of different ST packages raises clinical concern and hampers the widespread dissemination of the ST technology. In part, this is due to the lack of a solid commonly accepted quality assurance pipeline for ST packages. Recently, we have developed a framework to benchmark ST algorithms for 3-D echocardiography by using realistic simulated volumetric echocardiographic recordings. Yet, 3-D echocardiography remains an emerging technology, whereas the compelling clinical concern is, so far, directed to the standardization of 2-D ST only. Therefore, by building upon our previous work, we present in this paper a pipeline to generate realistic synthetic sequences for 2-D ST algorithms. Hereto, the synthetic cardiac motion is obtained from a complex electromechanical heart model, whereas realistic vendor-specific texture is obtained by sampling a real clinical ultrasound recording. By modifying the parameters in our pipeline, we generated an open-access library of 105 synthetic sequences encompassing: 1) healthy and ischemic motion patterns; 2) the most common apical probe orientations; and 3) vendor-specific image quality from seven different systems. Ground truth deformation is also provided to allow performance analysis. The application of the provided data set is also demonstrated in the benchmarking of a recent academic ST algorithm.

Serial assessment of left ventricular morphology and function in a rodent model of ischemic cardiomyopathy.

Left ventricular remodelling (LVr) occurs post myocardial infarction (MI), predisposing people to heart failure (HF). LV mechanics and morphology are important in this process. We hence sort to characterize LV mechanics and geometry in a post-MI rodent model. Thirty-two male Sprague-Dawley rats (150-200 g) sustained MI (n = 24) or sham (Sham; n = 8) surgery. In another six sham rats invasive blood pressure measurements were performed. Ultrasound imaging was done at baseline, and 1, 3, 7, 14, 30 and 60 days following surgery, and LV mechanics and morphology assessed. LV volumes increased with time (p < 0.01), at a greater rate in the MI group than the Sham group (p < 0.01). Strain was impaired in MI rats at day 1 (13.50 ± 6.64 vs. 25.71 ± 4.94%, p < 0.01) and remained impaired at day 60 (14.07 ± 5.37 vs. 22.98 ± 5.87%, p < 0.01). Strain rate was lower at day 1 (4.11 ± 1.29 vs. 8.10 ± 2.18%/s, p < 0.01), remained lower throughout follow-up (p < 0.01), and decreased at a greater rate in MI rats (p < 0.01). Mean systolic (204 ± 43 vs. 322 ± 75 1/m, p < 0.01) and diastolic (167 ± 21 vs. 192 ± 11 1/m, p < 0.01) curvature was lower in the MI rats at day 1 post surgery and throughout follow-up (p < 0.01). Maximum principal curvature decreased throughout time (p < 0.01), while minimum principal curvature did not (p = 0.86). Wall stress increased significantly after infarction in MI rats (p < 0.01). ST-elevation myocardial infarction (STEMI) changed LV shape and contractile function. The assessment of these indices may prove useful in understanding LVr and the development of HF.

Standardized Delineation of Endocardial Boundaries in Three-Dimensional Left Ventricular Echocardiograms.

BACKGROUND: Three-dimensional (3D) echocardiography is fundamental for left ventricular (LV) assessment. The aim of this study was to determine discrepancies in 3D LV endocardial tracings and suggest tracing guidance. METHODS: Forty-five 3D LV echocardiographic data sets were traced by three experienced operators, from different centers, according to predefined guidelines. The 3D meshes were compared with one another, and the endocardial areas of discrepancies were identified. A discussion and retracing protocol was used to reduce discrepancies. For each data set, an average 3D mesh was produced (reference mesh). Subsequently, four novice operators, divided into two groups, traced 20 of the data sets. Two operators followed the tracing protocol and two did not. RESULTS: The intraclass correlation coefficients among the three experienced operators for end-diastolic volume, end-systolic volume, and ejection fraction were 0.952, 0.955, and 0.932. The absolute distances between tracings were 1.11 ± 0.45 mm. The highest tracing discrepancies were at the apical cap and anterior and anterolateral walls in end-diastole and end-systole and also at the basal anteroseptum in end-systole. Agreement with the reference meshes was better for the novice operators who followed the guidance (10.9 ± 17.3 mL, 10.2 ± 14.7 mL, and -2.2 ± 4.1% for end-diastolic volume, end-systolic volume, and ejection fraction) compared with those who did not (16.3 ± 16.4 mL, 17.0 ± 16.0 mL, and -4.2 ± 4.1%, respectively). CONCLUSIONS: Comparing 3D LV tracings, the endocardial areas that are the most difficult to delineate were identified. The suggested protocol for LV tracing resulted in very good agreement among operators. The reference 3D meshes are available for online testing and ranking of LV tracing algorithms.

Fully automatic left ventricular myocardial strain estimation in 2D short-axis tagged magnetic resonance imaging.

Cardiovascular diseases are among the leading causes of death and frequently result in local myocardial dysfunction. Among the numerous imaging modalities available to detect these dysfunctional regions, cardiac deformation imaging through tagged magnetic resonance imaging (t-MRI) has been an attractive approach. Nevertheless, fully automatic analysis of these data sets is still challenging. In this work, we present a fully automatic framework to estimate left ventricular myocardial deformation from t-MRI. This strategy performs automatic myocardial segmentation based on B-spline explicit active surfaces, which are initialized using an annular model. A non-rigid image-registration technique is then used to assess myocardial deformation. Three experiments were set up to validate the proposed framework using a clinical database of 75 patients. First, automatic segmentation accuracy was evaluated by comparing against manual delineations at one specific cardiac phase. The proposed solution showed an average perpendicular distance error of 2.35 ± 1.21 mm and 2.27 ± 1.02 mm for the endo- and epicardium, respectively. Second, starting from either manual or automatic segmentation, myocardial tracking was performed and the resulting strain curves were compared. It is shown that the automatic segmentation adds negligible differences during the strain-estimation stage, corroborating its accuracy. Finally, segmental strain was compared with scar tissue extent determined by delay-enhanced MRI. The results proved that both strain components were able to distinguish between normal and infarct regions. Overall, the proposed framework was shown to be accurate, robust, and attractive for clinical practice, as it overcomes several limitations of a manual analysis.

heartBEATS: A hybrid energy approach for real-time B-spline explicit active tracking of surfaces.

In this manuscript a novel method is presented for left ventricle (LV) tracking in three-dimensional ultrasound data using a hybrid approach combining segmentation and tracking-based clues. This is accomplished by coupling an affine motion model to an existing LV segmentation framework and introducing an energy term that penalizes the deviation to the affine motion estimated using a global Lucas-Kanade algorithm. The hybrid nature of the proposed solution can be seen as using the estimated affine motion to enhance the temporal coherence of the segmented surfaces, by enforcing the tracking of consistent patterns, while the underlying segmentation algorithm allows to locally refine the estimated global motion. The proposed method was tested on a dataset composed of 24 4D ultrasound sequences from both healthy volunteers and diseased patients. The proposed hybrid tracking platform offers a competitive solution for fast assessment of relevant LV volumetric indices, by combining the robustness of affine motion tracking with the low computational burden of the underlying segmentation algorithm.

Cardiovascular magnetic resonance myocardial feature tracking using a non-rigid, elastic image registration algorithm: assessment of variability in a real-life clinical setting.

BACKGROUND: Cardiovascular magnetic resonance myocardial feature tracking (CMR-FT) is a promising technique for quantification of myocardial strain from steady-state free precession (SSFP) cine images. We sought to determine the variability of CMR-FT using a non-rigid elastic registration algorithm recently available in a commercial software package (Segment, Medviso) in a real-life clinical setting. METHODS: Firstly, we studied the variability in a healthy volunteer who underwent 10 CMR studies over five consecutive days. Secondly, 10 patients were selected from our CMR database yielding normal findings (normal group). Finally, we prospectively studied 10 patients with known or suspected myocardial pathology referred for further investigation to CMR (patient group). In the patient group a second study was performed respecting an interval of 30 min between studies. All studies were manually segmented at the end-diastolic phase by three observers. In all subjects left ventricular (LV) circumferential and radial strain were calculated in the short-axis direction (EccSAX and ErrSAX, respectively) and longitudinal strain in the long-axis direction (EllLAX). The level of CMR experience of the observers was 2 weeks, 6 months and >20 years. RESULTS: Mean contouring time was 7 ± 1 min, mean FT calculation time 13 ± 2 min. Intra- and inter-observer variability was good to excellent with an coefficient of reproducibility (CR) ranging 1.6% to 11.5%, and 1.7% to 16.0%, respectively and an intraclass correlation coefficient (ICC) ranging 0.89 to 1.00 and 0.74 to 0.99, respectively. Variability considerably increased in the test-retest setting with a CR ranging 4.2% to 29.1% and an ICC ranging 0.66 to 0.95 in the patient group. Variability was not influenced by level of expertise of the observers. Neither did the presence of myocardial pathology at CMR negatively impact variability. However, compared to global myocardial strain, segmental myocardial strain variability increased with a factor 2-3, in particular for the basal and apical short-axis slices. CONCLUSIONS: CMR-FT using non-rigid, elastic registration is a reproducible approach for strain analysis in patients routinely scheduled for CMR, and is not influenced by the level of training. However, further improvement is needed to reliably depict small variations in segmental myocardial strain.

Left Ventricular Myocardial Segmentation in 3-D Ultrasound Recordings: Effect of Different Endocardial and Epicardial Coupling Strategies.

Cardiac volume/function assessment remains a critical step in daily cardiology, and 3-D ultrasound plays an increasingly important role. Though development of automatic endocardial segmentation methods has received much attention, the same cannot be said about epicardial segmentation, in spite of the importance of full myocardial segmentation. In this paper, different ways of coupling the endocardial and epicardial segmentations are contrasted and compared with uncoupled segmentation. For this purpose, the B-spline explicit active surfaces framework was used; 27 3-D echocardiographic images were used to validate the different coupling strategies, which were compared with manual contouring of the endocardial and epicardial borders performed by an expert. It is shown that an independent segmentation of the endocardium followed by an epicardial segmentation coupled to the endocardium is the most advantageous. In this way, a framework for fully automatic 3-D myocardial segmentation is proposed using a novel coupling strategy.

Evaluation of the Transverse Oscillation Technique for Cardiac Phased Array Imaging: A Theoretical Study.

The transverse oscillation (TO) technique can improve the estimation of tissue motion perpendicular to the ultrasound beam direction. TOs can be introduced using plane wave (PW) insonification and bilobed Gaussian apodization (BA) on receive (abbreviated as PWTO). Furthermore, the TO frequency of PWTO can be doubled after a heterodyning demodulation process is performed (abbreviated as PWTO*). This paper is concerned with identifying the limitations of the PWTO technique in the specific context of myocardial deformation imaging with phased arrays and investigating the conditions in which it remains advantageous over traditional focused (FOC) beamforming. For this purpose, several tissue phantoms were simulated using Field II, undergoing a wide range of displacement magnitudes and modes (lateral, axial, and rotational motions). The Cramer-Rao lower bound was used to optimize TO beamforming parameters and theoretically predict the fundamental tracking performance limits associated with the FOC, PWTO, and PWTO* beamforming scenarios. This framework was extended to also predict the performance for BA functions that are windowed by the physical aperture of the transducer, leading to higher lateral oscillations. It was found that windowed BA functions resulted in lower jitter errors compared with traditional BA functions. PWTO* outperformed FOC at all investigated signal-to-noise ratio (SNR) levels but only up to a certain displacement, with the advantage rapidly decreasing when the SNR increased. These results suggest that PWTO* improves lateral tracking performance, but only when interframe displacements remain relatively low. This paper concludes by translating these findings into a clinical environment by suggesting optimal scanner settings.

A Comparison of the Performance of Different Multiline Transmit Setups for Fast Volumetric Cardiac Ultrasound.

It was previously demonstrated in 2-D echocardiography that a proper multiline transmit (MLT) implementation can be used to increase frame rate while preserving image quality. Initial findings for extending MLT to 3-D showed that it might address the low spatiotemporal resolution of current volumetric ultrasound systems. However, to date, it remains unclear how much transmit/receive parallelization would be possible using a 3-D MLT system. Therefore, the aim of this paper was to contrast different MLT setups for 3-D imaging by computer simulation in order to determine an optimal tradeoff between the amount of parallelization of an MLT system and the corresponding signal-to-noise ratio of the resulting images. Hereto, the image quality of several MLT setups was estimated by quantifying their crosstalk energy level. The results showed that for the tested setups, 4MLT broad beams and 9MLT narrow beams with Tukey ( α = 0.5 ) apodization in transmit and receive give the highest frame rate gain while maintaining an acceptable interbeam interference level. Moreover, although 16MLT narrow beams with Tukey/Tukey ( α = 0.5 ) apodization did show more pronounced interbeam interference, its gain in frame rate might outweigh its predicted loss in image quality. As such both 9MLT and 16MLT narrow beams were tested experimentally. For both systems, four receive lines were reconstructed from each transmit beam. The contrast-to-noise ratio of these imaging strategies was quantified and compared with the image quality obtained with line-by-line scanning. Despite some expected loss in image quality, the resulting images of the parallelized systems were very competitive to the benchmark, while speeding up the acquisition process by a factor of 36 and 64, respectively.

Detailed Evaluation of Five 3D Speckle Tracking Algorithms Using Synthetic Echocardiographic Recordings.

A plethora of techniques for cardiac deformation imaging with 3D ultrasound, typically referred to as 3D speckle tracking techniques, are available from academia and industry. Although the benefits of single methods over alternative ones have been reported in separate publications, the intrinsic differences in the data and definitions used makes it hard to compare the relative performance of different solutions. To address this issue, we have recently proposed a framework to simulate realistic 3D echocardiographic recordings and used it to generate a common set of ground-truth data for 3D speckle tracking algorithms, which was made available online. The aim of this study was therefore to use the newly developed database to contrast non-commercial speckle tracking solutions from research groups with leading expertise in the field. The five techniques involved cover the most representative families of existing approaches, namely block-matching, radio-frequency tracking, optical flow and elastic image registration. The techniques were contrasted in terms of tracking and strain accuracy. The feasibility of the obtained strain measurements to diagnose pathology was also tested for ischemia and dyssynchrony.

Anatomical Image Registration Using Volume Conservation to Assess Cardiac Deformation From 3D Ultrasound Recordings.

Myocardial deformation imaging can provide valuable insights in myocardial mechanics and help in the diagnosis, prognosis and follow-up of cardiac diseases. However, extracting these indices in 3D is challenging due to the limitations in spatial and temporal resolution of the current volumetric ultrasound systems. For this purpose, we developed an anatomical free-form deformation image registration framework which is locally adapted to the anatomy of the heart. In this work we explored whether incorporating a myocardial volume conservation regularizer would improve strain estimates. We evaluated our technique on in silico echo sequences featuring realistic speckle textures and showed the volume conservation regularizer to be beneficial in reducing strain errors further when used in combination with a smoothness penalty. This combination led to more physiological boundary conditions. It also made distinguishing ischemic from normal segments easier in clinical images.

Standardized Evaluation System for Left Ventricular Segmentation Algorithms in 3D Echocardiography.

Real-time 3D Echocardiography (RT3DE) has been proven to be an accurate tool for left ventricular (LV) volume assessment. However, identification of the LV endocardium remains a challenging task, mainly because of the low tissue/blood contrast of the images combined with typical artifacts. Several semi and fully automatic algorithms have been proposed for segmenting the endocardium in RT3DE data in order to extract relevant clinical indices, but a systematic and fair comparison between such methods has so far been impossible due to the lack of a publicly available common database. Here, we introduce a standardized evaluation framework to reliably evaluate and compare the performance of the algorithms developed to segment the LV border in RT3DE. A database consisting of 45 multivendor cardiac ultrasound recordings acquired at different centers with corresponding reference measurements from three experts are made available. The algorithms from nine research groups were quantitatively evaluated and compared using the proposed online platform. The results showed that the best methods produce promising results with respect to the experts' measurements for the extraction of clinical indices, and that they offer good segmentation precision in terms of mean distance error in the context of the experts' variability range. The platform remains open for new submissions.

Strain assessment in the carotid artery wall using ultrasound speckle tracking: validation in a sheep model.

The aim of this study was to validate carotid artery strain assessment in-vivo using ultrasound speckle tracking. The left carotid artery of five sheep was exposed and sonomicrometry crystals were sutured onto the artery wall to obtain reference strain. Ultrasound imaging was performed at baseline and stress, followed by strain estimation using an in-house speckle tracking algorithm tuned for vascular applications. The correlation between estimated and reference strain was r = 0.95 (p < 0.001) and r = 0.87 (p < 0.01) for longitudinal and circumferential strain, respectively. Moreover, acceptable limits of agreement were found in Bland-Altman analysis (longitudinally: -0.15 to 0.42%, circumferentially: -0.54 to 0.50%), which demonstrates the feasibility of estimating carotid artery strain using ultrasound speckle tracking. However, further studies are needed to test the algorithm on human in-vivo data and to investigate its potential to detect subclinical cardiovascular disease and characterize atherosclerotic plaques.

Automatic mitral annulus tracking in volumetric ultrasound using non-rigid image registration.

Analysis of mitral annular dynamics plays an important role in the diagnosis and selection of optimal valve repair strategies, but remains cumbersome and time-consuming if performed manually. In this paper we propose non-rigid image registration to automatically track the annulus in 3D ultrasound images for both normal and pathological valves, and compare the performance against manual tracing. Relevant clinical properties such as annular area, circumference and excursion could be extracted reliably by the tracking algorithm. The root-mean-square error, calculated as the difference between the manually traced landmarks (18 in total) and the automatic tracking, was 1.96 ± 0.46 mm over 10 valves (5 healthy and 5 diseased) which is within the clinically acceptable error range.

Ultrasound speckle tracking for radial, longitudinal and circumferential strain estimation of the carotid artery--an in vitro validation via sonomicrometry using clinical and high-frequency ultrasound.

Ultrasound speckle tracking for carotid strain assessment has in the past decade gained interest in studies of arterial stiffness and cardiovascular diseases. The aim of this study was to validate and directly contrast carotid strain assessment by speckle tracking applied on clinical and high-frequency ultrasound images in vitro. Four polyvinyl alcohol phantoms mimicking the carotid artery were constructed with different mechanical properties and connected to a pump generating carotid flow profiles. Gray-scale ultrasound long- and short-axis images of the phantoms were obtained using a standard clinical ultrasound system, Vivid 7 (GE Healthcare, Horten, Norway) and a high-frequency ultrasound system, Vevo 2100 (FUJIFILM, VisualSonics, Toronto, Canada) with linear-array transducers (12L/MS250). Radial, longitudinal and circumferential strains were estimated using an in-house speckle tracking algorithm and compared with reference strain acquired by sonomicrometry. Overall, the estimated strain corresponded well with the reference strain. The correlation between estimated peak strain in clinical ultrasound images and reference strain was 0.91 (p<0.001) for radial strain, 0.73 (p<0.001) for longitudinal strain and 0.90 (p<0.001) for circumferential strain and for high-frequency ultrasound images 0.95 (p<0.001) for radial strain, 0.93 (p<0.001) for longitudinal strain and 0.90 (p<0.001) for circumferential strain. A significant larger bias and root mean square error was found for circumferential strain estimation on clinical ultrasound images compared to high frequency ultrasound images, but no significant difference in bias and root mean square error was found for radial and longitudinal strain when comparing estimation on clinical and high-frequency ultrasound images. The agreement between sonomicrometry and speckle tracking demonstrates that carotid strain assessment by ultrasound speckle tracking is feasible.

Ultrasound speckle tracking strain estimation of in vivo carotid artery plaque with in vitro sonomicrometry validation.

Our objective was to validate a previously developed speckle tracking (ST) algorithm to assess strain in common carotid artery plaques. Radial and longitudinal strain was measured in common carotid artery gel phantoms with a plaque-mimicking inclusion using an in-house ST algorithm and sonomicrometry. Moreover, plaque strain by ST for seven patients (77 ± 6 y) with carotid atherosclerosis was compared with a quantitative visual assessment by two experienced physicians. In vitro, good correlation existed between ST and sonomicrometry peak strains, both radially (r = 0.96, p < 0.001) and longitudinally (r = 0.75, p < 0.01). In vivo, greater pulse pressure-adjusted radial and longitudinal strains were found in echolucent plaques than in echogenic plaques. This illustrates the feasibility of ultrasound ST strain estimation in plaques and the possibility of characterizing plaques using ST strain in vivo.

Fast automatic myocardial segmentation in 4D cine CMR datasets.

A novel automatic 3D+time left ventricle (LV) segmentation framework is proposed for cardiac magnetic resonance (CMR) datasets. The proposed framework consists of three conceptual blocks to delineate both endo and epicardial contours throughout the cardiac cycle: (1) an automatic 2D mid-ventricular initialization and segmentation; (2) an automatic stack initialization followed by a 3D segmentation at the end-diastolic phase; and (3) a tracking procedure. Hereto, we propose to adapt the recent B-spline Explicit Active Surfaces (BEAS) framework to the properties of CMR images by integrating dedicated energy terms. Moreover, we extend the coupled BEAS formalism towards its application in 3D MR data by adapting it to a cylindrical space suited to deal with the topology of the image data. Furthermore, a fast stack initialization method is presented for efficient initialization and to enforce consistent cylindrical topology. Finally, we make use of an anatomically constrained optical flow method for temporal tracking of the LV surface. The proposed framework has been validated on 45 CMR datasets taken from the 2009 MICCAI LV segmentation challenge. Results show the robustness, efficiency and competitiveness of the proposed method both in terms of accuracy and computational load.

Real-time 3D interactive segmentation of echocardiographic data through user-based deformation of B-spline explicit active surfaces.

Image segmentation is an ubiquitous task in medical image analysis, which is required to estimate morphological or functional properties of given anatomical targets. While automatic processing is highly desirable, image segmentation remains to date a supervised process in daily clinical practice. Indeed, challenging data often requires user interaction to capture the required level of anatomical detail. To optimize the analysis of 3D images, the user should be able to efficiently interact with the result of any segmentation algorithm to correct any possible disagreement. Building on a previously developed real-time 3D segmentation algorithm, we propose in the present work an extension towards an interactive application where user information can be used online to steer the segmentation result. This enables a synergistic collaboration between the operator and the underlying segmentation algorithm, thus contributing to higher segmentation accuracy, while keeping total analysis time competitive. To this end, we formalize the user interaction paradigm using a geometrical approach, where the user input is mapped to a non-cartesian space while this information is used to drive the boundary towards the position provided by the user. Additionally, we propose a shape regularization term which improves the interaction with the segmented surface, thereby making the interactive segmentation process less cumbersome. The resulting algorithm offers competitive performance both in terms of segmentation accuracy, as well as in terms of total analysis time. This contributes to a more efficient use of the existing segmentation tools in daily clinical practice. Furthermore, it compares favorably to state-of-the-art interactive segmentation software based on a 3D livewire-based algorithm.