The watershed transform has interesting properties that make it useful for many different image segmentation applications: it is simple and intuitive, can be parallelized, and always produces a complete division of the image. However, when applied to medical image analysis, it has important drawbacks (oversegmentation, sensitivity to noise, poor detection of thin or low signal to noise ratio structures). We present an improvement to the watershed transform that enables the introduction of prior information in its calculation. We propose to introduce this information via the use of a previous probability calculation. Furthermore, we introduce a method to combine the watershed transform and atlas registration, through the use of markers. We have applied our new algorithm to two challenging applications: knee cartilage and gray matter/white matter segmentation in MR images. Numerical validation of the results is provided, demonstrating the strength of the algorithm for medical image segmentation.

Purpose To characterize the trabeculated connective tissue microarchitecture of the lamina cribrosa (LC) in terms of total connective tissue volume (CTV), connective tissue volume fraction (CTVF), predominant beam orientation, and material anisotropy in monkeys with early experimental glaucoma (EG). Methods The optic nerve heads from three monkeys with unilateral EG and four bilaterally normal monkeys were three dimensionally reconstructed from tissues perfusion fixed at an intraocular pressure of 10 mm Hg. A three-dimensional segmentation algorithm was used to extract a binary, voxel-based representation of the porous LC connective tissue microstructure that was regionalized into 45 subvolumes, and the following quantities were calculated: total CTV within the LC, mean and regional CTVF, regional predominant beam orientation, and mean and regional material anisotropy. Results Regional variation within the laminar microstructure was considerable within the normal eyes of all monkeys. The laminar connective tissue was generally most dense in the central and superior regions for the paired normal eyes, and laminar beams were radially oriented at the periphery for all eyes considered. CTV increased substantially in EG eyes compared with contralateral normal eyes (82%, 44%, 45% increases; P < 0.05), but average CTVF changed little (−7%, 1%, and −2% in the EG eyes). There were more laminar beams through the thickness of the LC in the EG eyes than in the normal controls (46%, 18%, 17% increases). Conclusions The substantial increase in laminar CTV with little change in CTVF suggests that significant alterations in connective and nonconnective tissue components in the laminar region occur in the early stages of glaucomatous damage.

Recent advances in magnetic resonance (MR) imaging technology have unveiled a wealth of information regarding cardiac histoanatomical complexity. However, methods to faithfully translate this level of fine-scale structural detail into computational whole ventricular models are still in their infancy, and, thus, the relevance of this additional complexity for simulations of cardiac function has yet to be elucidated. Here, we describe the development of a highly detailed finite-element computational model (resolution: ∼125 μm) of rabbit ventricles constructed from high-resolution MR data (raw data resolution: 43 × 43 × 36 μm), including the processes of segmentation (using a combination of level-set approaches), identification of relevant anatomical features, mesh generation, and myocyte orientation representation (using a rule-based approach). Full access is provided to the completed model and MR data. Simulation results were compared with those from a simplified model built from the same images but excluding finer anatomical features (vessels/endocardial structures). Initial simulations showed that the presence of trabeculations can provide shortcut paths for excitation, causing regional differences in activation after pacing between models. Endocardial structures gave rise to small-scale virtual electrodes upon the application of external field stimulation, which appeared to protect parts of the endocardium in the complex model from strong polarizations, whereas intramural virtual electrodes caused by blood vessels and extracellular cleft spaces appeared to reduce polarization of the epicardium. Postshock, these differences resulted in the genesis of new excitation wavefronts that were not observed in more simplified models. Furthermore, global differences in the stimulus recovery rates of apex/base regions were observed, causing differences in the ensuing arrhythmogenic episodes. In conclusion, structurally simplified models are well suited for a large range of cardiac modeling applications. However, important differences are seen when behavior at microscales is relevant, particularly when examining the effects of external electrical stimulation on tissue electrophysiology and arrhythmia induction. This highlights the utility of histoanatomically detailed models for investigations of cardiac function, in particular for future patient-specific modeling.