Mixed integer evolution strategies for parameter optimization.

Evolution strategies (ESs) are powerful probabilistic search and optimization algorithms gleaned from biological evolution theory. They have been successfully applied to a wide range of real world applications. The modern ESs are mainly designed for solving continuous parameter optimization problems. Their ability to adapt the parameters of the multivariate normal distribution used for mutation during the optimization run makes them well suited for this domain. In this article we describe and study mixed integer evolution strategies (MIES), which are natural extensions of ES for mixed integer optimization problems. MIES can deal with parameter vectors consisting not only of continuous variables but also with nominal discrete and integer variables. Following the design principles of the canonical evolution strategies, they use specialized mutation operators tailored for the aforementioned mixed parameter classes. For each type of variable, the choice of mutation operators is governed by a natural metric for this variable type, maximal entropy, and symmetry considerations. All distributions used for mutation can be controlled in their shape by means of scaling parameters, allowing self-adaptation to be implemented. After introducing and motivating the conceptual design of the MIES, we study the optimality of the self-adaptation of step sizes and mutation rates on a generalized (weighted) sphere model. Moreover, we prove global convergence of the MIES on a very general class of problems. The remainder of the article is devoted to performance studies on artificial landscapes (barrier functions and mixed integer NK landscapes), and a case study in the optimization of medical image analysis systems. In addition, we show that with proper constraint handling techniques, MIES can also be applied to classical mixed integer nonlinear programming problems.

Integrated contour detection and pose estimation for fluoroscopic analysis of knee implants.

With fluoroscopic analysis of knee implant kinematics the implant contour must be detected in each image frame, followed by estimation of the implant pose. With a large number of possibly low-quality images, the contour detection is a time-consuming bottleneck. The present paper proposes an automated contour detection method, which is integrated in the pose estimation. In a phantom experiment the automated method was compared with a standard method, which uses manual selection of correct contour parts. Both methods demonstrated comparable precision, with a minor difference in the Y-position (0.08 mm versus 0.06 mm). The precision of each method was so small (below 0.2 mm and 0.3 degrees) that both are sufficiently accurate for clinical research purposes. The efficiency of both methods was assessed on six clinical datasets. With the automated method the observer spent 1.5 min per image, significantly less than 3.9 min with the standard method. A Bland-Altman analysis between the methods demonstrated no discernible trends in the relative femoral poses. The threefold increase in efficiency demonstrates that a pose estimation approach with integrated contour detection is more intuitive than a standard method. It eliminates most of the manual work in fluoroscopic analysis, with sufficient precision for clinical research purposes.

Neurological signs in relation to white matter hyperintensity volumes in memory clinic patients.

PURPOSE: To determine the frequency of neurological signs in a memory clinic population and to explore their associations with white matter hyperintensity (WMH). METHODS: We included patients with Alzheimer disease (AD; n = 210), vascular dementia (VaD; n = 34), mild cognitive impairment (MCI; n = 86) and subjective complaints (n = 153). The presence of extrapyramidal and unilateral signs was assessed from medical charts. On MRI, WMH volumes were extracted automatically. RESULTS: Extrapyramidal signs were found in 10% and unilateral signs in 12% of the patients. Age- and sex-adjusted extrapyramidal signs occurred more often in VaD compared to patients with subjective complaints. Unilateral signs were more prevalent in all groups compared to patients with subjective complaints. Two-way analysis of variance (ANOVA) with WMH as the dependent variable showed a main effect of diagnosis (p < 0.001), but not of extrapyramidal signs (p = 0.62). In contrast, 2-way ANOVA showed main effects of diagnosis (p < 0.001) and unilateral signs (p = 0.001). Furthermore, there was an interaction between these factors (p = 0.04); if unilateral signs were present, patients with subjective complaints and VaD showed more WMH, whereas there was no relation in AD and MCI. CONCLUSION: Extrapyramidal and unilateral signs are common in memory clinic patients, but are only modestly related to WMH.

Handling modular hip implants in model-based RSA: combined stem-head models.

