Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Statistical, Morphometric, Anatomical Shape Model (Atlas) of Calcaneus

  • Aleksandra U. Melinska ,

    aleksandra.melinska@pwr.edu.pl

    Current address: Wybrzeze Wyspianskiego 27, 50-370 Wroclaw, Poland

    Affiliation Department of Biomedical Engineering, Wroclaw University of Technology, Wroclaw, Poland

  • Patryk Romaszkiewicz,

    Affiliation Regional Specialist Hospital, Research and Development Centre, Wroclaw, Poland

  • Justyna Wagel,

    Affiliation Department of General Radiology, Interventional Radiology and Neuroradiology, Wroclaw Medical University, Wroclaw, Poland

  • Marek Sasiadek,

    Affiliation Department of General Radiology, Interventional Radiology and Neuroradiology, Wroclaw Medical University, Wroclaw, Poland

  • D. Robert Iskander

    Affiliation Department of Biomedical Engineering, Wroclaw University of Technology, Wroclaw, Poland

Abstract

The aim was to develop a morphometric and anatomically accurate atlas (statistical shape model) of calcaneus. The model is based on 18 left foot and 18 right foot computed tomography studies of 28 male individuals aged from 17 to 62 years, with no known foot pathology. A procedure for automatic atlas included extraction and identification of common features, averaging feature position, obtaining mean geometry, mathematical shape description and variability analysis. Expert manual assistance was included for the model to fulfil the accuracy sought by medical professionals. The proposed for the first time statistical shape model of the calcaneus could be of value in many orthopaedic applications including providing support in diagnosing pathological lesions, pre-operative planning, classification and treatment of calcaneus fractures as well as for the development of future implant procedures.

Introduction

Improving the diagnosis and supporting the therapy is a fundamental task of medical imaging [1, 2]. Planning the surgery and choosing the method of surgical treatment often requires deep knowledge of the shape and morphological characteristics of an object at hand. In case of traumatology the detailed assessment of the morphological characteristics of bone is important for the choice of stabilization method and subsequent healing process. Approximately 10% of all fractures occur in foot bones [3]. The calcaneus (lat.os calcis) fracture is the most common fracture of tarsal bones, with an incidence of about 60% [4, 5]. Calcaneus fractures are estimated at 2% of all fractures of human bones and can have long-term consequences for patients in terms of comfort and disability [6, 7]. This largest tarsal bone plays stabilizing function for human body and pathologies of this bone influence human mobility and quality of life [810]. The evaluation of calcaneus fractures and determination of the effect of treatment essentially uses three types of imaging techniques: X-ray CT, and MRI [1115]. Although X-ray study (radiography) remains the method of choice for the initial assessment of fracture [16, 17], CT imaging has become the current state-of-art in diagnostics of calcaneus fractures [1820]. An invaluable advantage of CT is the possibility of reconstructing the scanned object. The location and configuration at a joint can be provided having 3D morphological and architectural information about the individual foot bones [21]. Several studies considered characterizing the calcaneal bone shape. Early works of Stindel et at. [21, 22] focused on classification of foot type in terms of pathological deformations using MRI. Stephan et al. [23] obtained information on the joint surfaces and 3D orientations of calcaneus using CT to quantify the integrity of calcaneus joint surfaces. Guterkunst et al. [24] aimed at assessing the measurement precision of landmark-based 3D bone-to-bone orientations of hind foot and lesser tarsal bones to compare an atlas-based automated method to the expert rating. Qiang et al. [7] described morphological characteristics of the calcaneus based on CT images. The survey in the later the paper has an epidemiological value. They described calcaneus by size and distances between characteristic points on the surface that were manually outlined by an experienced operator. Hence, there are possibilities of bone shape description from which a model of calcaneus could be constructed. This subsequently can be used for building statistical anatomical atlases (statistical shape models, SSMs) [25, 26]. Many tasks in medical imaging involve automatic systems which use prior knowledge about an object. Model based methods have become popular because they have the potential to resolve possible confusion associated with structural complexity, to improve tolerance for noisy data and deliver more accurate results in case of missing data [2732]. In medical imaging, SSMs are mainly used in segmentation and recognition tasks. An atlas guide/model-based segmentation approach is recognized as one of the most successful methods for image analysis [3336]. So far, applications of an atlas based model include brain structures [3741], soft tissue structures in the abdominal and pelvic area [32, 42], cardiac structures [29, 43, 44], and several bone structures [4547], but not of a calcaneus. Shape representation is one of the most important issues in SSM based methods of medical image processing. Several shape representation forms are known [48]. Zhang and Liu [49] classify shape representation and description techniques into region-based methods and contour-based methods. Region-based methods consist of global methods such as geometric moment invariants, algebraic moment invariants, orthogonal moments, generic Fourier descriptors, grid based methods, and shape matrix as well as structural methods such as convex hull and medial axis. On the other hand, contour-based techniques exploit shape boundary including global methods such as simple shape descriptors (area, eccentricity, major axis orientation, blending energy), correspondence-based shape matching in space domain, shape signature (centroidal profile, complex coordinates, centroid distance, tangent angle, cumulative angle, curvature, area), boundary moments, elastic matching, stochastic method (autoregressive modeling), scale space method, and spectral transform as well as structural methods such as chain code representation, polygon decomposition, smooth curve decomposition, scale space method, syntactic analysis, and shape invariants. Shape representation forms are not limited to the classification proposed by Zhang and Liu [49]. Cottes et al. [25, 27, 28, 50] classified them as “Hand Crafted” models (built from specific scalable objects like lines, ellipses and arcs) [30], articulated models [51], active contour models [52], Fourier-based shape models [53], statistical models of shape, and final elements models [54, 55]. Also, for the 3D shape representation, Heiman and Meizer [56] considered medial models (pioneered by Pizer et al. [57]) that allow rendering 3D solids and spherical harmonics models [5860].

