Next Article in Journal
Enabling SDN in VANETs: What is the Impact on Security?
Next Article in Special Issue
Evaluation of Laser-Assisted Trans-Nail Drug Delivery with Optical Coherence Tomography
Previous Article in Journal
Cramer-Rao Lower Bound Evaluation for Linear Frequency Modulation Based Active Radar Networks Operating in a Rice Fading Environment
Previous Article in Special Issue
Neuron Stimulation Device Integrated with Silicon Nanowire-Based Photodetection Circuit on a Flexible Substrate
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bio-Photonic Detection and Quantitative Evaluation Method for the Progression of Dental Caries Using Optical Frequency-Domain Imaging Method

1
School of Electronics Engineering, College of IT Engineering, Kyungpook National University, 80, Daehak-ro, Buk-gu, Daegu 41566, Korea
2
Eaton-Peabody Laboratories, Massachusetts Eye and Ear Infirmary(MEEI), 243, Charles Street, Boston, MA 02114, USA
3
Department of Otology and Laryngology, Harvard Medical School, 243, Charles Street, Boston, MA 02114, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Sensors 2016, 16(12), 2076; https://doi.org/10.3390/s16122076
Submission received: 28 September 2016 / Revised: 26 November 2016 / Accepted: 2 December 2016 / Published: 6 December 2016

Abstract

:
The initial detection of dental caries is an essential biomedical requirement to barricade the progression of caries and tooth demineralization. The objective of this study is to introduce an optical frequency-domain imaging technique based quantitative evaluation method to calculate the volume and thickness of enamel residual, and a quantification method was developed to evaluate the total intensity fluctuation in depth direction owing to carious lesions, which can be favorable to identify the progression of dental caries in advance. The cross-sectional images of the ex vivo tooth samples were acquired using 1.3 μm spectral domain optical coherence tomography system (SD-OCT). Moreover, the advantages of the proposed method over the conventional dental inspection methods were compared to highlight the potential capability of OCT. As a consequence, the threshold parameters obtained through the developed method can be used as an efficient investigating technique for the initial detection of demineralization.

1. Introduction

Precise structural imaging and quantitative evaluation are critical for the diagnosis of dental caries and research to prevent the formation of initial demineralized regions [1,2]. The depth imaging of enamel, dentin, cavities, pits, and fissures with a high-resolution is particularly important for the study of anatomical and pathological changes of the dental structure [3]. Optical coherence tomography (OCT) is a rapidly advancing optical frequency-domain imaging modality, which can provide non-invasive high-resolution cross-sectional images of dental tissues and various biological tissues [4]. This near-infrared (NIR) biomedical imaging method provides images with high axial and lateral resolutions (i.e., below 8 μm and 15 μm, respectively) [5,6], and, furthermore, OCT has been widely used in different medical applications such as ophthalmology [7,8], dermatology [9], and otolaryngology [10,11]. The methods currently in use for the detection of dental caries such as radiography, microradiography, and X-rays do not provide sufficient resolution, sensitivity, and contrast compared to OCT. Radiography is the most frequently applied inspection method in dentistry with a resolution of 50 μm, which is comparatively lower than the resolution of OCT. Furthermore, radiography is not quantitative, and it is relatively difficult to apply for the initial detection of dental caries [12,13]. Microradiography is another method that can be used to analyze caries quantitatively, but it is hard to apply this method for clinical applications [14]. Other conventional diagnostic methods such as infrared (IR) imaging, dental explorer, and visual inspection are unable to provide more accurate cross-sectional images [15]. Several research groups have determined the mineral loss and the depth of enamel caries using a histology analysis method called transversal microradiography (TMR). However, due to the requirement of a thin sectioning process, applications of the method in dentistry have been scarce [16]. The main drawback of this method is that caries can be detected only at a relatively advanced stage when remineralization is no longer possible, and due to the incapability of obtaining precise quantitative measurements, it is hard to barricade the progression of caries. Thus, owing to the non-invasive and non-destructive imaging capability, and the capability of acquiring precise quantifications such as accurate thickness and volumetric measurements, OCT has gained a significant demand in the medical field as an early diagnosis method. Although OCT has been extensively used as a powerful dental imaging technique, quantitative evaluation of enamel thickness variation, depth dependent intensity fluctuation, and volumetric analysis of enamel residual has not been broadly studied for the initial diagnosis of demineralization.
Especially in dentistry, OCT has been used to produce longitudinal images of dental tissues and caries of an orientation similar to that of the B-scan ultrasound images [17,18]. In these studies, a reduction in enamel reflectivity was observed in areas of dental caries [19]. It is considered that the decrease in reflectivity during demineralization is related to the amount of mineral loss. Few studies have demonstrated that there is two to threefold increase in the scattering coefficient at a wavelength of 1.3 μm [20]. Furthermore, a polarization sensitive OCT (PS-OCT) endoscopic system using a swept source has been implemented as a compact system for dentistry application [21]. Similarly, surface demineralization can be detected using linearly polarized light and measured backscattered signal in two orthogonal axes [22]. Some other OCT techniques have been applied to diagnose dental caries as a result of changes in the optical properties of enamel after undergoing demineralization [23]. The obtained images were quantitatively evaluated by the identification of structures, dimensions, and properties [24]. Moreover, in several review reports and research studies, OCT based dental experiments were demonstrated to verify the stronger optical backscattering signals acquired from the demineralized enamel regions, oral tissue images, caries, periodontal diseases, and oral cancers [25,26]. In addition, infrared light with long wavelengths were used in OCT for clinical applications owing to the high depth penetration [27].
In this paper, we performed an initial ex vivo study using a 1.3 μm wavelength laser utilized spectral domain OCT (SD-OCT) system to introduce a quantitative method to calculate the thickness and volume of remaining enamel region (enamel residual). Furthermore, we developed an algorithm to analyze the total intensity fluctuation in depth direction of OCT images, which can be useful to identify the progression of initial caries. As a result, the proposed quantification can be implemented to identify the reduction of enamel region along with the progression of the demineralization. Moreover, the capability of our SD-OCT system to perform the proposed method was validated by obtaining images of ex vivo caries with a high resolution and a high depth penetration.

