Sagittal Craniosynostosis: Comparing Surgical Techniques Using 3D Photogrammetry

BACKGROUND

The aim of this study was to compare three surgical interventions for correction of sagittal synostosis-frontobiparietal remodeling (FBR), extended strip craniotomy (ESC), and spring-assisted correction (SAC)-based on three-dimensional (3D) photogrammetry and operation characteristics.

METHODS

Patients who were born between 1991 and 2019 and diagnosed with nonsyndromic sagittal synostosis who underwent FBR, ESC, or SAC and had at least one postoperative 3D photogrammetry image taken during one of six follow-up appointments until age 6 were considered for this study. Operative characteristics, postoperative complications, reinterventions, and presence of intracranial hypertension were collected. To assess cranial growth, orthogonal cranial slices and 3D photocephalometric measurements were extracted automatically and evaluated from 3D photogrammetry images.

RESULTS

A total of 322 postoperative 3D images from 218 patients were included. After correcting for age and sex, no significant differences were observed in 3D photocephalometric measurements. Mean cranial shapes suggested that postoperative growth and shape gradually normalized with higher occipitofrontal head circumference and intracranial volume values compared with normal values, regardless of type of surgery. Flattening of the vertex seems to persist after surgical correction. The authors' cranial 3D mesh processing tool has been made publicly available as a part of this study.

CONCLUSIONS

The findings suggest that until age 6, there are no significant differences among the FBR, ESC, and SAC in their ability to correct sagittal synostosis with regard to 3D photocephalometric measurements. Therefore, efforts should be made to ensure early diagnosis so that minimally invasive surgery is a viable treatment option.

CLINICAL QUESTION/LEVEL OF EVIDENCE

Therapeutic, III.

Reliability and Agreement of Automated Head Measurements From 3-Dimensional Photogrammetry in Young Children

