Biomedical subjects

P Suetens

Publications and source records attributed to P Suetens.

At least 19 recordsLinked to original sources

PubMed2001

Automated segmentation of multiple sclerosis lesions by model outlier detection.

This paper presents a fully automated algorithm for segmentation of multiple sclerosis (MS) lesions from multispectral magnetic resonance (MR) images. The method performs intensity-based tissue classification using a stochastic model for normal brain images and simultaneously detects MS lesions as outliers that are not well explained by the model. It corrects for MR field inhomogeneities, estimates tissue-specific intensity models from the data itself, and incorporates contextual information in the classification using a Markov random field. The results of the automated method are compared with lesion delineations by human experts, showing a high total lesion load correlation. When the degree of spatial correspondence between segmentations is taken into account, considerable disagreement is found, both between expert segmentations, and between expert and automatic measurements.

Algorithms↗

PubMed2001

An iterative maximum-likelihood polychromatic algorithm for CT.

A new iterative maximum-likelihood reconstruction algorithm for X-ray computed tomography is presented. The algorithm prevents beam hardening artifacts by incorporating a polychromatic acquisition model. The continuous spectrum of the X-ray tube is modeled as a number of discrete energies. The energy dependence of the attenuation is taken into account by decomposing the linear attenuation coefficient into a photoelectric component and a Compton scatter component. The relative weight of these components is constrained based on prior material assumptions. Excellent results are obtained for simulations and for phantom measurements. Beam-hardening artifacts are effectively eliminated. The relation with existing algorithms is discussed. The results confirm that improving the acquisition model assumed by the reconstruction algorithm results in reduced artifacts. Preliminary results indicate that metal artifact reduction is a very promising application for this new algorithm.

Algorithms↗

PubMed2001

Parametric study of the peak negative acoustic pressure distribution within the image plane of a phased array transducer.

A study is presented in which the influence of pulse and transducer characteristics on the distribution of the peak negative acoustic pressure within the image plane of a phased array transducer is shown. First, the influence of the center frequency, the bandwidth of the transmitted pulse, the position of the focus point, and the properties of the phased array are investigated by computer simulation. These simulations show that both the maximal amplitude and the homogeneity of the distribution of the peak negative pressure within a sector image depend on these parameters. Most remarkably, they show that the maximal peak pressure does not necessarily occur on the geometrical symmetry axis of the array transducer but can occur at the outermost lateral lines of a sector image. Second, this effect is demonstrated experimentally by hydrophone measurements of the sound field produced by a standard, clinically used, 2.5-MHz phased array transducer. The experimental results agree well with the theoretical model. The relevance of this study for contrast echography and the safety measurements of ultrasound is discussed.

Acoustics↗

PubMed2000

An image-based 3D planning environment for hemicricolaryngectomy and reconstruction by tracheal autotransplantation.

OBJECTIVE: A new conservation method for unilateral glottal cancer with significant subglottal extension consists of a hemicricolaryngectomy and subsequent reconstruction of the laryngeal defect with a transferable patch of revascularized cervical trachea. In order to restore the three crucial functions of the larynx - airway patency, speech, and swallowing--the reconstructed larynx should resemble a situation with one paralyzed and one intact vocal fold at the glottal level. We hypothesize that this result can be achieved when the tracheal patch meets a typical surgical constraint--the so-called paramedian position at the glottal level - and when the patch is rigidly applied. This hypothesis was tested by developing and using an image-based planning system in a number of situations. MATERIALS AND METHODS: An image-based surgery simulator was developed according to the working hypothesis. To validate the transfer from planning towards actual surgery, a cadaver study was set up. Based on a patient case-study, the sufficiency of the hypothesis for obtaining good functional results was evaluated. Finally, post-operative images of patients who were doing well after being operated without pre-operative planning were compared with retrospective planning to check the necessity of the hypothesis in obtaining good functionality of the neolarynx. RESULTS: We were able to design an efficient surgery simulator. The transfer from planning towards actual surgery is accurate, and satisfactory functional results were obtained from the patient case-study. The results from the retrospective planning showed that the hypothesis is not necessary to obtain good functional results, and that the smallest possible resection is preferred from the oncological point of view. CONCLUSION: The working hypothesis is sufficient for good functional results, but is unnecessary and often undesirable from an oncological point of view. Important lessons are drawn concerning our long-term goal of pre-operative determination of the shape and position of a tracheal patch suitable for optimizing the (neo-) laryngeal functions.