Migration measurements of hip prostheses using marker-based Roentgen stereophotogrammetric analysis (RSA) require the attachment of markers to the prostheses. The model-based approach, which does not require these markers, is, however, less precise. One of the reasons may be the fact that the spherical head has not been modelled. Therefore, we added a 3D surface model of the spherical head and estimated the position and orientation of the combined stem-head model. The new method using a combined stem-head model was compared in a phantom study on five prostheses (of different types) and in a clinical study using double examinations of implanted hip prostheses, with two existing methods: a standard model-based approach and one using elementary geometrical shapes. The combined model showed the highest precision for the rotation about the longitudinal axis in the phantom experiments. With a standard deviation of 0.69 degrees it showed a significant improvement (p=0.02) over the model-based approach (0.96 degrees ) on the phantom data, but no improvement on the clinical data. Overall, the use of elementary geometrical shapes was worse with respect to the model-based approach, with a standard deviation of 1.02 degrees on the phantom data and 0.79 degrees on the clinical data. This decrease in precision was significant (p<0.01) on the clinical data. With relatively small differences in the other migration directions, these results demonstrate that the new method with a combined stem-head model can be a useful alternative to the standard model-based approach.

Image-based RSA: Roentgen stereophotogrammetric analysis based on 2D-3D image registration.

Image-based Roentgen stereophotogrammetric analysis (IBRSA) integrates 2D-3D image registration and conventional RSA. Instead of radiopaque RSA bone markers, IBRSA uses 3D CT data, from which digitally reconstructed radiographs (DRRs) are generated. Using 2D-3D image registration, the 3D pose of the CT is iteratively adjusted such that the generated DRRs resemble the 2D RSA images as closely as possible, according to an image matching metric. Effectively, by registering all 2D follow-up moments to the same 3D CT, the CT volume functions as common ground. In two experiments, using RSA and using a micromanipulator as gold standard, IBRSA has been validated on cadaveric and sawbone scapula radiographs, and good matching results have been achieved. The accuracy was: |mu |< 0.083 mm for translations and |mu| < 0.023 degrees for rotations. The precision sigma in x-, y-, and z-direction was 0.090, 0.077, and 0.220 mm for translations and 0.155 degrees , 0.243 degrees , and 0.074 degrees for rotations. Our results show that the accuracy and precision of in vitro IBRSA, performed under ideal laboratory conditions, are lower than in vitro standard RSA but higher than in vivo standard RSA. Because IBRSA does not require radiopaque markers, it adds functionality to the RSA method by opening new directions and possibilities for research, such as dynamic analyses using fluoroscopy on subjects without markers and computer navigation applications.

Time continuous detection of the left ventricular long axis and the mitral valve plane in 3-D echocardiography.

Automated segmentation approaches for the left ventricle (LV) in 3-D echocardiography (3DE) often rely on manual initialization. So far, little effort has been put into automating the initialization procedure to get to a fully automatic segmentation approach. We propose a fully automatic method for the detection of the LV long axis (LAX) and the mitral valve plane (MVP) over the full cardiac cycle, for the initialization of segmentation algorithms in 3DE. Our method exploits the cyclic motion of the LV and therefore detects salient structures in a time-continuous way. Probabilities to candidate LV center points are assigned through a Hough transform for circles. The LV LAX is detected by combining dynamic programming detections on these probabilities in 3-D and 2D + time to obtain a time continuous solution. Subsequently, the mitral valve plane is detected in a projection of the data on a plane through the previously detected LAX. The method easily adjusts to different acquisition routines and combines robustness with good accuracy and low computational costs. Automatic detection was evaluated using patient data acquired with the fast rotating ultrasound (FRU) transducer (n=11 patients) and with the Philips Sonos 7500 ultrasound system (Philips Medical Systems, Andover, MA, USA), with the X4 matrix transducer (n=14 patients). For the FRU-transducer data, the LAX was estimated with a distance error of 2.85+/-1.70 mm (mean+/-SD) and an angle of 5.25+/-3.17 degrees; the mitral valve plane was estimated with a distance of -1.54+/-4.31 mm. For the matrix data, these distances were 1.96+/-1.30 mm with an angle error of 5.95+/-2.11 and -1.66+/-5.27 mm for the mitral valve plane. These results confirm that the method is very suitable for automatic detection of the LV LAX and MVP. It provides a basis for further automatic exploration of the LV and could therefore serve as a replacement of manual initialization of 3-D segmentation approaches.

