Skip to main content

Automatic method of analysis of OCT images in the assessment of the tooth enamel surface after orthodontic treatment with fixed braces

Abstract

Introduction

Fixed orthodontic appliances, despite years of research and development, still raise a lot of controversy because of its potentially destructive influence on enamel. Therefore, it is necessary to quantitatively assess the condition and therein the thickness of tooth enamel in order to select the appropriate orthodontic bonding and debonding methodology as well as to assess the quality of enamel after treatment and clean-up procedure in order to choose the most advantageous course of treatment. One of the assessment methods is computed tomography where the measurement of enamel thickness and the 3D reconstruction of image sequences can be performed fully automatically.

Material and method

OCT images of 180 teeth were obtained from the Topcon 3D OCT-2000 camera. The images were obtained in vitro by performing sequentially 7 stages of treatment on all the teeth: before any interference into enamel, polishing with orthodontic paste, etching and application of a bonding system, orthodontic bracket bonding, orthodontic bracket removal, cleaning off adhesive residue. A dedicated method for the analysis and processing of images involving median filtering, mathematical morphology, binarization, polynomial approximation and the active contour method has been proposed.

Results

The obtained results enable automatic measurement of tooth enamel thickness in 5 seconds using the Core i5 CPU M460 @ 2.5GHz 4GB RAM. For one patient, the proposed method of analysis confirms enamel thickness loss of 80 μm (from 730 ± 165 μm to 650 ± 129 μm) after polishing with paste, enamel thickness loss of 435 μm (from 730 ± 165 μm to 295 ± 55 μm) after etching and bonding resin application, growth of a layer having a thickness of 265 μm (from 295 ± 55 μm to 560 ± 98 μm after etching) which is the adhesive system. After removing an orthodontic bracket, the adhesive residue was 105 μm and after cleaning it off, the enamel thickness was 605 μm. The enamel thickness before and after the whole treatment decreased by about 125 μm.

Conclusions

This paper presents an automatic quantitative method for the assessment of tooth enamel thickness. This method has proven to be an effective diagnostic tool that allows evaluation of the surface and cross section of tooth enamel after orthodontic treatment with fixed thin-arched braces and proper selection of the methodology and course of treatment.

Introduction

Enamel is a tissue covering the dentine within the tooth crown. Its thickness in specific areas may exceed 1400 μm. Enamel is composed of about 98% of inorganic matter in the form of di-hydroxyapatite crystals. This chemical composition gives it unusual hardness. Organic compounds and water constitute the remaining 2% of enamel [1, 2]. Orthodontic treatment concerning a significant part of the population is related to enamel when it comes to installing fixed braces. In this case, there occur changes in the structure of enamel after treatment with fixed orthodontic appliances which depend on the used orthodontic hooks. There are different methods of attaching fixed braces and treating the tissue after hook removal. There is also different efficiency of cleaning a tooth surface off the composite material residue after the removal of orthodontic hooks when using selected methods of enamel treatment [3].

For example, there are techniques of adhesive attachment of orthodontic hooks to the enamel surface which enable rapid development of orthodontic treatment with fixed braces [46]. This method is continuously being developed and adhesive systems as well as new types of brackets are being improved. Unfortunately, fixed orthodontic appliances cause a lot of doubts, because of their potentially destructive influence on the enamel [7, 8] and others [912]. In [7], there are no statistically significant differences in volumetric changes after polishing between the different clean-up methods. However, sufficient clean-up without enamel loss was difficult to achieve. Enamel breakouts after debonding were detectable in 27% of all cases (mean depth of 44.9 mm). In [8], uncoated and precoated brackets exhibited similar debonding patterns. The debonding method tested in this study did not restore the original enamel surface, although there was no clinically relevant enamel damage. The use of coherence tomography in the areas of the mentioned enamel thickness or in the assessment of the strength of bonding is presented both in [9], with the use of a dental OCT handpiece, and in [10, 12]. Therefore, there is a need for quantitative and not only qualitative assessment of tooth enamel in order to select the appropriate methodology and course of treatment as well as to assess the quality of enamel after treatment [1316].

Evaluation of enamel quality can be done in many different ways. The known methods use an atomic force microscope (AFM) and a scanning electron microscope (SEM). It can also be evaluated using optical coherence tomography (OCT) [12]. This imaging method enables non-invasive structure visualization on the maximum enamel thickness and the desired tooth surface. The use of an optical tomograph having a wavelength of 840 nm to evaluate the enamel thickness is also interesting. For example, in [9] a central wavelength of 1310 nm and a spectral bandwidth of 47 nm were used for this purpose.

Unfortunately, at present it is relatively difficult to compare images taken at different time intervals and at different stages of treatment (diagnosis) [1719]. The difficulty lies in automatic finding of the same areas or points on two or more images. Therefore, there is a need to propose a new automatic algorithm able to analyse OCT images acquired at different time intervals and quantitatively assess the treatment progress [20, 21]. A profiled automatic algorithm for the analysis and processing of OCT images of tooth enamel, in particular, should enable to:

 quantify the changes in the enamel structure after treatment with fixed braces depending on the used orthodontic brackets, methods of attachment and tissue development after bracket removal,

 quantify the efficiency of cleaning a tooth surface off the composite material residue after the removal of orthodontic brackets when using selected ways of enamel development.