Cadaver↗

PubMed2000

Remote myocardial dysfunction after acute anterior myocardial infarction: impact of left ventricular shape on regional function: a magnetic resonance myocardial tagging study.

OBJECTIVES: We sought to evaluate regional morphology and function in patients in their first week after having a reperfused anterior myocardial infarction (MI) using magnetic resonance (MR) myocardial tagging. BACKGROUND: The mechanism of myocardial dysfunction in the remote, noninfarct-related regions is an unresolved issue to date. METHODS: Sixteen patients with a first reperfused transmural anterior MI were studied with MR tagging at 5 +/- 2 days after the event, and the results were compared with those of an age-matched control group regions. The left ventricle (LV) was divided into infarct, adjacent and remote regions. Magnetic resonance tagging provided information on the regional ventricular morphology and function. RESULTS: Morphologically, an increase of the circumferential radius of curvature was found in the remote myocardium, whereas the longitudinal radius of curvature was increased in all regions of the LV. A significant increase in apical sphericity was also found. A significant reduction in strain and function was found not only in the infarct region, but also in the adjacent and remote myocardium. The loss in regional ejection fraction in the remote myocardium (61.4 +/- 11.7% in patients vs. 68.7 +/- 10.0% in control subjects, p < 0.0001) was related to a significant reduction of the longitudinal and circumferential strain, whereas systolic wall thickening was preserved. CONCLUSIONS: Remote myocardial dysfunction contributes significantly to the loss in global ventricular function. This could be secondary to morphologic changes in the infarct region, leading to an increased systolic longitudinal wall stress without loss of intrinsic contractility in the remote regions.

Adult↗

PubMed2000

High frame rate myocardial integrated backscatter. Does this change our understanding of this acoustic parameter?

AIMS: Integrated backscatter (IB) and its cyclic variation (CV) derived from radio-frequency (RF) data have been used as parameters to attempt myocardial tissue characterization. Prior imaging systems used to measure IB and its CV typically acquired data at frame rates of 20-30 Hz and at a resolution of 6-8 bits. If changes in IB levels are in part related to specific short-lived events, occurring within the cardiac cycle, this frame rate and resolution could have been too low to resolve adequately what might be a more complex data set. METHODS AND RESULTS: To investigate this possibility, we acquired real time two-dimensional (2D) myocardial IQ data (the 'in-phase quadrature' sampled RF data) at high frame rate (> 100 Hz), high dynamic resolution (theoretical 19-bit) and a sector angle of 20 degrees. Several consecutive heart cycles of myocardial data were acquired from individual cardiac walls in five closed chest dogs and 10 healthy, young volunteers at normal heart rates. On the reconstructed RF data regions of interest were indicated, and IB and its CV were calculated. The extracted high frame rate curves showed that the CV of IB is not a smooth sinusoidal-like curve, but is made up of multiple reproducible peaks and troughs with local minima and maxima which are temporally related to active or passive mechanical events, i.e. systolic contraction, early ventricular relaxation and ventricular filling due to atrial contraction. CONCLUSIONS: This study shows that increasing the rate of real-time RF data acquisition results in a more complex, reproducible IB curve. The resolved maxima and minima in IB levels are related to specific phases of the myocardial contraction. Furthermore, spectral analysis showed that IB curves acquired at normal heart rates contain information up to 40 Hz. Hence, cardiac imaging data sets used to analyse regional myocardial function obtained at frequencies lower than 80 frames per second can contain aliased information.

Adult↗

PubMed2000

Regional strain and strain rate measurements by cardiac ultrasound: principles, implementation and limitations.