Clinical validation of an automated vessel-segmentation software of the extracranial-carotid arteries based on 3D-MRA: a prospective study.

OBJECTIVES: To determine the accuracy of automated vessel-segmentation software for vessel-diameter measurements based on three-dimensional contrast-enhanced magnetic resonance angiography (3D-MRA). METHOD: In 10 patients with high-grade carotid stenosis, automated measurements of both carotid arteries were obtained with 3D-MRA by two independent investigators and compared with manual measurements obtained by digital subtraction angiography (DSA) and 2D maximum-intensity projection (2D-MIP) based on MRA and duplex ultrasonography (US). In 42 patients undergoing carotid endarterectomy (CEA), intraoperative measurements (IOP) were compared with postoperative 3D-MRA and US. RESULTS: Mean interoperator variability was 8% for measurements by DSA and 11% by 2D-MIP, but there was no interoperator variability with the automated 3D-MRA analysis. Good correlations were found between DSA (standard of reference), manual 2D-MIP (rP=0.6) and automated 3D-MRA (rP=0.8). Excellent correlations were found between IOP, 3D-MRA (rP=0.93) and US (rP=0.83). CONCLUSION: Automated 3D-MRA-based vessel segmentation and quantification result in accurate measurements of extracerebral-vessel dimensions.

Accuracy and reproducibility of fast fractional flow reserve computation from invasive coronary angiography.

Fractional flow reserve (FFR) guided percutaneous coronary intervention (PCI) is associated with favourable outcome compared with revascularization based on angiographic stenosis severity alone. The feasibility of the new image-based quantitative flow ratio (QFR) assessed from 3D quantitative coronary angiography (QCA) and thrombolysis in myocardial infarction (TIMI) frame count using three different flow models has been reported recently. The aim of the current study was to assess the accuracy, and in particular, the reproducibility of these three QFR techniques when compared with invasive FFR. QFR was derived (1) from adenosine induced hyperaemic coronary angiography images (adenosine-flow QFR [aQFR]), (2) from non-hyperemic images (contrast-flow QFR [cQFR]) and (3) using a fixed empiric hyperaemic flow [fixed-flow QFR (fQFR)]. The three QFR values were calculated in 17 patients who prospectively underwent invasive FFR measurement in 20 vessels. Two independent observers performed the QFR analyses. Mean difference, standard deviation and 95% limits of agreement (LOA) between invasive FFR and aQFR, cQFR and fQFR for observer 1 were: 0.01 ± 0.04 (95% LOA: -0.07; 0.10), 0.01 ± 0.05 (95% LOA: -0.08; 0.10), 0.01 ± 0.04 (95% LOA: -0.06; 0.08) and for observer 2: 0.00 ± 0.03 (95% LOA: -0.06; 0.07), -0.01 ± 0.03 (95% LOA: -0.07; 0.05), 0.00 ± 0.03 (95% LOA: -0.06; 0.05). Values between the 2 observers were (to assess reproducibility) for aQFR: 0.01 ± 0.04 (95% LOA: -0.07; 0.09), for cQFR: 0.02 ± 0.04 (95% LOA: -0.06; 0.09) and for fQFR: 0.01 ± 0.05 (95% LOA: -0.07; 0.10). In a small number of patients we showed good accuracy of three QFR techniques (aQFR, cQFR and fQFR) to predict invasive FFR. Furthermore, good inter-observer agreement of the QFR values was observed between two independent observers.

Patients with type 3 severe von Willebrand disease are not protected against atherosclerosis: results from a multicenter study in 47 patients.

BACKGROUND: The results of a number of studies in pigs and mice suggest that absence of von Willebrand factor (vWF) protects against the development of atherosclerosis. We studied whether patients with a complete deficiency of vWF (type 3 von Willebrand disease [vWD]) develop fewer atherosclerotic vessel wall changes than healthy controls. METHODS AND RESULTS: This study included 47 individuals with type 3 vWD and 84 healthy controls. Early atherosclerotic changes were assessed by measuring the thickness of the intima-media in the carotid and femoral arteries by B-mode ultrasonography. Advanced atherosclerotic changes were quantified by summing the maximal thickness of atherosclerotic plaques in the carotid and femoral arteries and were expressed as a plaque score. Established risk factors were determined to adjust for possible differences between the groups. We found no substantial difference in intima-media thickness between vWD patients and controls (adjusted difference for carotid artery 0.007 mm, 95% CI -0.022 to 0.036 mm; femoral artery 0.069 mm, 95% CI -0.056 to 0.19 mm). Similar proportions of patients and controls had atherosclerotic plaques (19% and 17%, respectively). No difference was found in the plaque score between groups (adjusted difference -0.22 mm, 95% CI -0.69 to 0.26). Among vWD patients, we found no effect of treatment with vWF concentrates on intima-media thickness or plaque score. CONCLUSIONS: The results of this study indicate that vWF does not play a substantial role in human atherogenesis.

Optimizing the automatic segmentation of the left ventricle in magnetic resonance images.

Automatic segmentation of the left ventricular (LV) myocardial borders in cardiovascular MR (CMR) images allows a significant speed-up of the procedure of quantifying LV function, and improves its reproducibility. The automated boundary delineation is usually based on a set of parameters that define the algorithms. Since the automatic segmentation algorithms are usually sensitive to the image quality and frequently depend heavily on the acquisition protocol, optimizing the parameters of the algorithm for such different protocols may be necessary to obtain optimal results. In other words, using a default set of parameters may be far from optimal for different scanners or protocols. For the MASS-software, for example, this means that a total of 14 parameters need to be optimized. This optimization is a difficult and labor-intensive process. To be able to more consistently and rapidly tune the parameters, an automated optimization system would be extremely desirable. In this paper we propose such an approach, which is based on genetic algorithms (GAs). The GA is an unsupervised iterative tool that generates new sets of parameters and converges toward an optimal set. We implemented and compared two different types of the genetic algorithms: a simple GA (SGA) and a steady state GA (2SGA). The difference between these two algorithms lies in the characteristics of the generated populations: "nonoverlapping populations" and "overlapping populations," respectively "nonoverlapping" population means that the two populations are disjoint, and "overlapping" means that the best parameters found in the previous generation are included in the present population. The performance of both algorithms was evaluated on twenty routinely obtained short-axis examinations (eleven examinations acquired with a steady-state free precession pulse sequence, and nine examinations with a gradient echo pulse sequence). The optimal parameters obtained with the GAs were used for the LV myocardial border delineation. Finally, the automatically outlined contours were compared to the gold standard--manually drawn contours by experts. The result of the comparison was expressed as a degree of similarity after a processing time of less than 72 h to a 59.5% of degree of similarity for SGA and a 66.7% of degree of similarity for 2SGA. In conclusion, genetic algorithms are very suitable to automatically tune the parameters of a border detection algorithm. Based on our data, the 2SGA was more suitable than the SGA method. This approach can be generalized to other optimization problems in medical image processing.

A new type of model-based Roentgen stereophotogrammetric analysis for solving the occluded marker problem.

Roentgen stereophotogrammetric analysis (RSA) measures micromotion of an orthopaedic implant with respect to its surrounding bone. A problem in RSA is that the markers are sometimes overprojected by the implant itself. This study describes the so-called Marker Configuration Model-based RSA (MCM-based RSA) that is able to measure the pose of a rigid body in situations where less than three markers could be detected in both images of an RSA radiograph. MCM-based RSA is based on fitting a Marker Configuration model (MC-model) to the projection lines from the marker projection positions in the image to their corresponding Roentgen foci. An MC-model describes the positions of markers relative to each other and is obtained using conventional RSA. We used data from 15 double examinations of a clinical study of total knee prostheses and removed projections of the three tibial component markers, simulating occlusion of markers. The migration of the tibial component with respect to the bone, which should be zero, for the double examination is a measure of the accuracy of algorithm. With the new algorithm, it is possible to estimate the pose of a rigid body of which one or two markers are occluded in one of the images of the RSA radiograph with high accuracy as long as a proper MC-model of the markers in the rigid body is available. The new algorithm makes RSA more robust for occlusion of markers. This improves the results of clinical RSA studies because the number of lost RSA follow-up moments is reduced.

Lipoprotein(a) and ultrasonographically determined early atherosclerotic changes in the carotid and femoral artery.

Recent studies suggest that high lipoprotein(a) [Lp(a)] plasma levels are associated with symptomatic ischemic cardiovascular disease. We examined whether Lp(a) plasma levels are associated with early atherosclerotic vessel wall changes in a group of asymptomatic subjects. In a group of 142 asymptomatic men, the intima-media thickness (IMT) in the common carotid artery, the carotid bifurcation and the common femoral artery was determined by B-mode ultrasonography. In addition to Lp(a) levels, established risk factors, such as blood pressure and cholesterol levels were determined. Lipoprotein(a) values ranged from 2 mg L(-1) to 900 mg L(-1) (median 145 mg L(-1)). Linear regression analysis showed a clear association of IMT with the established risk factors but not with Lp(a) [regression coefficient carotid artery -0.0003, 95% confidence interval (CI) -0.002-0.001; regression coefficient femoral artery -0.0003, 95% CI -0.004-0.003]. We found no increased intima-media thickness in the carotid or femoral artery at high levels of Lp(a). Lipoprotein(a) levels are not associated with early atherosclerotic vessel wall changes in the carotid or femoral artery.

Arterial compliance adds to conventional risk factors for prediction of angiographic coronary artery disease.

BACKGROUND: Arterial compliance is related to left ventricular hypertrophy and risk for cardiovascular disease events; however, its association with coronary artery stenosis remains uncertain. We sought to assess the relation between lower extremity arterial compliance and presence of angiographically defined coronary artery disease. METHODS: Lower extremity arterial compliance was measured with the use of a noninvasive air plethysmography technique in 376 subjects undergoing routine diagnostic coronary angiography. RESULTS: Measures of calf arterial compliance were significantly associated with the presence of one or more stenoses > or =50% compared with no stenoses, even after adjustment for age, sex, smoking, diabetes, hypertension, hypercholesterolemia, and obesity (P =.03). Measures of thigh arterial compliance were also lower in subjects with disease, although this association did not reach statistical significance (P =.07). Receiver operator curves illustrate the incremental predictive ability of calf arterial compliance over and above age, sex, and conventional risk factors. CONCLUSIONS: Lower extremity arterial compliance is associated with presence of significant coronary stenoses in a cardiac catheterization laboratory referral population. This observation lends support for additional efforts to determine the utility of vascular stiffness measures in both clinical and pre-clinical populations to guide treatment and prevention efforts.

Arterial intima-media thickness and its relationship with cardiovascular disease and atherosclerosis: a possible contribution of medium-sized arteries.