2. Materials and Methods

2.1. Optical Frequency Domain Imaging (OFDI) Technique

The implemented optical frequency domain imaging technique was a customized 1.3 μm SD-OCT system. The speed of the SD-OCT system was 120 frames/s when the image size was 1024 × 500 pixels, and the average output power of the system was 16 mW. In Figure 1, the broadband light source that was used for light emission is a superluminescent diode (SLED) (Denselight Semiconductors, Singapore) with 1.3 μm central wavelength and 135 nm bandwidth. The axial resolution of the system was 6 μm (in air) and 3.61 μm (in tissue). The transverse resolution of the system was 25 μm. The detector was a 14-bit complementary metal-oxide semiconductor (CMOS) line scan camera (SU-1024LDM Compact; Goodrich, Charlotte, NC, USA) with 1024 pixels. A 50:50 optical fiber coupler was used to split the broadband light beam into the sample and reference arms. A galvanoscanner (GVS002, Thorlabs, Newton, NJ, USA) connected to the sample arm was used to scan the tooth samples. All the samples were scanned with a sufficient cross-sectional scanning range of 1 mm × 1 mm × 1 mm dimensions. A compact spectrometer was designed and contained a collimator, a diffraction grating, an achromatic doublet lens, and a line scan camera. The spectrometer was calibrated to compensate the distortion of the point spread function (PSF), and to improve the signal-to-noise ratio (SNR) up to 110 dB using previous literature reports [28,29]. Further details about the system configuration can be found in Table 1.

2.2. Specimen Preparation

For the proposed preliminary study, four types of ex vivo tooth samples, including partially demineralized canine tooth sample, partially demineralized pre-molar tooth sample, partially demineralized molar tooth sample, completely demineralized (carious) molar tooth sample, and in vivo healthy molar tooth sample were involved in the experiment. All tooth specimens were examined in patients before extraction. The experimented ex vivo tooth specimens were extracted after performing early childhood caries (ECC) surgeries for four orthodontic patients at different age groups (10–12 years) of the dental clinic of the Faculty of Dentistry, Kyungpook National University, Daegu, Korea. The details of the experimented volunteers and tooth specimens are illustrated in Table 2. Prior to the OCT inspection, all ex vivo tooth specimens were preserved in sterile filtered de-ionized water solution for 24 h at 30 °C to eliminate any possible superficial enamel cracks and maintain a standard smooth surface after the extraction. The experiments were performed in accordance with the guidelines of the Institutional Animal Care and Use Committee of Kyungpook National University (Daegu, Korea) and approvals from the human ethics committees of the Institute for Bio-diagnostics, Kyungpook National University (Daegu, Korea) were obtained prior to sample collection.

2.3. Intensity Fluctuation Analysis

To evaluate the proposed method precisely, the obtained cross-sectional images were involved in an amplitude scan (A-scan) depth profile analysis to verify the microstructural comparison between healthy, partially demineralized, and carious molar tooth samples. For the A-scan profile analysis, a software-based program was coded using Matlab (Mathworks, Natick, MA, USA) to search the intensity peaks in the depth direction. The acquired 2D OCT image was loaded and a peak search algorithm-based cropped window with 15 intensity signals (A-scan lines) was applied. The developed algorithm detects the maximum intensity in each individual A-scan line to search the peak position, and all the peak positions in all 15 A-scan lines were rearranged while matching the peak intensity index in the A-scans to flatten the region of interest. Owing to the non-flattened region of interest, the maximum intensity index positions vary, and therefore, the index positions with higher intensity values should be rearranged and matched linearly to obtain a flattened image. Finally, all the rearranged and flattened A-scan lines were summed up, averaged, and normalized to obtain a single A-scan depth profile of the region of interest. The applied refractive index of a tooth structure, which affects the depth scale of 2D OCT images was 1.63 [30,31]. Moreover, we performed an additional quantification method to analyze the total intensity fluctuation in deep microstructures according to imaging depth. Thus, an additional automated program was coded using Matlab to analyze the total pixel intensity of each depth range of demodulated 2D OCT images. The analysis was performed for the entire visible depth range of 1 mm, and the total pixel intensity was evaluated for each 250 μm depth range of the cross-sectional image to identify the depth dependent total intensity fluctuation of each tooth specimen category. Then, the entire total pixel intensity of each depth range was summed and averaged for each 2D OCT image. This study was a preliminary observational study, and the data analysis was primarily descriptive. A continuous variation of the optical laser source power was observed, which was compensated afterward. Due to the instability of the laser optical power, the entire intensity of 2D OCT images was compensated by multiplying ±5%.

2.4. Volumetric Analysis