This study aimed to assess the reliability and agreement of automated head measurements using 3-dimensional (3D) photogrammetry in young children. Specifically, the study evaluated the agreement between manual and automated occipitofrontal circumference (OFC) measurements (n = 264) obtained from 3D images of 188 patients diagnosed with sagittal synostosis using a novel automated method proposed in this study. In addition, the study aimed to determine the interrater and intrarater reliability of the automatically extracted OFC, cephalic index, and volume. The results of the study showed that the automated OFC measurements had an excellent agreement with manual measurements, with a very strong regression score ( R2 = 0.969) and a small mean difference of -0.1 cm (-0.2%). The limits of agreement ranged from -0.93 to 0.74 cm, falling within the reported limits of agreement for manual OFC measurements. High interrater and intrarater reliability of OFC, cephalic index, and volume measurements were also demonstrated. The proposed method for automated OFC measurements was found to be a reliable alternative to manual measurements, which may be particularly beneficial in young children who undergo 3D imaging in craniofacial centers as part of their treatment protocol and in research settings that require a reproducible and transparent pipeline for anthropometric measurements. The method has been incorporated into CraniumPy, an open-source tool for 3D image visualization, registration, and optimization, which is publicly available on GitHub ( https://github.com/T-AbdelAlim/CraniumPy ).

Efficiently compressing 3D medical images for teleinterventions via CNNs and anisotropic diffusion

PURPOSE

Efficient compression of images while preserving image quality has the potential to be a major enabler of effective remote clinical diagnosis and treatment, since poor Internet connection conditions are often the primary constraint in such services. This paper presents a framework for organ-specific image compression for teleinterventions based on a deep learning approach and anisotropic diffusion filter.

METHODS

The proposed method, deep learning and anisotropic diffusion (DLAD), uses a convolutional neural network architecture to extract a probability map for the organ of interest; this probability map guides an anisotropic diffusion filter that smooths the image except at the location of the organ of interest. Subsequently, a compression method, such as BZ2 and HEVC-visually lossless, is applied to compress the image. We demonstrate the proposed method on three-dimensional (3D) CT images acquired for radio frequency ablation (RFA) of liver lesions. We quantitatively evaluate the proposed method on 151 CT images using peak-signal-to-noise ratio ( PSNR ), structural similarity ( SSIM ), and compression ratio ( CR ) metrics. Finally, we compare the assessments of two radiologists on the liver lesion detection and the liver lesion center annotation using 33 sets of the original images and the compressed images.

RESULTS

The results show that the method can significantly improve CR of most well-known compression methods. DLAD combined with HEVC-visually lossless achieves the highest average CR of 6.45, which is 36% higher than that of the original HEVC and outperforms other state-of-the-art lossless medical image compression methods. The means of PSNR and SSIM are 70 dB and 0.95, respectively. In addition, the compression effects do not statistically significantly affect the assessments of the radiologists on the liver lesion detection and the lesion center annotation.

CONCLUSIONS

We thus conclude that the method has a high potential to be applied in teleintervention applications.

Multi-atlas image registration of clinical data with automated quality assessment using ventricle segmentation

Registration is a core component of many imaging pipelines. In case of clinical scans, with lower resolution and sometimes substantial motion artifacts, registration can produce poor results. Visual assessment of registration quality in large clinical datasets is inefficient. In this work, we propose to automatically assess the quality of registration to an atlas in clinical FLAIR MRI scans of the brain. The method consists of automatically segmenting the ventricles of a given scan using a neural network, and comparing the segmentation to the atlas ventricles propagated to image space. We used the proposed method to improve clinical image registration to a general atlas by computing multiple registrations - one directly to the general atlas and others via different age-specific atlases - and then selecting the registration that yielded the highest ventricle overlap. Finally, as an example application of the complete pipeline, a voxelwise map of white matter hyperintensity burden was computed using only the scans with registration quality above a predefined threshold. Methods were evaluated in a single-site dataset of more than 1000 scans, as well as a multi-center dataset comprising 142 clinical scans from 12 sites. The automated ventricle segmentation reached a Dice coefficient with manual annotations of 0.89 in the single-site dataset, and 0.83 in the multi-center dataset. Registration via age-specific atlases could improve ventricle overlap compared to a direct registration to the general atlas (Dice similarity coefficient increase up to 0.15). Experiments also showed that selecting scans with the registration quality assessment method could improve the quality of average maps of white matter hyperintensity burden, instead of using all scans for the computation of the white matter hyperintensity map. In this work, we demonstrated the utility of an automated tool for assessing image registration quality in clinical scans. This image quality assessment step could ultimately assist in the translation of automated neuroimaging pipelines to the clinic.

Contribution of general practitioners and sexual health centres to sexually transmitted infection consultations in five Dutch regions using laboratory data of Chlamydia trachomatis testing

Effective sexually transmitted infection (STI) control requires opportunities for appropriate testing, counselling and treatment. In the Netherlands, people may attend general practitioners (GPs) and sexual health centres (SHCs; also known as STI clinics) for STI consultations. We assessed the contribution of GPs and SHCs to STI consultations in five Dutch regions with different urbanization levels, using data of urogenital Chlamydia trachomatis (CT) testing. Data (2011–2016) were retrieved from laboratories, aggregated by gender and age group (15–24 and 25–64 years). Results show that test rates and GP contribution varied widely between regions. GP contribution decreased over time in Amsterdam (60–48%), Twente (79–61%), Maastricht (60–50%) and Northeast-Netherlands (82–77%), but not in Rotterdam (65–67%). Decreases resulted from increases in SHC test rates and slight decreases in GP test rates. GPs performed more tests for women and those aged 25–64 years compared to SHCs (relative risks ranging from 1.49 to 4.76 and 1.58 to 7.43, respectively). The average yearly urogenital CT positivity rate was 9.2% at GPs and 10.7% at SHCs. Overall, GPs accounted for most STI consultations, yet SHC contribution increased. Continued focus on good quality STI care at GPs is essential, as increasing demands for care can not be entirely covered by SHCs.

Integrating artificial intelligence into the clinical practice of radiology: challenges and recommendations

Artificial intelligence (AI) has the potential to significantly disrupt the way radiology will be practiced in the near future, but several issues need to be resolved before AI can be widely implemented in daily practice. These include the role of the different stakeholders in the development of AI for imaging, the ethical development and use of AI in healthcare, the appropriate validation of each developed AI algorithm, the development of effective data sharing mechanisms, regulatory hurdles for the clearance of AI algorithms, and the development of AI educational resources for both practicing radiologists and radiology trainees. This paper details these issues and presents possible solutions based on discussions held at the 2019 meeting of the International Society for Strategic Studies in Radiology. KEY POINTS: • Radiologists should be aware of the different types of bias commonly encountered in AI studies, and understand their possible effects. • Methods for effective data sharing to train, validate, and test AI algorithms need to be developed. • It is essential for all radiologists to gain an understanding of the basic principles, potentials, and limits of AI.

Why rankings of biomedical image analysis competitions should be interpreted with care

International challenges have become the standard for validation of biomedical image analysis methods. Given their scientific impact, it is surprising that a critical analysis of common practices related to the organization of challenges has not yet been performed. In this paper, we present a comprehensive analysis of biomedical image analysis challenges conducted up to now. We demonstrate the importance of challenges and show that the lack of quality control has critical consequences. First, reproducibility and interpretation of the results is often hampered as only a fraction of relevant information is typically provided. Second, the rank of an algorithm is generally not robust to a number of variables such as the test data used for validation, the ranking scheme applied and the observers that make the reference annotations. To overcome these problems, we recommend best practice guidelines and define open research questions to be addressed in the future.

Contour segmentation of the intima, media, and adventitia layers in intracoronary OCT images: application to fully automatic detection of healthy wall regions

PURPOSE

Quantitative and automatic analysis of intracoronary optical coherence tomography images is useful and time-saving to assess cardiovascular risk in the clinical arena.

METHODS

First, the interfaces of the intima, media, and adventitia layers are segmented, by means of an original front propagation scheme, running in a 4D multi-parametric space, to simultaneously extract three non-crossing contours in the initial cross-sectional image. Second, information resulting from the tentative contours is exploited by a machine learning approach to identify healthy and diseased regions of the arterial wall. The framework is fully automatic.

RESULTS

The method was applied to 40 patients from two different medical centers. The framework was trained on 140 images and validated on 260 other images. For the contour segmentation method, the average segmentation errors were [Formula: see text] for the intima-media interface, [Formula: see text] for the media-adventitia interface, and [Formula: see text] for the adventitia-periadventitia interface. The classification method demonstrated a good accuracy, with a median Dice coefficient equal to 0.93 and an interquartile range of (0.78-0.98).

CONCLUSION

The proposed framework demonstrated promising offline performances and could potentially be translated into a reliable tool for various clinical applications, such as quantification of tissue layer thickness and global summarization of healthy regions in entire pullbacks.

Non-Rigid Registration of Liver CT Images for CT-Guided Ablation of Liver Tumors

CT-guided percutaneous ablation for liver cancer treatment is a relevant technique for patients not eligible for surgery and with tumors that are inconspicuous on US imaging. The lack of real-time imaging and the use of a limited amount of CT contrast agent make targeting the tumor with the needle challenging. In this study, we evaluate a registration framework that allows the integration of diagnostic pre-operative contrast enhanced CT images and intra-operative non-contrast enhanced CT images to improve image guidance in the intervention. The liver and tumor are segmented in the pre-operative contrast enhanced CT images. Next, the contrast enhanced image is registered to the intra-operative CT images in a two-stage approach. First, the contrast-enhanced diagnostic image is non-rigidly registered to a non-contrast enhanced image that is conventionally acquired at the start of the intervention. In case the initial registration is not sufficiently accurate, a refinement step is applied using non-rigid registration method with a local rigidity term. In the second stage, the intra-operative CT-images that are used to check the needle position, which often consist of only a few slices, are registered rigidly to the intra-operative image that was acquired at the start of the intervention. Subsequently, the diagnostic image is registered to the current intra-operative image, using both transformations, this allows the visualization of the tumor region extracted from pre-operative data in the intra-operative CT images containing needle. The method is evaluated on imaging data of 19 patients at the Erasmus MC. Quantitative evaluation is performed using the Dice metric, mean surface distance of the liver border and corresponding landmarks in the diagnostic and the intra-operative images. The registration of the diagnostic CT image to the initial intra-operative CT image did not require a refinement step in 13 cases. For those cases, the resulting registration had a Dice coefficient for the livers of 91.4%, a mean surface distance of 4.4 mm and a mean distance between corresponding landmarks of 4.7 mm. For the three cases with a refinement step, the registration result significantly improved (p<0.05) compared to the result of the initial non rigid registration method (DICE of 90.3% vs 71.3% and mean surface distance of 5.1 mm vs 11.3 mm and mean distance between corresponding landmark of 6.4 mm vs 10.2 mm). The registration of the preoperative data with the needle image in 16 cases yielded a DICE of 90.1% and a mean surface distance of 5.2 mm. The remaining three cases with DICE smaller than 80% were classified as unsuccessful registration. The results show that this is promising tool for liver image registration in interventional radiology.

An automatic registration method for pre- and post-interventional CT images for assessing treatment success in liver RFA treatment

PURPOSE

In image-guided radio frequency ablation for liver cancer treatment, pre- and post-interventional CT images are typically used to verify the treatment success of the therapy. In current clinical practice, the tumor zone in the diagnostic, preinterventional images is mentally or manually mapped to the ablation zone in the post-interventional images to decide success of the treatment. However, liver deformation and differences in image quality as well as in texture of the ablation zone and the tumor area make the mental or manual registration a challenging task. Purpose of this paper is to develop an automatic framework to register the pre-interventional image to the post-interventional image.

METHODS

The authors propose a registration approach enabling a nonrigid deformation of the tumor to the ablation zone, while keeping locally rigid deformation of the tumor area. The method was evaluated on CT images of 38 patient datasets from Erasmus MC. The evaluation is based on Dice coefficients of the liver segmentation on both the pre-interventional and post-interventional images, and mean distances between the liver segmentations. Additionally, residual distances after registration between corresponding landmarks and local mean surface distance in the images were computed.

RESULTS

The results show that rigid registration gives a Dice coefficient of 87.9%, a mean distance of the liver surfaces of 5.53 mm, and a landmark error of 5.38 mm, while non-rigid registration with local rigid deformation has a Dice coefficient of 92.2%, a mean distance between the liver segmentation boundaries near the tumor area of 3.83 mm, and a landmark error of 2.91 mm, where a part of this error can be attributed to the slice spacing in the authors' CT images.

CONCLUSIONS

This method is thus a promising tool to assess the success of RFA liver cancer treatment.

Fusion of fibrous cap thickness and wall shear stress to assess plaque vulnerability in coronary arteries: a pilot study

Purpose

Identification of rupture-prone plaques in coronary arteries is a major clinical challenge. Fibrous cap thickness and wall shear stress are two relevant image-based risk factors, but these two parameters are generally computed and analyzed separately. Accordingly, combining these two parameters can potentially improve the identification of at-risk regions. Therefore, the purpose of this study is to investigate the feasibility of the fusion of wall shear stress and fibrous cap thickness of coronary arteries in patient data.

Methods

Fourteen patients were included in this pilot study. Imaging of the coronary arteries was performed with optical coherence tomography and with angiography. Fibrous cap thickness was automatically quantified from optical coherence tomography pullbacks using a contour segmentation approach based on fast marching. Wall shear stress was computed by applying computational fluid dynamics on the 3D volume reconstructed from two angiograms. The two parameters then were co-registered using anatomical landmarks such as side branches.

Results

The two image modalities were successfully co-registered, with a mean (±SD) error corresponding to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$8.6\,\pm \,6.7\,\%$$\end{document}8.6±6.7% of the length of the analyzed region. For all the analyzed participants, the average thinnest portion of each fibrous cap was \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$129\,\pm \,69\,\upmu \text {m}$$\end{document}129±69μm, and the average WSS value at the location of the fibrous cap was \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1.46\,\pm \,1.16\,\text {Pa}$$\end{document}1.46±1.16Pa. A unique index was finally generated for each patient via the fusion of fibrous cap thickness and wall shear stress measurements, to translate all the measured parameters into a single risk map.

Conclusion

The introduced risk map integrates two complementary parameters and has potential to provide valuable information about plaque vulnerability.

An Automatic 3D Facial Landmarking Algorithm Using 2D Gabor Wavelets

In this paper, we present a novel approach to automatic 3D facial landmarking using 2D Gabor wavelets. Our algorithm considers the face to be a surface and uses map projections to derive 2D features from raw data. Extracted features include texture, relief map, and transformations thereof. We extend an established 2D landmarking method for simultaneous evaluation of these data. The method is validated by performing landmarking experiments on two data sets using 21 landmarks and compared with an active shape model implementation. On average, landmarking error for our method was 1.9 mm, whereas the active shape model resulted in an average landmarking error of 2.3 mm. A second study investigating facial shape heritability in related individuals concludes that automatic landmarking is on par with manual landmarking for some landmarks. Our algorithm can be trained in 30 min to automatically landmark 3D facial data sets of any size, and allows for fast and robust landmarking of 3D faces.

The use of atlas registration and graph cuts for prostate segmentation in magnetic resonance images

PURPOSE

An automatic method for 3D prostate segmentation in magnetic resonance (MR) images is presented for planning image-guided radiotherapy treatment of prostate cancer.

METHODS

A spatial prior based on intersubject atlas registration is combined with organ-specific intensity information in a graph cut segmentation framework. The segmentation is tested on 67 axial T2-weighted MR images in a leave-one-out cross validation experiment and compared with both manual reference segmentations and with multiatlas-based segmentations using majority voting atlas fusion. The impact of atlas selection is investigated in both the traditional atlas-based segmentation and the new graph cut method that combines atlas and intensity information in order to improve the segmentation accuracy. Best results were achieved using the method that combines intensity information, shape information, and atlas selection in the graph cut framework.

RESULTS

A mean Dice similarity coefficient (DSC) of 0.88 and a mean surface distance (MSD) of 1.45 mm with respect to the manual delineation were achieved.

CONCLUSIONS

This approaches the interobserver DSC of 0.90 and interobserver MSD 0f 1.15 mm and is comparable to other studies performing prostate segmentation in MR.

Quantification of fibrous cap thickness in intracoronary optical coherence tomography with a contour segmentation method based on dynamic programming

OBJECTIVES

Fibrous cap thickness is the most critical component of plaque stability. Therefore, in vivo quantification of cap thickness could yield valuable information for estimating the risk of plaque rupture. In the context of preoperative planning and perioperative decision making, intracoronary optical coherence tomography imaging can provide a very detailed characterization of the arterial wall structure. However, visual interpretation of the images is laborious, subject to variability, and therefore not always sufficiently reliable for immediate decision of treatment.

METHODS

A novel semiautomatic segmentation method to quantify coronary fibrous cap thickness in optical coherence tomography is introduced. To cope with the most challenging issue when estimating cap thickness (namely the diffuse appearance of the anatomical abluminal interface to be detected), the proposed method is based on a robust dynamic programming framework using a geometrical a priori. To determine the optimal parameter settings, a training phase was conducted on 10 patients.

RESULTS

Validated on a dataset of 179 images from 21 patients, the present framework could successfully extract the fibrous cap contours. When assessing minimal cap thickness, segmentation results from the proposed method were in good agreement with the reference tracings performed by a medical expert (mean absolute error and standard deviation of 22 ± 18 μm) and were similar to inter-observer reproducibility (21 ± 19 μm, R = .74), while being significantly faster and fully reproducible.

CONCLUSION

The proposed framework demonstrated promising performances and could potentially be used for online identification of high-risk plaques.

Validation of renal artery dimensions measured by magnetic resonance angiography in patients referred for renal sympathetic denervation

RATIONALE AND OBJECTIVES

Magnetic resonance angiography (MRA) is a well-established modality for the assessment of renal artery stenosis. Using dedicated quantitative analyses, MRA can become a useful tool for assessing renal artery dimensions in patients referred for renal sympathetic denervation (RDN) and for providing accurate measurements of vascular response after RDN. The purpose of this study was to test the reproducibility of a novel MRA quantitative imaging tool and to validate these measurements against intravascular ultrasound (IVUS).

MATERIALS AND METHODS

In nine patients referred for renal denervation, renal artery dimensions were measured. Bland-Altman analysis was used to assess the intraobserver and interobserver reproducibility.

RESULTS

Mean lumen diameter was 5.8 ± 0.7 mm, with a very good intraobserver and interobserver variability of 0.7% (reproducibility: bias, 0 mm; standard deviation [SD], 0.1 mm) and 1.2% (bias, 0 mm; SD, 0.1 mm), respectively. Mean total lumen volume was 1035.3 ± 403.6 mm(3) with good intraobserver and interobserver variability of 2.9% (bias, -9.7 mm(3); SD, 34.0 mm(3)) and 2.8% (bias, -11.4 mm(3); SD, 42.4 mm(3)). The correlation (Pearson R) between mean lumen diameter measured with MRA and IVUS was 0.750 (P = .002).

CONCLUSIONS

Using a novel MRA quantitative imaging tool, renal artery dimensions can be measured with good reproducibility and accuracy. MRA-derived diameters and volumes correlated well with IVUS measurements.

Quantitative evaluation of noise reduction and vesselness filters for liver vessel segmentation on abdominal CTA images

Liver vessel segmentation in CTA images is a challenging task, especially in the case of noisy images. This paper investigates whether pre-filtering improves liver vessel segmentation in 3D CTA images. We introduce a quantitative evaluation of several well-known filters based on a proposed liver vessel segmentation method on CTA images. We compare the effect of different diffusion techniques i.e. Regularized Perona-Malik, Hybrid Diffusion with Continuous Switch and Vessel Enhancing Diffusion as well as the vesselness approaches proposed by Sato, Frangi and Erdt. Liver vessel segmentation of the pre-processed images is performed using a histogram-based region grown with local maxima as seed points. Quantitative measurements (sensitivity, specificity and accuracy) are determined based on manual landmarks inside and outside the vessels, followed by T-tests for statistic comparisons on 51 clinical CTA images. The evaluation demonstrates that all the filters make liver vessel segmentation have a significantly higher accuracy than without using a filter (p  <  0.05); Hybrid Diffusion with Continuous Switch achieves the best performance. Compared to the diffusion filters, vesselness filters have a greater sensitivity but less specificity. In addition, the proposed liver vessel segmentation method with pre-filtering is shown to perform robustly on a clinical dataset having a low contrast-to-noise of up to 3 (dB). The results indicate that the pre-filtering step significantly improves liver vessel segmentation on 3D CTA images.

Cardiac MR perfusion image processing techniques: a survey

First-pass cardiac MR perfusion (CMRP) imaging has undergone rapid technical advancements in recent years. Although the efficacy of CMRP imaging in the assessment of coronary artery diseases (CAD) has been proven, its clinical use is still limited. This limitation stems, in part, from manual interaction required to quantitatively analyze the large amount of data. This process is tedious, time-consuming, and prone to operator bias. Furthermore, acquisition and patient related image artifacts reduce the accuracy of quantitative perfusion assessment. With the advent of semi- and fully automatic image processing methods, not only the challenges posed by these artifacts have been overcome to a large extent, but a significant reduction has also been achieved in analysis time and operator bias. Despite an extensive literature on such image processing methods, to date, no survey has been performed to discuss this dynamic field. The purpose of this article is to provide an overview of the current state of the field with a categorical study, along with a future perspective on the clinical acceptance of image processing methods in the diagnosis of CAD.

Robust shape regression for supervised vessel segmentation and its application to coronary segmentation in CTA

This paper presents a vessel segmentation method which learns the geometry and appearance of vessels in medical images from annotated data and uses this knowledge to segment vessels in unseen images. Vessels are segmented in a coarse-to-fine fashion. First, the vessel boundaries are estimated with multivariate linear regression using image intensities sampled in a region of interest around an initialization curve. Subsequently, the position of the vessel boundary is refined with a robust nonlinear regression technique using intensity profiles sampled across the boundary of the rough segmentation and using information about plausible cross-sectional vessel shapes. The method was evaluated by quantitatively comparing segmentation results to manual annotations of 229 coronary arteries. On average the difference between the automatically obtained segmentations and manual contours was smaller than the inter-observer variability, which is an indicator that the method outperforms manual annotation. The method was also evaluated by using it for centerline refinement on 24 publicly available datasets of the Rotterdam Coronary Artery Evaluation Framework. Centerlines are extracted with an existing method and refined with the proposed method. This combination is currently ranked second out of 10 evaluated interactive centerline extraction methods. An additional qualitative expert evaluation in which 250 automatic segmentations were compared to manual segmentations showed that the automatically obtained contours were rated on average better than manual contours.

Reversible jump MCMC methods for fully automatic motion analysis in tagged MRI

Tagged magnetic resonance imaging (tMRI) is a well-known noninvasive method for studying regional heart dynamics. It offers great potential for quantitative analysis of a variety of kine(ma)tic parameters, but its clinical use has so far been limited, in part due to the lack of robustness and accuracy of existing tag tracking algorithms in dealing with low (and intrinsically time-varying) image quality. In this paper, we evaluate the performance of four frequently used concepts found in the literature (optical flow, harmonic phase (HARP) magnetic resonance imaging, active contour fitting, and non-rigid image registration) for cardiac motion analysis in 2D tMRI image sequences, using both synthetic image data (with ground truth) and real data from preclinical (small animal) and clinical (human) studies. In addition we propose a new probabilistic method for tag tracking that serves as a complementary step to existing methods. The new method is based on a Bayesian estimation framework, implemented by means of reversible jump Markov chain Monte Carlo (MCMC) methods, and combines information about the heart dynamics, the imaging process, and tag appearance. The experimental results demonstrate that the new method improves the performance of even the best of the four previous methods. Yielding higher consistency, accuracy, and intrinsic tag reliability assessment, the proposed method allows for improved analysis of cardiac motion.

Model-based approach for tracking embryogenesis in Caenorhabditis elegans fluorescence microscopy data

The nematode Caenorhabditis elegans (C. elegans) is a widely used model organism in biological investigations. Due to its well-known and invariant cell lineage tree, it can be used to study the effects of mutations and various disease processes. Effective and efficient analysis of the wealth of time-lapse fluorescence microscopy image data acquired in such studies requires automation of the cell segmentation and tracking tasks involved. This is hampered by many factors, including autofluorescence effects, low and uneven contrast throughout the images, high noise levels, large numbers of possibly simultaneous cell divisions, and touching or clustering cells. In this paper, we present a new algorithm for segmentation and tracking of cells in C. elegans embryogenesis image data. It is based on the model evolution framework for image segmentation and uses a novel multi-object tracking scheme based on energy minimization via graph cuts. Preliminary experiments on publicly available test data demonstrate the potential of the algorithm compared to existing approaches.

Quantitative comparison of spot detection methods in fluorescence microscopy

Quantitative analysis of biological image data generally involves the detection of many subresolution spots. Especially in live cell imaging, for which fluorescence microscopy is often used, the signal-to-noise ratio (SNR) can be extremely low, making automated spot detection a very challenging task. In the past, many methods have been proposed to perform this task, but a thorough quantitative evaluation and comparison of these methods is lacking in the literature. In this paper, we evaluate the performance of the most frequently used detection methods for this purpose. These include seven unsupervised and two supervised methods. We perform experiments on synthetic images of three different types, for which the ground truth was available, as well as on real image data sets acquired for two different biological studies, for which we obtained expert manual annotations to compare with. The results from both types of experiments suggest that for very low SNRs ( approximately 2), the supervised (machine learning) methods perform best overall. Of the unsupervised methods, the detectors based on the so-called h -dome transform from mathematical morphology or the multiscale variance-stabilizing transform perform comparably, and have the advantage that they do not require a cumbersome learning stage. At high SNRs ( > 5), the difference in performance of all considered detectors becomes negligible.

Conditional shape models for cardiac motion estimation

We propose a conditional statistical shape model to predict patient specific cardiac motion from the 3D end-diastolic CTA scan. The model is built from 4D CTA sequences by combining atlas based segmentation and 4D registration. Cardiac motion estimation is, for example, relevant in the dynamic alignment of pre-operative CTA data with intra-operative X-ray imaging. Due to a trend towards prospective electrocardiogram gating techniques, 4D imaging data, from which motion information could be extracted, is not commonly available. The prediction of motion from shape information is thus relevant for this purpose. Evaluation of the accuracy of the predicted motion was performed using CTA scans of 50 patients, showing an average accuracy of 1.1 mm.