The Fourier transform infrared (FT-IR) spectrum of a rock contains information about its constituent minerals. Using the wavelet transform, we roughly separate the mineralogical information in the FT-IR spectrum from the noise, using an extensive set of training data for which the true mineralogy is known. We ignore wavelet coefficients that vary too much among repeated measurements on rocks with the same mineralogy, since these are likely to reflect analytical noise. We also ignore those that vary too little across the entire training set, since they do not help to discriminate among minerals. We use the remaining wavelet coefficients as the data for the problem of estimating mineralogy from FT-IR data. For each mineral of interest, we construct an affine estimator <i>◯</i> of the mass fraction <i>x</i> of the mineral of the form <i>◯</i> = <i>a</i>·<i>w</i> + <i>b</i>, where <i>a</i> is a vector, <i>w</i> is the vector of retained wavelet coefficients, and <i>b</i> is a scalar. We find <i>a</i> and <i>b</i> by minimizing the maximum error of the estimator over the training set. When applying the estimator, we "truncate" to keep the estimated mineralogy between 0 and 1. The estimators typically perform better than weighted nonnegative least-squares.
You do not have subscription access to this journal. Cited by links are available to subscribers only. You may subscribe either as an OSA member, or as an authorized user of your institution.
Contact your librarian or system administrator
Login to access OSA Member Subscription