The volumetric measurements of enamel residual was obtained using 2D OCT image based 3D OCT volumetric images by implementing a pixel intensity based automated calculation method, which was performed through a software program coded in Matlab. Figure 2a shows the volumetric calculation algorithm of enamel residual of a single 2D OCT image along with the obtained 3D OCT volumetric image containing 500 2D OCT images (Figure 2c). In the developed algorithm based program, demodulated raw data was loaded and the intensity of the entire cross-sectional region was analyzed. An image window was applied to select the enamel residual region. For the precise selection (filter the region of interest) of the enamel residual region, we approximated the pixel intensity difference between the enamel residual region and other cross-sectional regions, and provided a pre-determined intensity threshold range for a separate evaluation of the enamel residual region. The entire cross-sectional intensity varies from 0 to 255, and the pre-determined enamel threshold (TH(en)) range can be expressed as, 45 ≤ TH(en) ≤ 255. The selected threshold range contains the intensity range of the enamel region and excludes the dark black region of the selected image window. The area of a single pixel, which belongs to the selected image window in Figure 2b can be expressed as
l x × l y = A p i x ,
where lx is the pixel size in the x-direction, ly is the pixel size in the y-direction, and Apix is the area of a single pixel. The pixel sizes in the x-, y-, and z-directions can be expressed as
l x = l y = 12   μ m ,   and   l z = 7   μ m .
The total number of z-direction pixels (B-mode images) is 500 owing to the composition of the 3D OCT image. Therefore, the enamel residual volume (volume of the remaining enamel) of 500 2D OCT images (3D volumetric image) can be calculated as
( N 1 × l x × l y ) × l z + ( N 2 × l x × l y ) × l z + ( N 3 × l x × l y ) × l z + ( N 500 × l x × l y ) × l z = V t o t ,
where Ni (i = 1, 2, 3, …, 500) is the number of image window pixels in each respective window, which satisfies the pre-determined threshold value range, and Vtot is the evaluated volume of the enamel residual. N1, N2, …, N500 represent the sequential number of respective image window pixels. The accuracy of the developed algorithm can be enhanced by providing a precise pre-determined threshold value range to detect the gradual changes of enamel structure as a result of demineralization.

3. Results and Discussion

3.1. Morphological Analysis of Dental Caries along with Quantitative Evaluations

Figure 3a–c show the cross-sectional comparison between healthy, partially demineralized, and completely demineralized (carious) molar tooth samples. Figure 3d,e present the three-dimensional volumetric images of partially demineralized and carious molar tooth samples for a better view. In Figure 3a, the structural layers such as, enamel, dentino-enamel junction, and dentin can be visualized along with the depth ranges of 250 μm, 600 μm, and 800 μm, respectively. The progression of the demineralization can be identified in Figure 3b due to the demineralized enamel region and the formation of pits region in the depth range of 500 μm below the enamel range, and the remaining dentino-junction and dentin layer thickness was about 300 μm. Moreover, a clearly distinguishable demineralized enamel and dentino-junction (in the depth range of 800 μm), which leads to a formation of a carious region including pits and fissures were identified in Figure 3c. In addition, the remaining dentin thickness was about 100 μm. The three-dimensional images (Figure 3d,e) emphasize the top view of the partially demineralized and carious molar teeth samples. The infected microstructures, formation of pits and fissures, as well as the enamel loss, can be clearly visualized through the obtained volumetric images, owing to the high depth penetration. The desired morphological changes owing to demineralization mostly occurred in the enamel and dentin regions, which belong to a depth range that can be sufficiently achieved from the developed OFDI technique, and therefore, the necessity of histological images could be minimized by using the two-dimensional cross-sectional images acquired from the developed non-invasive OCT system [32,33,34,35]. Although the customized SD-OCT system implemented in this preliminary study improves the anatomical evaluation of the enamel and dentin, due to the low image acquisition speed and low sensitivity, the image quality has a limitation in SD-OCT compared to optical frequency domain imaging based high-speed swept source OCT (SS-OCT) [36]. Though the implementation of high-speed 1.3 μm SS-OCT can be beneficial, the customized cost effective OFDI system of this preliminary study was capable of obtaining a sufficient depth visibility to confirm the desired morphological results.
The blue color box regions of Figure 3a–c depict the averaged A-scan depth profile regions. In Figure 4a, the blue solid line, red dotted line, and black dashed line represent the A-scan depth profiles of healthy, partially demineralized, and carious samples, respectively. Owing to the high scattering coefficient of the healthy structure region, a high backscattered signal intensity can be identified compared to partially demineralized and carious samples. In addition, it is difficult to detect dentin enamel junction and dentin regions of the partially demineralized and carious samples (within depth range of 500 μm to 1000 μm) due to low signal intensity from the structure. The demineralization mainly affects the hard calcium based structural surface of the tooth enamel region by decreasing the calcium constitution and simultaneously increase the progression of caries, which contain soft tissues and blood vessels. Owing to the aforementioned increase of the soft tissue region containing blood vessels, light scattering coefficient decreases. Therefore, the degree of demineralization affects the OCT signal in depth direction. As a consequence, low signal intensity could be detected from partially demineralized samples compared to healthy samples.
The aforementioned depth dependent intensity fluctuation according to each sample within each depth range is shown in the Figure 4b graph. The total intensity fluctuation of the healthy sample, partially demineralized sample, and carious sample are illustrated in blue, red, and gray color plots, respectively. Owing to the high optical scattering coefficient, the healthy sample performs the highest total intensity in all depth ranges, and the least intensity values could be identified in the carious sample. Furthermore, the partially demineralized sample performed less total intensity than the healthy sample and a higher intensity than the carious sample, which confirms the progression of dental caries. All the obtained intensity values are shown along with the corresponding graph bars. Furthermore, the enamel thickness was evaluated, and the next bar graph shown in Figure 4c represents the depth direction enamel thickness values of the aforementioned three samples using two-dimensional cross-sectional image based A-scan depth profile analysis. The thickness of the healthy enamel region was measured as 255.45 ± 15.03 μm, a partially demineralized molar specimen was measured as 150.30 ± 10.02 μm (with a reduction of 41.1% compared to the healthy sample), and the enamel residual of the carious region was measured as 100.20 ± 6.68 μm (with a reduction of 60.8% compared to the healthy sample), respectively.
To gain a better understanding about the enamel demineralization, we repeated the aforementioned similar experimental method using partially demineralized canine and pre-molar tooth samples, which have a healthy appearance externally. Figure 5a,c emphasize the 2D OCT images along with the acquired 3D OCT volumetric images (Figure 5b,d) and sample photographs. The scanned positions of the samples are shown by red dashed lines. Although the samples have a healthy appearance, the formation of the demineralized enamel and root cavity region were identified in cross-sectional images and 3D OCT images of both canine and pre-molar tooth samples. Moreover, the enamel of the pre-molar tooth sample was demineralized towards a 750 μm depth to form the cavity region, which can be quantitatively evaluated using A-scan depth profile analysis. To acquire the most accurate quantification of enamel residual and the demineralized enamel region, we repeated the previous quantitative methods under the same experimental conditions.
The obtained A-scans of partially demineralized canine tooth sample and partially demineralized pre-molar tooth sample were compared with a healthy sample and shown in Figure 6a. In the A-scan profile, the blue solid line, black dashed line, and green dotted line represent the A-scan depth profiles of healthy sample, canine sample, and pre-molar sample, respectively. Similar signal intensity behavior was confirmed as before, owing to the low scattering coefficient of partially demineralized samples, which confirms the changes of the sample composition compared to the healthy sample. Then, we analyzed and compared the total intensity fluctuation of both partially demineralized samples according to each imaging depth range (Figure 6b).
All the obtained intensity values are shown along with the corresponding graph bars. Although a reduction of the total intensity was noticed in partially demineralized samples (compared to the healthy sample), all the intensity levels were comparatively higher than the intensity levels of the carious sample. Therefore, the results confirm that the experimentally tested samples were neither healthy nor carious and verify the progression of demineralization. Next, we analyzed and compared enamel thickness of all three samples, and the bar graph shown in Figure 6c represents the obtained depth direction enamel thickness. The thickness of the healthy enamel region was measured as 255.45 ± 15.03 μm, partially demineralized pre-molar tooth sample was measured as 180.36 ± 12.02 μm (reduction of 29.4% compared to the healthy sample), and partially demineralized canine tooth sample was measured as 140.28 ± 9.35 μm (reduction of 45.1% compared to the healthy sample), respectively. Therefore, the obtained quantitative evaluations confirm the progression of early caries, and, moreover, the obtained quantitative evaluations can be utilized as threshold parameters to detect the progression of early caries. To gain a better understanding about the quantified thickness values and total intensity fluctuations of the experimented tooth specimens, the summarized quantifications are illustrated in Table 3.

