 Methodology
 Open Access
 Published:
Landmark precision and reliability and accuracy of linear distances estimated by using 3D computed microtomography and the opensource TINA Manual Landmarking Tool software
Frontiers in Zoologyvolume 12, Article number: 12 (2015)
Abstract
Introduction
The wider availability of nondestructive and highresolution methods, such as microcomputed tomography (microCT), has prompted its use in anatomical and morphometric studies. Yet, because of the actual scanning procedure and the processing of CT data by software that renders 3D surfaces or volumes, systematic errors might be introduced in placing landmarks as well as in estimating linear distances. Here we assess landmark precision and measurement reliability and accuracy of using microCT images of toad skulls and the TINA Manual Landmarking Tool software to place 20 landmarks and extract 24 linear distances. Landmark precision and linear distances calculated from 3D images were compared to the same landmarks and distances obtained with a 3D digitizer in the same skulls. We also compared landmarks and linear distances in 3D images of the same individuals scanned with distinct filters, since we detected variation in bone thickness or density among the individuals used.
Results
We show that landmark precision is higher for microCT than for the 3D digitizer. Distance reliability was very high withinmethods, but decreased in 20 % when 3D digitizer and microCT data were joined together. Still, we did not find any systematic bias in estimating linear distances with the microCT data and the betweenmethods errors were similar for all distances (around 0.25 mm). Absolute errors correspond to about 6.5 % of the distance’s means for microCT resolutions and 3D digitizer comparisons, and to 3 % for the filter type analysis.
Conclusions
We conclude that using microCT data for morphometric analysis results in acceptable landmark precision and similar estimates of most linear distances compared to 3D digitizer, although some distances are more prone to discrepancies betweenmethods. Yet, caution in relation to the scale of the measurements needs to be taken, since the proportional betweenmethod error is higher for smaller distances. Scanning with distinct filters does not introduce a high level of error and is recommended when individuals differ in bone density.
Introduction
The increasing use of nondestructive and highresolution data acquisition methods, such as microcomputed tomography (microCT), have provided researchers with the opportunity to study the anatomy and morphology of organisms with more detail and at a wider phylogenetic spectrum (e.g. [1–5]). Accordingly, 3D image processing software has been developed (e.g. OsiriX: [6], Amira: www.amira.com), with some designed to place 3D landmarks for shape or morphometric analysis (e.g. TINA Manual Landmarking Tool: [7]). However, there is no guarantee that the scanning procedure and the software used to process CT data, by creating 3D surface or 3D volume renderings, do not introduce systematic errors in the data [8, 9]. In addition, the landmark positioning process in the 3D images might also introduce systematic and random errors in the estimation of linear distances (measurements). Thus, the precision of placing landmarks in 3D images with software, as well as the accuracy and reliability of the distances taken by the use of CT data must be validated [8–10].
In this study, we evaluated the precision of placing the same landmarks in the same individuals of a toad species with two distinct methods: (1) the real skulls and a 3D digitizer to place the landmarks, and (2) 3D images of these skulls, obtained by microCT scans at two distinct resolutions, and software to place the landmarks (TINA Manual Landmarking Tool ([7]; hereafter called TINALandmark). The 3D digitizer is an articulated arm that creates a 3D coordinate system in which any point of an object can be identified in relation to a reference point. TINALandmark is recently developed opensource software created to enhance the precision of landmark positioning in 3D images by using volume rendering instead of surface rendering and by showing the crosssection images connected to the 3D volume [7]).We chose to compare the landmarks obtained from the 3D images with the ones obtained by the 3D digitizer because, for the latter, the landmarks are taken in the actual skulls, with no processing of data, not the case for the construction of 3D images. Also, in zoological studies, several authors measure the specimens with 3D digitizers (e.g. [11–14], just to cite a few), being a widely accepted technique in the morphometry field.
In addition to comparing the landmarks between the two methods, it is also of interest to evaluate the consequence of potential biases in the estimation of linear distances extracted from the landmarks. Although measurement error is intrinsic to the landmarks and independent of the linear distances computed from them, the proportion of error varies with the distance length. That is, if the error in placing landmarks is the same for all landmarks, the proportion of error will be higher for smaller distances than for longer distances [15]. On the other hand, if the error in placing landmarks varies with landmark type or position in the material, linear distances extracted from these more variable landmarks are expected to have greater error when comparing distinct methods, although an association between distance length and error proportion will still exist. Therefore, we also calculated reliability and accuracy of linear distances extracted from the landmarks placed with the 3D digitizer and in the 3D images.
Finally, we also compared the same toad specimens scanned with two distinct filters, since we discovered that some individuals differ in bone thickness and/or density. Filters are thin sheets of metal set in front of the material being scanned and can have different thickness. Varying the filter thickness has an effect on the mean Xray energy irradiating on the material being scanned. For thinner bones, a lower Xray energy is necessary to achieve the best 3D volumes and more precise positioning of the landmarks (see Fig. 1). Thus, changing filter type might be an additional source of error when placing landmarks and taking linear measurements in 3D images. We consider this last comparison quite relevant in zoological studies because other organisms might present the same variation in bone density and scanning with distinct filters will be indispensable.
Although there is no specific theory relating the potential effects on landmarks and on distances when using distinct resolutions, filters and reconstruction algorithms, some expectations based on landmark position and bone thickness can be made. Some of the landmarks that we selected in the toad skulls (Table 1 and Fig. 2) were more difficult to visualize in the 3D images than in the real skulls because of their position (landmarks 4, 10, 11, 17 and 20) or because they were placed at thinner bones (landmarks 6, 7, 8 and 13). Thus, we expect more variation in the positioning of these landmarks among methods, and as a consequence, less reliability and accuracy of the linear distances extracted from them (Table 2 and Fig. 2). We report that placing landmarks in 3D images obtained by microCT scanning is more precise than placing the same landmarks with the 3D digitizer. Yet, average differences in linear distances among methods are acceptable and represent a low error proportion in relation to the distances lengths. Scanning with distinct resolutions and distinct filters do not introduce high errors.
Results and discussion
Landmark precision with distinct methods
We compared landmark precision betweenmethods by calculating the mean distance of each landmark in the individuals from the same landmarks of a mean shape estimated with a slightly modified General Procrustes Analysis, for all methods (see Methods section). Mean individual landmark distances from the mean sample landmarks were higher when using the 3D digitizer (DIG) to place the landmarks than when using the 3D images, regardless of the resolutions (medium [MED] or high [HIGH]). This result holds for the shape space (without scale, Table A1 in Additional file 1), as well as for retaining landmark scale information in mm (after multiplying the Procrustes configuration by mean centroid size, Table 3), indicating that placing landmarks in 3D images obtained by the microCT is more precise. Even though the precisions of both equipments used in this study are similar (3D digitizer: 0.01 mm; microCT resolutions MED: 0.018 mm and HIGH: 0.009 mm), this result may be due to the fact that the 3D images of the skulls are much bigger in the screen than the real ones (around tenfold increase), facilitating visualization of several bone structures. The range of deviations of individual landmarks from the mean shape landmarks for all methods (DIG: 0.27 to 0.52 mm; HIGH: 0.14 to 0.33 mm and MED: 0.16 to 0.34 mm) are comparable to the error found by Richtsmeier et al. (1995) [16] when placing the same landmarks in two different 3D images of the same individual (scanned twice with a 1.5 mm slice thickness resolution), which was 0.15 to 0.48 mm. Thus, the precision of both methods in placing landmarks is acceptable for morphometric studies, at least for distances as large as the ones obtained in this study.
The two landmarks with highest mean deviations for DIG (14: posterior tip of parasphenoid process and 19: posterior tip of pterygoid process) were not expected to have higher error, however they both are type II landmarks, suggesting that locating the exact maximum bone curvature was more difficult than locating bone sutures (type I landmarks). The higher error for type II landmarks compared to type I landmarks is expected when working with the real skulls [17]. Landmarks 3, 5, 6, 7, 9, 10, 11, 12, 15, 16 and 20 presented higher discrepancies in the landmarks for DIG when compared to the microCT resolutions (boldface in Table 3). For the landmarks expected to have higher error in the 3D images because of their positions in the skull or because they were located in thinner bones, only landmarks 4, 8 and 12 corresponded to the expectation. All of them are type I landmarks, agreeing with the findings of other authors working with 3D images of more error in placing landmarks at some particular bone sutures [18, 19].
MED and HIGH presented very similar results, indicating that there is no difference in switching the resolutions to visualize the landmarks. The fact that we had to downsample the HIGH data by a factor of two to load the crosssection sequences in TINALandmark probably did not interfere in landmark precision, since switching the compacting factors with the MED data also produced similar results (comparisons between stride = 2 [S2] and stride = 2 and downsample = 2 [S2/D2]; and stride = 3 [S3] and stride = 3 and downsample = 2 [S3/D2]; Table 3). This result is consistent with the fact that the crosssections are maintained in the original resolution even when downsampling the data (which affects only the resolution of the 3D volume), being possible to refine the landmarks positions in them. Finally, for the filter type analysis, landmarks 6, 8 and 13 presented higher mean individual landmark distances for aluminum 0.5 mm (AL0.5) than for aluminum 1.0 mm (AL1.0), coinciding with the ones expected as more variable for being located at thinner bones. The only exception was landmark 5, with high deviation for AL 0.5. This landmark is located in the nasal and frontoparietal lateral suture, coinciding with a neural crest. It is possible that the presence of the crest caused higher error in the positioning of the landmark.
Linear distances reliability and accuracy
The reliability of the linear distances was obtained by using distance repeatability, which indicates the reliability of multiple measurements in the same individuals. It describes the proportion of variance due to differences among individuals (betweengroup variance) in relation to the residual variance plus the betweengroup variance [20]. Withinmethod mean distance repeatabilities (i.e. considering only replicates measured with the same method) were very high, above 0.9 for all methods, as can be seen in Table 4 for microCT resolutions and 3D digitizer; and in Table 5 for image compacting factors and distinct filters. This result shows that measuring procedure is very reliable inside each method. However, betweenmethods mean distance repeatability is considerably lower (20 % lower) when replicates of MED or HIGH were joined with replicates of DIG, but not when replicates of distinct microCT resolutions were joined together (Table 4). When replicates of AL0.5 and AL1.0 were joined together, the mean repeatability also reduced (Table 5), yet the drop was less steep (13 % lower) than MED + DIG and HIGH + DIG. For the image compacting factors there was no drop in mean distance repeatability when replicates of different values of the compacting factors were considered together (Table 5). Looking at repeatabilities separately for each distance in the betweenmethods analysis (Fig. 3b,c), it can be noticed that distances 3, 5, 6, 9, 10, 11, 13, 14, 16, 23 and 24 all have low repeatabilities (below 0.8) and are composed of landmarks detected with higher mean deviations for DIG (landmarks 3, 5, 8, 9, 10, 11, 12 and 19). These distances comprise small as well as large distances, showing that the error is independent of distance length and related to the quality of the landmark. Although we did not make expectations of higher error for landmarks placed with DIG, the landmarks with higher deviations had as a consequence higher discrepancies in the betweenmethods repeatabilities. However, the distances 17 and 22, composed of the landmarks with the highest error for DIG (14 and 19), did not present low repeatabilities. This is a consequence of their longer mean lengths (17: 7.14 mm and 22: 5.06 mm) and a smaller proportional error as a consequence. On the other hand, distances 5, 10, 11, 13 and 24 are small distances (less than 3.0 mm of length) composed of highly discrepant landmarks, and therefore, have a larger proportional error. These differences in proportional error depending on the mean distance length is shown by a significant positive association between mean distance repeatability and mean distance length for the joined data sets MED + DIG and HIGH + DIG (r = 0.45, d.f. = 22, P = 0.02 and r = 0.44, d.f. = 22, P = 0.03, respectively; Fig. 3b,c), as well as for S2 + S2/D2 and AL0.5 + AL1.0 (r = 0.44, d.f. = 24, P = 0.03 and r = 0.52, d.f. = 24, P = 0.007, respectively; Fig. 3e,i). Although a significant relation between repeatability and distance length exists, practically all distances composed of landmarks with higher differences betweenmethods presented lower repeatability when compared to withinmethod repeatability.
Additionally to the distance repeatability analysis, we also compared mean raw and absolute differences in the distances within and between methods. The mean raw differences within and between methods were close to zero for all data sets, being sometimes positive and sometimes negative (Tables 4 and 5), indicating that in average there is no consistent bias when measuring the same individual two times or measuring the same individuals with different methods. Yet, when looking at the raw differences between MEDDIG and HIGHDIG for the distances separately, we can see that some distances have higher differences, ranging from −0.38 mm to 0.36 mm (Table A2 in the Additional file 1). These distances correspond to different bones in the toad skulls, indicating that the deviations in the relative landmark positions are not localized in a few bones. We can see that the higher differences (above 0.17 mm in magnitude, the highest difference found for MEDHIGH, values in boldface in Table A2) have approximately the same magnitude in mm for small and long distances. For instance, small distances 5, 8 and 13 (distance means below 3.0 mm) have a mean betweenmethods difference of 0.23 mm, while long distances have a mean of 0.26 mm. However, as mentioned above, the same magnitude of betweenmethods error for small and long distances results in a higher proportion of error in relation to the distances means when considering the smaller ones (Table A3 in Additional file 1). When looking at the betweenmethods mean percentage error, the error in relation to the distances mean lengths, we can see that the cases around 10 % correspond to three small distances: 5 (frontoparietal bone), 13 (nasal bone) and 25 (parasphenoid bone), and just one longer distance: 16 (squamosal bone), which is composed of two landmarks detected with high deviations for DIG (landmarks 8 and 12). The magnitudes of differences that we found for DIG and the microCT resolutions are similar to differences reported by other authors when comparing measurements taken with CT data and digital calipers [18, 21], although the last authors did find a systematic bias in the CT data (all distances underestimated). Yet, when comparing our results with other authors that measured the same specimens with a 3D digitizer (Polhemus 3Space) and CT [15, 19], our error betweenmethods is much smaller (around five to ten times smaller). This is probably due to the fact that we have used a much higher resolution than these authors in both the microCT as well as the digitizer since current equipments are an order of magnitude more accurate than 20 years ago.
The only distances estimated with more error when scanning at different resolutions were distances17 and 19 (0.17 mm of difference), both distances in the parasphenoid bone, one of the thinner bones in the toad skull. It is possible that in this case scanning with HIGH enhanced the visualization of the sutures in the parasphenoid bone compared to MED. Scanning with HIGH instead of with MED increases the scanning time per skull over 4 times (20 min versus 1.5 h), so the advantage in enhancing the accuracy of some distances by using higher resolution depends on the question being asked and the time available. The betweenmethods absolute differences for each distance (Fig. 4) also indicate that some distances had higher magnitudes of difference. The worst case was for distance 16 (squamosal bone), the only long distance that presented percentage error around 10 % and an error around 0.5 mm. This magnitude of error is very close to what we considered a gross error when measuring the individuals inside each method. When looking at the distinct filter results (AL1.0AL0.5), we can see that both the raw and absolute differences between individuals scanned with distinct filters are much lower than differences between microCT resolutions and 3D digitizer. Similarly, in table A3 we can notice that the highest mean percentage error for filter type comparisons was around 7 %, and all related to small distances.
Conclusions
In order to use 3D images of organisms in morphometric studies, one needs to know if the relative positions among landmarks are kept the same for different scanning procedures or rendering algorithms. By comparing landmark precision and linear distances in toad skulls scanned with two resolutions and measured with the software TINALandmark, with the same skulls measured with a 3D digitizer, we conclude that the degree of discrepancy is acceptable in general, although several distances had betweenmethods discrepancies above 5 %. Yet, we must consider that 30 % of the distances are below 3.0 mm in length, being quite small distances, thus for several other bigger organisms the error proportion will probably be lower than what we report. Attention needs to be paid in relation to the scale of the distances, as smaller distances might be estimated with proportionally more error. Using distinct microCT resolutions, distinct compacting factors for loading 2D crosssection sequences in TINALandmark, as well as using distinct filters to scan individuals that differ in bone density do not introduce high errors. We recommend the switching of filters when individuals to be compared differ in bone density, at least with the resolutions that we have tested.
Methods
Species and scanning procedures
For the microCT and 3D digitizer comparisons, as well as for the 3D image compacting factor analysis, we used 20 adult individuals of the toad species Rhinella granulosa, collected in January 2012 at a site near the municipality of Angicos (5° 39′ S, 36° 36′W), in the state of Rio Grande do Norte, NE Brazil. The toads were sacrificed by peritoneal application of an anesthetic in excess and then were fixed in 70 % alcohol. We scanned the toad skulls with an Xray microCT system (SkyScan 1176, Konitch, Belgium) placed at the Instituto de Biociências, Universidade de São Paulo. All individuals were scanned with a 1.0 mm AL filter at two different resolutions: medium (MED: 18 μm, 70 kV, 356 μA) and high (HIGH: 9 μm, 65 kV, 380 μA). These resolutions correspond to different voxel sizes, the smallest volume unit in the 3D volumes: 3.85 10^{−6} mm^{3} for MED and 1.92 10^{−6} mm^{3} for HIGH, which are much smaller than conventional CTs (voxel sizes ranging from 0.1 to 5 mm^{3}). Before scanning, toads were wrapped with Parafilm to avoid too much alcohol evaporation and consequent dehydration, which could lead to blurred images by sample contraction. The scanning time per skull for MED was 20 min. and for HIGH was 1.5 h. After scanning, the skulls were reconstructed using NRecon software (SkyScan, Konitch, Belgium). In this process, 2D projection images are reconstructed to crosssection images by use of a mathematical algorithm (Feldkamp). The first step in the reconstruction process was to choose the lower and upper limits of the threshold for the linear attenuation coefficient (AC). The AC measures how much the intensity of the Xray beam is reduced as it passes through the material being scanned, and is related to the density of the materials. Based on the histograms of soft tissue and bone density of the skulls, we chose an AC threshold of 0.0 and 0.05 for all specimens scanned. Afterwards, we applied different types of corrections in the reconstruction process to soften some undesirable effects in the images: postalignment = −5.0 to 2.0; ring artifact reduction = 2.0 to 4.0; beamhardening correction = 30 % and smoothing = 2.0. These corrections are important to avoid blurring and artifacts in the 3D images and for all corrections the values used were small compared to maximum values. Therefore, we do not expect these corrections to interfere in landmark visualization, but instead to improve its identification. The sequences of crosssection images varied from 400 to 600 images for MED (500 MB to 1 GB per skull file size) and from 800 to 1,200 images for HIGH (1 to 5 GB per skull file size) and were all in BMP extension.
For the filter type analysis, we used 20 adult museum specimens of Rhinella pygmaea (MZUSP, São Paulo, Brazil). This species was chosen because several specimens presented too much transparency or even holes in their skulls, especially in the squamosal, sphenethmoid and parasphenoid bones, preventing the placement of some landmarks and also suggesting that there were differences in bone density among individuals (Fig. 1). We scanned the 20 toads at medium resolution using two distinct filters: 1.0 mm AL (AL1.0: 70 kV, 356 μA) and 0.5 mm AL (AL0.5: 50 kV, 500 μA). Filters are thin metal sheets that are set in front of the Xray source and can have different thickness. Filtration retains a part of the low energy photons of the Xray, thus increasing mean Xray energy. The reduction in AL filter thickness from 1.0 mm to 0.5 mm results in a lower mean energy of the Xray because photons with lower energy traverse the filter and achieve the thinner bones, being retained on them [22]. The individuals used in this analysis presented low transparency and no holes since the objective was to evaluate whether the use of distinct filters changed the landmarks positions and the linear distances obtained in the same individuals. In these toads all landmarks could be placed, which would not be the case if we had used specimens with holes or high transparency. CT system, scanning and reconstruction procedures were the same as described above.
Landmarking procedure and linear distances in 3D images
We placed 20 landmarks at bone sutures (type I landmarks) or bone processes (type II landmarks; 35 landmarks if counting both sides of the skull, left and right) in all views (dorsal, ventral and lateral; Table 1 and Fig. 1) using TINALandmark software [7]. This software provides different views of the data: a 3D volume of the whole skull and three 2D views of the orthogonal crosssections (axial, sagittal and transversal) throughout the image sequence, enabling landmarks to be precisely placed given that the views are linked (i.e. when setting a landmark in the 3D volume, it also appears in all 2D views allowing a refinement of the landmark position). The skull’s crosssection sequences were loaded in TINALandmark by using the “Sequence Tool” after the BMP files were converted to DICOM files. The Sequence Tool presents two types of image compacting factors for loading heavy files: “Stride”, which reduces the sequence in the interslice direction (i.e., stride = 2.0 and “Stride average” = ON means that an average image will be loaded at every two crosssections from the sequence, resulting in half of the original size of the sequence file); and “Downsample”, which reduces the sequence along the x and y directions (i.e., downsamples the size of the pixels; [23]). MED resolution data were loaded using stride = 2.0 whereas HIGH resolution data were loaded using stride = 3.0 and downsample = 2.0, since they were heavier and failed to load with the same compacting factor than MED data. The downsampling procedure only affects the visualization of the 3D volume of the skulls, not changing the visualization of the crosssection images (axial, sagittal and transversal views).
Appropriate visualization of 3D volumes and determination of bone threshold value (the average density of soft tissue and bone to guarantee that landmarks are placed exactly where the mouse is pointing in the 3D image) were obtained by following the recommendations contained in the TINA Geometrics Morphometrics Toolkit manual [23]. Finally, after all steps of image adjustments, a landmark list was loaded with the 35 points and each was placed two times (two replicates per individual) in the skull’s images of R. granulosa and R. pygmea individuals, so we could assess measuring reliability. The corresponding x, y and z landmark coordinates were saved in a TXT file, which was later loaded in the R programming environment [24] where distances were calculated.
We determined 24 linear distances in the skull of the toads, all distances representing individual bone dimensions (Table 2, except distance 7 which corresponds to orbit size). To obtain distances in mm, the coordinates of MED and HIGH were multiplied by 0.01742 and 0.00871, respectively, which are the size of the pixels in mm for each resolution. In the case of the R. pygmea specimens that were scanned with different filters, landmark coordinates were multiplied by 0.01742 because they were only scanned at medium resolution. Distances from replicates of the same individual (within each method) were inspected for gross measurement error (difference between replicates above 0.5 mm for most distances, except for distances 5, 8, 10, 11, 13, 17 and 24, which have means smaller than 3.0 mm, and were controlled for error above 0.3 mm), and when detected, the landmarks correspondent were placed again in both replicates and corrected. The reference to consider the magnitude of a gross error was based on the precision of the 3D digitizer, which is 0.01 mm. This procedure was adopted because these gross errors could lead to misleading conclusions about the betweenmethods analysis, since they are actually referred to gross human error inside each method. In order to analyze the possible effect on the landmarks and distances of switching the compacting factors (stride and downsample) when loading the crosssection sequences, we placed the landmarks in the R. granulosa skulls scanned in medium resolution using the following values in TINALandmark: (S2) stride = 2.0; (S3) stride = 3.0; (S2/D2) stride = 2.0 and downsample = 2.0; and (S3/D2) stride = 3.0 and downsample = 2.0. Again, each individual at each of these four compacting values were measured twice and gross measurement error between replicates inside each compacting factor was controlled. Thus, total sample size for the compacting factor analysis was 20 individuals for each of the four situations, and 80 in total.
Landmarking procedure and linear distances using 3D digitizer
The same R. granulosa adults that we scanned with MED and HIGH resolutions were cleaned and had their skin removed. In this process, we lost four of the skulls by crushing, 16 skulls remaining. We placed the same landmarks described above in the cleaned skulls by using a 3D digitizer (Microscribe 3DX, IL). Since the toad skulls are quite small (total length around 20.0 mm), we used a binocular loupe to mark the landmarks with a pencil. Then, we placed the digitizer pen in the graphite marks to digitize the landmarks. We placed landmarks in all individuals twice and gross error in the distances from one replicate to the other was also controlled as described before at the time of landmark digitalization. A TXT file was created for the coordinates obtained from DIG and was loaded in R environment in order to calculate the same distances as MED and HIGH. The landmarks and distances obtained by using DIG were considered as the reference for comparisons.
Withinmethod mean individual landmark distances
To assess the differences of placing landmarks in the skulls 3D volumes compared to placing them in the real skulls with the 3D digitizer, we used a superimposition method, the Generalized Procrustes Analysis (GPA), to estimate the sample mean landmarks (mean shape) and to calculate individual distances from the mean landmarks within each method. The GPA sumperimposes landmark configurations of several individuals, scaling them to unit centroid size, and uses leastsquare estimates for translation and rotation parameters [25]. These transformations of the landmarks are needed because there is no natural coordinate system that is common to all individuals that are digitized and a common shape space is achieved when doing this procedure [26, 27]. However, this transformation process confounds variation at different homologous landmarks, i.e., the most variable landmarks among individuals (nonisotropic variation) have their variation spread across other landmarks [16, 26, 27]. Thus, to avoid this problem in our landmark variation analysis, we followed the idea presented in van der Linde and Houle (2009) [27] and excluded from the GPA one landmark at a time to estimate the sample mean shape and the rotation matrices for the 19 remaining landmarks for each individual. By doing that, we precluded that the variation in the landmark of interest got spread through the other landmarks. Afterwards, we multiplied the excluded landmark of each individual by its corresponding rotation matrix (all landmarks of an individual are rotated in the same angle) to have all landmarks in the same new coordinate system, the so called shape space. The last step was to calculate the individual distance of the landmark of interest from the mean landmark (in the shape space, therefore with no scale). This procedure was repeated 20 times for each method, so that all 20 landmarks were excluded at each time and the distances from the individuals and the mean landmark could be calculated. To compare the different methods, we also calculated the mean individual landmark distance from the mean landmarks. Finally, in order to compare the methods of placing landmarks in a more intuitive scale, we multiplied the mean deviations and SD values by its corresponding mean sample centroid size in mm. Only the medial and left side landmarks were used in this analysis and the GPA was done with the “geomorph” [28] and “shapes” [29] packages in R environment.
Distance repeatability within and between methods
We calculated repeatability values for all distances in every data subset using the calculations described in [20]. Each individual is a group with two replicates and the betweengroup variance is the sum of squares of the deviations of the group means from the total mean, whereas the withingroup variance corresponds to the sum of squares of the deviations of each replicate from its own group mean. The index is calculated as follows:
s^{2} _{A} being the variance among groups and s^{2} the residual variance, both calculated from the sum of squares in an ANOVA.
We calculated repeatabilities for MED, HIGH and DIG separately, and also for joined data sets, (MED + HIGH), (MED + DIG), (HIGH + DIG), with the three betweenmethod data sets made up of four replicates per individual. In the same manner, we calculated repeatability for replicates of the different values of stride and downsample, S2, S3, S2/D2 and S3/D2 separately; and also of joined data sets, (S2 + S3), (S2 + S2/D2), (S2 + S3/D2), (S3 + S2/D2), (S3 + S3/D2) and (S2/D2 + S3/D2). Finally, we did the same calculations for the different filters data sets separately, AL1.0 and AL0.5, and also joining all the replicates (AL1.0 + AL0.5). The distance repeatabilities for separated data sets indicate the withinmethod reliability of measuring the skulls, whereas the repeatabilities for joined data sets comprise both within and betweenmethod reliability.
Raw and absolute differences within and between methods
In addition to the repeatability analysis, we also calculated raw and absolute differences of the distances between replicates, within and between methods. In the first case, we just subtracted the distances between replicates of the same individual and calculated the mean raw and absolute differences across all individuals for each distance. These calculations were made for each data set separately: MED, HIGH, DIG, S2, S3, S2/D2, S3/D3, AL1.0 and AL0.5. The betweenmethod analysis was done by subtracting the distances between the same individuals measured with different methods (e.g., in the case of MEDHIGH, distances from individual 1 of MED were subtracted from the distances of individual 1 of HIGH, and so on for all 20 individuals) and calculating the raw and absolute mean differences across all individuals for each distance. All the differences that involved DIG (MEDDIG and HIGHDIG) had 16 (4 skulls were lost in the cleaning procedure), while the rest of the differences between methods had 20 individuals. The raw differences show whether any method presents a consistent bias (e.g. under or overestimating the distances) in the data. The absolute differences indicate the magnitude of the error within and between methods, independent if the errors in the distances are in one direction or another. To get an idea of the error in relation to the distances means, we calculated the percentage error by dividing the absolute mean differences by their corresponding means and multiplying by 100.
Statistical analysis
Correlation tests were all done using Pearson product moment correlation and significance level of 0.05. All ANOVAs in the repeatability analysis, correlation tests and graphics were done in the R environment [23].
Abbreviations
 AL0.5:

Aluminum filter 0.5 mm
 AL1.0:

Aluminum filter 1.0 mm
 DIG:

3 D digitizer
 HIGH:

High resolution
 MED:

Medium resolution
 D2:

Downsample = 2.0
 S2:

Stride = 2.0
 S3:

Stride = 3.0
 s^{2} :

Withingroup variance
 s^{2} _{A} :

Amonggroup variance
References
 1.
Costantini D, Alonso ML, Moazen M, Bruner E. The Relationship Between Cephalic Scales and Bones in Lizards: A Preliminary Microtomographic Survey on Three Lacertid Species. Anat Rec Adv Integr Anat Evol Biol. 2010;293:183–94.
 2.
Ekdale EG. Ontogenetic Variation in the Bony Labyrinth of Monodelphis domestica (Mammalia: Marsupialia) Following Ossification of the Inner Ear Cavities. Anat Rec Adv Integr Anat Evol Biol. 2010;293:1896–912.
 3.
Wilkinson M, San Mauro D, Sherratt E, Gower DJ. A ninefamily classification of caecilians (Amphibia: Gymnophiona). Zootaxa. 2011;2874:41–64.
 4.
Cuff AR, Rayfield EJ. Feeding Mechanics in Spinosaurid Theropods and Extant Crocodilians. PLoS One. 2013;8:e65295.
 5.
Gignac PM, Kley NJ. Iodineenhanced microCT imaging: Methodological refinements for the study of the softtissue anatomy of postembryonic vertebrates. J Exp Zool B Mol Dev Evol. 2014;322:166–76.
 6.
Rosset A, Spadola L, Ratib O. OsiriX: An OpenSource Software for Navigating in Multidimensional DICOM Images. J Digit Imaging. 2004;17:205–16.
 7.
Schunke AC, Bromiley PA, Tautz D, Thacker NA. TINA manual landmarking tool: software for the precise digitization of 3D landmarks. Front Zool. 2012;9:6.
 8.
Kohn LA, Cheverud JM. Issues in evaluating repeatability of an imaging system for use in anthropometry, Proceedings of the Electronic imaging of the Human bodyworking group. OH: Crew System Ergonomics Information Analysis Center, Dayton; 1992.
 9.
Kim G, Jung HJ, Lee HJ, Lee JS, Koo S, Chang SH. Accuracy and Reliability of Length Measurements on ThreeDimensional Computed Tomography Using OpenSource OsiriX Software. J Digit Imaging. 2012;25:486–91.
 10.
HalperinSternfeld M, Machtei E, Horwitz J. Diagnostic Accuracy of Cone Beam Computed Tomography for Dimensional Linear Measurements in the Mandible. Int J Oral Maxillofac Implants. 2014;29:593–9.
 11.
Zelditch ML. Ontogenetic variation in patterns of phenotypic integration in the laboratory rat. Evolution. 1988;42:28–41.
 12.
Cheverud JM. Morphological Integration in the Saddleback tamarin (Saguinus fuscicollis) cranium. Am Nat. 1995;145:63–89.
 13.
Young NM, Hallgrímsson B. Serial homology and the evolution of mammalian limb covariation structure. Evolution. 2005;59:2691–704.
 14.
Porto A, de Oliveira FB, Shirai LT, De Conto V, Marroig G. The Evolution of Modularity in the Mammalian Skull I: Morphological Integration Patterns and Magnitudes. Evol Biol. 2009;36:118–35.
 15.
Corner BD, Lele S, Richtsmeier JT. Measuring Precision of ThreeDimensional Landmark Data. J Quant Antrophology. 1992;3:347–59.
 16.
Richtsmeier JT, Paik CH, Elfert PC, Cole TM, Dahlman HR. Precision, Repeatability and Validation of the Localization of Cranial Landmarks Using Computed Tomography Scans. Cleft Palate Craniofac J. 1995;32:217–27.
 17.
Zeldicth ML, Swiderski DL, Sheets HD, Fink WL. Geometric Morphometrics for Biologists: A Primer. New York and London: Elsevier Academic Press; 2004. p. 237.
 18.
Richard AH, Parks CL, Monson KL. Accuracy of standard craniometric measurements using multiple data formats. Forensic Sci Int. 2014;242:177–85.
 19.
Stull KE, Tise ML, Ali Z, Fowler DR. Accuracy and Reliability of measurements obtained from computed tomography 3D volume rendered images. Forensic Sci Int. 2014;238:133–40.
 20.
Lessells CM, Boag PT. Unrepeatable Repeatabilities: A Common Mistake. Auk. 1987;104:116–21.
 21.
Fernandes, TMF, Adamczyk, J, Poleti, ML, Henriques, JFC, Friedland, B, Garib, DG: Comparison between 3D volumetric rendering and multiplanar slices on the reliability on linear measurements on CBCT images: an in vitro study. Journal of Applied Oral Science 2014.
 22.
SkyScan. SkyScan 1176. In vivo XRay Microtomograph Instruction Manual. 2011.
 23.
Bromiley PA, Ragheb H, Thacker NA. The TINA Morphometrics Geometric Toolkit, Tina Memo No 2010–007. 2012.
 24.
R Core Team: R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria 2013. URL: http://www.Rproject.org/.
 25.
Bookstein, FL: Morphometrics tools for landmark data: Geometry and Biology. Cambridge University Press. 1997. p. 425.
 26.
Lele S. Euclidean distance matrix analysis (EDMA): Estimation of mean form and mean form difference. Math Geol. 1993;25:573–602.
 27.
van der Linde K, Houle D. Inferring the Nature fo Allometry from Geometric Data. Evol Biol. 2009;36:311–22.
 28.
Adams DC, OtarolaCastillo E. geomorph: an R package for the collection and analysis of geometric morphometric shape data. Methods Ecol Evol. 2013;4:393–9.
 29.
Dryden, IL: shapes: Statistical Shape Analysis. URL: http://cran.rproject.org/web/packages/shapes/index.html.
Acknowledgments
We thank Dr. Taran Grant from MZUSP for the loan of the R. pygmaea specimens used in this study. MNS also thanks Francisco de Assis and family for providing house in the field, and Pedro Leite Ribeiro for helping to collect and preserve R. granulosa specimens used in the study. We thank the technician Marco Antônio Puodizius for helping cleaning the toad skulls and the Departamento de Zoologia at Instituto de Biociências, USP, for allowing the use of their facilities. We also thank Guilherme Garcia for assistance with the landmark analysis and the two anonymous reviewers for their very careful revision of our manuscript and for their important criticisms. This study was supported by Fundação de Amparo à Pesquisa de São Paulo (FAPESP), Ph.D. scholarship 2011/075842, multiuser project 2009/542034 and thematic project2011/142957.
Author information
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
MNS conceived the general idea of the study, collected and analyzed the data and drafted the manuscript. GM conceived the design of the study, analyzed and discussed the data and reviewed the manuscript. All authors read and approved the final manuscript.
Monique Nouailhetas Simon and Gabriel Marroig contributed equally to this work.
Additional file
Additional file 1: Table A1.
Withinmethods mean individual landmark distance without scale. Table A2. Between methods mean raw differences (in mm) in all the distances for joined data sets of microCT resolutions and 3D digitizer and of different filter type. Table A3. Mean percentage error in all distances for joined data sets of microCT resolutions and 3D digitizer and for distinct filter types.
Rights and permissions
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Received
Accepted
Published
DOI
Keywords
 Bone density
 Distance repeatability
 Geometric morphometrics
 Traditional morphometry
 Toad skulls
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate. Please note that comments may be removed without notice if they are flagged by another user or do not comply with our community guidelines.