The purpose of this work was to develop an anatomically accurate SSM of calcaneus which includes its morphological characteristics and provides mathematical representation of its shape.

Methods

Measurement setup

Retrospective data from regular hospital records have been used in this study. All patient records were anonymized and de-identified prior to processing according to the standard data release procedures of the two hospitals involved in the study. The Review Board of the Department of Radiology, Wroclaw Medical University, Wroclaw has approved the study. Poland Volume data have been acquired with three commercially available CT scanners (Siemens Somatom Definition Flash (2013), Forchheim, Germany; GE Medical Systems Discovery CT 750 HD (2010), Buckinghamshire, UK; and GE Medical Systems Light Speed CT VCT (2007), Buckinghamshire, UK). The tube voltage was 120 kVp in all three instruments. Data were retrospectively collected from the Department of Radiology at Regional Specialist Hospital, Research and Development Center and University Clinical Hospital in Wroclaw, Poland, saved in DICOM 3.0 format. All patients were scanned in feet-first, supine (FFS) position. Each volume image I, size [m, n, o] is made of different number of slices o, with each slice having size of 512×512 pixels. The parameters such as Pixel Spacing (0028, 0030)—xt, yt and Slice Thickness (0018, 0050)—zt were obtained from DICOM metadata [61] to extract physical dimensions of calcaneus. Volume data of 18 left and 18 right feet of 28 male subjects were collected. Subjects, aged from 17 to 62 years (mean 36.8 ± std 21.2 years), had no known foot pathology in the scanned foot. Subjects were aged matched, so they were no statistically significant differences between the age group averages of the left and the right foot subjects.

Data processing

A scheme for automatically building a morphometric and anatomical atlas was described by Subsol et al. [2]. It consists of the following stages: feature extraction, common feature identification, averaging feature position to obtain mean geometry and variability analysis. The process for building a morphometric and anatomical atlas of calcaneus is shown in Fig 1. In the following, we describe in detail the particular steps of the process. All procedures have been developed in Matlab (MathWorks, Inc., Natick, MA, USA) and are available from authors for free upon request.

thumbnail
Fig 1. The process of building the morphometric and anatomical atlas of calcaneus.

Hand icon denotes manual or semi-manual assisted processing lead by an experienced operator. Arrows icon denotes two isometric operations (translation and rotation) of the shape. Steps A to H are described in detail Section 2B.