A description of the dedicated automatic algorithm is discussed in detail later in the paper.

Material

OCT images of 180 teeth were obtained from the Topcon 3D OCT-2000 camera - Figure 1. The images were obtained in vitro by performing sequentially 7 stages of treatment on each of the 180 teeth: before any interference into enamel, polishing with orthodontic paste, etching and application of a bonding system, orthodontic bracket bonding, orthodontic bracket removal, cleaning off adhesive residue. 7 OCT images were obtained at each of the 7 stages (the first image refers to a tooth without any interference).

Figure 1
figure 1

Image showing the method of acquisition of OCT images of the teeth. The following items are depicted: 1 - OCT tomograph, 2 - digital camera for taking images in visible light, 3 - screen of the tomograph, 4 - joystick enabling changes object position, 5 - method of attachment of the tooth in the device.

All the teeth were extracted and then stored for 30 days in distilled water with 0.1% crystal thymol at room temperature, and then cleaned with pumice using a eraser (TopDental), rinsed with distilled water and dried with compressed air for 15 seconds. The composite material BluGloo (ORMCO, USA) was used to attach orthodontic brackets. The vestibular tooth surface was etched for 30 seconds with 36% phosphoric acid - Blue-Etch (CERKAMED, Poland), rinsed with distilled water and dried with compressed air. The adhesive system Ortho Solo (ORMCO, USA) was applied to the etched surface. A bracket was attached with the use of clamping tweezers for tooth enamel so that its centre was 3.5 mm below the edge of the occlusal surface, in the middle of the mesial-distal axis of the tooth. The test teeth were stored in water at 37˚C for 24 hours. Adhesive material residue was removed from enamel in two ways: by the use of a carbide cutter and by means of cup-shaped elastics OneGloss 0181 (Shofu Dental Corp., USA).

The acquired images (B-scans) had a resolution of M × N = 884 × 512 pixels (where M - number of rows, N - number of columns) where a single pixel covered the area of 5 × 11.7 μm at the colour resolution of 8 bits/pixel. For each patient I = 128 B-scans (every 47 μm) were obtained, which allowed for full 3D tooth image reconstruction. All data were analysed in a source format *.fda, *.fds. The data were anonymised and tooth extraction in all the patients was performed in accordance with the Declaration of Helsinki.

Method

Pre-processing

Pre-processing of images is related to automatic reading of a sequence of OCT images from a source file with the extension *.fds, thus obtaining a matrix of individual images. The read data also include information on the location of subsequent B-scans, their size, date of test, patient data and others. The sequence of images L GRAY (m,n,i) (Figure 2) thus obtained is further subjected to median filtering with a mask h 1 sized M h1  × N h1 × I h1  = 3 × 3 × 3. The mask size was chosen arbitrarily taking into account the size of noise and artefacts entering the optical path. The filtered image L M (m,n,i) (where m,n are the coordinates of the rows and columns of the matrix m(1,M) and n(1,N)) is further subjected to preliminary transformations. These include determination of the enamel layer, measurement of its thickness in the entire area, comparison with other images of the same tooth and determination of the value of the proposed enamel quality coefficient.

Figure 2
figure 2

Reconstruction of a sequence of images L GRAY ( m,n,i ) for M×N×I= 884 × 512 × 128 pixels. The image shows the result of reconstruction of an image sequence with a sample B-scan obtained for i = 60. The analysis aims at automatic determination of the edge on a sequence of OCT images (B-scans) in order to determine enamel thickness.

Processing

The filtered image L M (m,n,i) free of noise and artefacts is the basis for further transformations. These include determination of enamel thickness and imposition of two 2D images (C-scans) in order to compare changes in enamel thickness.

Determination of enamel thickness

Determination of enamel thickness is related to determining the position of two enamel layers: outer L en (n,i) and inner L in (n,i) directly in the image L M (m,n,i). Determination of the outer layer L en (n,i) started from setting the upper limit of L ef (n,i) between the object (objects) and background:

L ef n , i = min m L efm m , n , i
(1)

where:

L efm n , i = m if L M m , n , i > p r N other
(2)

and p r – binarization threshold determined automatically – Otsu [22].

The upper limit of the object L ef (n,i) is shown in green in Figure 3a. The presented method according to the equations (1) and (2) gives correct results in 100% of cases (for more than 20,000 images). This limit is used to determine the correct upper and lower boundaries of enamel (only a portion of the upper limit of objects L ef (n,i) is the appropriate outer boundary of the enamel L en (n,i)).

Figure 3
figure 3

Example of the result of the upper boundary analysis and selected parts of the analysis. The image a) shows the result (in green) of automatic detection of the upper boundary L ef (n,i). The image b) shows changes in brightness for the selected A-scan (L M (m,n = 240, i = 60)) and the image c) shows the histogram on the basis of which the binarization threshold is selected in order to detect the enamel inner boundary.

