Advanced Search

Journal Navigation

Journal Home

Subscriptions

Archive

Contact Us

Table of Contents

CiteULike is a free service for managing and discovering scholarly references - click here to get started.

Sign In to gain access to subscriptions and/or personal tools.
Journal of Dental Research
This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to Saved Citations
Right arrow Download to citation manager
Right arrowRequest Permissions
Right arrow Request Reprints
Right arrow Add to My Marked Citations
Citing Articles
Right arrow Citing Articles via Google Scholar
Right arrow Citing Articles via Scopus
Google Scholar
Right arrow Articles by Tohnak, S.
Right arrow Articles by Crozier, S.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Tohnak, S.
Right arrow Articles by Crozier, S.
Social Bookmarking
 Add to CiteULike   Add to Complore   Add to Connotea   Add to Del.icio.us   Add to Digg   Add to Reddit   Add to Technorati   Add to Twitter  
What's this?

Biomaterials & Bioengineering

Synthesizing Dental Radiographs for Human Identification

S. Tohnak1, A.J.H. Mehnert1,*, M. Mahoney2 and S. Crozier1

1 School of Information Technology and Electrical Engineering, The University of Queensland, Brisbane QLD 4072, Australia; and
2 School of Dentistry, The University of Queensland, Brisbane QLD 4072, Australia

Correspondence: * corresponding author, mehnert{at}itee.uq.edu.au


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS & METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
The task of identifying human remains based on dental comparisons of post mortem (PM) and ante mortem (AM) radiographs is labor-intensive, subjective, and has several drawbacks, including: inherently poor image quality, difficulty matching the viewing angles in PM radiographs to those taken AM, and the fact that the state of the dental remains may entirely preclude the possibility of obtaining certain types of radiographs PM. The aim of the present study was to investigate the feasibility of using radiograph-like images reconstructed from PM x-ray computed tomography (CT) data to overcome the shortcomings of conventional radiographic comparison. Algorithms for computer synthesis of panoramic, periapical, and bitewing images are presented. The algorithms were evaluated with data from clinical examinations of two persons. The results demonstrate the efficacy of the CT-based approach and that, in comparison with conventional radiographs, the synthesized images exhibit minimal geometric distortion, reduced blurring, and reduced superimposition of oral structures.

Key Words: computed tomography • radiograph • Radon transform • dental comparison • forensic identification


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS & METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Dental comparison provides one of the best avenues for the forensic identification of human remains, because teeth, unlike other body parts, are highly resistant to decomposition and are not easily destroyed. Post mortem (PM) dental information is compared with ante mortem (AM) information in an attempt to establish the identity of the remains. Dental information typically includes x-ray-based radiographs—including periapical film and bitewing radiographs, and panoramic radiographs or orthopantomographs (OPGs)—and written notes and charts from patient records.

Dental comparison is based on characteristics of the teeth, jaws, and bucco-lingual structures, as well as on the characteristics of any restorations or dental work (Adams, 2003). The classic approach to dental identification involves manually summarizing dental records and radiographs. A PM summary can then be compared with AM summaries to establish a short list of possible identities. Manual comparisons of the candidate AM radiographs with the PM radiographs are necessary before a positive identification can be made.

The classic approach to dental comparison has several shortcomings. First, it is very labor-intensive (Kogon et al., 1974). Second, it is error-prone, because mistakes can be made when information is collected and summaries are compiled (Kogon et al., 1974). Third, it is very subjective (Kogon et al., 1974). Finally, and perhaps most significantly, the short list of possible matches is based on summary information rather than on features, e.g., tooth and root morphology, derived directly from the PM and AM radiographs. These shortcomings have motivated research into, and development of, computer-assisted methods to improve the accuracy and efficiency of dental identification. These fall into two categories: (i) methods that aim to improve the efficiency and accuracy of the classic approach; and (ii) methods that aim to extend the classic approach.

Software packages such as CAPMI (Lorton et al., 1988), WinID (McGivney, 2006), IDIS (Chomdej et al., 2006), and OdontoSearch (Adams, 2003) belong to the first category. These packages essentially permit the creation of AM/PM computer databases using standardized summary forms, and text-based searching of records. With regard to the second category, several approaches to content-based matching of digitized radiographs have been proposed (Jain et al., 2003; Jain and Chen, 2004; Chen and Jain, 2004; Said et al., 2004; Zhou and Abdel-Mottaleb, 2004, 2005; Mahoor and Abdel-Mottaleb, 2005; Nomir and Abdel-Mottaleb, 2005). These approaches utilize image processing and statistical pattern recognition techniques to perform automated or semi-automated matching of digitized radiographs based on tooth shape.

The majority of existing approaches to computer-assisted dental comparison have focused on the use of conventional x-ray radiographs. The reason for this is that conventional radiographs are the most abundant source of AM image data. However, this approach has several drawbacks, including: (i) reduced image quality in both AM and PM images, as a result of improper film processing and geometric distortion; (ii) difficulty discerning tooth structure because of overlapping adjacent teeth and/or other oral structures; (iii) difficulty matching views in corresponding PM and AM radiographs, because of differing projection geometries at the time of acquisition; and (iv) the fact that the quality of the remains may preclude the possibility of acquiring certain types of radiographs PM. An alternative approach, which overcomes these drawbacks, is to use PM x-ray computed tomography (CT) for matching against AM radiographs (Jackowski et al., 2006; Thali et al., 2006).

CT is an imaging modality that reveals the three-dimensional (3D) relationships between teeth and their supporting tissues in a series of two-dimensional (2D) digital images without geometric distortion. Unfortunately, CT is not an abundant source of AM dental image data. This is because of the higher cost and the higher dose of radiation to which a person is exposed compared with conventional radiography. Nevertheless, CT data can be easily acquired PM and, given that CT is a true 3D modality, used to generate 2D images that mimic conventional radiographs. Moreover, these images can be generated with arbitrary projection geometry, reduced blurring, and reduced superimposition of other oral structures. The aim of the present study was to investigate the feasibility of this approach. Algorithms for computer synthesis of panoramic, periapical, and bitewing radiograph-like images from CT data are presented. In addition, an evaluation of the algorithms, using real clinical data, is presented.


    MATERIALS & METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS & METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Overview of Conventional Dental Radiography
Periapical and bitewing radiographs are the most commonly acquired radiographs in clinical practice. These are referred to as ’intra-oral’ radiographs because the film is placed in the person’s mouth. They are used to assess dental caries. A periapical radiograph is used to examine individual teeth and the tissue surrounding the root apices (Whaites, 1996). Each film typically presents from 1 to 4 teeth. Bitewing radiographs are named after the original technique, which required the person to bite on a small wing attached to the intra-oral film. This technique is used to examine the crowns of both the premolar and molar teeth on one side of the jaw (Whaites, 1996). Three or four upper and lower tooth crowns are typically displayed in each film. Both periapical and bitewing radiographs aim to record a select range of teeth with the least magnification error. The ideal imaging geometry requires the film and the teeth to be parallel to one another, and for parallel x-ray beams to be transmitted perpendicular to the teeth and film (Whaites, 2002). The orientation of the person, film, and x-ray tube head are fixed during the exposure of the film. Consequently, given the limits in film positioning because of oral anatomy, in some cases the film may not be in the ideal position, leading to geometric distortion in the resulting image.

A panoramic radiograph, or OPG, is an extra-oral radiograph used to supplement bitewing and periapical radiographs. In panoramic radiography, both the x-ray tube head and the film are simultaneously moved around the person’s head during the exposure. The process of image formation is known as tomography. The x-ray source emits a narrow vertical beam, and as the system rotates, a 2D image is exposed on film. The method permits a single curved layer through the jaws to be imaged, while at the same time blurring out shadows of oral structures superimposed from other layers. This curved layer is known as the ‘focal trough’ or ‘image plane’. Manufacturers design the size and shape of the focal trough to accommodate the average jaw.

Overview of X-ray Computed Tomography
X-ray CT is a more recent, but less widely used, modality for dental imaging (Vannier et al., 1997). In contrast to conventional dental radiography, very sensitive crystal or gas detectors rather than film are used to acquire an image. A dental CT scan yields a sequence of cross-sectional images (slices) axially through the mandible and maxilla. In principle, a single axial slice image is acquired when the x-ray source and detectors are rotated around the person’s head, and successive slices are obtained by moving the person progressively through the scanner. Modern CT scanners use a spiral scanning method that permits data to be acquired as a volume of contiguous slices. Each CT slice image is a 2D matrix of grey values. Each element corresponds to a volume element (voxel). A voxel’s grey value represents the average attenuation for all of the x-rays that passed through that voxel.

Image Projection and the Radon Transform
A conventional dental radiograph fundamentally records the cumulative x-ray attenuation along each ray transmitted from the x-ray source to the film. For this reason, it is possible to mimic conventional dental radiographs by taking 2D projections through the 3D CT data. An important mathematical tool here is the 2D Radon transform (RT). The RT is the projection or line integral of a 2D function over the set of all possible lines (Deans, 1983). However, given that CT data are discrete, and, in particular, that each slice is a matrix (2D discrete function), a discrete version of the transform is of interest where integrals are replaced by finite pixel sums. The discrete Radon transform (DRT) is a mapping from a 2D discrete function to a set of 1D projections or sums (Kingston and Svalbe, 2003).

Algorithms for Reconstructing Radiograph-like Images from CT Data
During a dental CT scan, the person’s maxillary and mandibular teeth are usually kept apart by a separating device of low radiographic density, characteristically an empty syringe. In addition, to help reduce metal artifacts from dental restorations, the imaging of the mandible and maxilla may be performed in two separate, but sequential, scans, and in each scan the angle between the axial plane and the occlusal plane may be different. Consequently, it is necessary to align the mandible and maxilla volumes spatially before attempting to reconstruct radiograph-like images from the CT data (Table 1Go, Steps 1 to 6).


View this table:
[in this window]
[in a new window]

 
Table 1. Algorithm for Reconstructing Periapical- and Bitewing-like Images
 
Our proposed algorithm for reconstructing periapical- and bitewing-like images from CT data (Table 1Go) requires that the user manually define a 3D region of interest (rectangular prism) around the tooth or teeth that are to appear in the projection.

Our proposed algorithm for reconstructing OPG-like images from CT data (Table 2Go) requires that the user manually define a mask on the maximum intensity projection through the CT data in the axial plane. The aim is to define a mask that isolates the voxels of the dental arch, and whose medial axis passes through the center of the dental arch. The algorithm tracks along this medial axis line, projecting through each axial slice of the CT data normal to the tangent. This strategy aims to minimize magnification of the teeth. Projections for the left and right sides of the jaw are performed in isolation, to avoid superimposition of dental structures.


View this table:
[in this window]
[in a new window]

 
Table 2. Algorithm for Reconstructing an OPG-Like Image
 
Evaluation Using Routine Clinical Data
Retrospective anonymized image data for two persons were obtained from a local dental imaging practice. The data were chosen to include a person for whom a spiral CT dataset and an OPG film had been acquired (Person A), and a person with no metallic restorations (Person B). Approval from the Human Research Ethics Committee of the University of Queensland was obtained for this study.

The CT data for Person A consisted of 109 axial slices through the mandible (voxel size of 0.24 mm x 0.24 mm x 0.50 mm) and 115 axial slices through the maxilla (voxel size of 0.18 mm x 0.18 mm x 0.50 mm). The CT data for Person B consisted of 113 axial slices through the mandible (voxel size of 0.22 mm x 0.22 mm x 0.50 mm).

The reconstruction algorithms were implemented in MATLAB 7.1 (The MathWorks Inc., Natick, MA, USA, 2005). The algorithm for reconstructing periapical- and bitewing-like images was applied to the CT data for Person A with a variety of user-defined regions of interest. The algorithm for reconstructing an OPG-like image was applied to the combined mandible and maxilla CT data for Person A, and to the mandible-only CT data for Person B.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS & METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
The reconstructed periapical- and bitewing-like images (Fig. 2aGo) revealed tooth composition, dental caries, dental restorations, and, in the case of the periapical images, also the root shape and supporting bone in the anterior teeth. Negatively, the reconstructed images contained intensity artifacts as the result of metal artifacts in the original CT data, originating from metallic dental restorations. The presence of these artifacts can make interpretation difficult (e.g., Fig. 2aGo, bottom-left).


Figure 2
View larger version (77K):
[in this window]
[in a new window]

 
Figure 2. Radiograph-like images generated by the proposed reconstruction algorithms. (a) Examples of reconstructions of periapical images (top two rows) and bitewing images (bottom row) from the CT data for Person A. (b) Reconstructed OPG-like image from the CT data for Person A. (c) Real OPG for Person A. (d) Reconstructed OPG-like image from the CT data for Person B. (e) Volume rendering of the spiral CT data for Person B.

 
In the case of Person A, a visual comparison between the reconstructed OPG (Fig. 2bGo) and the real OPG (Fig. 2cGo) revealed that the former exhibited reduced geometric distortion, reduced blurring, and reduced overlapping of dental structures. Of particular note was that in the reconstructed OPG: (i) there was less horizontal magnification in the premolar area; (ii) there was minimal overlapping of the posterior teeth; (iii) the root curvature of the permanent lower molars was preserved; (iv) the dental restorations in both the crown and root portions of restored teeth could be seen clearly; and (v) there was no superimposition of other oral tissues. Negatively, the reconstructed OPG contained intensity artifacts stemming from metal artifacts manifest in the CT data. Similarly, in the case of Person B, a visual comparison between the reconstructed OPG (Fig. 2dGo) and the CT volume rendering (Fig. 2eGo) revealed that the former exhibited little geometric distortion, blurring, and overlapping of dental structures. The reconstruction, in this case, was free from intensity artifacts, because the person had no metallic restorations.


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS & METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
The results demonstrated that it is possible to reconstruct 2D images from spiral CT data that mimic conventional periapical-, bitewing-, and panoramic radiographs. In addition, they demonstrated the efficacy of the proposed computer algorithms for reconstructing these images. Collectively, the results demonstrated that it is indeed feasible to use these images in lieu of conventional radiographs for matching against real AM radiographs, to establish the identity of human remains. This approach to dental comparison opens up the possibility for automated, objective comparison by computer. The approach overcomes several limitations of conventional radiographic comparison, including: inherently poor image quality, difficulty matching the viewing angles in PM radiographs to those AM, and the fact that the state of the dental remains may entirely preclude the possibility of obtaining certain types of radiographs PM. Metallic dental restorations lead to streak-like intensity artifacts in the CT image data. These occur because x-rays that pass through metal objects are much more attenuated than those that pass through soft tissues and bone. These artifacts, in turn, lead to intensity artifacts in the reconstructed radiograph-like images produced by our algorithms. They occur, in particular, on and around teeth with restorations. They affect not only the visual quality of the reconstructions, but also their interpretability. For this reason, we intend to incorporate metal artifact reduction as a pre-processing step in future versions of our reconstruction algorithms.

This research represents the first step in a larger initiative to develop a computerized image analysis system to assist the forensic dentist with the task of identifying human remains. In particular, the system will permit radiograph-like images to be reconstructed from PM CT data dynamically ("on the fly"), and will enable these to be matched against AM digitized conventional radiographs. Matching will be performed on the basis of quantitative features extracted from these images that characterize both the morphology of the teeth and jaws and their spatial interrelationships.


Figure 1
View larger version (47K):
[in this window]
[in a new window]

 
Figure 1. Illustration of selected steps in the proposed reconstruction algorithms. First reconstruction algorithm: (a) location of the occlusal plane in the sagittal maximum intensity projections (MIPs) of the interpolated mandible and maxilla CT data; (b) definition of the rectangular mask around the teeth of interest on the axial MIP of the combined CT data; and (c) extension of the mask in (b) to define a rectangular prism needed to reconstruct the desired periapical- or bitewing-like image projection. Second reconstruction algorithm: (d) axial MIP of the combined CT data with the jaw mask superimposed (dashed line); (e) medial axis of the jaw mask superimposed with the midline (dashed line), and the tangent and normal (dotted lines) at one of the points on the medial axis; and (f) distribution of tangent angles along the medial axis, proceeding from left to right, with a polynomial fit superimposed.

 

    ACKNOWLEDGMENTS
 
Ms. Tohnak acknowledges financial support from Naresuan University and the Commission on Higher Education, Ministry of Education, Thailand.

Received for publication October 17, 2006. Revision received June 28, 2007. Accepted for publication July 9, 2007.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS & METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 

  • Adams BJ (2003). Establishing personal identification based on specific patterns of missing, filled, and unrestored teeth. J Forensic Sci 48(Pt 3):487–496.[Medline] [Order article via Infotrieve]
  • Chen H, Jain AK (2004). Tooth contour extraction for matching dental radiographs. In: Proceedings of the 17th International Conference on Pattern Recognition, ICPR 2004, Aug 23–26, 2004. Kittler J, Petrou M, Nixon M, editors. Cambridge, UK: IEEE Computer Society, pp. 522–525.
  • Chomdej T, Pankaow W, Choychumroon S (2006). Intelligent dental identification system (IDIS) in forensic medicine. Forensic Sci Int 158(Pt 1):27–38.[Medline] [Order article via Infotrieve]
  • Deans SR (1983). The Radon transform and some of its applications. New York: Wiley.
  • Jackowski C, Aghayev E, Sonnenschein M, Dirnhofer R, Thali MJ (2006). Maximum intensity projection of cranial computed tomography data for dental identification. Int J Legal Med 120(Pt 3):165–167.[CrossRef][Medline] [Order article via Infotrieve]
  • Jain AK, Chen H (2004). Matching of dental x-ray images for human identification. Pattern Recognition 37(Pt 7):1519–1532.
  • Jain AK, Chen H, Minut S (2003). Dental biometrics: human identification using dental radiographs. In: Audio- and Video-Based Biometric Person Authentication. Proceedings, 4th International Conference, AVBPA 2003, Guilford, UK, Jun 9–11, 2003. Kittler J, Nixon MS, editors. Guilford, UK: Springer, pp. 429–437.
  • Kingston A, Svalbe I (2003). Mapping between digital and continuous projections via the discrete Radon transform in Fourier space. In: Proceedings of the Seventh International Conference on Digital Image Computing: Techniques and Applications, DICTA 2003, Dec 10–12, 2003. Sun C, Talbot H, Ourselin S, Adriaansen T, editors. Sydney, NSW, Australia: CSIRO Publishing, pp. 263–272.
  • Kogon SL, Petersen KB, Locke JW, Petersen NO, Ball RG (1974). A computerized aid to dental identification in mass disasters. Forensic Sci 3:151–162.[Medline] [Order article via Infotrieve]
  • Lam L, Lee S-W, Suen CY (1992). Thinning methodologies—a comprehensive survey. IEEE Trans Pattern Anal Mach Intell 14(Pt 9):869–885.[CrossRef]
  • Lorton L, Rethman M, Frideman R (1988). The computer-assisted postmortem identification (CAPMI) system: a computer-based identification program. J Forensic Sci 33(Pt 4):977–984.[Medline] [Order article via Infotrieve]
  • Mahoor MH, Abdel-Mottaleb M (2005). Classification and numbering of teeth in dental bitewing images. Pattern Recognition 38(Pt 4):577–586.
  • McGivney J (2006). WinID3 dental identification system. www.winid.com
  • Nomir O, Abdel-Mottaleb M (2005). A system for human identification from x-ray dental radiographs. Pattern Recognition 38(Pt 8):1295–1305.
  • Said E, Fahmy GF, Nassar D, Ammar H (2004). Dental x-ray image segmentation. In: Proceedings of the SPIE-The International Society for Optical Engineering, Biometric Technology for Human Identification, Apr 12–13, 2004. Jain AK, Ratha NK, editors. Orlando, FL, USA: SPIE-International Society for Optical Engineering, pp. 409–417.
  • Thali MJ, Markwalder T, Jackowski C, Sonnenschein M, Dirnhofer R (2006). Dental CT imaging as a screening tool for dental profiling: advantages and limitations. J Forensic Sci 51:113–119.[CrossRef][Medline] [Order article via Infotrieve]
  • The MathWorks Inc. (2005). MATLAB (computer program). Natick, MA, USA: The MathWorks, Inc.
  • Vannier MW, Hildebolt CF, Conover G, Knapp RH, Yokoyama-Crothers N, Wang G (1997). Three-dimensional dental imaging by spiral CT: a progress report. Oral Surg Oral Med Oral Pathol Oral Radiol Endod 84:561–570.[CrossRef][Medline] [Order article via Infotrieve]
  • Whaites E (1996). Essentials of dental radiography and radiology. 2nd ed. New York: Churchill Livingstone.
  • Whaites E (2002). Essentials of dental radiography and radiology. 3rd ed. New York: Churchill Livingstone.
  • Zhou J, Abdel-Mottaleb M (2004). Automatic human identification based on dental x-ray images. In: Proceedings of SPIE-The International Society for Optical Engineering, Biometric Technology for Human Identification, Apr 12–13, 2004. Vol. 5404. Jain AK, Ratha NK, editors. Orlando, FL, USA: SPIE-International Society for Optical Engineering, pp. 373–380.
  • Zhou J, Abdel-Mottaleb M (2005). A content-based system for human identification based on bitewing dental x-ray images. Pattern Recognition 38(Pt 11):2132–2142.

Journal of Dental Research, Vol. 86, No. 11, 1057-1062 (2007)
DOI: 10.1177/154405910708601107


Add to CiteULike CiteULike   Add to Complore Complore   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us   Add to Digg Digg   Add to Reddit Reddit   Add to Technorati Technorati   Add to Twitter Twitter    What's this?



This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to Saved Citations
Right arrow Download to citation manager
Right arrowRequest Permissions
Right arrow Request Reprints
Right arrow Add to My Marked Citations
Citing Articles
Right arrow Citing Articles via Google Scholar
Right arrow Citing Articles via Scopus
Google Scholar
Right arrow Articles by Tohnak, S.
Right arrow Articles by Crozier, S.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Tohnak, S.
Right arrow Articles by Crozier, S.
Social Bookmarking
 Add to CiteULike   Add to Complore   Add to Connotea   Add to Del.icio.us   Add to Digg   Add to Reddit   Add to Technorati   Add to Twitter  
What's this?