 Research article
 Open Access
 Published:
Standardization of threedimensional pose of cylindrical implants from intraoral radiographs: a preliminary study
BMC Oral Health volume 21, Article number: 100 (2021)
Abstract
Background
To introduce a theoretical solution to a posteriori describe the pose of a cylindrical dental fixture as appearing on radiographs; to experimentally validate the method described.
Methods
The pose of a conventional dental implant was described by a triplet of angles (phipitch, thetaroll, and psiyaw) which was calculated throughout vector analysis. Radiographic and simulatedimage obtained with an algorithm were compared to test effectiveness, reproducibility, and accuracy of the method. The length of the dental implant as appearing on the simulated image was calculated by the trigonometric function and then compared with real length as it appeared on a twodimensional radiograph.
Results
Twenty radiographs were analyzed for the present in silico and retrospective study. Among 40 fittings, 37 resulted as resolved with residuals ≤ 1 mm. Similar results were obtained for radiographic and simulated implants with absolute errors of − 1.1° ± 3.9° for phi; − 0.9° ± 4.1° for theta; 0° ± 1.1° for psi. The real and simulated length of the implants appeared to be heavily correlated. Linear dependence was verified by the results of the robust linear regression: 0.9757 (slope), + 0.1344 mm (intercept), and an adjusted coefficient of determination of 0.9054.
Conclusions
The method allowed clinicians to calculate, a posteriori, a single real triplet of angles (phi, theta, psi) by analyzing a twodimensional radiograph and to identify cases where standardization of repeated intraoral radiographies was not achieved. The a posteriori standardization of twodimensional radiographs could allowed the clinicians to minimize the patient’s exposure to ionizing radiations for the measurement of marginal bone levels around dental implants.
Background
Osseointegrated dental implants supporting fixed and removable prostheses were found to be highly effective rehabilitation strategies in the management of edentulous patients.
The radiographic assessment of the ongoing periimplant vertical bone loss after implant placement was an essential issue for clinicians to address [1].
Periapical and panoramic radiographs and computerized tomography allowed the clinician to measure marginal bone levels around a dental implant to evaluate its longterm stability, but with a varying degree of accuracy and reliability [1].
In twodimensional radiography, different types of intraoral and extraoral equipment met the requirements for standardized evaluation of periimplant marginal bone loss [2]. These types of equipment took into account the analytic basis of projective geometry which consisted of specific features used for describing the spatial relationships between xray source and object, and the complex relationships existing between the various layers (object plane image and detector plane image) [3].
The determination of the threedimensional position and orientation of known objects from single perspective projections were initially proposed by Hoffmann and Esthappan with a priori knowledge of positions of the specific markers [4]. Generally, external markers (alone or in combination with a holding device) were and are employed to allow for an accurate estimate a posteriori of the object’s position and orientation, thereby providing both the referral points in space (set of three radiopaque ball bearings or spheres) and the known dimensions of the object being analyzed [5, 6].
Others suggested the possibilities of estimating fixture angulation, dimension, and position through distortion on its shadow as it appeared in the projective plane [7,8,9,10,11].
In computer vision, the pose of an object was the combination of position and orientation relative to some external and independent coordinate system. Here, pose referred to the orientation of the cylindrical implant concerning the imagecoordinate frame. The “withinobject features” was a direct consequence of the abovementioned footsteps, proposing to use the dental implant as a reference object of known dimensions in the image without any external markers [12]. That's was made to minimize information about the position and to focus solely on the orientation. The procedure consisted of making the most of the more complex radiographic shadow of a conventional dental implant, which had, at least in part, the shape of a right circular cylinder with known diameter, brought.
In previous studies, rotational poses of the nonconventional implants, with or without effects of image distortion, had been described; the components of the lengthwise and crosswisevector allowed clinicians to estimate sets of angles describing the position of an implant based on its appearance on a single twodimensional radiograph [13].
However, the abovementioned methods worked only with a noncylindrical object, such as a blade implant, in which lengthwise and crosswisevector were orthogonal to one another and easy to find using an implant projection. At first glance, it might seem that rotation of the implant around its main axis does not provide meaningful information since the implant had rotational symmetry. Rotating the implant about the long axis did not roughly influence its radiographic shadow but any degree turn about this axis resulted in a rotated image projection of the bone surrounding the implant by the same degree. So the calculus of this angle from a radiographic shadow was extremely important to ensure consistent reproducibility in marginal bone assessment, which was measured along the line of the implant.
The meaning of the paper was that the authors took patient radiographs with implantimages (implant dimensions were known from CT or optical scans) and calculated their “pose” relative to a referring point using an algorithm [3]. Then the authors simulated this position again using another algorithm [3] and compared computed differences between real measurements and the ones from the simulations.
The present article will:

(1)
introduce a theoretical solution to obtain a triplet of angles (φ, θ, and ψ) describing the pose of the cylindrical fixture as shown, a posteriori, in a single twodimensional radiograph;

(2)
present in silico and experimental application of the method for accuracy;

(3)
verify the effectiveness of the method, attempting to relate reallength data to those calculated by the angular correction factor simulation.
Methods
Theory
Especially in the measurement of the marginal bone level around dental implants in repeated intraoral radiographs of the same patient, a very similar acquisition process, for which it's hard to quantify the number of errors attributed to a series of characteristics, was needed. To properly fulfill standardized twodimensional periapical radiography of a cylindrical dental implant in the alveolar bone, the following general requirements had to be respected:

Threedimensional position: angulations and translations in threedimensional space of the dental implant as equal to each other as possible could be achieved with a customized bite block often formed of a rigid material combined with an extension cone paralleling device (an XCPORA, the eXtention Cone ParallelingOne Ring and Arm positioning system, Dentsply International, Elgin, Illinois, USA);

Magnification: an image’s size could be either enlarged or reduced depending on the distance from the xray source to the detector, that is directly proportional to the length of the XCPORA’s arm;

Distortion: it could arise from an angle between the implant longitudinal axis and the imaging plane or a very oblique projection. An error in the misalignment could most likely be generated because the Indicating Device (PID) of the Dental Xray Tube (DXT) head was manually and poorly placed by the clinician without the help of the XCP device or because the patient moved during the measurement.
The effects of both the image magnification level and the relative translation between dental implant and frame of the image were not significant because the dental implant itself was the object of the investigation and acted as a benchmark for further analyses. After the components related to translation and magnification had been deleted, just three rotation components might be determined for describing implant orientation. The simplest way to do it is to describe the initial position of the cylindrical dental implant, given by the position of the apices described as the vector v_{0} = (x_{0}, y_{0}, z_{0}). which is perpendicular to the ground (that is, perpendicular to the XYplane along the implant direction Z_{i} and parallel to both the XZ and YZplane, with implant directions Y_{i} and X_{i}, respectively as appearing in Fig. 1). After applying a rotational matrix to the vector, its final position will be v_{1} = (x_{1}, y_{1}, z_{1}) which sets the viewpoint after the rotationalsteps have been completed with three angles of rotation, pitchφ, rollθ, and yawψ, those correspond to three separate and successive rotations along the implant main direction X_{i}, Z_{i}, and Y_{i}.
The vector is rotated in space and projected on the detector plane (YZplane) as described by the following components which depended on the triad of angles (φ, θ, ψ) and the magnitude of the vector’s components. For each main axis a Rodrigues rotational matrix, \({R}_{\widehat{\omega }}\left(\alpha \right)={e}^{\widehat{\omega }\alpha }\) corresponding to a rotation by an angle α about a fixed axis specified by the unit vector \(\widehat{\omega }\) with the following explicit formula \({R}_{\widehat{\omega }}\left(\alpha \right)=\widehat{i}+\widehat{\omega }\mathrm{sin}\alpha +{\widehat{\omega }}^{2}\left(1\mathrm{cos}\alpha \right)\) has been applied. T_{1} is the rotational matrix obtained by rotating the implant along the Y_{i}direction and angle φ with its explicit form from Rodrigues’s second equation:
T_{2} is the second rotational matrix obtained by rotating the implant along the Z_{i}direction (after rotation of φ) with angle θ:
T_{3} is the third rotational matrix obtained by rotating the implant along the X_{i}direction (after φ and θ rotations) with angle ψ:
When the three rotational matrices were combined in the following order T1 ⊗ T2 ⊗ T3, the matrix T_{4} simplifies to become:
When the transformation is expressed in the form of vectors, v_{1} = T_{4} ⋅ v_{0}
the vector v_{1} can be written in the following explicit form:
Equations 9, 10 and 11 can be used to describe any vector in the space (like how a combination of arbitrary x_{0}, y_{0}, and z_{0} vectors).
To resolve the threedimensional rotation of a cylindrical object two vectors are required. So lengthwise and crosswise vectors are introduced to simplifies the equation system. As depicted in Fig. 1, the lengthwise vector, \(\overrightarrow{\mathrm{OA}}\), goes from the point located at the center of the implant shoulder, that is point \(\mathrm{O}=\left(\mathrm{0,0},0\right)\), to the point in the center of the implant apex, that is the point \(\mathrm{A}=\left({\mathrm{A}}_{\mathrm{X}},{\mathrm{A}}_{\mathrm{Y}},{\mathrm{A}}_{\mathrm{Z}}\right)\).
The crosswise vector, \(\overrightarrow{\mathrm{OR}}\), goes from the point \(\mathrm{O}=\left(\mathrm{0,0},0\right)\) to point \(\mathrm{R}=\left({\mathrm{R}}_{\mathrm{X}},{\mathrm{R}}_{\mathrm{Y}},{\mathrm{R}}_{\mathrm{Z}}\right)\).
The method allows clinicians to calculate the implant pose with a known reference object. In this perspective, the center of another implant shoulder (point R of implant 2 in Fig. 1), which has a different pose from the currently analyzed implant (implant 1 in Fig. 1), can be used to disclose the pose rotation along the long axis of the implant 1. Any other fixed reference point can be used in place of the shoulder of implant 2, however (referring object of Fig. 1).
Since the depth in a single twodimensional image is not taken into account, it's been arbitrarily established that the dimension along Xaxis is disregarded. So the system of equations obtained with the known magnitude of \(\overrightarrow{\mathrm{OA}}\) and \(\overrightarrow{\mathrm{OR}}\) are:
The components viewed from the detector plane and measured within a single twodimensional radiograph are \({\overrightarrow{\mathrm{OA}}}_{\mathrm{Y}}\),\({\overrightarrow{\mathrm{OA}}}_{\mathrm{Z}}\) and \({\overrightarrow{\mathrm{OB}}}_{\mathrm{Y}}\),\({\overrightarrow{\mathrm{OB}}}_{\mathrm{Z}}\), and depend on the triad of angles (φ, θ, ψ) and on the vectors’ components \({\overrightarrow{\mathrm{OA}}}_{0}=\left(\begin{array}{c}0\\ 0\\ l\end{array}\right)\) and \({\overrightarrow{\mathrm{OR}}}_{0}=\left(\begin{array}{c}0\\ {\mathrm{R}}_{\mathrm{Y}}\\ {\mathrm{R}}_{\mathrm{Z}}\end{array}\right)\)The four vectorcomponent equations undergo simplification using the known magnitude of the lengthwise (\(\overrightarrow{\mathrm{OA}}\)) and crosswise vector (\(\overrightarrow{\mathrm{OR}}\)).
in which R_{Y} and R_{Z} can be directly measured by a threedimensional viewing software (dentascan or CAD/CAM ComputerAided Design/ComputerAided Manufacturing technology) or explicitly calculated by Eqs. 20 and 21 as reported in the “Appendix”.
The R_{X} parameter is 0, whereas R_{Y} and R_{Z} are described as the following:
with n from 1 to 3, and m from 1 to 3 in the interest of brevity of the equations which depend on the three dimensional coordinates of \(\mathrm{O}=\left({o}_{1},{o}_{2},{o}_{3}\right)\), \(\mathrm{A}=\left({a}_{1},{a}_{2},{a}_{3}\right)\) and \(\mathrm{R}=\left({r}_{1},{r}_{2},{r}_{3}\right)\) obtained from either CT or optical scans.
If the marginal bone level (MBL) has to be measured along the \(\overrightarrow{\mathrm{OA}}\), a correction factor by using the dental implant itself as internal standard must be employed. The correction factor can be obtained by considering the length of the implant (len) expressed by the equation of the magnitude of \(\overrightarrow{\mathrm{OA}}\) along the plane of the detector (YZplane).
The value of the correction factor CF is given by Eq. 23:
where MBL and MBL^ are, respectively, the measured and corrected values of the marginal bone level.
The calculated correction factor is identical to that calculated in a previous study [13].
Experimental evaluation
CT scans and radiographs of patients who had undergone dental implant placement were retrospectively collected from the case sheets between the years 2014 and 2016. Patients were informed about the type of the present in silico study and signed an informed consent form before data selection and analysis and approval for this retrospective analysis was obtained from the local Ethical Committee. Implants were usually placed according to the instructions of the manufacturers employing an electric handpiece, electromagnetic (Magnetic Mallet, Osseotouch www.osseotouch.com, Turbigo, Milano, Italy) and piezoelectric devices (Piezosurgery, Mectron, Carasco, Italy). Ethics approval was not required for this in vitro study. The present paper was prepared according to the SQUIRE2.0 checklist.
The conebeam computerized tomography scan (Gendex GXCB500, Gendex Dental Systems) provides a finelydetailed threedimensional overview of the bony architecture around the implants inserted in it. The following set of default parameters have been applied in the acquisition of scans: 120 kV, 30.89 mAs, 0.2 mm × 0.2 mm × 0.2 mm isotropic voxel size, and 8.72 mm diameter of Field Of View. Digital intraoral periapical radiographs were taken (70 kVp, 7 mA) using a CMOS (Complementary Metal Oxide Semiconductor) digital sensor (Schick CDR Elite, Schick Technologies Inc., Long Island City, NY, USA) and imaging software (FONA Computed Dental Radiography DICOM 4.5, Schick Technologies Inc., Long Island City, NY, USA) with a pixel pitch of 100 μm (0.1 mm × 0.1 mm; bit depth, 8bit grayscale).
Experimental simulation
Experimental evidence that supported the validity of the abovementioned vector equations and verified the model's accuracy was given by an algorithm using components \({\overrightarrow{\mathrm{OA}}}_{\mathrm{Y}}\), \({\overrightarrow{\mathrm{OA}}}_{\mathrm{Z}}\), \({\overrightarrow{\mathrm{OR}}}_{\mathrm{Y}}\), \({\overrightarrow{\mathrm{OR}}}_{\mathrm{Z}}\) directly measured on the radiographs as per in Fig. 1 and “Appendix”. Digital prototypes of cylinder screwed, rootform, 1 mmtread pitch dental implants were obtained by a Cone Beam Computerized Tomography scanner (Gendex GXCB500, Gendex Dental Systems) with the following setting: 120 kV, 30.89 mAs, isotropic voxel size of 200 μm. The virtual 3D phantomimplant was voxelized and interpolated from the original.dcm file with the following setting: 100 μm × 100 μm × 100 μm, bit depth 8bit grayscale.
A simulated radiograph could be generated through the overlapping of all the phantomlayers depicting implant in threedimensions along the direction of the xaxis, so obtaining a projected phantomimplant in the YZplane (like it is the detector plane) by a subroutine described in a preceding article [3].
Variables
Determination of numerical variables
For each dental implant, a dentascan program was used to acquire the threedimensional coordinates of the three reference points A, O, and R to calculate l, R_{Y}, and R_{Z}, that is the magnitude of the vectors (\({\overrightarrow{\mathrm{OA}}}_{0}\) and \({\overrightarrow{\mathrm{OR}}}_{0}\)), in other words, input parameters to be inserted in the righthand side of equations from 16 to 19. A free standalone software (Osiris 4.19 the University of Genève. Switzerland) applying a mousedriven measurement tool was used to acquire the twodimensional coordinates of the points A, O, and R to calculate the components \({\overrightarrow{\mathrm{OA}}}_{\mathrm{Y}}\),\({\overrightarrow{\mathrm{OA}}}_{\mathrm{Z}}\) and \({\overrightarrow{\mathrm{OR}}}_{\mathrm{Y}}\),\({\overrightarrow{\mathrm{OR}}}_{\mathrm{Z}}\), that is, input parameters to be inserted in the lefthand side of equations from 16 to 19. Measurements of coordinates and components were taken twice by two independent investigators (CS, PDM). The diameter of the dental implant at healing cap/implantabutment connection (∅ from 4 to 5 mm), which was one of the cleanest and more consistent features to be certainly identified in a cylindrical shaped object, allowed the investigators to calibrate all the images. Numbers in vectors analysis provided a resolution of a hundredth of a millimeter.
Pose angle triplet
One set of angles (triad) was given by a triplet of rotational angles φ, θ and ψ (pitch, roll, and yaw, respectively). The position of each radiographic cylindrical dental implant could be described by a radiographic triad of angles (_{r}φ, _{r}θ, _{r}ψ) which was estimated by numerically solving the system of equations from 16 to 19 by using the algorithm suggested by Toti et al. [13], radiographic components \({{}_{r}\overrightarrow{\mathrm{OA}}}_{\mathrm{Y}}\), \({}_{r}{\overrightarrow{\mathrm{OA}}}_{\mathrm{Z}}\), \({}_{r}{\overrightarrow{\mathrm{OR}}}_{\mathrm{Y}}\), \({}_{r}{\overrightarrow{\mathrm{OR}}}_{\mathrm{Z}}\), and following setting for the algorithms running on a matrix elaborator (MatLab 7.13, The MathWorks, Natick, MA, USA): “Levenberg–Marquardt” ON.
For each radiographic implant (with input: length, diameter triad of angles, and referring points) a projected phantom of a dental implant was generated; analyzing simulated images, simulated components \({{}_{s}\overrightarrow{\mathrm{OA}}}_{\mathrm{Y}}\),\({}_{s}{\overrightarrow{\mathrm{OA}}}_{\mathrm{Z}}\),\({}_{s}{\overrightarrow{\mathrm{OR}}}_{\mathrm{Y}}\),\({}_{s}{\overrightarrow{\mathrm{OR}}}_{\mathrm{Z}}\) were obtained, and the simulated triad of angles (_{s}φ, _{s}θ, _{s}ψ) was calculated as shown in Fig. 2. So radiographic and simulated triads could be statistically compared.
Statistical analysis
The solution of the system of equations was computed with a specific algorithm (MatLab7.13, The MathWorks, Natick, MA, USA). The correlation between the radiographic length of implants and that simulated by using the angular correction factor of implants was tested by robust linear regression. Inter and intraobserver agreements were tested by Bland–Altman analysis. Absolute differences between the radiographic and virtual triad of angles (test of accuracy) were calculated.
Results
Twenty intraoral radiographs (10 in the upper and 10 in the lowerjaw) were selected and analyzed for the present in silico retrospective study. Y and Zcomponents of the twodirectional vectors \(\overrightarrow{\mathrm{OA}}\) and \(\overrightarrow{\mathrm{OR}}\) were measured for each selected cylindrical dental implant. The components were used in a nonlinear fitting procedure, and triads of angles—describing the rotations—were estimated. Table 1 shows the values of the components of the estimated angles with the results of statistical analysis.
Among all the fitting results (40), 37 resulted as resolved with residuals ≤ 1 mm, whereas 32 of them reported on the possible inaccuracies. Results for the 3 nonlinear fittings indicated unsolved equations with residuals > 1 mm. However, the system was well over determined (more equations, four, than unknowns, three).
As was to be expected, how angles and rotations had been set out in the theoretical section indicated that in the upperjaw implants ψ might be of values close to 180° (downwards, in the range between 148.2° and 193.4° with a mean of 169.2°), whereas in the lowerjaw implants ψ ranged from − 33.9° to 52.3° with a mean of 2.2 (upwards) (Fig. 3).
Reproducibility of the angular measurements
Intraexaminer and interexaminer reproducibility were described in detail with Bland and Altman plots in Fig. 2. The intraexaminer differences ranged from − 0.18 to + 0.2 mm for observer 1, and from − 0.25 to + 0.25 mm for observer 2, respectively. The interexaminer differences ranged from − 0.18 to + 0.33 mm for session 1, and from − 0.28 to + 0.16 mm for session 2, respectively.
Accuracy of the angular measurements
The description of the output triad of angles for real radiographs was of − 7.8° ± 29.3°, 30.9° ± 17.1°, and 2.2° ± 27.6° respectively for angles φ, θ, and ψ in the lower jaw, whereas was of 6.5° ± 21.2°, 4.3° ± 29.5°, and 169.2° ± 15.6° respectively for angles φ, θ, and ψ in the upper jaw. Very similar values were obtained for the simulated triad of angles, as they appear in the virtual phantoms: − 4.8° ± 31.9° and 5.7° ± 20.7° for angle φ; 30.6° ± 16.9° and 6.3° ± 25.9° for angle θ; 2.6° ± 28.0° and 168.8° ± 15.2° for angle ψ, for the lower and the upperjaws, respectively. Absolute error, as shown in Fig. 4, was: − 1.1° ± 3.9° for angle φ; − 0.9° ± 4.1° for angle θ; 0° ± 1.1° for angle ψ.
Correction factors for marginal bone level
Radiographs with the related phantoms of cylindrical dental implants, simulated by using the results of triadfitting, are shown in Fig. 2. As shown in Fig. 4, radiographic lengths of the implant as measured on the radiographs and simulated lengths calculated with the angular correction factors appeared to be heavily correlated. Such a linear dependence was verified by the results of robust linear regression, which employed an iteratively reweighted leastsquares algorithm less sensitive to outliers. All collected data were included without any data exclusion and the following coefficients resulted 0.9757 (slope) and + 0.1344 (intercept), with an adjusted coefficient of determination of 0.9054.
Discussion
The number of degrees of freedom in intraoral radiograph describing the projection of an object on the detector’s surface could be reduced by the use of devices for image standardization, that is, bite block with fixedlength extension cone paralleling instruments coupled with equipment indicating the direction of the dental xray tube. In this way, the degrees of freedom went from nine to six; moreover, those related to the object’s translation in space had a negligible effect because the dental implant was the reference object for the radiograph.
The basic idea of the present study was that a cylindrical object did not have enough information, at this stage, to resolve all degrees of freedom about the position and direction of the dental implant.
Several authors attempted to reduce the number of degrees of freedom to simplify the analysis of the physical system. Some authors let the coordinate system be located external to the relevant structures using at least three radiopaque spheres of known dimension [4,5,6]; one of them later had been quite successful in achieving a solution to the descriptions of the geometric structure about the position of the cylindrical dental implant in space [12]; he designed the system of the rigid object to be as compact as possible while describing translational and rotational vector graphics. However, the degree of freedom related to the rotation along the main axis of the dental implant remained unknown. The author suggested that lost rotation might be solved by simply adding an extra reference point not collinear with the dental implant axis. Theoretically, a known pose of the cylindrical dental implant through three reference points, two inside the implant body, and the third in a mesial or distal position (further dental implant, tooth, or identifiable anatomic landmark) reduced all degrees of freedom to one.
In a previous paper, Toti and coworkers suggested that “vector analysis” attempted to estimate the pose of a blade implant [13]. It should be noted that components of the lengthwise vector were the same as those reported in Eqs. 6 and 7; the only adjustment was the value of the parameter l which represented the length of the cylindrical dental implant. On the contrary, the y and z components of the crosswise vector had more complex equations in which the two terms (R_{Y} and R_{Z}) depended on both the threedimensional position of the third reference point (R) and the triad of angles (φ, θ, ψ).
The system given by equations from 6 to 9 admitted explicit solutions; thereby disregarding solutions could be doubled or quadrupled if inverse trigonometric functions were introduced. An equation solver allowed clinicians to obtain very easily an explicit solution given in terms of a single triad of angles. The measurements of fixed parameters (l, R_{Y}, R_{Z}) and input variables (\({\overrightarrow{\mathrm{OA}}}_{\mathrm{Y}}\),\({\overrightarrow{\mathrm{OA}}}_{\mathrm{Z}}\),\({\overrightarrow{\mathrm{OR}}}_{\mathrm{Y}}\),\({\overrightarrow{\mathrm{OR}}}_{\mathrm{Z}}\)) appeared to give reproducible and reliable results when analyses had been performed with Bland and Altman test (inter and intraobserver differences ranging from − 0.28 mm to + 0.33 mm).
Unfortunately, the present formulation did not take into account the distorting effects produced by the misalignment between the detector ring radius and the xray beam caused by neglect of radiological guidelines or involuntary patient movement [14]. Even if a theoretical approach to assess the quality of the image standardization had not yet been developed for a cylindrical object, data points on implant lengths radiographically measured and then calculated using the angular correction factors (Eq. 11) well matched the straight line, so the assumption of the linear dependence was justified and the linear regression model could be applied with acceptable accuracy (value of the adjusted coefficient of determination of 0.9054).
Results of nonlinear fittings accounting for differences between radiographic (_{r}φ, _{r}θ, _{r}ψ) and a simulated triad of angles (_{s}φ, _{s}θ, _{s}ψ) showed that absolute errors had means ranging from − 1.1° to 0° with low measures of dispersion (e.g. standard deviations ≤ 4.1°). The results were in agreement with those described in a previous paper (absolute error from − 0.29 to + 0.64 where 4.35° was the maximum standard deviation) [13].
And if anyone asked where experimental errors arose, it’s supposed to be caused by the cylindrical shape of the dental implant for which an exact calculation of the coordinates of three referring points (O, A, and R) could difficulty result. Moreover, virtual phantom obtained from twodimensional radiographic simulation had not high enough resolution for such application (pixel resolution of 0.05 mm). Another source of error in the calculation of components could arise from the divergence between the xray beam and the normal vector to the detector surface. The problem of the misalignment could not yet be resolved a posteriori for a cylindrical system, at least so far.
Cliniciancs need to be especially mindful to understand that all the processes implemented using algorithms [13], such as measurements of the components of the lengthwise and crosswise vectors and angles calculation of the pose or virtual simulation of the pose of the cylindrical implant are completely independent from each other; differences between true and simulated could only be due to measurement errors (and maybe discretization errors due to the pixel values).
Instead of making a radiograph under known geometric conditions, then computing the pose of the implant from it, and then comparing the pose to the known pose in the imaging geometry, the authors decided to directly measure the pose of the implant on a real radiograph and to simulate obtained angular results. This is because it is very difficult to experimentally check the position of a cylindrical object through three angles one of which including the component describing a rotation of the object along the main axis (hidden component). Moreover, there is also a slight misassumption included in the definition of the projection of the realworld points O, A, and R: due to the projection geometry, the center of the elliptical shadow cast by the (upper) end of the cylinder does not coincide with the realworld centerpoint along the longitudinal axis (A). However, the influence of this misassumption on the accuracy of the method is expected to be negligible, as for common geometries describing the difference could be of subpixel magnitude.
Externalization of the present results (system of Eqs. 16–19) could be used for pose determination of any conical/cylindrical object observed in twodimensional radiographs. For example, the present method to evaluate the angulations of teeth, knowing their length and through the use of only one referring point, offered an interesting perspective in the fields of orthodontics.
There is no doubt that any improvement in the resolution of the detector and automated detection of parameters and components could help refine the method. The authors are trying to solve equations involving parameters related to errors in alignment.
There are limitations in the present study. The main concern is the small sample size with only twenty intraoral radiographs (ten each from upper and lower jaw) that may be not enough to have confirmative conclusions. Another potential bias is that the intraoral radiographs were taken from anterior or posterior regions with possible different results due to different bone thickness.
Conclusions
Theoretical and experimental methods, estimating the pose of a cylindrical dental implant in threedimensional space allowed clinicians to calculate a posteriori single triad of angles (φ, θ, ψ) by analyzing twodimensional images. The value of the correction factor calculated based on data obtained by analysis of angles supplies the clinician with additional information for the real measurement of the implant marginal bone level. The a posteriori standardization of twodimensional radiographs could allow the clinicians to minimize the patient's exposure to ionizing radiations for the measurement of success or failure rates. Finally, information about the rotation of a cylindrical like structure (tooth) could be useful in various fields of dentistry.
Availability of data and materials
All data generated and analysed during this study are published in Table 1 and figures, plus all data generated are available from the corresponding author upon reasonable request.
Abbreviations
 CT:

Computerized tomography
 XCPORA:

EXtention Cone ParallelingOne Ring and Arm positioning system
 PID:

Portable indicating device
 DXT:

Dental Xray tube
 CAD/CAM:

Computeraided design/computeraided manufacturing
 CF:

Correction factor
 MBL:

Marginal bone level
 CMOS:

Complementary metal oxide semiconductor
References
Cassetta M, Di Giorgio R, Barbato E. Are intraoral radiographs reliable in determining periimplant marginal bone level changes? The correlation between open surgical measurements and periapical radiographs. Int J Oral Maxillofac Surg. 2018;47(10):1358–64. https://doi.org/10.1016/j.ijom.2018.05.018.
Meijer HJ, Steen WH, Bosman F. Standardized radiographs of the alveolar crest around implants in the mandible. J Prosthet Dent. 1992;68:318–21. https://doi.org/10.1016/00223913(92)90337a.
MenchiniFabris GB, Covani U, Toti P, Marconcini S, Barone A, Martuscelli R. A methodological approach to standardize and control the quality of the position and alignment of lamina implants on twodimensional radiographs. Oral Radiol. 2020;36:288–306. https://doi.org/10.1007/s11282019003745.
Hoffmann KR, Esthappan J. Determination of threedimensional positions of known sparse objects from a single projection. Med Phys. 1997;24:555–64. https://doi.org/10.1118/1.597938.
Schulze R, Bruellmann DD, Roeder F, d’Hoedt B. Determination of projection geometry from quantitative assessment of the distortion of spherical references in singleview projection radiography. Med Phys. 2004;31:2849–54. https://doi.org/10.1118/1.1796951.
Schulze RK, Weinheimer O, Brullmann DD, Roder F, d’Hoedt B, Schoemer E. Software for automated application of a referencebased method for a posteriori determination of the effective radiographic imaging geometry. Dentomaxillofac Radiol. 2005;34:205–11. https://doi.org/10.1259/dmfr/56357032.
Sewerin IP. Estimation of angulation of Brånemark titanium fixtures from radiographic thread images. Clin Oral Implants Res. 1991;2:20–3. https://doi.org/10.1034/j.16000501.1991.020102.x.
Roeder F, Brüllmann D, d’Hoedt B, Schulze R. Ex vivo radiographic tooth length measurements with the reference sphere method (RSM). Clin Oral Investig. 2010;14:645–51. https://doi.org/10.1007/s0078400903509.
Stoetzer M, Wagner MEH, Wenzel D, Lindhorst D, Gellrich NC, von See C. Nonradiological method for 3 dimensional implant position assessment using an intraoral scan: new method for postoperative implant control. Implant Dent. 2014;23:612–6. https://doi.org/10.1097/ID.0000000000000118.
Vazquez L, Nizamaldin Y, Combescure C, Nedir R, Bischof M, Ehrenfest DM, Carrel JP, Belser UC. Accuracy of vertical height measurements on direct digital panoramic radiographs using posterior mandibular implants and metal balls as reference objects. Dentomaxillofac Radiol. 2013;42:20110429. https://doi.org/10.1259/dmfr.20110429.
von See C, Wagner MEH, Schumann P, Lindhorst D, Gellrich NC, Stoetzer M. Nonradiological method for threedimensional implant position evaluation using an intraoral scan method. Clin Oral Implants Res. 2014;25:1091–3. https://doi.org/10.1111/clr.12214.
Schulze RKW. Pose determination of a cylindrical (dental) implant in threedimensions from a single twodimensional radiograph. Dentomaxillofac Radiol. 2010;39:33–41. https://doi.org/10.1259/dmfr/12523158.
Toti P, Barone A, Marconcini S, MenchiniFabris GB, Martuscelli R, Covani U. Pose determination of a blade implant in three dimensions from a single twodimensional radiograph. Dentomaxillofac Radiol. 2018;47:20170258. https://doi.org/10.1259/dmfr.20170258.
FernandezFormoso N, Rilo B, Mora MJ, MartinezSilva I, Santana U. A paralleling technique modification to determine the bone crest level around dental implants. Dentomaxillofac Radiol. 2011;40:385–9. https://doi.org/10.1259/dmfr/45365752.
Acknowledgements
None.
Funding
No funding was obtained for this study.
Author information
Authors and Affiliations
Contributions
Conceptualization, PT, SC; Methodology, PT, DPO; Software, SC, PT; Validation, MPD, UC, BCB; Formal Analysis, SC, MPD; Investigation, UC, BCBB; Resources, SC; Data Curation, PT, UC, BCB; Writing—Original Draft Preparation, PT, SC, DPO; Writing—Review and Editing, MPD, SC, BCB; Visualization, PT, UC; Supervision, DPO; Project Administration, UC, BCB; Funding Acquisition, UC”. All authors have read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
This retrospective chart review study involving human participants was in accordance with the ethical standards of the institutional and national research committee and with the 1964 Helsinki Declaration and its later amendments or comparable ethical standards. The Human Investigation Committee (IRB) of University of Pisa approved present data analysis (Ethical Approval Form 2626/2008 Protocol Number 58183).
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interest.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Appendix
Appendix
The crosswire vector
A line thorough \(O=\left({o}_{1},{o}_{2},{o}_{3}\right)\) and \(A=\left({a}_{1},{a}_{2},{a}_{3}\right)\) is formed by points \(P=\left(x,y,z\right)\) for which vector \(\overrightarrow{OP}\) has the same direction of the vector \(\overrightarrow{AO}\) as per
So
And the system of coordinates for P is
The plane perpendicular to the vector \(\overrightarrow{AO}\) and passing through \(R=\left({r}_{1},{r}_{2},{r}_{3}\right)\) is
from which Cartesian equation of the plane is the following:
Then equation of the straight line intersecting the piano in the point H can be obtained by the following system:
in which t has the following equation:
and the point \(H=\left({h}_{1},{h}_{2},{h}_{3}\right)\) has the following coordinates:
The Y coordinate is the distance between H and R that is:
And the coordinate Z is:
or
and explicit magnitude equations of the initial crosswise vector \({\overrightarrow{OR}}_{0}\) from the point O (0, 0, 0) to point R (R_{X}, R_{Y}, R_{Z}) are:
Four components of the lengthwise \(\overrightarrow{\mathrm{OA}}\) and the crosswise vector \(\overrightarrow{\mathrm{OR}}\) are:
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.
About this article
Cite this article
Cosola, S., Toti, P., PeñarrochaDiago, M. et al. Standardization of threedimensional pose of cylindrical implants from intraoral radiographs: a preliminary study. BMC Oral Health 21, 100 (2021). https://doi.org/10.1186/s12903021014489
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12903021014489
Keywords
 Dental informatics/bioinformatic
 Computer simulation
 Dental implant(s)
 Digital imaging/radiology
 Mathematical modeling