Detection of the enamel lower boundary L in (n,i) is related to the analysis of each column of the image L M (m,n,i) in the range from L ef (n,i) to the last row. Sample data (levels of brightness) are shown in Figure 3b and consist of two parts. Part I (Figure 3b concerns reduction of brightness of a fixed speed whereas part II refers to a constant brightness level. Part III is a shift of the mean (line II) by a standard deviation of the mean from this range to the top. The grey level in the range from p r2 to p r3 (II to III) is a place indicating tooth translucency reduction related to the end of the enamel layer and the start of another layer, namely dentine. The decision which of the threshold values should be taken as the correct one (in the range from p r2 to p r3 ) is connected in this case with a limit shift of about 70 pixels which corresponds to 70·5 = 350 μm. The value of this threshold can be automatically determined and clarified on the basis of the previous and subsequent columns and histogram analysis - Figure 3c. The threshold value p r2 (n,i) (simplified to p r2 ) was calculated automatically based on the histogram:

p r 2 n , i = arg b max b L hist n , b , i
(3)

where:

L hist n , b , i = m = 1 M L s m , n , b , i
(4)
L s m , n , b , i = 1 if L M m , n , i = b 0 other
(5)

The obtained result, image L B (m,n,i) for the adopted binarization threshold, is shown in Figure 4a:

L B m , n , i = 1 if L M m , n , i > p r 2 0 other
(6)
Figure 4
figure 4

Example of the binarization result of an OCT image and the results of analysis of inner and outer boundaries of the enamel. The image a) shows problems that occur when determining the enamel inner boundary on a binary image whereas the image b) shows an example of the result of analysis based on the proposed algorithm. The outer boundary is marked in green, the inner boundary before filtration in red and the inner boundary after filtration in blue.

The binary image L B (m,n,i) is further subjected to median filtering with a mask h 2 sized M h2  × N h2 × I h2  = 3 × 3 × 3. Then, the resulting image L MB (m,n,i) undergoes a morphological closing operation, i.e.:

L C m , n , i = min SE max SE L MB m , n , i
(7)

where SE – structural element sized 3 × 3.

The size of the structural element was chosen on the basis of measurements and average noise sizes (isolated ones) that make up clusters consisting of no more than 2 or 3 pieces. The next step involves filling holes and labelling the image L C (m,n,i). All further operations are carried out for each of the two-dimensional images for constant i, that is L C (m,n) (L C (m,n,i = const.)). Labelling enables to give a separate label (number) to each cluster of ones L W (m,n). In this case, it is possible to eliminate all of the clusters which do not have a maximum surface area. It results in the image L W (m,n) shown in Figure 4. In the next step, the inner enamel boundary L in (n) is determined as a 3rd order polynomial approximation of the course of the inner edge L inr (n) of an object in the image L W (m,n). An example of a polynomial takes the following form:

L in n = 0.0001 · n 3 + 0.018 · n 2 3.21 · n + 693.39
(8)

The outer edge of the enamel, L en (n) for short, is part of the previously defined edge of the object L ef (n) (for i = const.), i.e.:

L en n = L ef n if L ef n + 1 = 1 L ef n < N 0 other
(9)

Based on the determined edges L ef (n) and L in (n), their difference is calculated:

L d n = L in n L en n
(10)

This difference is calculated for every i (for each B-scan), i.e. L d (n,i) forms a matrix of differences. The obtained matrices of differences for individual tooth treatment steps are shown in Figure 5 in the following order: the tooth before interference L da (n,i) - Figure 5a, after polishing with paste L db (n,i) - Figure 5b, after etching L dc (n,i) - Figure 5c, after the application of the adhesive system L dd (n,i) - Figure 5d, the orthodontic bracket L de (n,i) - Figure 5e, after orthodontic bracket removal L df (n,i) - Figure 5f, after cleaning off adhesive residue L dg (n,i) - Figure 5g. These matrices (from L da (n,i) to L dg (n,i)) are considered in the subsequent stages of analysis.

Figure 5
figure 5

Examples of results of enamel thickness measurement for the subsequent stages of tooth treatment. It shows: a) enamel thickness before any interference L da (n,i), b) after polishing with paste L db (n,i), c) after etching L dc (n,i), d) after application of the adhesive system L dd (n,i), e) an orthodontic bracket L de (n,i), f) after bracket removal L df (n,i), g) after cleaning off adhesive residue L dg (n,i). These images then form the basis for the analysis of the impact of treatment methods and choice of orthodontic systems on changes in enamel thickness. The thickness is expressed in pixels for the artificial colour palette (ranges from blue to red, and passes through the colours cyan, yellow, and orange).

Imposition of two 2D images

The images from L da (n,i) to L dg (n,i) show enamel thickness in various stages of tooth treatment. From a practical point of view, a change in thickness is important and, to be more exact, a change in enamel volume for each step. Therefore, it is necessary to propose a method of automatic imposition of individual images from L da (n,i) to L dg (n,i). Imposition of images, finding homogeneous points, is connected in this case with finding three values, namely a shift in the axis Δn, a shift in the axis Δi and zooming Δz. In practice, due to the specificity of attachment used, the errors resulting from rotation did not exceed 1°. Because of this, it was assumed that there is no error in the images caused by rotation. The searched shifts in the axes Δn and Δi are further expressed in pixels whereas zooming Δz is expressed as a percentage of the original image, i.e. 10% means a 10% zoom with respect to the reference image.

