This clinical retrospective investigation was conducted from January 2022 to December 2023 in XXX Hospital, XXX University. For this investigation, the sample comprised two imaging modalities: Orthopantomography (PAN, n1 = 290) and with corresponding Cone Beam Computed Tomography (CBCT, n2 =30). The research was in conformity with the guidance of the World Medical Association Helsinki Declaration for biomedical research involving human subjects and the imaging data involved was approved by the institutional review boards (protocol number: XXXXXXX-X-XXXX–XX). Informed consent was waived because this retrospective study was based on existing data. All LCRs were anonymized. The marking and the measurement of key points on two types of images were completed by two senior dentists after unified training.
2.1 Clinical diagnosis and eligibility criteria
This retrospective study utilized clinical investigation records of CAL that were assessed via direct intraoral examination. The radiograph either PAN or CBCT captured with Guidelines on radiology standards for primary dental care.
The inclusion criteria for this retrospective study required that patients have permanent teeth, as per the FDI tooth numbering system, with at least seven teeth in each quadrant and clear presentation of both teeth and periodontal tissues.
The exclusion criteria targeted severe maxillofacial anomalies such as facial or jaw deformities, and patients with a history of dental implants or removable dentures. Also excluded were cases involving significant gingival swelling, unclear tooth or tissue boundaries, abnormal tooth morphology, severe dentin defects. Additional exclusions covered ABL unrelated to periodontal diseases, traumatic gingival recession, cervical caries, distal clinical attachment level retraction, pus discharge from dental pulp, and vertical root fractures.
For reliability testing specifically, the study selected 30 patients who had undergone CBCT within a three-month window surrounding their Panoramic radiography acquisition.
2.2 Panoramic radiography measurements and annotation
The Labelme (version 5.2.1, Anaconda, Austin, TX) annotation tool was used to label panoramic images of 8120 teeth from 290 patients, as displayed in Table 1. It was imperative to assign a distinct identifier to each physiological and anatomical feature of the teeth in question. In accordance with the guidelines set forth by the Fédération Dentaire Internationale (FDI), orientation descriptors included the distal (D), mesial (M), occlusal (O), and gingival (G) surfaces of the teeth on 2D diagnostic images. To streamline the calculation of anatomical teeth height and their attachments, we used specific landmarks depicted in Fig. 2, abbreviated according to these teeth surfaces and regions. The labeling process was halted at the distal aspect of the second molars. In situations where the third molar was missing and recent extraction records confirmed its removal, the analysis exclude the distal aspect of the second molar. This exclusion was taken because of the common presence of localized ABL following an extraction. This exclusion was particularly relevant when annotating the alveolar crest (C) and cemento-enamel junction (J) groups. Conversely, the pulp (T) and root (R) groups extended up to the second molar in each quadrant.
Table 1: Artificial discrimination of the ABL and detailed grouping according to the segmented regions a
|
Characteristic
|
n(%)
|
Age, median
|
|
13~77, 39
|
Sex, n (%)
|
|
|
Male
|
|
116, 40.00%
|
Female
|
|
174, 60.00%
|
Alveolar bone non-resorption group (Normal)
|
131, 45.17%
|
Alveolar bone resorption group a (Abnormal)
|
159, 54.83%
|
Maxillary regions diagnosed ABL b
|
137, 47.24%
|
Maxillary anterior region
|
b-F c
|
121, 41.73%
|
Maxillary posterior region
|
|
129, 44.48%
|
|
a-R
|
123, 42.41%
|
|
c-L
|
112, 38.62%
|
Maxillary regions diagnosed ABL
|
139, 47.93%
|
Mandibular anterior region
|
e-F
|
137, 47.24%
|
Mandibular posterior region
|
|
120, 41.38%
|
|
d-L
|
116, 40.00%
|
|
f-R
|
101, 34.83%
|
a. Clinical examination for periodontitis / alveolar bone resorption and overall diagnosis: Based on the retained medical records, confirm whether the patient is diagnosed as periodontitis (or periodontitis has been cured), and record the site of alveolar bone resorption (the smallest unit in the research is the segment, for example, maxillary anterior area, mandibular posterior area)
b. Clinical examination for segmented diagnosis
c. F: front, R: right, L-left
|
2.3 CBCT measurements and annotation
The CBCT data, with a slice thickness of 0.25 mm, were browsed and processed using Mimics (version 21.0, Materialise, Belgium) set to a slice thickness and interval of 0.2 mm. Under the guidance of a professional periodontologist, two experienced dentists conducted annotations and measurements following a clockwise chart divided into six regions according to the distal surface of the canine: starting with the upper right posterior (a) and anterior (b), moving to the upper left (c), then to the lower right (f), the lower anterior (e), and finally the lower left (d).
Each tooth was initially examined by establishing a cross-sectional view, which was then adjusted to the optimal plane for analysis. The selected image aligned the Y-axis with the long axis of the tooth under examination in both the sagittal and coronal planes, while ensuring that the X-axis of the sagittal plane ran parallel to the tooth's mesiodistal orientation. By utilizing the coronal plane for observation, the ABL for each tooth was meticulously quantified. This process revealed discernible differences between the anterior and posterior teeth, as depicted in Fig. 3. The ABL for each tooth was measured, also revealing notable differences between anterior and posterior teeth. To ensure consistent measurements across the same patient, each measurement plane was refined to create three distinct sections, with the central section selected to mitigate edge blurring that could compromise accuracy. The root length, ABL, and crown length (in millimeters) were each measured four times on both the mesiodistal and buccal-lingual aspects of the teeth. Across all four planes, the root and crown lengths were determined by averaging the two higher values (for posterior teeth, the index derived from the buccal-lingual side was considered acceptable), while ABL was quantified using two absolute linear measurements: the mesiodistal average and the buccal-lingual average, as illustrated in Fig. 3. The absolute linear distances of the buccal-lingual or mesiodistal ABL were compared with the average root length, and the resulting ratios were averaged to calculate the relative distance of alveolar bone resorption.
It is crucial to note that, when the molars with multiple long roots, C-shaped roots, and curved roots, each root was considered as a singular unit. The length of each root was individually measured, and corresponding absolute and relative ABL values were determined. Ultimately, for the target molar, the average ABL of two or three roots was recorded to provide a comprehensive assessment of the tooth's periodontal condition.
2.4 Model architecture and Computer processing
a. ABL analysis in PAN-POL model
The PAN-POL model was processed by PyCharm 2023.3.3(JetBrains, Prague, Czech Republic). For all PANs that completed the labeling of landmarks, the polynomial function “polyfit” was employed to draw a curve that accurately fit all labeled landmarks within each segment.
Considering the morphological characteristics in clinical diagnosis, the curve fitting method is categorized into full-mouth and sectional curve fitting, with the fitting function model limited to a quadratic function as shown in Fig. 4. Each curve was meticulously drawn through points labeled in various groups as described above: the cuspid tip (group T), the cemento-enamel junction (group J), the alveolar crest (group C), and the root apex (group R). By considering the connecting line on PAN, we aimed to minimize the overall fitting error of the coefficient of determination (R2). This approach enhances the accuracy and intuitive visualization of both full-mouth and localized ABL. For full-mouth fitting, tooth point coordinate data and a polynomial order are inputted into the polyfit function. The least squares method is then used to calculate the polynomial coefficients, ensuring the polynomial curve aligns as closely as possible with these data points. From this, we derive three inherent curves for the cuspid tip (rT), cementoenamel junction (rJ), and root (rR), along with one ABL description curve (rC) relatively. The final index of these curves is displayed in Table 2. The abscissa at the end of each curve does not necessarily align; in such cases, the remaining three curves are truncated by the abscissa closest to the median line on each curve. Subsequently, the average distance between the curves is calculated.
Table 2 The marker processing on PAN
|
Indicators
|
Paraphrase
|
Measurement reference
|
Absolute linear distance
|
|
|
h1
|
The vertical distance between the rT and the rJ in the corresponding section (Crown length).
|
Fitting curve rJ
|
h2
|
The vertical distance between the rR and the rJ in the corresponding section (Root length).
|
h3
|
The vertical distance between the rc and the rJ in the corresponding section (ABL).
|
Relative linear distance
|
|
|
n0
|
Crown/Root Ratio, n0=h1/h2
|
The vertical distance between the rR and the rJ in the corresponding section (Root length).
|
n1
|
Relative ABL, n1=h3/h2
|
Goodness of fit R²
|
In the case of the same abscissa, the consistency between the data points and the approximate fitting curves
|
According to the arrangement of the six quadrants, which include the maxillary and mandibular posterior and anterior regions, segmented curve fitting was employed to accurately reflect local alveolar bone resorption. The indices derived showed no differences to those from the full-mouth fit. Additionally, the values n0 and n1 were calculated through computer processing of the marked labels.
b. CBCT ABL analysis
Each absolute linear distance index on CBCT (h1~3) covered multiple tooth positions within each tooth area, with each tooth position encompassing several measurement sites—each site corresponding to one tooth root in both the buccal-lingual and mesiodistal directions. Patient data collection was directly managed using Microsoft Excel. For each dental region, the mean value of the relative linear distance indices, derived from various measurement planes as noted earlier at each tooth position, was calculated. These mean values were then aggregated with corresponding indices from other teeth within the same region. The software then calculated the average of these combined values, providing a comprehensive assessment for that particular dental region.
Similarly, the full-mouth assessment, which can be viewed as an expanded local region, employed the same methodology. This systematic and standardized approach to processing the required indicators ensures a thorough analysis of ABL across the entire dentition. It allows for a detailed evaluation of periodontal health and the effectiveness of various treatment modalities on bone structure.
c. Reliability and validity test
To assess the consistency between different observers, the annotation work was carried out by two senior dentists. Initially, CBCT data from four patients not involved in the study—two with alveolar bone resorption and two without—were measured. Additionally, the PAN data from 39 not involved randomly selected patients were marked and the data output according to established rules to ensure the randomness of image selection and uniformity of annotation standards. The measurements of ABL (n1) and the crown-root ratio (n0) by the two dentists were then compared.
d. Intraclass correlation efficient test
We employed the Intraclass Correlation Coefficient (ICC) to analyze the two sets of measured or calculated values, using the absolute agreement definition in the two-way random model. The results indicated an ICC range of 0.69 to 1.00 (P<0.05, Supplement 1 Table 3). Except for the ICC of 0.69 for the ABL in c-segment on the PAN, all other ICC values were greater than 0.75, demonstrating a high level of agreement between the observers for most measurements.
2.5 Statistical analysis
Statistical analysis was conducted using data from each patient record, which were simultaneously entered into Microsoft® Excel® 2019 MSO (Microsoft Corp., Redmond, WA) and analyzed using IBM SPSS Statistics 24.0 (IBM, Chicago, IL, USA).
To demonstrate the distinction between the two sets of model-processed data, we compared the full-mouth and segmented ABL(n1) between the ABL-abnormal and ABL-normal groups, which had been initially classified based on the presence of any ABL as well as ABL in specific segmented regions. And then the proportion indices from 290 manually labeled panoramic images were collected and subjected to a normality test. For groups that demonstrated normal distribution and homogeneity of variance, a single-sample t-test was conducted. Conversely, nonparametric tests were applied to groups exhibiting either non-normal distributions or non-homogeneous variances. After performing a Test of Normality, we used a single-sample t-test or a Mann-Whitney Test to analyze differences in ABL between the ABL-normal and ABL-abnormal groups, as well as the difference between two groups.
Following this, we compared ABL and crown-root ratio (n0) values from CBCT measurements with those from the curve-fitted PAN-POL model in 30 patients, aiming to verify the model’s reliability. The medians and interquartile ranges were reported, facilitating the evaluation of reference ranges for normal and abnormal values. Subsequent to a further Test of Normality, both the Intraclass Correlation Coefficient (ICC) and the Wilcoxon Signed Ranks Test were employed to assess the reliability of the PAN-POL model in comparison to CBCT measurements.
In addition, R2 was utilized to gauge how well the curves fit the points and to determine if the curves effectively represented the key point positions for each group. Validation of the model relied on clinical records and was rigorously tested using the R method to ensure its accuracy.