3.2. Volumetric Evaluation Technique to Identify Initial Caries

We quantified the volume of the enamel residual by applying the described volumetric algorithm in Section 2.4. Thus, we calculated the enamel residual of a selected particular position determined by the expert orthodontist for the healthy tooth specimen, three partially demineralized but healthy appearing tooth specimens, and the carious tooth specimen. All the performed calculations were based on the refractive index of 1.63, and the calculated enamel residual volume evaluations are illustrated along with the parameters in Table 4.
Hence, the obtained results confirmed that the proposed volumetric evaluation method will be more useful to detect the progression of caries, since the gradual reduction of the enamel volume owing to the gradual growth of caries can be detected quantitatively in advance. Therefore, medical treatments can be initiated immediately in order to obstruct the progression of caries, once the volume reduction of teeth is identified.

3.3. Structural Comparison between OCT and Conventional Methods

Figure 7 shows the structural analysis of a carious molar tooth sample and a comparison between imaging results obtained from various inspection methods. Figure 7a shows the in vivo radiographic image of the carious tooth, which was captured before early childhood caries (ECC) surgery performed on a 10-year-old male volunteer. The images were acquired to inspect the dental caries and cavity filling portions using a system that is currently applied in standard clinical practice: ultra speed (D-speed) film (Kodak, Rochester, NY, USA), 150 kVp, 15 mA, and 20 impulses. In this radiographic method, resolving a sub-millimeter tissue structure proves to be difficult, and only the surface structures of the cavity fillings, carious region, partially demineralized regions, pulps, and root canals along with healthy tooth were visualized. Thus, the obtained results were neither quantitative nor sensitive, and the cavity depths could not be imaged as well. Hence, precise radiographic detection of demineralization is a challenging task, since minimally demineralized regions are unable to reach the threshold of resolution. Figure 7b shows the photograph of the carious tooth, which was obtained after early childhood caries (ECC) surgery. Figure 7c,d represent the ex vivo 3D OCT images of the same sample and show the top and the side views of the sample. A precise enhancement could be identified, compared to radiographic images. Both 3D OCT figures give a clear view of the distinguishable anatomical structures e.g., dentin tubules, pulp, root canals, and cement owing to the high axial and lateral resolutions. Therefore, the applicability and the reliability of our system were verified because the dental caries, demineralization, and the inner microstructures of dentin were confirmed simultaneously through the obtained results.

4. Conclusions

