Deformable triangular surfaces using fast 1-D radial Lagrangian dynamics-segmentation of 3-D MR and CT images of the wrist

• Published in: IEEE transactions on medical imaging 2002-08, Vol.21 (8), p.888-903
• Main Authors: Snel, J.G., Venema, H.W., Grimbergen, C.A.
• Format: Article
• Language: English
• Subjects: Algorithms; Biological and medical sciences; Bone and Bones - anatomy & histology; Bone and Bones - diagnostic imaging; Bones; Computed tomography; Deformable models; Elasticity; Filters; Humans; Image Enhancement - methods; Image segmentation; Imaging, Three-Dimensional - methods; Lagrangian functions; Magnetic resonance; Magnetic Resonance Imaging - methods; Medical sciences; Models, Biological; Motion; Reproducibility of Results; Robustness; Sensitivity and Specificity; Solid modeling; Stochastic Processes; Stress, Mechanical; Studies; Surface Properties; Tomography, X-Ray Computed - methods; Wrist; Wrist Joint - anatomy & histology; Wrist Joint - diagnostic imaging
• ISSN: 0278-0062; 1558-254X
• DOI: 10.1109/TMI.2002.803127

We developed a new triangulated deformable surface model, which is used to detect the boundary of the bones in three-dimensional magnetic resonance (MR) and computed tomography (CT) images of the wrist. This surface model is robust to initialization and provides wide geometrical coverage and quantitative power. The surface is deformed by applying one-dimensional (1-D) radial Lagrangian dynamics. For initialization a tetrahedron is placed within the bone to be segmented. This initial surface is inflated to a binary approximation of the boundary. During inflation, the surface is refined by the addition of vertices. After the surface is fully inflated, a detailed, accurate boundary detection is obtained by the application of radial scale-space relaxation. In this optimization stage, the image intensity is filtered with a series of 1-D second-order Gaussian filters. The resolution of the triangulated mesh is adapted to the width of the Gaussian filter. To maintain the coherence between the vertices, a resampling technique is applied which is based on collapsing and splitting of edges. We regularized the triangulated mesh by a combination of volume-preserving vertex averaging and equi-angulation of edges. In this paper, we present both qualitative and quantitative results of the surface segmentations in eight MR and ten CT images.