The intima-media thickness (IMT) of the arterial wall, measured by B-mode ultrasonography, has been related to cardiovascular disease and atherosclerosis in many studies, most of which have considered carotid and common femoral arteries. No study has related the IMT of medium-sized arteries (e.g. brachial) to the atherosclerotic process. The aim of this study was to evaluate whether more peripheral arteries (e.g. brachial, superficial femoral and popliteal) represent a good model in studies of arterial wall IMT, as well as carotid artery. Twenty-six male patients with demonstrated coronary artery disease (median age 63 yrs; range 41-70) and twenty-four male controls (median age 62 yrs; range 53-74) were studied. The far-wall IMT of the common carotid, brachial, superficial femoral and popliteal arteries was measured by high-resolution B-mode ultrasonography. Its relationship with ischemic cardiovascular disease and the presence of atherosclerotic plaques in all the arteries was also evaluated. Mean IMT values were higher in patients than controls, more so in superficial femoral artery (IMT difference 0.08 mm [95% confidence interval 0.02 - 0.14]) and popliteal artery (0.08 [0.01-0.15]) than in common carotid artery (0.02 [-0.07-0.11]) or brachial artery (0.01[-0.01-0.03]). The difference did not change after adjustment for age. The prevalence of plaques was consistently higher in patients than controls, except for brachial artery, in which no plaque was found either in patients and controls. In all the arteries except brachial we found an association between increase in IMT and number of plaques. Age was strongly related to the presence of plaques in the carotid artery, and less markedly in superficial femoral and popliteal arteries. The brachial artery does not seem a good model to study atherosclerosis by ultrasound measurements of arterial IMT, whereas superficial femoral and popliteal arteries might be chosen for these studies besides carotid artery.

A new model-based RSA method validated using CAD models and models from reversed engineering.

Roentgen stereophotogrammetric analysis (RSA) was developed to measure micromotion of an orthopaedic implant with respect to its surrounding bone. A disadvantage of conventional RSA is that it requires the implant to be marked with tantalum beads. This disadvantage can potentially be resolved with model-based RSA, whereby a 3D model of the implant is used for matching with the actual images and the assessment of position and rotation of the implant. In this study, a model-based RSA algorithm is presented and validated in phantom experiments. To investigate the influence of the accuracy of the implant models that were used for model-based RSA, we studied both computer aided design (CAD) models as well as models obtained by means of reversed engineering (RE) of the actual implant. The results demonstrate that the RE models provide more accurate results than the CAD models. If these RE models are derived from the very same implant, it is possible to achieve a maximum standard deviation of the error in the migration calculation of 0.06 mm for translations in x- and y-direction and 0.14 mm for the out of plane z-direction, respectively. For rotations about the y-axis, the standard deviation was about 0.1 degrees and for rotations about the x- and z-axis 0.05 degrees. Studies with clinical RSA-radiographs must prove that these results can also be reached in a clinical setting, making model-based RSA a possible alternative for marker-based RSA.

Accurate object localization in gray level images using the center of gravity measure: accuracy versus precision.

A widely used subpixel precision estimate of an object center is the weighted center of gravity (COG). We derive three maximum-likelihood estimators for the variance of the two-dimensional (2-D) COG as a function of the noise in the image. We assume that the noise is additive, Gaussian distributed and independent between neighboring pixels. Repeated experiments using 2500 generated 2-D bell-shaped markers superimposed with an increasing amount of Gaussian noise were performed, to compare the three approximations. The error of the most exact approximative variance estimate with respect to true variance was always less than 5% of the latter. This deviation decreases with increasing signal-to-noise ratio. Our second approximation to the variance estimate performed better than the third approximation, which was originally presented by Oron et al. by up to a factor approximately 10. The difference in performance between these two approximations increased with an increasing misplacement of the window in which the COG was calculated with respect to the real COG.

Relationship of peripheral arterial compliance and standard cardiovascular risk factors.

Abnormalities of peripheral arterial compliance are clinically useful markers of atherosclerosis and risk of vascular events. Local peripheral arterial compliance can be easily and accurately assessed in the clinic by computer-controlled pulse volume recordings (air plethysmography). The purpose of this study was to investigate the relationship between clinical cardiovascular risk factors, a surrogate of atherosclerotic burden, and peripheral arterial compliance in the thigh and calf determined by quantification of local pulse volume recordings in patients undergoing coronary angiography. Peripheral arterial compliance in the thigh and calf was measured in 346 patients undergoing diagnostic cardiac catheterization at 4 centers. Demographic and cardiovascular risk factor data were collected, and their relationship to local arterial compliance examined using a new device that assesses maximal local arterial volume change in an extremity segment. Pulse volume recordings detected decreased local arterial compliance in the thigh associated with a history of hypertension (p < 0.0001), diabetes mellitus (p = 0.0001), and hyperlipidemia (p = 0.0007). In the calf, this arterial compliance measure was associated with a history of hypertension (p < 0.0001) and diabetes mellitus (p = 0.002). Females had lower arterial compliance than males in the thigh (p = 0.003) and calf (p < 0.0001). Limited evidence of lower arterial compliance in the thigh was found for those with obesity (p = 0.07). This procedure also demonstrated that subjects with multiple cardiovascular risk factors had lower arterial compliance in the thigh than subjects with no or 1 risk factor (p = 0.0001). Peripheral arterial compliance determined by air plethysmography is strongly associated with standard cardiovascular risk factors. The noninvasive measurement of local arterial compliance by regional pulse volume recording may be a useful adjunct for cardiovascular risk stratification early in the course of the disease as well as for monitoring vascular response to therapy.