We have demonstrated an optical frequency-domain imaging technique based quantitative evaluation method as an initial ex vivo study to detect the progression of dental caries by comparing partially and completely demineralized tooth samples with healthy tooth specimens. The quantification techniques were carried out by evaluating precise volume and thickness of enamel residual. Next, the total intensity fluctuation in each imaging depth range of all the specimens was quantified to confirm the changes that occurred in the internal composition of partially and completely demineralized tooth samples compared to a healthy sample. The performed study was a preliminary descriptive observational study, which was performed to confirm the feasibility of the three developed quantification techniques. The representative ex vivo tooth specimens as well as the experimental procedure was conducted according to the guidelines provided by an expert orthodontist. The results obtained using our high-resolution OFDI system revealed anatomic and quantitative information in a relatively nondestructive manner. The threshold parameters to detect the progression of early caries were determined on the basis of the quantitative results obtained from partially demineralized samples. Therefore, the physicians were able to diagnose the tooth volumetric and thickness changes at an initial stage by considering the obtained results as promising threshold parameters, which will be useful to barricade the progression of caries. To enhance the accuracy of the threshold parameters, quantitative (thickness and volumetric) information of multiple in vivo specimens will be evaluated, averaged, and normalized along with clinical trials in future studies.

Acknowledgments

This research was supported by the Industrial Strategic Technology Development Program, Grant No. 10047943; the “Development of Micro-Surgical Apparatus Based on 3D Tomographic Operating Microscope” program, funded by the Ministry of Trade, Industry and Energy (MI, Korea, No. 10047943). This study was also supported by the BK21 Plus project funded by the Ministry of Education, Korea (21A20131600011).

Author Contributions

