  minimally invasive dentistry is an approach that seeks to maintain the patient 's oral health with preventive measures and to treat possible disturbances of health as early as possible and with as little intervention as possible .
this requires that caries is detected at an early stage of development and that its status can be monitored frequently .
however , the current methods for diagnosing caries are able to detect caries only at a relatively advanced stage .
accordingly , methods for early detection of caries have been researched for the past twenty years .
many of these methods still require extensive research before they can be used in clinical practice .
optical caries diagnosis methods are based on the fact that caries cause changes in the tooth 's optical properties at an early stage of development .
this was a pilot study to investigate whether diffuse reflectance visible / near - infrared spectroscopy ( vis / nir - s ) can be used to detect dental caries lesions .
reflectance spectroscopy measures the intensity of light at several different wavelengths , that is , its spectra , after the light has reflected from the studied object .
diffuse reflectance refers to light that has been reflected from the inside of the object , rather than from its surface . in this
study the intensity was measured at wavelengths in the visible range and at wavelengths in the near - infrared range , covering wavelengths in the range 4201000  nanometers . within this range ,
the intensity was measured at 2305 different wavelengths , so that the difference between consecutive wavelengths was approximately 0.25  nm .
this study was limited to studying natural caries lesions that could be diagnosed with fiber - optic illumination , on smooth surfaces of extracted tooth .
a theory of caries diagnosis using near - infrared spectroscopy emerges from the previous studies of detecting caries lesions with near - infrared light [ 27 ] .
according to this theory , the development of a caries lesion increases the porosity of the affected tissue , which in turn leads to an increased scattering of light in the lesion .
wavelengths in the near - infrared range are considered better than the wavelengths in the visible range , because the former can penetrate deeper into the tissue and are less affected by stains on the tooth .
the purpose of this study was to provide additional evidence in support of this theory .
the dental services of the city of vaasa provided extracted human teeth for the study . the teeth were stored immersed in denatured alcohol in order to disinfect them and to keep them hydrated . before inspection and measurements ,
the teeth were inspected by the first author with fiber - optic illumination , after the technique was introduced to him by the second author , in order to detect healthy areas of enamel and areas of enamel that contained caries lesions . in total 21 teeth were used in the study .
a total of 109 points of enamel were measured on the teeth , consisting of 69 points which were thought to represent healthy enamel and 40 points which were thought to represent caries lesions .
each measurement point produced a spectra , a sample for the rest of the analysis . in pattern recognition terminology the diagnosis of a given sample , as either healthy or carious , is called the label of the sample .
the analysis of the samples tries to create a method which estimates the diagnosis , the label , of the sample based only on the measurements .
an optical fiber , placed in contact with the sample , conveys light from a light source to the sample .
another optical fiber is placed in contact with the sample at a small distance from the first fiber .
some fraction of the light which scatters inside the sample will eventually exit the sample so that it enters the second optical fiber .
it then gets conveyed to a spectrometer , which measures the spectra of the reflected light .
photonics describes the key properties with the absorption coefficient and the scattering coefficient of the material .
the measured spectra is analyzed in order to deduce information about the sample material .  
the measurements were made with a spectrometer hr4000 ( ocean optics inc . , dunedin , fl , usa ) and with a general purpose transmission dip probe model t300-rt - vis / nir ( ocean optics inc . ,
the probe contains two optical fibers , both with a diameter of 300  m , housed in a stainless steel assembly with a diameter of 3.175  mm .
one of the fibers is connected to a light source and brings light to the sample .
the light source used in this study was a tungsten halogen lamp hl-2000 ( ocean optics inc . , dunedin , fl , usa ) .
thus , the study had to be carried out with a probe that was readily available in our laboratory .
however , it was considered to be suitable for this study when the ferrule enclosing the inner assembly was removed , exposing the stainless steel assembly that houses the fiber optics .
the period of time for which the spectrometer collects light when it is making one measurement is called the integration time . in this study integration time
a longer integration time produces better measurement results than a short one , because the intensity of the collected light increases at all wavelengths , yielding a better signal - to - noise ratio ( snr ) .
therefore , the integration time is typically set as long as possible . however , if the intensity of the collected light at a given wavelength exceeds the measurement range of the spectrometer , the spectrometer saturates . in that case
the intensity can not be measured , and we know only that it exceeded the maximum measurable value . in order to make the measurement results comparable to results that would have been obtained with the same spectroscope using another light source or another integration time
this is done by measuring the smallest and the greatest intensity value that a measurement can produce with the given integration time for each wavelength and by scaling all other measurements to that range .
the lowest possible intensity values are obtained by measuring the so - called dark current , which is caused by thermal noise .
measuring a white reference sample produces the greatest possible intensity values . in this study , the integration time was set so that the white reference sample ( a white reference tile ws-2 , avantes inc . ,
a spectrometer must also be calibrated for its detector , so that its measurement results are comparable to those obtained by other spectroscopes .
this is done by measuring the spectra of a sample whose spectra is known . in this study
the used spectroscope was calibrated for its detector by the manufacturer as part of its construction .
a spectroscopic measurement result contains many small random errors , which are collectively called ( thermal ) noise .
when the dark current is subtracted from the spectra , the mean value of the effect caused by the noise is shifted to zero , and thus the effect of noise is observed as errors which have a normal distribution with a zero mean . in order to minimize the effect of noise in the samples , each point was measured one hundred times consecutively , and the resulting spectra were averaged .
this meant that the probe needed to stay as motionless as possible for two seconds .
however , a far shorter time period would have probably been sufficient . as a further measure against noise ,
the samples were smoothed by using the savitzky - golay method with a window length of 61 and sixth degree polynomials .
this method selects the coefficients of a sixth degree polynomial so that the polynomial is the best possible approximation for the measurement result , that is , the spectra , for the 30 wavelengths before a given wavelength and for the 30 wavelengths after it .
the value of the polynomial at the given wavelength replaces the measured intensity at that wavelength .
this removes , or smoothens , fast and small changes in the spectra , which are mainly caused by noise . a simple computational algorithm , based on exhaustive search ,
was then used to find a set of rules that could be used for detecting caries lesions . at this point ,
the goal was to classify the samples into two classes : points on healthy enamel ( healthy samples ) and points on caries lesions ( carious samples ) .
for this , a set of rules was searched for , so that every rule had the following format : if the sample 's normalized intensity at a given wavelength  is greater than ( or smaller than ) a given threshold i , the sample is classified as carious ; otherwise , the sample is classified as healthy .
thus , each rule had three parameters : the wavelength  , the intensity threshold i , and whether or not the threshold is an upper or lower limit for the intensity .
if , and only if , one or more of the rules classified the sample as carious , the sample was classified as carious . if none of the rules considered the sample as carious , it was classified as healthy .
it was hoped that the algorithm would select a set of rules which resembles the results found in earlier studies on this subject .
a number of wavelengths were selected from the range of available wavelengths ( 4201000  nm ) as options for parameter  in the search , so that the intervals between the wavelengths were equal and the first and the last wavelength were always selected as options .
the search was done with different numbers of wavelengths . for each of the selected wavelengths , the algorithm sorted the samples ' intensities at that wavelength and considered the midpoint between each two consecutive intensities as a possible threshold i in a rule .
the algorithm calculated the classification accuracy on the training set for each of the pairs  and i described above , using the threshold i first as an upper limit for classifying the sample as carious and then using it as a lower limit , and chose the values of  and i and the type of threshold , which gave the best accuracy ( see pseudo - code at pseudocode 4 ) .
after a rule had been selected this way , the algorithm selected another rule with the same method , so that the new rule gave the best possible accuracy when used together with the previously selected rule(s ) .
this was continued until the maximum allowed number of rules , here five rules , was reached , or until the classifier was unable to find a new rule which would improve the classification accuracy .
this algorithm , like every machine learning method , requires a set of samples which is used for searching for the rules and a separate set of samples which is used for evaluating the accuracy that is achieved with the resulting rules .
the former set of samples is called the training set and the latter set is called the validation set .
this may cause problems for the machine learning method when the samples are divided into a training set and a validation set , because some types of samples may become overrepresented in the training set , misleading the learning method as it tries to recognize what discerns the two classes from each other . in this study , this risk was alleviated by using a 4-fold cross - validation . in this method , the samples are divided into four groups , and one of them is used as the validation set while the other three groups form the training set .
each group in turn is used as the validation set , and the results from these four  folders  are averaged . this way
each sample is a part of the training set in three folders and a part of the validation set in one folder .
it is unlikely that the same types of samples would be overrepresented in all four training sets , unless the entire set of available samples has this problem .
a single training set which has this problem would stand out from the others , and the skewed learning results from it would be corrected by the results from the other training sets . while the small set of samples may still give a skewed representation of the kinds of samples which are being studied , the cross - validation seeks to minimize this problem . in this study
the averaging was done so , that a median rule set was constructed from the rules which the algorithm selected for the folders , and all of the samples were then classified with the median rule set .
median of the numbers of rules in the folders determined the number of rules in the median rule set .
some manual deliberation was used when constructing the rules of the set . for each rule in the final set , a temporary rule set was composed by selecting one rule from each folder 's rule set , so that the rules in the temporary set resembled each other , if that was possible given the available rules .
the median of the wavelengths used in the rules in the temporary rule set determined the wavelength for the rule in the median rule set .
the intensity threshold and the type of threshold were selected similarly for the rule in the median rule set .
each sample was diagnosed as either healthy or carious by the first author , and the selected rules estimated each sample to be either healthy or carious . based on these two properties ,
samples which were diagnosed as healthy and which were estimated to be healthy by the rules are called true negatives ( tns ) .
similarly , carious samples which were correctly estimated are called true positives ( tps ) .
a healthy sample which was estimated to be carious is called false positive ( fp ) and a carious sample which was estimated to be healthy is called a false negative ( fn ) .
these four values can be used to calculate the following five values which describe the accuracy of the selected rules .
positive predictive value ( ppv ) is the probability that the classifier , that is , the set of rules , is correct when it estimates a sample to be carious .
negative predictive value ( npv ) is the probability that a healthy estimate is correct .
accuracy is the fraction of the samples which were correctly estimated , that is , where the rules gave the correct answer .  
positive predictive value ( ppv ) is the probability that the classifier , that is , the set of rules , is correct when it estimates a sample to be carious .  
negative predictive value ( npv ) is the probability that a healthy estimate is correct .  
accuracy is the fraction of the samples which were correctly estimated , that is , where the rules gave the correct answer .  
after the classification rules had been selected and the samples had been classified according to them , there were fifteen samples which the author had diagnosed as carious but which were classified as healthy ( false negatives ) .
the spectra of these samples were virtually indistinguishable from the spectra of the healthy samples ( see figure 3(a ) ) , at least for the analysis methods used in this study .
thus , a hypothesis was made that these samples , the false negative cases , had been misdiagnosed by the author and subsequently mislabeled .
the rules that were selected by the algorithm suggested that a short wavelength , namely , 420  nm , was relatively useful in the diagnosis of caries .
therefore , another hypothesis was made , according to which a number of samples had been diagnosed by the author as carious while in fact the measured points were only stained and were thus false positive cases of the diagnosis , even if they had been classified correctly by the classifier .
notation i( ) refers to the normalized intensity of the spectra at wavelength . in other words , the sample was thought to represent a stain if it had a small scattering coefficient at both a long wavelength ( 815  nm , which is in the near - infrared range ) and a short wavelength ( 420  nm ) .
the samples , or the spectra of the measured points , are presented in figure 2 .
the number of wavelengths which were selected as options for the rule 's parameter  , that is , parameter wavelengthoptioncount , had only a small effect on the accuracy of the resulting median rule set .
when only the shortest wavelength ( 420  nm ) and the longest wavelength ( 1000  nm ) were available as options , the median rule set had an accuracy of 82% . with three wavelengths to choose from ,
when the number of options was between four and six , the accuracy was 85% . with greater numbers of wavelengths available , the accuracy was 84% .
this suggested that the rules depicted a phenomenon which was consistently present in all four folders .
when the number of options for the rules ' wavelengths was 15 , the median rule set indicated that a sample is carious if , and only if , i(  420 ) 
a confusion matrix of the classification accuracy that is achieved with these rules is presented in table 1 , showing that these rules reached an accuracy of 84% .   as can be seen in figure 3(a ) and in table 1 , there were fifteen carious samples which were classified as healthy ( false negatives ) , and whose spectra was virtually indistinguishable from the spectra of the healthy samples .
as explained in section 2.4 , this leads to a hypothesis that these samples had been misdiagnosed and subsequently mislabeled , and that they therefore represented healthy samples and were in fact classified correctly .  
according to the theory on optical caries diagnosis , an elevated intensity in the near - infrared range is the best indication of a dental caries lesion .
however , the rules selected by the search algorithm indicated that a short wavelength , namely , 420  nm , was relatively useful in the diagnosis of caries . as explained in section 2.4 ,
another hypothesis was thus made , according to which a number of samples had been diagnosed as carious while in fact they were only stained .
all samples that were identified as stained had been diagnosed and classified as carious , and thus appeared to be true positive cases .
these suspected misdiagnoses had not lowered the apparent accuracy of the classification , but they may have caused the rule set to erroneously consider stains as caries lesions .
when the search algorithm was run again , giving 15 options for the parameter  , after first relabeling the fifteen false negative cases as healthy samples ( first hypothesis ) and then relabeling the eight suspected stains as healthy samples ( second hypothesis ) , the algorithm selected only one rule in every cross - validation folder .
all rules set an upper limit for the normalized intensity at a wavelength in the near - infrared range .
the median of those rules was i(  791 )  0.3255 , which is consistent with the theory .
this study suffers from a small number of samples . although the study used 109 measurements , they were taken from only 21 individual teeth
this fact is significant , because it is probable that samples taken from a single tooth resemble each other more than samples taken from different teeth or from different patients .
fifteen of those measurements were considered to be misdiagnosed by the first hypothesis , and further eight measurements were considered to be misdiagnosed by the second hypothesis . therefore , further study is needed to increase the reliability of the accuracy estimate of this method .
the measurement results together with the theory on the topic suggest that many of the measurements which were supposedly made from a caries lesion are in fact made from healthy enamel , which was in some cases stained . when we make these suggested corrections to the labeling of the samples , the samples seem to fit well to the theory and the samples can easily be accurately classified
. these kinds of diagnostic mistakes , or false positive diagnoses , are a credible explanation , because the diagnoses were made by a novice on the subject .
however , such corrections also pose a risk that the measurement results are relabeled to make them fit the theory , which would inflate the accuracy of the method .
the composition of the dental tissues varies from tooth to tooth and between different sites of a given tooth .
as can be seen in figure 2 , the spectra of the different healthy samples vary quite a bit , especially at the visible wavelengths .
this suggests that the threshold intensity or intensities for diagnosing a suspected lesion as carious might also vary similarly . in order to compensate for the inter - tooth and intra - tooth variance
, we might consider measuring the average spectra for a given tooth by measuring several points on the tooth surface , that is , by scanning the surface and by evaluating how much the spectra of the suspected lesion differ from the tooth 's average .
unfortunately , this approach could potentially make this method less effective for its original purpose .
the method is being developed for the detection of caries lesions at an early stage of development .
thus , the dentist does not necessarily notice all of the lesions which are detected by the device .
if the inspected tooth surface contains several developing caries lesions , the average spectra of the surface could be something in between the healthy enamel and the carious enamel , making the lesions appear too similar to the average surface to be diagnosed as carious .
the scanning method would also make it rather awkward to inspect several teeth per patient .
quantitative light - induced fluorescence ( qlf ) and laser - induced fluorescence ( lf ) are two optical methods for the detection of caries lesions .
they are based on fluorescence , or the phenomenon that when the tooth sample is illuminated with a light source , some of the light is absorbed in the sample , after which the sample emits light at a longer wavelength . for both methods
study the sample was considered carious if the measured intensity was greater than a fixed threshold , i(  791 )  0.3255 .
the proposed explanation is that the increased scattering due to caries causes more light to be reflected to the measuring fiber optic .
 > 520  nm because increased scattering due to caries interferes with the detection of the fluorescence , and lf expects to find increased intensity at the near - infrared range caused by fluorescence from organic molecules in the sample .
since the samples in this study were stored in denatured alcohol , they were probably relatively free of organic molecules .
further study is required to determine whether the fluorescence from organic molecules , that is , the phenomenon measured by lf , interferes with the detection method outlined in this study , especially for in vivo measurements .
if it does interfere , it probably makes the method more eager to label a sample as carious , thus increasing its sensitivity and reducing its specificity .
this effect may be modified , at least in part , by selecting a new set of rules based on results from in vivo measurements . incidentally ,
in contrast , authors of this study felt that the method outlined in this paper helped them to increase specificity .
an ability to measure the amount of dental tissue lost to caries could be pursued by inducing caries in vitro to a tooth sample ( see [ 6 , 17 , 18 ] ) so that the amount of mineral dissolved from the tooth could be measured without destroying the sample , and by measuring the spectra of the sample at varying degrees of mineral loss .
one possible method for this would be to cycle the tooth sample in de- and remineralization solutions and to measure the amount of mineral dissolved to the solutions with a mass spectrometer .
finding a method to calculate the amount of the mineral loss from the spectra would be a regression problem .
it seems that spectroscopic measurements can help to reduce false positives at in vitro setting , including those caused by stains .
however , the work reported in this paper was a pilot study , and further research is required to evaluate the strength of the evidence for the method 's performance at in vitro setting and to extend the measurements to in vivo setting .
 
 background . 
 a caries lesion causes changes in the optical properties of the affected tissue . 
 currently a caries lesion can be detected only at a relatively late stage of development . 
 caries diagnosis also suffers from high interobserver variance 
 . methods . 
 this is a pilot study to test the suitability of an optical diffuse reflectance spectroscopy for caries diagnosis . 
 reflectance visible / near - infrared spectroscopy ( vis / nirs ) was used to measure caries lesions and healthy enamel on extracted human teeth . 
 the results were analysed with a computational algorithm in order to find a rule - based classification method to detect caries lesions . 
 results . 
 the classification indicated that the measured points of enamel could be assigned to one of three classes : healthy enamel , a caries lesion , and stained healthy enamel . 
 the features that enabled this were consistent with theory . conclusions . 
 it seems that spectroscopic measurements can help to reduce false positives at in vitro setting . 
 however , further research is required to evaluate the strength of the evidence for the method 's performance . 
1. Introduction
2. Methods
3. Results
4. Discussion
5. Conclusions