The non-invasive quantification of regional myocardial function is an important goal in clinical cardiology. Myocardial thickening/thinning indices is one method of attempting to define regional myocardial function. A new ultrasonic method of quantifying regional deformation has been introduced based on the principles of 'strain' and 'strain rate' imaging. These new imaging modes introduce concepts derived from mechanical engineering which most echocardiographers are not familiar with. In order to maximally exploit these new techniques, an understanding of what they measure is indispensable. This paper will define each of these modalities in terms of physical principles and will give an introduction to the principles of data acquisition and processing required to implement ultrasonic strain and strain rate imaging. In addition, the current status of development of the technique and its limitations will be discussed, together with examples of potential clinical applications.

Biomedical Engineering↗

PubMed2000

Are changes in myocardial integrated backscatter restricted to the ischemic zone in acute induced ischemia? An in vivo animal study.

Integrated backscatter (IB) from a myocardial region, calculated from radiofrequency echocardiographic data, has been proposed as a useful parameter for investigating changes in myocardial tissue induced by ischemia. In 10 closed-chest dogs, 5 minutes of myocardial ischemia was induced by either a proximal occlusion of the circumflex coronary artery (CX) (5 dogs), resulting in extensive ischemia in the posterior wall, or by occluding the distal CX vessel (5 dogs) to produce a small localized ischemic zone in the posterior wall. High-resolution digital radiofrequency data from the whole left ventricular myocardium, in the imaging plane during one complete heart cycle, were acquired with a whole-image real-time acquisition approach. Regions in the septum and posterior wall (both ischemic tissue and, in the case of distal occlusions, tissue surrounding the ischemic zone) were chosen for analysis, and IB and cyclic variation (CV) of IB were calculated. Post occlusion, an increase in mean IB values was found in the ischemic segment. However, an increase in CV was also observed in the peri-ischemic zone for the distal CX occlusion and in the septum after proximal CX occlusion. These findings show that changes in CV are not restricted to the ischemic zone but may also occur in distal myocardium. This may be explained by changes in the regional contractile state and loading conditions of the "normal" myocardium, which are altered in response to the distal ischemia.

Animals↗

PubMed1999

Robustness of integrated backscatter for myocardial tissue characterization.

Integrated backscatter (IB) has been used for ultrasonic tissue characterization. To assess the potential variables in IB measurements, we performed both theoretical simulations and in vitro phantom measurements. First, we simulated data in which the scatterer position randomly was varied. IB values for the resulting images were calculated. Second, RF data from a tissue-mimicking phantom were acquired. Third, an adapted imaging approach, based on phase insensitivity, was evaluated. For both the simulations and phantom measurements, IB showed a standard deviation of +/-20%. These large deviations can be explained by variations in interference of signals and are not related to the state of the tissue. Small deviations in position of the scatterers resulted in important variations in IB. They must be taken into account and may limit the use of IB in cardiological applications. An improvement potentially can be obtained using phase insensitivity in new ultrasound processing schemes.

Computer Simulation↗

PubMed1999

Comparative evaluation of multiresolution optimization strategies for multimodality image registration by maximization of mutual information.

Maximization of mutual information of voxel intensities has been demonstrated to be a very powerful criterion for three-dimensional medical image registration, allowing robust and accurate fully automated affine registration of multimodal images in a variety of applications, without the need for segmentation or other preprocessing of the images. In this paper, we investigate the performance of various optimization methods and multiresolution strategies for maximization of mutual information, aiming at increasing registration speed when matching large high-resolution images. We show that mutual information is a continuous function of the affine registration parameters when appropriate interpolation is used and we derive analytic expressions of its derivatives that allow numerically exact evaluation of its gradient. Various multiresolution gradient- and non-gradient-based optimization strategies, such as Powell, simplex, steepest-descent, conjugate-gradient, quasi-Newton and Levenberg-Marquardt methods, are evaluated for registration of computed tomography (CT) and magnetic resonance images of the brain. Speed-ups of a factor of 3 on average compared to Powell's method at full resolution are achieved with similar precision and without a loss of robustness with the simplex, conjugate-gradient and Levenberg-Marquardt method using a two-level multiresolution scheme. Large data sets such as 256(2) x 128 MR and 512(2) x 48 CT images can be registered with subvoxel precision in <5 min CPU time on current workstations.