https://doi.org/10.1371/journal.pone.0134603.g001

  1. Image pre-processing [Fig 1A and 1B] The volume image I is decomposed in the sagittal plane into a series of 2D images Ik, k = 1, 2, …, n. Each Ik is first normalized in two steps: (1) and then contrast is enhanced using sigmoid function [62]: (2) where g is a gain, and c is cut-off value, Ikmax is maximal and Ikmin is minimal pixel intensity value. Both parameters were set empirically to 0.4 and 0.7, respectively.
  2. Contour extraction [Fig 1C]
    Extracting the contour of calcaneus in each image Ik can be achieved in a variety of ways. For the purpose of this study we use the region growing algorithm [63, 64] with manually assisted assignment of a seed point by an expert operator.
  3. 3D point cloud to surface [Fig 1D]
    The result of the region growing algorithm is a set of coordinates W = [x, y, z], with x = [x1, x2, …, xn]′, y = [y1, y2, …, yn]′, and z = [z1, z2, …, zn]′ representing the contour of the calcaneus gathered respectively from each Ik, k = 1, 2, …, n. For such a point-cloud set the oriented normals are estimated and then using Poisson surface reconstruction method the surface mesh is generated [65]. Meshlab (Pisa, Italy) software was used to generate meshes [66]. Surface representation is needed to further land-marking of the characteristic points of the calcaneus.
  4. Land-marking [Fig 1E]
    Following Qiang et al. [7] the anatomic landmarks are manually marked on bone surface mesh (denoted by letters A to F in Fig 1E by an expert operator. They are: the highest point of the posterior articular facet (point A), the bottom of the posterior articular facet at the lateral surface intersecting the anterior process of the calcaneus (point B), the highest point of the calcaneocuboid joint (point C), the lowest point of the calcaneocuboid joint (point D), the most posterior point of the calcaneal tuberosity (point E), and the highest point on the superior edge of the calcaneal tuberosity (point F).
  5. Averaging feature position and orientation [Fig 1F]
    The objects of study are not scanned in the same position so their orientation and location need to be normalized for the statistical bone atlas to be built. For this, the point-cloud set W and the characteristic points C, D, E are used. First, the plane π1 that includes points C, D, E is obtained
    This is followed by calculating the angle α between the plane π1 and the plane π: z = 0. Then W is translated by the vector to set the most posterior point of the calcaneal tuberosity (point E) in the origin. Next, W is rotated about the x axis by the angle α.
    The rotation matrix: is applied and the angle β between x axis and vector is calculated. Finally, W is rotated about the x axis by the angle β using the rotation matrix:
  6. Spherical harmonics decomposition [Fig 1G] The mathematical shape description is achieved by using the spherical harmonics (SPHARM) [58, 6770]. Using linear-in-parameters least square procedure, data (cloud of extracted points) are first fitted with a best sphere: (3) and then translated to a new coordinate system centred at the estimated origin : (4)
    The center of the new origin is set as a center of the mass: (5)
    Then, transformation from the Cartesian to spherical coordinates is made: (6) (7) , where θd is the azimuth, ϕd is the elevation angle and Rd is the radius of sphere.
    Then the normalization to a unit sphere is performed: (8) , where Rs = max(Rd)
    The complete complex-valued spherical harmonics are solutions to the Laplace equation: (9) are given by: (10) where is the normalization factor defined as: (11)
    are the associated Legendre functions: (12)
    The real-valued spherical harmonics can be expressed as: (13)
    A spherical harmonics based decomposition of the (θd, ϕd, Rd/Rs) calcaneus data can be expressed as: (14) or alternatively: (15) where l = ⌈d⌉ -1 and m alternates between -l, l + 1, …, 0, l-1, l.
    The coefficient estimates in the spherical harmonic expansion , q = 1, 2, … Q can be easily evaluated using the method of linear least-squares by concatenating the data into a single column vector as described earlier. The estimate of the spherical harmonics based calcaneus surface is then (16)
    Transforming this result back to Cartesian coordinate system: (17) and then translating it back to the starting origin (see Fig 1G and 1H) creates an estimate of the calcaneus and allows to evaluate the quality of the approximation at the original points (xd, yd, zd), d = 1, 2, …, D.
  7. Model and model order selection [Fig 1H] The optimal model order of the SPHARM expansion can be evaluated using one of the standard information criteria such as Akaike Information Criterion or the Rissanen Minimum Description Length (MDL) criterion [71, 72]. For the purpose of this works the latter criterion was used. Given the optimal model order, further reduction of the number of coefficients in the SPHARM expansion is achieved by applying a multiple hypotheses procedure to test each of the coefficient values for zero.

Results

The MDL criterion estimated the optimal radial order of SPHARM expansion at 11 resulting in 144 characteristic coefficients . Fig 2 shows the box-plots for the first 25 SPHARM coefficients (corresponding to radial order of 4) for the left and right bones. The first, second and third quartile of were estimated and for them the model shape reconstructed. In Fig 3 we present reconstructed shape (as an example) for the right calcaneus. The mean values and the corresponding standard deviation (error bars) for the first 25 SPHARM coefficients are presented in Fig 4. Comparing the correspondent coefficients for the left and right bone we can see that some of them are similar and some of them have opposite sign. This is indicative of the symmetry of the objects. Fig 5 shows the correlation between mean values of SPHARM coefficients for left and right foot. When including the entire set of coefficients, the correlation of shapes is obviously very high (r2 = 0.99) and statistically significant (p ≪ 0.01). This correlation remains moderate when the three highest coefficients omitted amounting to r2 = 0.59, p ≪ 0.01 (see Fig 5). Application of a multiple hypotheses procedure to test each of the coefficient values for zero reduced the set of coefficient to 41 and 47 representative coefficients for the right and left calcaneus, respectively. Fig 6 shows an example of reconstructed SSM of right calcaneus using the whole set of 144 spherical harmonic coefficients and the reduced set of 41 representative coefficients.

thumbnail
Fig 2. The statistics of SPHARM estimates (box-plots) of calcaneus (the first 25 coefficients) for the group of 18 left (top) and 18 right calcanei (bottom).

Crosses indicate outliers.

https://doi.org/10.1371/journal.pone.0134603.g002

thumbnail
Fig 3. An example of reconstructed SSM of right calcaneus using (from left to right) the 25th, 50th and 75th quartile of the spherical harmonic coefficients.

https://doi.org/10.1371/journal.pone.0134603.g003

thumbnail
Fig 4. The mean values of SPHARM estimates of calcaneus (the first 25 coefficients) for the group of 18 right (dark gray) and 18 left (light grey) calcanei.

Error bars indicate one standard deviation.

https://doi.org/10.1371/journal.pone.0134603.g004

thumbnail
Fig 5. Correlation between the SPHARM coefficients of the left and the right calcaneus.

Linear correlation is performed on the reduced set of coefficients excluding the three highest coefficients (see inset).

https://doi.org/10.1371/journal.pone.0134603.g005

thumbnail
Fig 6. An example of reconstructed SSM of right calcaneus using the whole set of 144 spherical harmonic coefficients (left) and the reduced set of 41 representative coefficients (right).

https://doi.org/10.1371/journal.pone.0134603.g006

Table 1 shows the morphological characteristics of the calcaneus based on landmarks distances measurement as follows (see Methods): the length of the calcaneal axis (LCA): the distance from point G (the midpoint between C and D) to point E, the height of the posterior facet (HPF): the perpendicular distance from point A to the calcaneal axis, the length of the posterior facet (LPF): the distance between A and B, the length of the anterior process (LAP): the distance between B and C, the height of the anterior process (HAP): the distance between C and D, the critical angle of Gissane (GA): ∠ABC, the tuber angle of Böhler (BA): (180°-∠CAF). No statistically significant differences were found between the parameters of the left and right foot (t-test, p > 0.05).

thumbnail
Table 1. The group statistics (MEAN±SD) of the morphological characteristics of the calcaneus based on landmark distances.

https://doi.org/10.1371/journal.pone.0134603.t001

Discussion and Conclusions

Building an anatomically accurate statistical shape model of calcaneus requires close collaboration of image processing engineers, radiologists and orthopaedists. Following a general scheme for building such an atlas in an automatic function, as proposed by Subsol et al. [2], we were able to bring down the characterization of the mean geometry of calcaneus to several dozen representative coefficients of the spherical harmonic expansion. It should be noted that the scheme was not fully automatic and that the assistance of the experienced operator was necessary in two critical steps of building the model. That is, in the process of contour extraction where the selection of the seed point was manually assisted and in the difficult and tedious characterization of the representative anatomical features of the calcaneus. It is worth noting that the later was in a close agreement with the results of Qiang et al. [7] who also manually identified the anatomic landmarks of calcaneus. Shapes may be represented in many different ways and an argument can be made whether the choice of spherical harmonic representation was the most appropriate. Given the many advantages of spherical harmonics [69] such as their orthogonality and completeness, as well as the ability to estimate the coefficients using linear-in-parameters least-squares method provides a sufficient motivation for using this specific representation. It is important to acknowledge that other forms of 3D basis functions could also be used in place [73]. However, comparing different representations in terms of goodness-of-fit and the minimum number of representative coefficients, although interesting, was not within the scope of this work. One of the limitations of the study was the lack of complete CT data for both feet. It is expected that otherwise the symmetry of the model, as estimated via the two sets of spherical harmonic expansions, would be more evident. Another drawback was the acquisition of a single CT scanning per subjects (to limit the dose of radiation and cost), which prevented assessing the reproducibility of the measurement and modelling procedure. Nevertheless, other studies related to foot bone imaging with CT show relatively high test-retest reliability [24] that can be anticipated to hold for imaging calcaneus. Finally our study examined a set of male calcanei but the exact procedure can be readily applied for the female counterparts. To the best of our knowledge, the proposed statistical anatomical atlas of the calcaneus is presented for the first time. Tools of image processing, 3D shape modelling, and expert manual assistance were necessary to achieve the goal. The manual assistance in building a statistical shape model may appear to be a hindrance, but could be attempted, after verification, with robust shape registration methods [74, 75].

However, at the same time it provides an assurance of anatomical accuracy sought by medical professionals. In this view, the proposed statistical model of calcaneus including the quantitative mathematical shape description could be of value to diagnosing pathological changes by matching a model which contains information about the expected shape [76, 77]. It can be applied to monitoring growth process in childhood, bone and surrounding tissue age degeneration and osteoarthritis progression, as well as to detection of lesion (cysts, tumour-like changes). The model of calcaneus finds the application in pre-operative planning [78, 79], classifications [6] and follow-up of treatment of calcaneus fractures [8, 80], assessing the implant procedures [81, 82], and 3D reconstruction of bones [83]. Consequently, the model could also be used in image processing and computer vision software supporting those medical imaging tasks like segmentation or recognition.

As the summary demonstrates, there is already a considerable amount of 3D SSMs employed in medical image analysis. Extending this area to models of calcaneus could further support medicine in better diagnosis and treatment.

Supporting Information

S1 Data Set. The full data set for all 144 SPHARM coefficients (Fig 2).

https://doi.org/10.1371/journal.pone.0134603.s001

(XLS)

S2 Data Set. Data used to reconstruction of 3D shape for the right calcaneus (Fig 3).

https://doi.org/10.1371/journal.pone.0134603.s002

(XLS)

Acknowledgments

This research was supported by National Science Centre (NCN) Poland, Preludium UMO-2013/09/N/ST7/03418.

Author Contributions

Conceived and designed the experiments: AM DRI. Performed the experiments: AM PR JW. Analyzed the data: AM DRI. Contributed reagents/materials/analysis tools: PR JW. Wrote the paper: AM JW PR MS DRI.

References

  1. 1. Brief J, Hassfeld S, Wom H. Statistical analysis of the morphology of three-dimensional objects arid pathologic structures using spherical harmonics. Stud Health Technol Inform. 2001;81:103. pmid:11317720
  2. 2. Subsol G, Thirion JP, Ayache N. A scheme for automatically building three-dimensional morphometric anatomical atlases: application to a skull atlas. Med Image Anal. 1998;2(1):37–60. pmid:10638852
  3. 3. Brockwell J, Yeung Y, Griffith JF. Stress fractures of the foot and ankle. Sports Med Arthrosc. 2009;17(3):149–159. pmid:19680111
  4. 4. Freedman D, Zhang T. Interactive graph cut based segmentation with shape priors. In: CVPR. vol. 1. IEEE; 2005. p. 755–762.
  5. 5. Potter MQ, Nunley JA. Long-term functional outcomes after operative treatment for intra-articular fractures of the calcaneus. Bone Joint J. 2009;91(8):1854–1860.
  6. 6. Daftary A, Haims AH, Baumgaertner MR. Fractures of the Calcaneus: A Review with Emphasis on CT 1. Radiographics. 2005;25(5):1215–1226. pmid:16160107
  7. 7. Qiang M, Chen Y, Zhang K, Li H, Dai H. Measurement of three-dimensional morphological characteristics of the calcaneus using CT image post-processing. J Foot Ankle Res. 2014;7(1):19. pmid:24625107
  8. 8. Rammelt S, Amlang M, Barthel S, Gavlik JM, Zwipp H. Percutaneous treatment of less severe intraarticular calcaneal fractures. Clin Orthop Relat Res. 2010;468(4):983–990. pmid:19582524
  9. 9. Sanders R. Current Concepts Review-Displaced Intra-Articular Fractures of the Calcaneus*. Bone Joint J. 2000;82(2):225–50.
  10. 10. Stagni R, Leardini A, O’Connor JJ, Giannini S. Role of passive structures in the mobility and stability of the human subtalar joint: a literature review. Foot Ankle Int. 2003;24(5):402–409. pmid:12801196
  11. 11. Kang C, Speller R. Comparison of ultrasound and dual energy X-ray absorptiometry measurements in the calcaneus. Br J Radiol. 1998;71(848):861–867. pmid:9828799
  12. 12. Liang J, Pan AC, Huang YH, Fan XY. Fracture identification of X-ray image. In: ICWAPR. IEEE; 2010. p. 67–73.
  13. 13. Revol-Muller C, Benoit-Cattin H, Carillon Y, Odet C, Briguet A, Peyrin F. Bone MRI segmentation assessment based on synchrotron radiation computed microtomography. In: IEEE Nucl Sci Symp Conf Rec. vol. 3. IEEE; 2000. p. 19–7.
  14. 14. Saha PK, Wehrli FW. Measurement of trabecular bone thickness in the limited resolution regime of in vivo MRI by fuzzy distance transform. IEEE Trans Med Imaging. 2004;23(1):53–62. pmid:14719687
  15. 15. Udupa JK, Herman GT. 3D imaging in medicine. CRC press; 1999.
  16. 16. Bernstein J, Ahn J. In brief: Fractures in brief: calcaneal fractures. Clin Orthop Relat Res. 2010;468(12):3432–3434. pmid:20387018
  17. 17. Schepers T, Ginai A, Mulder P, Patka P. Radiographic evaluation of calcaneal fractures: to measure or not to measure. Skeletal Radiol. 2007;36(9):847–852. pmid:17554537
  18. 18. Bandyopadhyay O, Chanda B, Bhattacharya BB. Entropy-based automatic segmentation of bones in digital X-ray images. In: Pattern Rec Mach Intell. Springer; 2011. p. 122–129.
  19. 19. Cronier P, Pillet J, Talha A, Ghestem P, Liebault B, Rieux D, et al. Scanographic study of the calcaneus: normal anatomy and clinical applications. Surg Radiol Anat. 1988;10(4):303–310. pmid:3145572
  20. 20. Wu J, Tang Y, Davuluri P, Ward K, Belle A, Hobson R, et al. Fracture detection and quantitative measure of displacement in pelvic CT images. In: BIBMW. IEEE; 2011. p. 600–606.
  21. 21. Stindel E, Udupa JK, Hirsch BE, Odhner D. A characterization of the geometric architecture of the peritalar joint complex via MRI, an aid to classification of foot type. IEEE Trans Med Imaging. 1999;18(9):753–763. pmid:10571380
  22. 22. Stindel E, Udupa JK, Hirsch BE, Odhner D, Couture C. 3D MR image analysis of the morphology of the rear foot: application to classification of bones. Comput Med Imaging Graph. 1999;23(2):75–83. pmid:10227373
  23. 23. Stephan D, Panzer S, Göttlinger M, Augat P. Analysis of the intra-individual differences of the joint surfaces of the calcaneus. Comput Methods Biomech Biomed Engin. 2013;(ahead-of-print):1–7.
  24. 24. Gutekunst DJ, Liu L, Ju T, Prior FW, Sinacore DR. Reliability of clinically relevant 3D foot bone angles from quantitative computed tomography. J Foot Ankle Res. 2013;6(1):1–9.
  25. 25. Cootes TF, Taylor CJ. Statistical models of appearance for medical image analysis and computer vision. In: Proc. SPIE Medical Imaging. International Society for Optics and Photonics; 2001. p. 236–248.
  26. 26. Yushkevich PA, Zhang H, Gee JC. Statistical modeling of shape and appearance using the continuous medial representation. In: MICCAI. Springer; 2005. p. 725–732.
  27. 27. Cootes TF, Edwards GJ, Taylor CJ. Active appearance models. IEEE Trans Pattern Anal Mach Intell. 2001;23(6):681–685.
  28. 28. Cootes TF, Hill A, Taylor CJ, Haslam J. Use of active shape models for locating structures in medical images. Image Vis Comput. 1994;12(6):355–365.
  29. 29. Frangi AF, Rueckert D, Schnabel JA, Niessen WJ. Automatic construction of multiple-object three-dimensional statistical shape models: Application to cardiac modeling. IEEE Trans Med Imaging. 2002;21(9):1151–1166. pmid:12564883
  30. 30. Hill A, Thornham A, Taylor CJ. Model-Based Interpretation of 3D Medical Images. In: BMVC; 1993. p. 1–10.
  31. 31. Kelemen A, Székely G, Gerig G. Three-dimensional model-based segmentation of brain MRI. In: Proc. Workshop on Biomed Image Anal. IEEE; 1998. p. 4–13.
  32. 32. Lamecker H, Wenckebach TH, Hege HC. Atlas-based 3D-shape reconstruction from X-ray images. In: ICPR. vol. 1. IEEE; 2006. p. 371–374.
  33. 33. Becker M, Magnenat-Thalmann N. Deformable Models in Medical Image Segmentation. In: 3D Multiscale Physiological Human. Springer; 2014. p. 81–106.
  34. 34. Kelemen A, Szekely G, Gerig G. Elastic model-based segmentation of 3-D neuroradiological data sets. IEEE Trans Med Imaging. 1999;18(10):828–839. pmid:10628943
  35. 35. Park H, Bland PH, Meyer CR. Construction of an abdominal probabilistic atlas and its application in segmentation. IEEE Trans Med Imaging. 2003;22(4):483–492. pmid:12774894
  36. 36. Van Ginneken B, Frangi AF, Staal JJ, ter Haar Romeny BM, Viergever MA. Active shape model segmentation with optimal features. IEEE Trans Med Imaging. 2002;21(8):924–933. pmid:12472265
  37. 37. Guimond A, Meunier J, Thirion JP. Average brain models: A convergence study. Comput Vis Image Underst. 2000;77(2):192–210.
  38. 38. Leventon ME, Grimson WEL, Faugeras O. Statistical shape influence in geodesic active contours. In: CVPR. vol. 1. IEEE; 2000. p. 316–323.
  39. 39. Rueckert D, Frangi AF, Schnabel JA. Automatic construction of 3-D statistical deformation models of the brain using nonrigid registration. IEEE Trans Med Imaging. 2003;22(8):1014–1025. pmid:12906255
  40. 40. Styner M, Gerig G, Lieberman J, Jones D, Weinberger D. Statistical shape analysis of neuroanatomical structures based on medial models. Med Image Anal. 2003;7(3):207–220. pmid:12946464
  41. 41. Tao X, Prince JL, Davatzikos C. Using a statistical shape model to extract sulcal curves on the outer cortex of the human brain. IEEE Trans Med Imaging. 2002;21(5):513–524. pmid:12071622
  42. 42. Acosta O, Simon A, Monge F, Commandeur F, Bassirou C, Cazoulat G, et al. Evaluation of multi-atlas-based segmentation of CT scans in prostate cancer radiotherapy. In: Biom Imaging: From Nano to Macro. IEEE; 2011. p. 1966–1969.
  43. 43. Isgum I, Staring M, Rutten A, Prokop M, Viergever MA, van Ginneken B. Multi-atlas-based segmentation with local decision fusionApplication to cardiac and aortic segmentation in CT scans. IEEE Trans Med Imaging. 2009;28(7):1000–1010. pmid:19131298
  44. 44. Mitchell SC, Bosch JG, Lelieveldt BP, van der Geest RJ, Reiber JH, Sonka M. 3-D active appearance models: segmentation of cardiac MR and ultrasound images. IEEE Trans Med Imaging. 2002;21(9):1167–1178. pmid:12564884
  45. 45. Lefaix G, Riot X, Haigron P, Collorec R, Ramee A. 3D modeling and deformation analysis of the vertebra with spherical harmonics. In: EMBS. vol. 1. IEEE; 1997. p. 422–425.
  46. 46. Lu YC, Untaroiu CD. Statistical shape analysis of clavicular cortical bone with applications to the development of mean and boundary shape models. Comp Meth Prog Bio. 2013;111(3):613–628.
  47. 47. van de Giessen M, Foumani M, Streekstra GJ, Strackee SD, Maas M, van Vliet LJ, et al. Statistical descriptions of scaphoid and lunate bone shapes. J Biomech. 2010;43(8):1463–1469. pmid:20185138
  48. 48. Brechbühler C, Gerig G, Kübler O. Parametrization of closed surfaces for 3-D shape description. Computer vision and image understanding. 1995;61(2):154–170.
  49. 49. Zhang D, Lu G. Review of shape representation and description techniques. Pattern Recogn. 2004;37(1):1–19.
  50. 50. Cootes TF, Taylor CJ, Cooper DH, Graham J. Active shape models-their training and application. Comput Vis Image Underst. 1995;61(1):38–59.
  51. 51. Beinglass A, Wolfson HJ. Articulated object recognition, or: How to generalize the generalized Hough transform. In: Proc. CVPR. IEEE; 1991. p. 461–466.
  52. 52. Kass M, Witkin A, Terzopoulos D. Snakes: Active contour models. Int J Comput Vision. 1988;1(4):321–331.
  53. 53. Scott GL. The Alternative Snake-and Other Animals. In: Alvey Vision Conf. Citeseer; 1987. p. 1–8.
  54. 54. Cohen LD, Cohen I. Finite-element methods for active contour models and balloons for 2-D and 3-D images. IEEE Trans Pattern Anal Mach Intell. 1993;15(11):1131–1147.
  55. 55. Karaolani P, Sullivan GD, Baker KD, Baines M. A Finite Element Method for Deformable Models. In: Alvey Vision Conf. Citeseer; 1989. p. 1–6.
  56. 56. Heimann T, Meinzer HP. Statistical shape models for 3D medical image segmentation: A review. Med Image Anal. 2009;13(4):543–563. pmid:19525140
  57. 57. Pizer SM, Fletcher PT, Joshi S, Thall A, Chen JZ, Fridman Y, et al. Deformable m-reps for 3D medical image segmentation. Int J Comput Vision. 2003;55(2–3):85–106.
  58. 58. Iskander DR. Modeling videokeratoscopic height data with spherical harmonics. Optom Vis Sci. 2009 May;86(5):542–547. JID: 8904931; ppublish. pmid:19319007
  59. 59. Nikou C, Bueno S, Heitz F, Armspach JP. A joint physics-based statistical deformable model for multimodal brain image analysis. IEEE Trans Med Imaging. 2001;20(10):1026–1037. pmid:11686438
  60. 60. Styner M, Oguz I, Xu S, Brechbühler C, Pantazis D, Levitt JJ, et al. Framework for the statistical shape analysis of brain structures using SPHARM-PDM. The insight journal. 2006;(1071):242. pmid:21941375
  61. 61. The DICOM Standard 2015a;. Available from: http://dicom.nema.org/standard.html.
  62. 62. Tumblin J, Hodgins JK, Guenter BK. Two methods for display of high contrast images. ACM Trans Graph. 1999;18(1):56–94.
  63. 63. Chang YL, Li X. Adaptive image region-growing. IEEE Trans Image Process. 1994;3(6):868–872. pmid:18296257
  64. 64. Pham DL, Xu C, Prince JL. Current methods in medical image segmentation 1. Annu Rev Biomed Eng. 2000;2(1):315–337. pmid:11701515
  65. 65. Kazhdan M, Bolitho M, Hoppe H. Poisson surface reconstruction. In: Proc. ACM SPG. vol. 7; 2006.
  66. 66. CNR VCLI. MeshLab;. Http://meshlab.sourceforge.net/.
  67. 67. Bookstein FL. Shape and the information in medical images: A decade of the morphometric synthesis. In: Proc MMBIA. IEEE; 1996. p. 2–12.
  68. 68. Gerig G, Styner M, Jones D, Weinberger D, Lieberman J. Shape analysis of brain ventricles using spharm. In: MMBIA. IEEE; 2001. p. 171–178.
  69. 69. Matheny A, Goldgof DB. The use of three-and four-dimensional surface harmonics for rigid and nonrigid shape recovery and representation. IEEE Trans Pattern Anal Mach Intell. 1995;17(10):967–981.
  70. 70. Shen L, Farid H, McPeek MA. Modeling three-dimensional morphological structures using spherical harmonics. Evolution. 2009;63(4):1003–1016. pmid:19154365
  71. 71. Rissanen J. Modeling by shortest data description. Automatica. 1978;14(5):465–471.
  72. 72. Davies RH, Twining CJ, Cootes TF, Waterton JC, Taylor CJ. A minimum description length approach to statistical shape modeling. Medical Imaging, IEEE Transactions on. 2002;21(5):525–537.
  73. 73. Tangelder JW, Veltkamp RC. A survey of content based 3D shape retrieval methods. Multimed Tools Appl. 2008;39(3):441–471.
  74. 74. Trouvé A. Diffeomorphisms Groups and Pattern Matching in Image Analysis. International Journal of Computer Vision. 1998;28(3):213–221.
  75. 75. Durrleman S, Prastawa M, Charon N, Korenberg JR, Joshi S, Gerig G, et al. Morphometry of anatomical shape complexes with dense deformations and sparse parameters. NeuroImage. 2014;101:35–49. pmid:24973601
  76. 76. Ebraheim NA, Zeiss J, Skie MC, Jackson WT. Radiological evaluation of peroneal tendon pathology associated with calcaneal fractures. J Orthop Trauma. 1991;5(3):365–369. pmid:1941321
  77. 77. Guyer B, Levinsohn E, Fredrickson B, Bailey G, Formikell M. Computed tomography of calcaneal fractures: anatomy, pathology, dosimetry, and clinical relevance. AJR Am J Roentgenol. 1985;145(5):911–919. pmid:3876749
  78. 78. Chen Y, Qiang M, Zhang K, Li H, Dai H. Novel computer-assisted preoperative planning system for humeral shaft fractures: report of 43 cases. Int J Med Robotics Comput Assist Surg. 2014;.
  79. 79. Chen Y, Zhang K, Qiang M, Li H, Dai H. Computer-assisted preoperative planning for proximal humeral fractures by minimally invasive plate osteosynthesis. Chin Med J (Engl). 2014;127(18):3278–3285.
  80. 80. Rammelt S. An Update on the Treatment of Calcaneal Fractures. J Orthop Trauma. 2014;28(10):549–550. pmid:25243848
  81. 81. Cao L, Weng W, Song S, Mao N, Li H, Cai Y, et al. Surgical Treatment of Calcaneal Fractures of Sanders Type II and III by a Minimally Invasive Technique Using a Locking Plate. J Foot Ankle Surg. 2015;54(1):76–81. pmid:25441282
  82. 82. Görres J, Brehler M, Franke J, Barth K, Vetter SY, Córdova A, et al. Intraoperative detection and localization of cylindrical implants in cone-beam CT image data. Int J Comput Assist Radiol Surg. 2014;9(6):1045–1057.
  83. 83. Cash C, Treece G, Berman L, Gee A, Prager R. 3D reconstruction of the skeletal anatomy of the normal neonatal foot using 3D ultrasound. 2005;931(78):587–595.