It is most convenient to find homogeneous points in standardized input images after median filtering and not in images showing different enamel thickness, i.e. (assuming that the denominator is different from zero):

L T n , i = m = 1 M L M m , n , i min m , i m = 1 M L M m , n , i max m , i m = 1 M L M m , n , i min m , i m = 1 M L M m , n , i
(11)

The first image L Ta (n,i) is here treated as the reference one, symbols of the other ones are analogous to L d (n,i), these are subscripts from ‘b’ to ‘g’ i.e. from L Tb (n,i) to L Tg (n,i). For each new shift and zoom relative to the reference image L Ta (n,i), the value of the match criterion J was calculated. The criterion J is the sum of absolute differences between the reference image L Ta (n,i) and the image being compared, e.g. L Tb (n,i). Therefore a point reaching the minimum and the values of Δn, Δi and Δz for which it occurred are looked for. It is expressed as a percentage of changes with respect to the reference image, i.e., e.g. -10% means 90% of the reference image for a given magnification of the image (superscript is expressed as a percentage, e.g. for 120% relative to a standard, L Tb (100)(n + Δn,i + Δi). Therefore, the criterion J calculated for a pair of images L Ta and L Tb takes the following form:

J = 1 I · N i = 1 I n = 1 N L Ta n , i L Tb 100 n + Δ n , i + Δ i
(12)

for Δn(−40,40) pixels and Δi(−40,40) pixels.

The missing or redundant rows and columns in relation to L Ta (n,i) were removed or zero was written so that the numbers of rows and columns are identical. The position in both axes was changed every pixel, whereas zooming was carried out in the range from −20% to 20% every 10% relative to the reference image. The zoom for which calculations were made was marked as a superscript, e.g. if there is no image magnification then L Tb (0)(n + Δn,i + Δi). The obtained results, changes in the criterion J for different positions of the image L Tb (0)(n + Δn,i + Δi), in the absence of magnification, are shown Figure 6. After imposition of all the images and appropriate adjustment of the position of images from L db (n,i) to L dg (n,i) (after adjustment starting from L db *(n,i) to L dg *(n,i)), a region of interest (ROI) is determined. It has to meet two criteria: its total area must be within the matrix area and it must coincide with the contour of the tooth. The second criterion requires determination of the tooth area.

Figure 6
figure 6

Graph of changes in the criterion J for a sample pair of images L Tb (0) ( n + Δn,i + Δi ) and L Ta ( n,i ) and images after correction. The value looked for in the graph is the minimum for which the values of shifts in both axes (Δn, Δi) and the zoom which gives the smallest value are read. In this case, the graph refers to a constant zoom of the image L Tb (0) (0%) and for the example shown the following values were read from the graph: Δn = 9 and Δi = 4 pixels. In the left hand side of the image, there are the results obtained after shifting: images L Ta and L Tb *.

Determination of the tooth contour

Determination of the tooth contour is related to the analysis of one selected image L da (n,i) or L Ta (n,i). The other images from L db *(n,i) to L dg *(n,i) have contours in the same place, which results from their earlier imposition (finding homogeneous points). Just like in the case of imposition of images, the image L Ta (n,i) is easier to analyse than its corresponding image of the enamel thickness L da (n,i). In the first step, the image L Ta (n,i) is subjected to median filtering with a mask h 3 sized N h3 × I h3  = 31 × 31 pixels. The mask size was selected on the basis of the average width of the tooth and the size of artefacts being the image noise. The output image L Ra (n,i) is the basis for determining the contour using the standard method and the active contour method. For the latter method, the following operating parameters were adopted: changes in the axes of rows and columns in the range of ±20 pixels and averaging the grey level in the areas sized 6 × 6 pixels on both sides of the analysed pixel. The algorithm stops work at a time when there is no evidence of differences between successive steps of bringing the contour closer to the object. The results are shown in Figure 7. The determined ROI is the basis for further analysis of the results. The resulting images from L ROIa to L ROIg contain information about the thickness of the layer only in the designated area of the tooth.

Figure 7
figure 7

Figure showing operation of the active contour method. The input image L Ta after median filtering (L Ra ) is analysed using the active contour method. The standard whose size is independent of the size of the object is highlighted in green in the image. The blue shows the final stage of approximation - the area of ROI analysis.

A block diagram of the proposed algorithm is shown in Figure 8.

Figure 8
figure 8

Block diagram of the proposed method for image analysis and processing. First, the input image is subjected to pre-processing (reading of data from the format *.fds, filtration) and then the enamel layer is recognized on subsequent B-scans. After that, C-scans obtained for subsequent phases of the enamel treatment are imposed on each other and the difference in enamel thickness in the ROI is measured. The last stage of processing involves comparison of changes in the thickness of enamel - its volume for the subsequent processing steps.

Results

The obtained images from L ROIa to L ROIg containing information about the thickness of the enamel layer are the basis for the calculation of mean thickness, standard deviation, and minimum and maximum values. The results are shown in Table 1, whereas Table 2 presents information about enamel thickness differences between the images. According to the adopted symbols from L ROIa to L ROIg , these are: tooth before interference (index ‘a’), after polishing with paste (index ‘b’), after etching (index ‘c’), after application of an adhesive system (index ‘d’), an orthodontic bracket (index ‘e’), after removal of the orthodontic bracket (index ‘f’) and after cleaning off adhesive residue (index ‘g’). The differences confirm enamel thickness loss of 80 μm (enamel thickness at this stage is 650 ± 129 μm) after polishing with paste, loss of 435 μm (enamel thickness at this stage is 295 ± 55 μm) after bonding, an increase in the layer thickness of 265 μm (enamel thickness at this stage is 560 ± 98 μm) which is the bonding system. After the removal of the orthodontic bracket, the adhesive residue is 105 μm (enamel thickness at this stage is 455 ± 91 μm) and after removing the adhesive residue, the enamel thickness is 605 μm. The thickness of the enamel before and after the whole treatment decreased by about 125 μm.

Table 1 Mean, minimum and maximum values of the enamel thickness in the ROI and standard deviation of the mean for the subsequent processing steps (according to the ROI designated in various stages - Figure 5)
Table 2 Changes in the enamel thickness in the ROI for the subsequent processing steps (according to the ROI designated in various stages - Figure 5 )

The results in Table 1 contain information about the minimum and maximum tooth enamel thickness. The differences in these values (minimum and maximum) result from actual variations in the enamel thickness in various locations of the tooth and the method errors. These errors are influenced to the greatest extent by the automatic selection of binarization thresholds p r2 and p r3 shown in detail in Figure 3 and the consequences in the binary image in Figure 4.

Accuracy verification of the proposed method

The accuracy of the proposed method was assessed on the basis of the comparison of the results of enamel thickness measurements obtained with the proposed algorithm with independent measurements carried out by an expert. They involved manual marking of the enamel boundary by an expert on a tomographic image for 20 teeth (20·128 = 256 B-scans). In the case of the manual method, an expert dentist indicated two boundaries: the enamel boundary L in (E)(n) and the outer tooth boundary L en (E)(n) - Figure 9. On this basis, the enamel thickness was determined (analogously to (10)):

L d E n = L in E n L en E n
(13)
Figure 9
figure 9

Results of the automatic operation of the algorithm and the manual method of determining the enamel thickness. The expert dentist marked two contours - the outer tooth contour L en (E)(n) (in magenta) and the enamel boundary L in (E)(n) (in yellow). On this basis and on the basis of the results of automatic analysis (L en (n) green, L in (n) blue) with the proposed algorithm, the measurement error was calculated per 20 teeth. This error did not exceed 15%.

The measurement error, adopted as the measurement error of the proposed algorithm, was designated as:

δ d = 1 N · n = 1 N L d n L d E n L d E n · 100 %
(14)

For the analysed 20 teeth, the error δ d did not exceed 15%. Individual results of the manual and automatic analyses are shown in Figure 9. In none of the analysed 20 teeth, the algorithm was mistaken by more than a few pixels when defining the outer tooth layer - L en (n). The contour of this layer was in each case very clearly visible and the differences between layer designation by an expert and the algorithm were visible only after zooming (in the lower part of the image Figure 9). In contrast, the enamel layer marked by the expert was only visible in selected ranges of columns. This situation, for layers L in (n) and L in (E)(n), is shown enlarged in Figure 9 at the bottom. Individual yellow points spaced apart by a few pixels indicated by the expert are visible. These points were further connected with segments. It was possible owing to information derived from anthropometric data, on whose basis the curvature of the contour and its possible variation per a single tooth are known. In summary, using the described algorithm, the measurement error δ d does not exceed 15%.

Comparison with other known methods

The method of analysis and processing of OCT images of tooth enamel developed by the authors can be compared to other known methods of analysis. These include:

 methods of analysis and processing of OCT images applicable in ophthalmology (e.g. recognition of retinal layers) and

 teeth image analysis carried out in visible light.

In the first group, there are methods which enable detection of individual layers in OCT images of the eye fundus [2325]. The class of images is very similar to OCT images of tooth enamel. These methods work in a time comparable to that of the algorithm (from 10 to 200 ms/2D image) [2628]. However, specificity of OCT images for the eye enables a significant simplification of the algorithm in the detection area of individual layers. This is due to a much better contrast between the layers. The advantage of OCT images of tooth enamel is the necessity of determining usually only one layer.

The other group of methods is closely related to dentistry but the ray spectrum range is significantly different (visible light). The acquired images do not allow for assessment of the thickness of any tissue, however, some of the methods of image imposition can be transferred to OCT imaging. An important difference is the possibility of using, in the case of visible light, RGB components which enable to distinguish a tooth from a handle and background. The applied active contour method brings correct results in both cases (visible light and OCT).

The forms of imaging, namely in vitro and in vivo, are also extremely important here. The differences are very interesting not only from the clinical point of view but also in terms of the algorithm of image analysis and processing. In the first case (in vitro), repeatable attachment of the tooth greatly facilitates image analysis - the segmentation of the tooth boundaries. In the other case (in vivo), the situation is more difficult because not only the background but also the mutual adjacent teeth arrangement change. In this case, there may also occur other problems with imaging such as placing the tomograph perfectly perpendicular to the tooth etc. However, preliminary tests confirm the usefulness of the proposed method, also in in vivo studies. The proposed segmentation process is resistant, in the area declared, to the impact of the background of the surrounding teeth (in the case of in vivo tests) on the results obtained. This is achieved owing to a segmentation method using the active contour. In this method, an important element is the contour of the tooth and not the background as may be the case in the method based on thresholding (binarization).

The presented method can also be compared to other methods of imaging of the enamel layer. As mentioned in the introduction, the known methods of tooth enamel analysis include assessment by means of an atomic force microscope (AFM) [29, 30] and a scanning electron microscope (SEM) [31]. The class of images and the depth of imaging in these cases are quite different from the ones in OCT and typically confined to the outer enamel surface analysis – like, e.g., in [29]. A different class of images determines a different type of a profiled algorithm for image analysis and processing.

The other known methods of enamel thickness analysis do not enable automatic, quantitative measurement of the enamel thickness present in the ROI and automatic comparison of image groups. This is the case in [32], where comparisons between specific areas of the tooth enamel were made manually in OCT images. Automatic measurement was presented only in [33]. However, it concerns polarization sensitive optical coherence tomography (PS-OCT) and is not related to the problem of overlap of individual images in the subsequent processing stages of the tooth as shown in this paper.

Summary

The paper presents the possibility of using the profiled algorithm proposed by the authors for the analysis of OCT images of tooth enamel. The algorithm of OCT image analysis and processing has the following features:

 works fully automatically - without any operator intervention,

 allows to quantify the changes in the structure of enamel,

 allows for quantitative assessment of the effectiveness of cleaning the tooth surface and the effectiveness of the use of selected methods of enamel development,

 analysis time of a sequence of 2D images does not exceed 5 seconds when using the Core i5 CPU M460 @ 2.5GHz 4GB RAM,

 the results of the mean thickness of the tooth enamel, minimum and maximum values as well as standard deviation are analysed automatically and saved to text files *.txt and Excel *.xls.

The presented algorithm is only one possibility of completing the task described in the paper. Automatic analysis of tooth enamel thickness provides a number of further possibilities. These include area analysis of enamel thickness (for each individual tooth area separately) and enamel texture analysis. Imaging and quantitative measurement of the enamel structure before installation of braces and after their removal enables to expose the tooth tissue damage extent depending on the used brackets and method of attachment. This makes it possible to deduce which brackets and what technique of their installation is the safest for tooth enamel.

The presented methodology does not fully exhaust the issue of automatic measurement of enamel thickness and the effect of the method of attachment of orthodontic brackets. The presented tool will be used in future work to assess the impact of multiple parameters. For example, the impact of the OCT settings on the results will be analysed. Moreover, the effect of time of in vitro studies on the result will be assessed. The analysis of reproducibility of results obtained in tomographic cameras at different wavelengths in the range from 600 to 1400 nm is equally interesting. An important aspect here will be confirmation of the universality of the proposed algorithm of image analysis and processing for different types of OCT devices or even the impact of the type of studies, in vitro or in vivo, on the process of segmentation.

It is also possible to modify the approach to the presented problem by using, for example, texture analysis [34], morphological approach [35], spectral analysis or other methods used for slightly different applications which, in practice, generate correct results [3641]. In the future, the authors intend to perform and describe these methods, especially area and morphological analysis of the tooth enamel, by comparing the results obtained and measurement errors resulting from various adopted methods of analysis.

Abbreviations

ROI:

Region of interest

AFM:

Atomic force microscope

SEM:

Scanning electron microscope.

References

  1. Jones RS, Darling CL, Featherstone JD, Fried D: Remineralization of in vitro dental caries assessed with polarization-sensitive optical coherence tomography. J Biomed Opt 2006,11(1):014016. 10.1117/1.2161192

    Article  Google Scholar 

  2. Bishara SE, Oonsombat C, Soliman MM, Warren JJ, Laffoon JF, Ajlouni R: Comparison of bonding time and shear bond strength between a conventional and a new integrated bonding system. Angle Orthod 2005,75(2):237–242.

    Google Scholar 

  3. Buyukyilmaz T, Usumez S, Karaman AI: Effect of self-etching primers on bond strength-are they reliable? Angle Orthod 2003, 73: 64–70.

    Google Scholar 

  4. Bishara SE, Von Wald L, Laffoon JF, Warren JJ: Effect of a self-etch primer/adhesive on the shear bond strength of orthodontic brackets. Am J Orthod Dentofacial Orthop 2001, 119: 621–624. 10.1067/mod.2001.113269

    Article  Google Scholar 

  5. Yamada R, Hayakawa T, Kasai K: Effect of using self-etching primer for bonding orthodontic brackets. Angle Orthod 2002,7(6):558–564.

    Google Scholar 

  6. Grubisa HS, Heo G, Raboud D, Glover KE, Major PW: An evaluation and comparison of orthodontic bracket bond strengths achieved with self-etching primer. Am J Orthod Dentofacial Orthop 2004, 126: 213–219. 10.1016/j.ajodo.2004.01.016

    Article  Google Scholar 

  7. Ryf S, Flury S, Palaniappan S, Lussi S, van Meerbeek S, Zimmerli B: Enamel loss and adhesive remnants following bracket removal and various clean-up procedures in vitro. Eur J Orthod 2012,34(1):25–32. 10.1093/ejo/cjq128

    Article  Google Scholar 

  8. Bonetti GA, Zanarini M, Incerti Parenti S, Lattuca M, Marchionni S, Gatto MR: Evaluation of enamel surfaces after bracket debonding: An in-vivo study with scanning electron microscopy. Am J Orthod Dentofacial Orthop 2011,140(5):6967–02.

    Google Scholar 

  9. Colston BW, Sathyam US, DaSilva LB, Everett MJ: Dental OCT. Opt Express 1998,3(6):230–238. 10.1364/OE.3.000230

    Article  Google Scholar 

  10. Otis LL, Everett MJ, Sathyam US, Colston BW: Optical coherence tomography: a new imaging technology for dentistry. J Am Dent Assoc 2000,131(4,5):11–14.

    Google Scholar 

  11. Tao YC, Fried DJ: Near-infrared image-guided laser ablation of dental decay. J Biomed Opt 2009,14(5):054045. 10.1117/1.3253390

    Article  Google Scholar 

  12. Huang D, Swanson EA, Lin CP, Schuman JS, Stinson WG, Chang W, Hee MR, Flotte T, Gregory K, Puliafito CA, Fujimoto JG: Optical coherence tomography. Science 1991,254(5035):1178–1181. 10.1126/science.1957169

    Article  Google Scholar 

  13. Sumi Y, Ozawa N, Nagaosa S, Minakuchi S, Umemura O: Application of optical coherence tomography (OCT) to nondestructive inspection of dentures. Arch Gerontol Geriatr 2011,53(2):237–241. 10.1016/j.archger.2010.11.022

    Article  Google Scholar 

  14. Baek JH, Na J, Lee BH, Choi E, Son WS: Optical approach to the periodontal ligament under orthodontic tooth movement: a preliminary study with optical coherence tomography. Am J Orthod Dentofacial Orthop 2009,135(2):252–259. 10.1016/j.ajodo.2007.10.037

    Article  Google Scholar 

  15. Gracez AS, Suzuki SS, Ribeiro MS, Mada EY, Freitas AZ, Suzuki H: Biofilm retention by 3 methods of ligation on orthodontic brackets: a microbiologic and optical coherence tomography analysis. Am J Orthod Dentofacial Orthop 2011,140(4):193–198. 10.1016/j.ajodo.2011.04.019

    Article  Google Scholar 

  16. Fried D, Featherstone JD, Darling CL, Jones RS, Ngaotheppitak P, Bühler CM: Early caries imaging and monitoring with near-infrared light. Dent Clin N Am 2005,49(4):771. 10.1016/j.cden.2005.05.008

    Article  Google Scholar 

  17. Wu J, Fried D: High contrast near-infrared polarized reflectance images of demineralization on tooth buccal and occlusal surfaces at lambda = 1310 nm. Lasers Surg Med 2009,41(3):208–213. 10.1002/lsm.20746

    Article  Google Scholar 

  18. Kang H, Darling CL, Fried D: Nondestructive monitoring of the repair of enamel artificial lesions by an acidic remineralization model using polarization-sensitive optical coherence tomography. Dent Mater 2012,28(5):48894.

    Article  Google Scholar 

  19. Lee C, Darling CL, Fried D: Polarization-sensitive optical coherence tomographic imaging of artificial demineralization on exposed surfaces of tooth roots. Dent Mater 2009,25(6):721–728. 10.1016/j.dental.2008.11.014

    Article  Google Scholar 

  20. Kang H, Jiao JJ, Lee C, Le MH, Darling CL, Fried DJ: Nondestructive assessment of early tooth demineralization using cross-polarization optical coherence tomography. IEEE J Sel Top Quantum Electron 2010,16(4):870–876.

    Article  Google Scholar 

  21. Van Waes H, Matter T, Krejci I: Three-dimensional measurement of enamel loss caused by bonding and debonding of orthodontic brackets. Am J Orthod Dentofacial Orthop 1997, 112: 666–669. 10.1016/S0889-5406(97)70232-4

    Article  Google Scholar 

  22. Otsu N: A threshold selection method from gray-level histograms. IEEE Trans Syst Man Cybern 1979,9(1):62–66.

    Article  MathSciNet  Google Scholar 

  23. Mishra A, Wong A, Bizheva K, Clausi DA: Intra-retinal layer segmentation in optical coherence tomography images. Opt Express 2009,17(26):23719–23728. 10.1364/OE.17.023719

    Article  Google Scholar 

  24. Lu S, Cheung C, Liu J, Lim S, Leung C, Wong T: Automated layer segmentation of optical coherence tomography images. IEEE Trans Biomed Eng 2010,57(10):2605–2608.

    Article  Google Scholar 

  25. DeBuc C, Somfai GM: Early detection of retinal thickness changes in diabetes using optical coherence tomography. Med Sci Monit 2010,16(3):MT15-MT21.

    Google Scholar 

  26. Koprowski R, Wrobel Z: Identification of layers in a tomographic image of an eye based on the canny edge detection. Inf Technol Biomed Adv Intell Soft Comput 2008, 47: 232–239. 10.1007/978-3-540-68168-7_26

    Article  Google Scholar 

  27. Koprowski R, Wróbel Z: Image processing in optical coherence tomography: using matlab. Katowice, Poland: University of Silesia; 2011. http://www.ncbi.nlm.nih.gov/books/NBK97169/

    Google Scholar 

  28. Koprowski R, Wróbel Z: Layers recognition in tomographic eye image based on random contour analysis. Computer recognition systems 3. Adv Intell Soft Comput 2009, 57: 471–478. 10.1007/978-3-540-93905-4_56

    Article  Google Scholar 

  29. Finke M, Parker DM, Jandt KD: Influence of soft drinks on the thickness and morphology of in situ acquired pellicle layer on enamel. J Colloid Interface Sci 2002,251(2):263–270. 10.1006/jcis.2002.8428

    Article  Google Scholar 

  30. Watari F: In situ quantitative analysis of etching process of human teeth by atomic force microscopy. J Electron Microsc (Tokyo) 2005,54(3):299–308. 10.1093/jmicro/dfi056

    Article  Google Scholar 

  31. Lorenzo MC, Portillo M, Moreno P, Montero J, García A, Santos-Del Riego SE, Albaladejo A: Ultrashort pulsed laser conditioning of human enamel: in vitro study of the influence of geometrical processing parameters on shear bond strength of orthodontic brackets. Lasers Med Sci 2013, 1: 10.

    Google Scholar 

  32. Wilder-Smith CH, Wilder-Smith P, Kawakami-Wong H, Voronets J, Osann K, Lussi A: Quantification of dental erosions in patients with GERD using optical coherence tomography before and after double-blind, randomized treatment with esomeprazole or placebo. Am J Gastroenterol 2009,104(11):2788. 10.1038/ajg.2009.441

    Article  Google Scholar 

  33. Le MH, Darling CL, Fried D: Automated analysis of lesion depth and integrated reflectivity in PS-OCT scans of tooth demineralization. Lasers Surg Med 2010,42(1):62. 10.1002/lsm.20862

    Article  Google Scholar 

  34. Korzynska A, Hoppe A, Strojny W: Investigation of a combined texture and contour method for segmentation of light microscopy cell images. In Proceedings of the Second IASTED International Conference on Biomedical Engineering. Anaheim, Calif. Calgary: Zurich ACTA Press; 2004:234–239.

    Google Scholar 

  35. Korzynska A, Iwanowski M: Multistage morphological segmentation of bright-field and fluorescent microscopy images. Opt-Electron Rev 2012,20(2):87–99.

    Google Scholar 

  36. Siedlecki D, Nowak J, Zajac M: Placement of a crystalline lens and intraocular lens: retinal image quality. J Biomed Opt 2006,11(5):054012. 10.1117/1.2358959

    Article  Google Scholar 

  37. Porwik P: Efficient spectral method of identification of linear Boolean function. Control Cybern 2004,33(4):663–678.

    MathSciNet  Google Scholar 

  38. Korus P, Dziech A: Efficient method for content reconstruction with self-embedding. IEEE Trans Image Process 2013,22(3):1134–1147.

    Article  MathSciNet  Google Scholar 

  39. Wilczyński S, Ptaszkiewicz M, Pierzchała E, Pilawa B, Swakoń J, Olko P: Electron paramagnetic resonance studies of gamma irradiated azithromycin. Curr Top in Biophysics 2008, 31: 1–4.

    Google Scholar 

  40. Tadeusiewicz R, Ogiela MR: Automatic understanding of medical images new achievements in syntactic analysis of selected medical images. Biocybern Biomed Eng 2002,22(4):17–29.

    Google Scholar 

  41. Wilczyński S, Pilawa B, Koprowski R, Wróbel Z, Ptaszkiewicz M, Swakoń J, Olko P: EPR studies of free radical decay and survival in gamma irradiated aminoglycoside antibiotics: sisomicin, tobramycin and paromomycin. Eur J Pharm Sci 2012, 45: 251–262. 10.1016/j.ejps.2011.11.009

    Article  Google Scholar 

Download references

Acknowledgements

No outside funding was received for this study.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Robert Koprowski.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

RK suggested the algorithm for image analysis and processing, implemented it and analysed the images. MM, KW, ZW performed the acquisition of the tomography images and consulted the obtained results. All authors have read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Koprowski, R., Machoy, M., Woźniak, K. et al. Automatic method of analysis of OCT images in the assessment of the tooth enamel surface after orthodontic treatment with fixed braces. BioMed Eng OnLine 13, 48 (2014). https://doi.org/10.1186/1475-925X-13-48

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1475-925X-13-48

Keywords