Algorithms↗

PubMed1999

Predictability of a three-dimensional planning system for oral implant surgery.

OBJECTIVES: To compare 2D CT alone with 2D + 3D reconstruction for pre-operative planning of implant placement. METHODS: Spiral CT scans of 33 consecutive patients were used for both reformatted 2D and 3D computer-assisted planning. The number, site and size of implants and the occurrence of anatomical complications during planning and implant placement were statistically compared using the percentage agreement and the Kendall's correlation coefficients (tau). Although planning was performed in 33 patients, implants were only placed in 21 patients. In 11 patients surgery was based on 2D + 3D imaging and in ten patients on 2D planning. RESULTS: Agreement between planning and placement of implants was highly significant for the implant sites selected. For 2D based planning and placement, agreement reached 68% (tau = 0.94). For 2D + 3D based planning and placement, agreement attained 73% (tau = 0.89). For planning and placement of implant size based on 2D images, agreement was 31% and not significant (tau = 0.23). When based on 2D + 3D images, agreement for implant size was 44% (tau = 0.5). Agreement was not significant for anatomical complications: 69% for 2D planning and 71% for 2D + 3D planning (tau = 0.24 for 2D and tau = 0.21 for 2D + 3D). CONCLUSIONS: The 3D planning system is a reliable tool for pre-operative assessment of implant placement. Both 2D and 2D + 3D planning have a good predictability for the number and site of the implants but less so for anatomical complications. However, the 2D + 3D planning provides a better pre-operative assessment of implant size.

Adult↗

PubMed1999

Simultaneous maximum a posteriori reconstruction of attenuation and activity distributions from emission sinograms.

In order to perform attenuation correction in emission tomography an attenuation map is required. We propose a new method to compute this map directly from the emission sinogram, eliminating the transmission scan from the acquisition protocol. The problem is formulated as an optimization task where the objective function is a combination of the likelihood and an a priori probability. The latter uses a Gibbs prior distribution to encourage local smoothness and a multimodal distribution for the attenuation coefficients. Since the attenuation process is different in positron emission tomography (PET) and single photon emission tomography (SPECT), a separate algorithm for each case is derived. The method has been tested on mathematical phantoms and on a few clinical studies. For PET, good agreement was found between the images obtained with transmission measurements and those produced by the new algorithm in an abdominal study. For SPECT, promising simulation results have been obtained for nonhomogeneous attenuation due to the presence of the lungs.

Abdomen↗

PubMed1999

Automated model-based bias field correction of MR images of the brain.

We propose a model-based method for fully automated bias field correction of MR brain images. The MR signal is modeled as a realization of a random process with a parametric probability distribution that is corrupted by a smooth polynomial inhomogeneity or bias field. The method we propose applies an iterative expectation-maximization (EM) strategy that interleaves pixel classification with estimation of class distribution and bias field parameters, improving the likelihood of the model parameters at each iteration. The algorithm, which can handle multichannel data and slice-by-slice constant intensity offsets, is initialized with information from a digital brain atlas about the a priori expected location of tissue classes. This allows full automation of the method without need for user interaction, yielding more objective and reproducible results. We have validated the bias correction algorithm on simulated data and we illustrate its performance on various MR images with important field inhomogeneities. We also relate the proposed algorithm to other bias correction algorithms.

Algorithms↗

PubMed1999

Automated model-based tissue classification of MR images of the brain.

We describe a fully automated method for model-based tissue classification of magnetic resonance (MR) images of the brain. The method interleaves classification with estimation of the model parameters, improving the classification at each iteration. The algorithm is able to segment single- and multispectral MR images, corrects for MR signal inhomogeneities, and incorporates contextual information by means of Markov random Fields (MRF's). A digital brain atlas containing prior expectations about the spatial location of tissue classes is used to initialize the algorithm. This makes the method fully automated and therefore it provides objective and reproducible segmentations. We have validated the technique on simulated as well as on real MR images of the brain.

Algorithms↗

PubMed1999

Nonlinear propagation effects on broadband attenuation measurements and its implications for ultrasonic tissue characterization.

A study is presented in which the influence of the pressure amplitude of the incident pulse on the estimated frequency dependency of the attenuation coefficient is shown. First, the effect is demonstrated with a simple theoretical model for both transmission and reflection measurements. Simulations show that for both measurement techniques a high-amplitude incident pulse results in a biased estimate of the attenuation coefficient due to nonlinear interaction of the different frequency components of the incident pulse. It is shown that in transmission and reflection measurements the biases have opposite signs. The effect of bandwidth, central frequency, and phase of the incident pulse on this bias is investigated. Second, the effect is demonstrated both in vitro, using a broadband through-transmission substitution technique on a tissue mimicking gelatine phantom, and in vivo, using reflection measurements with standard clinical equipment. The experimental results agree well with the theoretical model. The relevance of this study for ultrasonic tissue characterization is shown.

Culture Techniques↗

PubMed1998

Iterative reconstruction for helical CT: a simulation study.

Iterative reconstruction algorithms for helical CT are presented. The algorithms are derived from two-dimensional reconstruction algorithms, by adapting the projector/backprojector to the helical orbit of the source, and by constraining the axial frequencies with a Gaussian sieve. Simulations have been carried out and the performance of the iterative algorithms is compared to that of filtered backprojection of synthetic (interpolated) two-dimensional sinograms. The iterative algorithms produce superior bias-noise curves. Axial resolution is superior, but disturbing edge-artefacts are introduced.

Algorithms↗

PubMed1998

An image-guided planning system for endosseous oral implants.

A preoperative planning system for oral implant surgery was developed which takes as input computed tomographies (CT's) of the jaws. Two-dimensional (2-D) reslices of these axial CT slices orthogonal to a curve following the jaw arch are computed and shown together with three-dimensional (3-D) surface rendered models of the bone and computer-aided design (CAD)-like implant models. A technique is developed for scanning and visualizing an eventual existing removable prosthesis together with the bone structures. Evaluation of the planning done with the system shows a difference between 2-D and 3-D planning methods. Validation studies measure the benefits of the 3-D approach by comparing plans made in 2-D mode only with those further adjusted using the full 3-D visualization capabilities of the system. The benefits of a 3-D approach are then evident where a prosthesis is involved in the planning. For the majority of the patients, clinically important adjustments and optimizations to the 2-D plans are made once the 3-D visualization is enabled, effectively resulting in a better plan. The alterations are related to bone quality and quantity (p < 0.05), biomechanics (p < 0.005), and esthetics (p < 0.005), and are so obvious that the 3-D plan stands out clearly (p < 0.005). The improvements often avoid complications such as mandibular nerve damage, sinus perforations, fenestrations, or dehiscences.

Dental Implantation, Endosseous↗

PubMed1997

Evaluation of manual vs semi-automated delineation of liver lesions on CT images.

In this paper we compare a semi-automated delineation method with totally manual delineation for area quantification, with respect to efficiency, quality, and intra- and interobserver variability. Liver lesions on 28 CT images were delineated by three observers, twice using completely manual delineation and twice using a semi-automated method. Quantitative comparisons were performed with respect to delineated area and time required for the delineation tasks. Subjective comparisons were performed with respect to efficiency and perceived quality of the semi-automated method. The areas obtained using semi-automated delineation were significantly smaller (11 %) than those obtained using totally manual delineation. Intraobserver and interobserver variability with the semi-automated method were approximately three times lower than with manual delineation. Efficiency of the semi-automated method was subjectively rated favorable, although further improvements are possible. With respect to quality, the semi-automated method was ranked better than the manual method in 73 % of cases.

Humans↗