R.E.W. and N.H.C. conducted the experiment, and analyzed and drafted the manuscript. K.P. designed the experimental setup with software and hardware components. M.J. and J.K. designed the study and critically revised and finalized the intellectual contents of manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Selwitz, R.H.; Ismail, A.I.; Pitts, N.B. Dental caries. Lancet 2007, 369, 51–59. [Google Scholar] [CrossRef]
  2. Fejerskov, O. Changing paradigms in concepts on dental caries: Consequences for oral health care. Caries Res. 2004, 38, 182–191. [Google Scholar] [CrossRef] [PubMed]
  3. Dove, S.B.; McDavid, W. A comparison of conventional intra-oral radiography and computer imaging techniques for the detection of proximal surface dental caries. Dentomaxillofac. Radiol. 1992, 21, 127–134. [Google Scholar] [CrossRef] [PubMed]
  4. Fercher, A.F. Optical coherence tomography. J. Biomed. Opt. 1996, 1, 157–173. [Google Scholar] [CrossRef] [PubMed]
  5. Ding, Z.; Ren, H.; Zhao, Y.; Nelson, J.S.; Chen, Z. High-resolution optical coherence tomography over a large depth range with an axicon lens. Opt. Lett. 2002, 27, 243–245. [Google Scholar] [CrossRef] [PubMed]
  6. Boppart, S.A.; Herrmann, J.; Pitris, C.; Stamper, D.L.; Brezinski, M.E.; Fujimoto, J.G. High-resolution optical coherence tomography-guided laser ablation of surgical tissue. J. Surg. Res. 1999, 82, 275–284. [Google Scholar] [CrossRef] [PubMed]
  7. Nassif, N.; Cense, B.; Park, B.; Pierce, M.; Yun, S.; Bouma, B.; Tearney, G.; Chen, T.; de Boer, J. In vivo high-resolution video-rate spectral-domain optical coherence tomography of the human retina and optic nerve. Opt. Express 2004, 12, 367–376. [Google Scholar] [CrossRef] [PubMed]
  8. Wijesinghe, R.E.; Park, K.; Kim, P.; Oh, J.; Kim, S.-W.; Kim, K.; Kim, B.-M.; Jeon, M.; Kim, J. Optically deviated focusing method based high-speed sd-oct for in vivo retinal clinical applications. Opt. Rev. 2016, 23, 307–315. [Google Scholar] [CrossRef]
  9. Welzel, J. Optical coherence tomography in dermatology: A review. Skin Res. Technol. 2001, 7, 1–9. [Google Scholar] [CrossRef] [PubMed]
  10. Cho, N.H.; Lee, J.W.; Cho, J.-H.; Kim, J.; Jang, J.H.; Jung, W. Evaluation of the usefulness of three-dimensional optical coherence tomography in a guinea pig model of endolymphatic hydrops induced by surgical obliteration of the endolymphatic duct. J. Biomed. Opt. 2015, 20, 036009. [Google Scholar] [CrossRef] [PubMed]
  11. Lee, J.; Kim, K.; Wijesinghe, R.E.; Jeon, D.; Lee, S.H.; Jeon, M.; Jang, J.H. Decalcification using ethylenediaminetetraacetic acid for clear microstructure imaging of cochlea through optical coherence tomography. J. Biomed. Opt. 2016, 21, 081204. [Google Scholar] [CrossRef] [PubMed]
  12. Akkaya, N.; Kansu, Ö.; Kansu, H.; Çağirankaya, L.; Arslan, U. Comparing the accuracy of panoramic and intraoral radiography in the diagnosis of proximal caries. Dentomaxillofac. Radiol. 2014. [CrossRef] [PubMed]
  13. Seneadza, V.; Koob, A.; Kaltschmitt, J.; Staehle, H.; Duwenhoegger, J.; Eickholz, P. Digital enhancement of radiographs for assessment of interproximal dental caries. Dentomaxillofac. Radiol. 2014. [CrossRef] [PubMed]
  14. Featherstone, J.; Ten Cate, J.; Shariati, M.; Arends, J. Comparison of artificial caries-like lesions by quantitative microradiography and microhardness profiles. Caries Res. 1983, 17, 385–391. [Google Scholar] [CrossRef] [PubMed]
  15. Bühler, C.M.; Ngaotheppitak, P.; Fried, D. Imaging of occlusal dental caries (decay) with near-ir light at 1310-nm. Opt. Express 2005, 13, 573–582. [Google Scholar] [CrossRef] [PubMed]
  16. Cochrane, N.J.; Iijima, Y.; Shen, P.; Yuan, Y.; Walker, G.D.; Reynolds, C.; MacRae, C.M.; Wilson, N.C.; Adams, G.G.; Reynolds, E.C. Comparative study of the measurement of enamel demineralization and remineralization using transverse microradiography and electron probe microanalysis. Microsc. Microanal. 2014, 20, 937–945. [Google Scholar] [CrossRef] [PubMed]
  17. Fried, D.; Xie, J.; Shafi, S.; Featherstone, J.D.; Breunig, T.M.; Le, C. Imaging caries lesions and lesion progression with polarization sensitive optical coherence tomography. J. Biomed. Opt. 2002, 7, 618–627. [Google Scholar] [CrossRef] [PubMed]
  18. Ishibashi, K.; Ozawa, N.; Tagami, J.; Sumi, Y. Swept-source optical coherence tomography as a new tool to evaluate defects of resin-based composite restorations. J. Dent. 2011, 39, 543–548. [Google Scholar] [CrossRef] [PubMed]
  19. Baumgartner, A.; Dichtl, S.; Hitzenberger, C.; Sattmann, H.; Robl, B.; Moritz, A.; Fercher, A.; Sperr, W. Polarization–sensitive optical coherence tomography of dental structures. Caries Res. 1999, 34, 59–69. [Google Scholar] [CrossRef]
  20. Huynh, G.D.; Darling, C.L.; Fried, D. Changes in the optical properties of dental enamel at 1310 nm after demineralization. In Proceedings of the Biomedical Optics 2004, San Jose, CA, USA, 28 May 2004; pp. 118–124.
  21. Cahill, L.; Lee, A.M.; Pahlevaninezhad, H.; Ng, S.; MacAulay, C.E.; Poh, C.; Lane, P. Passive endoscopic polarization sensitive optical coherence tomography with completely fiber based optical components. In Proceedings of the SPIE BiOS, San Francisco, CA, USA, 2 March 2015; p. 930413.
  22. Jones, R.S.; Staninec, M.; Fried, D. Imaging artificial caries under composite sealants and restorations. J. Biomed. Opt. 2004, 9, 1297–1304. [Google Scholar] [CrossRef] [PubMed]
  23. Popescu, D.P.; Sowa, M.G.; Hewko, M.D. Assessment of early demineralization in teeth using the signal attenuation in optical coherence tomography images. J. Biomed. Opt. 2008, 13, 054053. [Google Scholar] [CrossRef] [PubMed]
  24. Shimada, Y.; Sadr, A.; Burrow, M.F.; Tagami, J.; Ozawa, N.; Sumi, Y. Validation of swept-source optical coherence tomography (ss-oct) for the diagnosis of occlusal caries. J. Dent. 2010, 38, 655–665. [Google Scholar] [CrossRef] [PubMed]
  25. Cara, A.C.; Zezell, D.M.; Ana, P.A.; Maldonado, E.P.; Freitas, A.Z. Evaluation of two quantitative analysis methods of optical coherence tomography for detection of enamel demineralization and comparison with microhardness. Lasers Surg. Med. 2014, 46, 666–671. [Google Scholar] [CrossRef] [PubMed]
  26. Hsieh, Y.-S.; Ho, Y.-C.; Lee, S.-Y.; Chuang, C.-C.; Tsai, J.-C.; Lin, K.-F.; Sun, C.-W. Dental optical coherence tomography. Sensors 2013, 13, 8928–8949. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Keane, P.A.; Ruiz-Garcia, H.; Sadda, S.R. Clinical applications of long-wavelength (1000-nm) optical coherence tomography. Ophthalmic Surg. Lasers Imaging Retin. 2011, 42, S67–S74. [Google Scholar] [CrossRef] [PubMed]
  28. Jeon, M.; Kim, J.; Jung, U.; Lee, C.; Jung, W.; Boppart, S.A. Full-range k-domain linearization in spectral-domain optical coherence tomography. Appl. Opt. 2011, 50, 1158–1163. [Google Scholar] [CrossRef] [PubMed]
  29. Jung, U.-S.; Cho, N.-H.; Kim, S.-H.; Jeong, H.-S.; Kim, J.-H.; Ahn, Y.-C. Simple spectral calibration method and its application using an index array for swept source optical coherence tomography. J. Opt. Soc. Korea 2011, 15, 386–393. [Google Scholar] [CrossRef]
  30. Lee, C.; Lee, S.-Y.; Kim, J.-Y.; Jung, H.-Y.; Kim, J. Optical sensing method for screening disease in melon seeds by using optical coherence tomography. Sensors 2011, 11, 9467–9477. [Google Scholar] [CrossRef] [PubMed]
  31. Wijesinghe, R.E.; Lee, S.-Y.; Kim, P.; Jung, H.-Y.; Jeon, M.; Kim, J. Optical inspection and morphological analysis of diospyros kaki plant leaves for the detection of circular leaf spot disease. Sensors 2016, 16, 1282. [Google Scholar] [CrossRef] [PubMed]
  32. Shellis, R. Relationship between human enamel structure and the formation of caries-like lesions in vitro. Arch. Oral Biol. 1984, 29, 975–981. [Google Scholar] [CrossRef]
  33. Silverstone, L.; Poole, D. The effect of saliva and calcifying solutions upon the histological appearance of enamel caries. Caries Res. 1968, 2, 87–96. [Google Scholar] [CrossRef] [PubMed]
  34. Ekstrand, K.; Kuzmina, I.; Bjørndal, L.; Thylstrup, A. Relationship between external and histologic features of progressive stages of caries in the occlusal fossa. Caries Res. 1995, 29, 243–250. [Google Scholar] [CrossRef] [PubMed]
  35. Go, E.-J.; Jung, H.-S.; Kim, E.-S.; Jung, I.-Y.; Lee, S.-J. Histology of dental pulp healing after tooth replantation in rats. J. Korean Acad. Conserv. Dent. 2010, 35, 273–284. [Google Scholar] [CrossRef]
  36. Yun, S.; Tearney, G.; de Boer, J.; Iftimia, N.; Bouma, B. High-speed optical frequency-domain imaging. Opt. Express 2003, 11, 2953–2963. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Schematic diagram of the spectral-domain optical coherence tomography (SD-OCT) system. Note the use of the following acronyms in the figure: BS: broadband source, C: collimator, DG: diffraction grating, FC: fiber coupler, FG: frame grabber, GS: galvanoscanner, L: lens, LSC: line scan camera, M: mirror, PC: polarization controller, and ST: sample tooth.