Evaluation of three pose estimation algorithms for model-based roentgen stereophotogrammetric analysis.

Model-based roentgen stereophotogrammetric analysis (RSA) uses a three-dimensional surface model of an implant in order to estimate accurately the pose of that implant from a stereo pair of roentgen images. The technique is based on minimization of the difference between the actually projected contour of an implant and the virtually projected contour of a model of that same implant. The advantage of model-based RSA over conventional marker-based RSA is that it is not necessary to attach markers to the implant. In this paper, three pose estimation algorithms for model-based RSA are evaluated. The algorithms were assessed on the basis of their sensitivities to noise in the actual contour, to the amount of drop-outs in the actual contour, to the number of points in the actual contour and to shrinkage or expansion of the actual contour. The algorithms that were studied are the iterative inverse perspective matching (IIPM) algorithm, an algorithm based on minimization of the difference (DIF) between the actual contour and the virtual contour, and an algorithm based on minimization of the non-overlapping area (NOA) between the actual and virtual contour. The results of the simulation and phantom experiments show that the NOA algorithm does not fulfil the high accuracy that is necessary for model-based RSA. The IIPM and DIF algorithms are robust to the different distortions, making model-based RSA a possible replacement for marker-based RSA.

Quantitative analysis of vascular images, in particular of abdominal aorta aneurysms from 3D CTA data sets.

This article presents a combination of well known image processing techniques to automatically segment CTA images of the Abdominal Aortic Aneurysm. Current results are that about 80% of the contours need no manual corrections. The remaining 20% fail due to calcified plaque close to the lumen border. After correction a 3D surface model is created from the 2D contours which is used as input for flow simulations and for parameter extraction of the AAA by clinicians for selecting the proper size and shape endograft, and to plan the placement procedure of this endograft in the patient.

Mitral valve regurgitation: accurate blood flow quantification with MRI.

BACKGROUND: The quantification of transvalvular blood flow through the mitral valve (MV) and regurgitant flow in particular is difficult with echocardiography, which is the method of choice to diagnose patients selected for valve repair or replacement. With magnetic resonance imaging, information on the intraventricular blood flow can be obtained. Several scanning techniques have attempted to assess the regurgitant flow. These techniques either do not directly assess the complete flow through the MV, or they do not measure the flow at the location of the valve. AIM: To investigate the accuracy of a novel method using three-directional velocity-encoded MRI to acquire the transvalvular blood flow directly from the intraventricular blood flow field, also representing the regurgitant flow during systole. METHODS: Ten volunteers without cardiac valvular disease were recruited. The transvalvular MV flow volume was measured with three-directional velocity-encoded MRI (3-dir MV flow). RESULTS: The transvalvular flow measurements correlate very well with the flow measured in the aorta (rp=0.92, p<0.01). The small differences (mean -5±7 ml) are insignificant (p=0.06) and demonstrate the high accuracy of the new method. Intra- and inter-observer studies showed non-significant mean differences of 0.9±5.1 ml and 1.3±5.6 ml, respectively, thereby proving the high reproducibility. CONCLUSION: Three-directional velocity-encoded MRI is a patient-friendly and easy-to-use method suitable for quantifying the regurgitant MV flow in clinical practice.