Figure 1. Schematic diagram of the spectral-domain optical coherence tomography (SD-OCT) system. Note the use of the following acronyms in the figure: BS: broadband source, C: collimator, DG: diffraction grating, FC: fiber coupler, FG: frame grabber, GS: galvanoscanner, L: lens, LSC: line scan camera, M: mirror, PC: polarization controller, and ST: sample tooth.
Sensors 16 02076 g001
Figure 2. Volumetric evaluation algorithm for enamel residual. (a) Sequential 2D OCT images along with the applied image window; (b) Detected pixels, which satisfy the applied pre-determined threshold range; (c) Acquired 3D OCT volumetric image.
Figure 2. Volumetric evaluation algorithm for enamel residual. (a) Sequential 2D OCT images along with the applied image window; (b) Detected pixels, which satisfy the applied pre-determined threshold range; (c) Acquired 3D OCT volumetric image.
Sensors 16 02076 g002
Figure 3. Two-dimensional OCT image comparison between healthy, partially demineralized, and completely demineralized (carious) molar tooth samples along with three-dimensional OCT images. (a) 2D OCT image of a healthy molar tooth region; (b) 2D OCT image of a partially demineralized molar tooth region; (c) 2D OCT image of a completely demineralized (carious) molar tooth region; (d) 3D OCT volumetric image of a partially demineralized molar tooth sample; and (e) 3D OCT volumetric image of a carious molar tooth sample.
Figure 3. Two-dimensional OCT image comparison between healthy, partially demineralized, and completely demineralized (carious) molar tooth samples along with three-dimensional OCT images. (a) 2D OCT image of a healthy molar tooth region; (b) 2D OCT image of a partially demineralized molar tooth region; (c) 2D OCT image of a completely demineralized (carious) molar tooth region; (d) 3D OCT volumetric image of a partially demineralized molar tooth sample; and (e) 3D OCT volumetric image of a carious molar tooth sample.
Sensors 16 02076 g003
Figure 4. The quantitative evaluation method for healthy, partially demineralized, and completely demineralized (carious) molar tooth samples. (a) A-scan depth profiles of healthy, partially demineralized, and carious samples; (b) The total intensity fluctuation of healthy, partially demineralized molar, and completely demineralized (carious) molar tooth samples for the entire visible depth range of 1 mm with a gap of 250 μm depth range; (c) The depth direction enamel thickness values of the healthy, partially carious, and carious samples.
Figure 4. The quantitative evaluation method for healthy, partially demineralized, and completely demineralized (carious) molar tooth samples. (a) A-scan depth profiles of healthy, partially demineralized, and carious samples; (b) The total intensity fluctuation of healthy, partially demineralized molar, and completely demineralized (carious) molar tooth samples for the entire visible depth range of 1 mm with a gap of 250 μm depth range; (c) The depth direction enamel thickness values of the healthy, partially carious, and carious samples.
Sensors 16 02076 g004
Figure 5. Two-dimensional OCT image comparison between a partially demineralized but healthy appearing canine tooth sample and a partially demineralized but healthy appearing pre-molar tooth sample with three-dimensional OCT images. (a) 2D OCT image of a partially demineralized but healthy appearing canine tooth sample; (b) 3D OCT volumetric image of a partially demineralized but healthy appearing canine tooth sample; (c) 2D OCT image of a partially demineralized but healthy appearing pre-molar tooth sample; (d) 3D OCT volumetric image of a partially demineralized but healthy appearing pre-molar tooth sample.
Figure 5. Two-dimensional OCT image comparison between a partially demineralized but healthy appearing canine tooth sample and a partially demineralized but healthy appearing pre-molar tooth sample with three-dimensional OCT images. (a) 2D OCT image of a partially demineralized but healthy appearing canine tooth sample; (b) 3D OCT volumetric image of a partially demineralized but healthy appearing canine tooth sample; (c) 2D OCT image of a partially demineralized but healthy appearing pre-molar tooth sample; (d) 3D OCT volumetric image of a partially demineralized but healthy appearing pre-molar tooth sample.
Sensors 16 02076 g005
Figure 6. The quantitative evaluation for a healthy sample and partially demineralized but healthy appearing canine and pre-molar tooth samples. (a) A-scan depth profiles of healthy and partially demineralized but healthy appearing canine and pre-molar tooth samples; (b) the total intensity fluctuation of healthy, canine, and pre-molar tooth samples for the entire visible depth range of 1 mm with a gap of 250 μm depth range; (c) the depth direction enamel thickness values of the healthy, partially demineralized but healthy appearing canine and pre-molar tooth samples.
Figure 6. The quantitative evaluation for a healthy sample and partially demineralized but healthy appearing canine and pre-molar tooth samples. (a) A-scan depth profiles of healthy and partially demineralized but healthy appearing canine and pre-molar tooth samples; (b) the total intensity fluctuation of healthy, canine, and pre-molar tooth samples for the entire visible depth range of 1 mm with a gap of 250 μm depth range; (c) the depth direction enamel thickness values of the healthy, partially demineralized but healthy appearing canine and pre-molar tooth samples.
Sensors 16 02076 g006
Figure 7. Structural analysis of the carious tooth sample using various inspection methods. (a) Radiographic image of the carious tooth sample before the surgery; (b) Photograph of the carious tooth sample after the surgery; (c,d) top and side views of 3D OCT images of the carious tooth sample.
Figure 7. Structural analysis of the carious tooth sample using various inspection methods. (a) Radiographic image of the carious tooth sample before the surgery; (b) Photograph of the carious tooth sample after the surgery; (c,d) top and side views of 3D OCT images of the carious tooth sample.
Sensors 16 02076 g007
Table 1. The details of the optical frequency domain imaging system.
Table 1. The details of the optical frequency domain imaging system.
System ParametersSpecification
Central Wavelength1310 nm
Spectral bandwidth135 nm
Axial resolution air/tissue6 μm/3.61 μm
Transverse resolution25 μm
Maximum imaging width8 mm
Maximum imaging depth>6 mm
Optical power variation±5%
Table 2. The details of the experimented tooth specimens.
Table 2. The details of the experimented tooth specimens.
Experimented VolunteerTooth ClassificationInspection Category
11-year-old maleMolar toothHealthy
11-year-old femaleMolar toothPartially demineralized
10-year-old maleMolar toothCarious
11-year-old maleCanine toothPartially demineralized
12-year-old femalePre-molar toothPartially demineralized
Table 3. The quantified enamel thickness and depth direction total intensity fluctuations.
Table 3. The quantified enamel thickness and depth direction total intensity fluctuations.
Specimen CategoryEnamel Thickness (μm)Total Intensity Fluctuation in Each Depth Range (a.u.)
0–250 μm250–500 μm500–750 μm750–1000 μm
Healthy molar255.45 ± 15.032.34 ± 0.23.10 ± 0.22.14 ± 0.21.21 ± 0.2
Dem. molar150.30 ± 10.022.09 ± 0.22.21 ± 0.21.24 ± 0.20.77 ± 0.1
Carious molar100.20 ± 6.681.24 ± 0.10.74 ± 0.050.41 ± 0.020.11 ± 0.01
Dem. canine140.28 ± 9.351.10 ± 0.11.21 ± 0.10.99 ± 0.10.65 ± 0.1
Dem. premolar180.36 ± 12.021.91 ± 0.11.65 ± 0.11.34 ± 0.10.72 ± 0.05
Table 4. The volumetric evaluation results of the enamel residual.
Table 4. The volumetric evaluation results of the enamel residual.
Tooth SpecimenTotal Number of 2D OCT ImagesTotal Number of Enamel Residual Pixels Enamel Residual Volume (mm3)
Healthy molar tooth5002.13 × 10728.72
Part.dem. molar5001.31 × 10717.70
Carious molar5000.91 × 10712.26
Part.dem. canine5001.28 × 10717.20
Part.dem. pre-molar5001.42 × 10719.15

Share and Cite

MDPI and ACS Style

Wijesinghe, R.E.; Cho, N.H.; Park, K.; Jeon, M.; Kim, J. Bio-Photonic Detection and Quantitative Evaluation Method for the Progression of Dental Caries Using Optical Frequency-Domain Imaging Method. Sensors 2016, 16, 2076. https://doi.org/10.3390/s16122076

AMA Style

Wijesinghe RE, Cho NH, Park K, Jeon M, Kim J. Bio-Photonic Detection and Quantitative Evaluation Method for the Progression of Dental Caries Using Optical Frequency-Domain Imaging Method. Sensors. 2016; 16(12):2076. https://doi.org/10.3390/s16122076

Chicago/Turabian Style

Wijesinghe, Ruchire Eranga, Nam Hyun Cho, Kibeom Park, Mansik Jeon, and Jeehyun Kim. 2016. "Bio-Photonic Detection and Quantitative Evaluation Method for the Progression of Dental Caries Using Optical Frequency-Domain Imaging Method" Sensors 16, no. 12: 2076. https://doi.org/10.3390/s16122076

APA Style

Wijesinghe, R. E., Cho, N. H., Park, K., Jeon, M., & Kim, J. (2016). Bio-Photonic Detection and Quantitative Evaluation Method for the Progression of Dental Caries Using Optical Frequency-Domain Imaging Method. Sensors, 16(12), 2076. https://doi.org/10.3390/s16122076

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop