Quantitative Analysis of Mixed Minerals with Finite Phase Using Thermal Infrared Hyperspectral Technology

It is crucial but challenging to detect intermediate or end products promptly. Traditional chemical detection methods are time-consuming and cannot detect mineral phase content. Thermal infrared hyperspectral (TIH) technology is an effective means of real-time imaging and can precisely capture the emissivity characteristics of objects. This study introduces TIH to estimate the content of potassium salts, with a model based on Competitive Adaptive Reweighted Sampling (CARS) and Partial Least Squares Regression (PLSR). The model takes the emissivity spectrum of potassium salt into account and accurately predicts the content of Mixing Potassium (MP), a mineral mixture produced in Lop Nur, Xinjiang. The main mineral content in MP was measured by Mineral Liberation Analyzer (MLA), mainly including picromerite, potassium chloride, magnesium sulfate, and less sodium chloride. 129 configured MP samples were divided into calibration (97 samples) and prediction (32 samples) sets. The CARS-PLSR method achieved good prediction results for MP mineral content (picromerite: correlation coefficient of correction set (Rp2) = 0.943, predicted root mean square error (RMSEP) = 2.72%, relative predictive deviation (RPD) = 4.24; potassium chloride: Rp2 = 0.948, RMSEP = 2.86%, RPD = 4.42). Experimental results convey that TIH technology can effectively identify the emissivity characteristics of MP minerals, facilitating quantitative detection of MP mineral content.


Introduction
During the mineral production process, the composition of intermediate and terminal products can change within a certain range due to factors such as alterations in the original ore's nature, production operations, or other variables. In this type of product, although the mineral content varies within a certain range, the types of minerals are basically unchanged and can be referred to as "finite phase minerals". Not adjusting production parameters in a timely manner based on actual production conditions can cause production instability, leading to a waste of resources and increased economic costs. As a result, it is essential to conduct intermittent sampling and tests to gather production information that can aid in the production process [1]. The chemical volumetric method is the primary analytical technique in production. In this method, minerals are dissolved, and specific chemical reagents are used to identify relevant elements. Although this method can meet production requirements to some extent, the low detection speed and incapable of detecting mineral phases cannot be neglected. Other mineral quantitative detection methods [2,3], such as hyperspectral to accurately detect salt mineral content [29], it is promising in detecting salt minerals that have a tendency to form dense mixtures. Based on this, this study utilized TIH technology to acquire and quantitatively analyze the emissivity data of potassium salt. The contributions of this study are as follows: (1) TIH technology was employed to obtain thermal infrared spectral characteristic information of potassium salt; (2) The thermal infrared spectrum characteristics of pure minerals in potassium salt were obtained and the potassium salt grade was detected by CARS-PLSR method.

Technical Route
In this study, we aim to figure out the feasibility and accuracy of the TIH technique in the application of the mineral phases of potassium salts detection. The potash samples used in the study came from Xinjiang Lop Nur Potash Co., China's largest potash sulfate producer. The main research process ware as follows: (1) Based on the actual MP samples, mixed potassium samples were designed and prepared by adding varying amounts of Pure minerals. A batch of potassium salt samples were obtained with sufficient and known mineral phase composition. (2) The samples were standardized by the pellet method, which stabilized the roughness of each sample and reduced the "noise" in the spectrum data. Subsequently, the spectral data of the potassium salt samples were collected using the Hyper-Cam hyperspectral imager in an open environment. (3) The TIH data for potassium salt was processed to obtain its emissivity data. Based on this emissivity data, prediction models were established for the mineral phase composition of potassium salt using the CARS and PLSR methods. (4) The prediction results of the model were evaluated and analyzed. The flowchart of the experimental and modeling process is shown in Figure 1.
contents in soil, there is currently limited literature available on the identification of various salt contents [14,17]. Since TIH technology could capture feature information with hyperspectral to accurately detect salt mineral content [29], it is promising in detecting salt minerals that have a tendency to form dense mixtures. Based on this, this study utilized TIH technology to acquire and quantitatively analyze the emissivity data of potassium salt. The contributions of this study are as follows: (1) TIH technology was employed to obtain thermal infrared spectral characteristic information of potassium salt; (2) The thermal infrared spectrum characteristics of pure minerals in potassium salt were obtained and the potassium salt grade was detected by CARS-PLSR method.

Technical Route
In this study, we aim to figure out the feasibility and accuracy of the TIH technique in the application of the mineral phases of potassium salts detection. The potash samples used in the study came from Xinjiang Lop Nur Potash Co., China's largest potash sulfate producer. The main research process ware as follows: (1) Based on the actual MP samples, mixed potassium samples were designed and prepared by adding varying amounts of Pure minerals. A batch of potassium salt samples were obtained with sufficient and known mineral phase composition. (2) The samples were standardized by the pellet method, which stabilized the roughness of each sample and reduced the "noise" in the spectrum data. Subsequently, the spectral data of the potassium salt samples were collected using the Hyper-Cam hyperspectral imager in an open environment. (3) The TIH data for potassium salt was processed to obtain its emissivity data. Based on this emissivity data, prediction models were established for the mineral phase composition of potassium salt using the CARS and PLSR methods. (4) The prediction results of the model were evaluated and analyzed. The flowchart of the experimental and modeling process is shown in Figure 1.

MP Samples
Based on a random survey of 738 MP mineral phase data from the Lop Nur Potash production site in the first half of 2021, 138 MP samples of 5 g each were designed and

MP Samples
Based on a random survey of 738 MP mineral phase data from the Lop Nur Potash production site in the first half of 2021, 138 MP samples of 5 g each were designed and configured. The 400 g MP sample from the production site of potassium salt (dried at 50 • C, mixed evenly by shifting cone method) was used as a master sample. The mineral phases and contents of the master sample were accurately analyzed with multiple detection methods. The 138 MP samples were then configured by adding one or more salts (contained in the MP) to the master sample.

Mineral Facies Analysis
By comprehensive analysis of X-ray diffraction analysis (XRD) [30], MLA [31], Chemical titration, and ASS [32], the types and specific contents of various minerals in the MP master sample were defined. The specific parameters of the test method are as follows: XRD: The D8 Advance X-ray diffractometer with Cu target from Bruker, Germany was used for XRD testing. The acceleration voltage is 40 kV and the anode current is 150 mA. Wide-angle diffraction patterns are collected in the range of 5-70 • (2θ) at a scanning speed of 5 • /min. MLA: The system consists of a ZEISS Sigma model 300 high-resolution field emission scanning electron microscopy (FESEM), a Bruker XFlash 6|60 type X-ray spectrometer (EDS), and a set of AMICS software (including 4 subroutines, AMICSTool, Investigator, MineralSTDManager, and AMICSProcess). Experimental conditions: accelerating voltage 20 kV, High vacuum mode, a working distance of about 8.5 mm, backscattered electron detector (HDBSD), High Current mode. The sample was sprayed before testing to eliminate the charge on the surface of the non-conductive sample. The samples were polished using saturated MP mother liquor as a polish.
Chemical Titration: Titration of K + : Sodium tetraphenylboron-quaternary amine salt volumetric method. Determination principle: Potassium ions react with sodium tetraphenolate boron to form potassium tetraphenolate boron deposits, and excess sodium tetraphenolate boron reacts with quaternary amine salt to form insoluble quaternary amine salt sodium tetraphenolate boron double salt.
Titration of Mg 2+ : EDTA complexometric titration. Determination principle: The cation in EDTA complexes with the calcium and magnesium cation in the test solution, so that the color of the solution changed from tartar red to azure blue, which was the endpoint.
Titration of SO 2− 4 : Moore method. Determination principle: Potassium chromate as indicator was titrated directly with silver nitrate standard solution in neutral weak alkaline solution.
Titration of Cl -: Barium chloride volumetric method. Determination principle: In an acidic solution (pH = 2~3), SO 2− 4 and BaCl 2 generated excessive barium salt precipitation and alizarin red~S indicator to generate a yellow complex, titrated with a known standard solution, adding ethanol to improve the sensitivity of the endpoint of the titration.
AAS: Titration of Na + : According to the general rules of atomic absorption spectrometry, the sodium ions in the solution were analyzed by using the German continuous light source atomic absorption spectrometer at the pressure of 0.3 MPa.

Mineral Composition of MP Samples
XRD result in Figure 2 shows that the main components of MP are picromerite (KMg(SO 4 ) 2 ·6H 2 O) and potassium chloride (KCl). In the MLA test results shown in Figure 3, the element Mg:S:K in the red region has a mole ratio of approximately 1:2.03:1.85, indicating that picromerite is the primary mineral, followed by a small amount of magnesium sulfate (MgSO 4 ·7H 2 O). The mole ratio of element K:Cl in the yellow region is about 1.01:1.00, which is expressed as KCl. The mole ratio of the element Na:Cl in blue is about 1.00:1.00, represented by sodium chloride (NaCl). The MLA detection result shows that the MP sample contains a large amount of picromerite and sodium chloride in addition to less potassium sulfate and potassium chloride.  However, during sample preparation for the MLA test, the MP sample was partially dissolved due to the requirement to wash and polish the sample simultaneously, despite the saturated MP mother liquor being used as a detergent. As shown in Table 1, 8.18% of the MLA test field is not uniform and could not be measured, which means, 5.37% of the MP sample could not be identified. Therefore, chemical volumetric methods were used to supplement the MLA results. The result of the chemical volumetric test in Figure 3 shows that the mass of potassium ion is 27.18%, a magnesium ion is 4.30%, sulfate ion is 33.67%, chloride ion is 13.57%, and sodium ion is 0.32%. Finally, the composition of the 400 g MP sample can be inferred as follows: picromerite, 65.70%, potassium chloride, 27.50%, magnesium sulfate, 5.97%, sodium chloride, 0.81%, combined with the ionic content data from the chemical titration in Table 2.  However, during sample preparation for the MLA test, the MP sample was partially dissolved due to the requirement to wash and polish the sample simultaneously, despite the saturated MP mother liquor being used as a detergent. As shown in Table 1, 8.18% of the MLA test field is not uniform and could not be measured, which means, 5.37% of the MP sample could not be identified. Therefore, chemical volumetric methods were used to supplement the MLA results. The result of the chemical volumetric test in Figure 3 shows that the mass of potassium ion is 27.18%, a magnesium ion is 4.30%, sulfate ion is 33.67%, chloride ion is 13.57%, and sodium ion is 0.32%. Finally, the composition of the 400 g MP sample can be inferred as follows: picromerite, 65.70%, potassium chloride, 27.50%, magnesium sulfate, 5.97%, sodium chloride, 0.81%, combined with the ionic content data from the chemical titration in Table 2. the red area is picromerite, the dark gray area is low counting rate, the yellow area is potassium chloride, the blue area is sodium chloride, and light gray area is an unknown mineral.
However, during sample preparation for the MLA test, the MP sample was partially dissolved due to the requirement to wash and polish the sample simultaneously, despite the saturated MP mother liquor being used as a detergent. As shown in Table 1, 8.18% of the MLA test field is not uniform and could not be measured, which means, 5.37% of the MP sample could not be identified. Therefore, chemical volumetric methods were used to supplement the MLA results. The result of the chemical volumetric test in Figure 3 shows that the mass of potassium ion is 27.18%, a magnesium ion is 4.30%, sulfate ion is 33.67%, chloride ion is 13.57%, and sodium ion is 0.32%. Finally, the composition of the 400 g MP sample can be inferred as follows: picromerite, 65.70%, potassium chloride, 27.50%, magnesium sulfate, 5.97%, sodium chloride, 0.81%, combined with the ionic content data from the chemical titration in Table 2.

Samples Library of MP
As shown in Figure 4a, among the 738 MP samples investigated, the content of picromerite ranged from 30.73% to 76.79%, with an average of 51.20% and a median of 51.19%. In 138 manually configured MP samples, the content parameters of picromerite ranged from 30.46% to 74.49%, with an average of 53.91% and a median of 53.43%. The content distribution of 138 MP samples was similar to that of 738 MP samples data. While Table 3 shows that potassium chloride, magnesium sulfate, and sodium chloride are also similar, which indicates 138 MP samples manually configured can completely meet the needs of field production.

Standardized Processing of MP Samples
The MP samples were then slightly ground and placed in 50 mL centrifuge tubes and mixed using a vortex oscillating mixer at 300 r/min for 6 min. Finally, the MP samples were successively put into a square mold with an inner diameter of 2 cm and pressed for 30 s under a pressure equivalent to 6 t to complete sample preparation.

Hyperspectral Imaging System and Image Acquisition
In this study, a Hyper-Cam hyperspectral imager (Telops, Canada) was used for data acquisition, and its imaging mode was Fourier interference imaging. The spectral range of the instrument is 7.7-11.8 μm, the spectral resolution range is 0.25-150 cm −1 , the Field of View (FOV) (°) is 40, and the imaging field is 320 × 256 array.

Standardized Processing of MP Samples
The MP samples were then slightly ground and placed in 50 mL centrifuge tubes and mixed using a vortex oscillating mixer at 300 r/min for 6 min. Finally, the MP samples were successively put into a square mold with an inner diameter of 2 cm and pressed for 30 s under a pressure equivalent to 6 t to complete sample preparation.

Hyperspectral Imaging System and Image Acquisition
In this study, a Hyper-Cam hyperspectral imager (Telops, Canada) was used for data acquisition, and its imaging mode was Fourier interference imaging. The spectral range of the instrument is 7.7-11.8 µm, the spectral resolution range is 0.25-150 cm −1 , the Field of View (FOV) ( • ) is 40, and the imaging field is 320 × 256 array.
As shown in Figure 5, the Hyper-Cam hyperspectral imaging system was used to collect TIH data from samples and a diffuse reflector (gold-plated plate) in a relatively open environment. The open environment can effectively avoid the interference of other objects in the surrounding environment and reduce the noise of acquiring data. At the same time, the sample was preheated with a heating table, and the heating temperature was set to 50°C. Fluke 54-IIB dual channel contact thermometer was used to record the temperature of the sample and gold foil before and after imaging, and the mean value was taken as the temperature of the sample and gold foil. Then, the absolute value of the difference between the measured temperature and the inverted temperature was used to evaluate the quality of the sample spectral data. It is generally believed that if the absolute value of the temperature difference is less than 1 K, the quality of collected sample data is acceptable and can be used for accurate quantitative analysis and study.

Emissivity Inversion
The method of obtaining sample emissivity in this study refers to the study [33]. Iterative Spectrally Smooth Temperature and Emissivity Separation (ISSTES) algorithm is widely used in remote sensing applications to accurately calculate the emissivity of solid objects from their spectral radiance measurements. The ISSTES takes advantage of the fact that the emissivity spectrum of solid objects is smoother than the characteristics of the atmosphere. ISSTES algorithm defines a group of emissivity curves indexed by surface temperature as follows: where L sensor (λ) is the total radiance reaching the sensor, λ is the central wavelength of each channel (µm), ε(λ) is the emissivity, and E(λ) is environmental radiance measured by the gold-plated plate. τ(λ) and Lu(λ) are atmospheric transmittance and upwelling radiation, respectively, which are simulated by MODTRAN ® 5.2. B(λ,T) is the Planck radiance at panel

Emissivity Inversion
The method of obtaining sample emissivity in this study refers to the study [33]. Iterative Spectrally Smooth Temperature and Emissivity Separation (ISSTES) algorithm is widely used in remote sensing applications to accurately calculate the emissivity of solid objects from their spectral radiance measurements. The ISSTES takes advantage of the fact that the emissivity spectrum of solid objects is smoother than the characteristics of the atmosphere. ISSTES algorithm defines a group of emissivity curves indexed by surface temperature as follows: where Lsensor(λ) is the total radiance reaching the sensor, λ is the central wavelength of each channel (μm), ε(λ) is the emissivity, and E(λ) is environmental radiance measured by the gold-plated plate. τ(λ) and Lu(λ) are atmospheric transmittance and upwelling radiation, respectively, which are simulated by MODTRAN ® 5.2. B(λ,T) is the Planck radiance at panel temperature T. The units for Lsensor(λ), Lu(λ), B(λ, T) and E(λ), are all W/(m²·sr·μm). ε(λ) and τ(λ) are dimensionless quantities.
As can be seen from Equation (1), given the surface temperature and known atmospheric parameters, the corresponding emissivity spectrum can be calculated. Then, if the smoothness criterion of the emissivity curve is determined, the smoothness of a series of emissivity curves at different temperatures near the real surface temperature can be calculated, and the temperature corresponding to the estimated emissivity spectrum is the smoothest result. The smoothness criterion used by ISSTES is designed as follows: where S is the calculated smoothness criterion value, and N is the frequency band number of hyperspectral data.
In this study, the pretreatment of TIH data as well as the CARS and PLSR algorithms were calculated using MATLAB 2019a. As can be seen from Equation (1), given the surface temperature and known atmospheric parameters, the corresponding emissivity spectrum can be calculated. Then, if the smoothness criterion of the emissivity curve is determined, the smoothness of a series of emissivity curves at different temperatures near the real surface temperature can be calculated, and the temperature corresponding to the estimated emissivity spectrum is the smoothest result. The smoothness criterion used by ISSTES is designed as follows: where S is the calculated smoothness criterion value, and N is the frequency band number of hyperspectral data. In this study, the pretreatment of TIH data as well as the CARS and PLSR algorithms were calculated using MATLAB 2019a.

Data Quality Evaluation by Temperature
In this study, the accuracy of the emissivity data was evaluated indirectly by comparing the measured temperature of the sample with the temperature obtained through inverted emissivity data using the algorithm. Since temperature and emissivity are coupled. When the emissivity error is about 0.015, the temperature error is about 1 k [34]. By comparing the measured and inverted temperatures, any discrepancies or errors in the emissivity data could be identified, thereby improving the accuracy of the spectral data and subsequent analysis. The difference between the measured and inverted temperatures was used as an indicator of the quality of the inverted data. This approach is crucial for the successful analysis and interpretation of the TIH data [35]. Figure 6a illustrates that the temperature difference of all samples is within ±2 K. In Figure 6b, 71.0% of the samples are in the range of absolute temperature difference ≤1 K, which indicates that the data collection and data processing process of MP samples is effective and the emissivity data obtained are of high quality. These emissivity results can be used for quantitative analysis and research of MP. and research of MP.
As can be seen from the temperature distribution of some MP samples in Figure 6c, the surface emissivity of MP samples is relatively uniform, indicating that MP samples are fully and evenly mixed during preparation. In addition, the mean emissivity of the 15 × 20 pixel at the center of the MP sample surface is taken as the emissivity of the corresponding sample to further reduce the effect of systematic errors in sample preparation.  As can be seen from the temperature distribution of some MP samples in Figure 6c, the surface emissivity of MP samples is relatively uniform, indicating that MP samples are fully and evenly mixed during preparation. In addition, the mean emissivity of the 15 × 20 pixel at the center of the MP sample surface is taken as the emissivity of the corresponding sample to further reduce the effect of systematic errors in sample preparation.

Prediction Model
Firstly, 138 samples were tested for outliers and removed. To avoid some unreliable samples in the modeling calculation [36], the 138 samples were sorted and classified according to their picromerite content and then divided into 5 groups by regular interval sampling. Then, 4 groups of samples were randomly selected as the correction set, and the remaining group of samples was the prediction set. PLSR calculations were performed 5 times. Therefore, each sample had 4 calibration residual errors (CRE) and 1 predicted residual error (PRE). Then, PLSR [37] was used for analysis and calculation, and the CARS [38] method was introduced to choose feature wavelength. The PLSR model and the CARS-PLSR model were constructed separately. CARS were employed to screen out the sensitive bands corresponding to the minerals, while PLSR is utilized to establish a linear relationship between the mineral-sensitive bands and mineral content. The resulting model was then employed to estimate the content of salt ore. The PLSR method, by projection, obtains orthogonal feature vectors for independent and dependent variables, respectively, and then establishes a unitary linear regression relation between the feature vectors. This approach not only resolves the issue of collinearity but also emphasizes the importance of independent variables in explaining and predicting dependent variables during the selection of feature vectors. By selecting several key variables, this method enables the model to achieve the minimum number of variables necessary for optimal performance. For the CARS method, the absolute value of the regression coefficient of the PLSR model is used as the index to evaluate the importance of each wavelength. In each sampling run, a fixed proportion of samples is randomly selected to establish the correction model. The number of samples selected in this study is 80% of the number of calibrations set. The wavelength for selection bands with larger absolute regression coefficients in the PLSR model is defined through adaptive reweighting sampling (ARS) technology. That is, CARS can often find the optimal combination of some key wavelengths that explain the chemical properties of interest.

Model Evaluation
To fully verify the performance of the model, the samples in the MP sample bank were divided into a calibration set and a testing set. The calculated parameters for the calibration set include the correlation coefficient (R 2 c ) and root mean square error (RMSEC). The calculated parameters for the prediction set include the correlation coefficient (R 2 p ), root mean square error (RMSEP), and relative prediction deviation (RPD). The higher R 2 c , R 2 p , and RPD, the lower RMSEC, and RMSEP, the better the predictive effect of the model. If R 2 c is close to R 2 p , or RMSEC is close to RMSEP, the model is robust. In particular, an RPD value between 1.8 and 2.0 indicates a good model, between 2.0 and 2.5 indicates a very good model, and above 2.5 indicates excellent model performance [39].
In Formulas (3)-(5), y i represents the predicted sample mineral content, y i represents the mean value of the actual sample mineral content, y i represents the actual sample mineral content, and n represents the number of samples involved in the calculation.

Outlier Detection
The outlier observation map of the MP sample was established using MCRE with abscissa and PRE with ordinate, as shown in Figure 7. The mean of corrected residual error (MPCE) and predicted residual error (PRE) were taken as the reference, and the larger samples were eliminated. Among them, the points with large MCRE and PRE values are removed. The samples with numbers 108, 84, and 131 are removed with picromerite; samples with numbers 21, 68, 84, 119, and 131 are removed with potassium chloride; samples with number 44 and 55 are removed with magnesium sulfate; samples with number 94 and 106 are removed with sodium chloride. That is, a total of 9 outliers are detected, and the remaining 129 points are used for subsequent model calculation.
larger samples were eliminated. Among them, the points with large MCRE and PRE values are removed. The samples with numbers 108, 84, and 131 are removed with picromerite; samples with numbers 21, 68, 84, 119, and 131 are removed with potassium chloride; samples with number 44 and 55 are removed with magnesium sulfate; samples with number 94 and 106 are removed with sodium chloride. That is, a total of 9 outliers are detected, and the remaining 129 points are used for subsequent model calculation.

Calibration Set and Prediction Set
As shown in Table 4 and Figure 8, a statistical analysis of the mineral composition of the four MP in the prediction and calibration sets is performed in terms of mean, standard deviation, minimum, median, and maximum. It illustrates that the range of sample parameters of the prediction set are all within the range of the calibration set. Moreover, from the two statistical terms mean and median, it can be found that the distribution of

Calibration Set and Prediction Set
As shown in Table 4 and Figure 8, a statistical analysis of the mineral composition of the four MP in the prediction and calibration sets is performed in terms of mean, standard deviation, minimum, median, and maximum. It illustrates that the range of sample parameters of the prediction set are all within the range of the calibration set. Moreover, from the two statistical terms mean and median, it can be found that the distribution of sample parameters of the prediction set and the calibration set is close, which means the results of the prediction set can fully illustrate the model performance of the calibration set.

MP Content Prediction
In the calibration set, PLSR conducts internal validation of the model based on leaveone-out validation, and the calculation results are shown in Table 5 and Figure 9. PLSR is the calculation method of multivariate linear relationship analysis, and the relationship between the percentage content of the four minerals in the MP is not clear. Therefore, univariate PLSR (U-PLSR) and multivariate PLSR (U-PLSR) are used for calculation. From the comparison of RPD values of the prediction set, the prediction accuracy of U-PLSR is higher than that of M-PLSR. The RPD values of four minerals of U-PLSR are 4.59, 4.06, 2.16, and 1.57, respectively, and those of four minerals of M-PLSR are 4.00, 4.06, 2.08, and 1.57, respectively. In the process of PLSR calculation, the components extracted from the joint composition matrix Y of the four minerals cannot represent the various minerals themselves, which deduced that there should be no correlation between the percentage contents of the four minerals. sample parameters of the prediction set and the calibration set is close, which means the results of the prediction set can fully illustrate the model performance of the calibration set.

MP Content Prediction
In the calibration set, PLSR conducts internal validation of the model based on leaveone-out validation, and the calculation results are shown in Table 5 and Figure 9. PLSR is the calculation method of multivariate linear relationship analysis, and the relationship

Sensitive Bands Selection of Pure Mineral
The CARS method is implemented through PLSR calculation of a single dependent variable, which can screen out the respective sensitive wavelengths of each mineral in MP samples. The CARS method could simplify the entire prediction method and enhance the feasibility of application in realistic scenarios of potassium salts. Table 6 shows that there In the calibration set of U-PLSR, the R 2 c of picromerite and potassium chloride is larger, while the R 2 c of magnesium sulfate and sodium chloride is smaller, which may be related to the ratio of four minerals in the MP sample. It means that the content of picromerite and potassium chloride in the calibration set is greater than that of magnesium sulfate and sodium chloride. In MP samples, the signal strength may be more pronounced when the content of a particular mineral is larger, and the final prediction accuracy is higher. In the U-PLSR, the RMSEP values of the four minerals are 2.52%, 3.12%, 1.17%, and 0.35%, respectively, which increase by 0.52%, 0.76%, 0.34%, and 0.14% compared with the RMSEC values. The results indicate that the prediction set is accurate and the PLSR model is robust to the mineral content prediction of MP.

Sensitive Bands Selection of Pure Mineral
The CARS method is implemented through PLSR calculation of a single dependent variable, which can screen out the respective sensitive wavelengths of each mineral in MP samples. The CARS method could simplify the entire prediction method and enhance the feasibility of application in realistic scenarios of potassium salts. Table 6 shows that there are 12 sensitive wavelengths for picromerite, 11 sensitive wavelengths for potassium chloride, 11 sensitive wavelengths for magnesium sulfate, and 8 sensitive wavelengths for sodium chloride. Figure 10 maps and compares the sensitive wavelengths of each mineral, pure mineral, and MP mean emissivity spectrum data.  9.13, 9.53, 9.62, 9.98, 10.17, 11.01 Potassium chloride 11 8.48, 8.72, 8.86, 9.21, 9.53, 9.58, 9.62, 9. 9.84, 9.98, 11.01, 11.30, 11.36 The sensitive wavelength distribution of picromerite is relatively uniform and distributed in the whole wavelength, among which there are four sensitive wavelengths concentrated in the trough of 8.5-9.0 µm. For the sensitive wavelengths of magnesium sulfate, the distribution of sensitive wavelengths is far away from the trough at 9.1 um, except for one near the trough at about 8.9 µm. The emissivity curves of potassium chloride and sodium chloride are similar in shape, with a broad peak between 8.5 and 9.5 µm, but the sensitive wavelengths do not overlap. The sensitive wavelength of potassium chloride is uniformly distributed over the entire thermal infrared band, while the sensitive wavelength of sodium chloride is mainly distributed on both sides of the thermal infrared band.
Since the sensitive wavelengths selected by the CARS method for four minerals in MP overlaps, U-PLSR and M-CARS-PLSR methods are used to establish prediction models for U-CARS-PLSR and M-CARS-PLSR (the input wavelength is the union of the sensitive wavelengths of four minerals, 31 wavelengths), as shown in Figure 11 and 9. 13, 9.53, 9.62, 9.98, 10.17, 11.01 Potassium chloride 11 8.48, 8.72, 8.86, 9.21, 9.53, 9.58, 9.62, 9.66, 10.22, 10.68, 11.07 Magnesium sulfate 11 8.48, 8.51, 8.94, 9.45, 9.62, 9.93, 10.12, 10.42, 10.68, 10.73, 10.84 Sodium chloride 8 8.34, 8.54, 8.58, 9.84, 9.98, 11.01, 11.30, 11.36  Since the sensitive wavelengths selected by the CARS method for four minerals in MP overlaps, U-PLSR and M-CARS-PLSR methods are used to establish prediction models for U-CARS-PLSR and M-CARS-PLSR (the input wavelength is the union of the sensitive wavelengths of four minerals, 31 wavelengths), as shown in Figure 11 and Table 7. In the calibration set, the values of R c 2 and RMSEC of the U-CARS-PLSR model are higher than those of the M-CARS-PLSR for picromerite and potassium chloride, while the results reverse for magnesium sulfate and sodium chloride. In the prediction set, the values of RMSEP and RPD for picromerite are almost the same (RMSEP: 2.72%, 2.78%, and RPD: 4.24, 4.15; RMSEP of sodium chloride: 0.40%, 0.40%, and RPD: 1.37, 1.36), but the results of potassium chloride and magnesium sulfate were significantly different (RMSEP of potassium chloride: 2.86%, 3.16%, and RPD: 4.42, 4.00; RMSEP of magnesium sulfate: 1.45%, 1.29%, and RPD: 1.73, 1.96).  Considering that the percentage of magnesium sulfate is much smaller than potassium chloride (standard deviation: potassium chloride, 12.66%; magnesium sulfate: 2.52%), it can be concluded that the performance of U-CARS-PLSR is better than that of M-CARS-PLSR. As Figure 12 shows, four minerals of RPD values of U-PLSR and U-CARS-PLSR are compared: picromerite 4.59, 4.24; potassium chloride, 4.06, 4.42; magnesium sulfate 2.16, 1.73; sodium chloride 1.37, 1.57. The results deduced that the prediction accuracy of U-CARS-PLSR has decreased to a certain extent except for potassium chloride. However, the CARS method reduced the number of selected calculated wavelengths (from 69 to 31 wavelengths), and the RPD values of major minerals in MP samples are both greater than 4.00, indicating that the prediction effect was still acceptable. In other words, CARS have a great role in reducing the cost of TIH technology and then being applied to the production and detection of potassium salt.   fate 2.16, 1.73; sodium chloride 1.37, 1.57. The results deduced that the prediction accuracy of U-CARS-PLSR has decreased to a certain extent except for potassium chloride. However, the CARS method reduced the number of selected calculated wavelengths (from 69 to 31 wavelengths), and the RPD values of major minerals in MP samples are both greater than 4.00, indicating that the prediction effect was still acceptable. In other words, CARS have a great role in reducing the cost of TIH technology and then being applied to the production and detection of potassium salt.

Conclusions
To solve the problem of rapid detection in the process of mineral production, this paper introduced the application of TIH technology in the detection of potassium salt samples as an example. The relationship between the potassium salt spectrum and mineral content was constructed, and the potassium salt content was successfully predicted under the PLSR method. The following conclusions are reached: (1) TIH imaging revealed highly prominent emissivity characteristics of the MP samples in the thermal infrared band. Furthermore, the temperature discrepancy between the inverted and actual temperatures of the samples was relatively small, with 71% of the samples exhibiting a temperature difference of less than 1 K. This suggests that the emissivity accuracy of the samples obtained in this experimental process is high and can be approximated as the emissivity spectrum corresponding to potassium salt samples. (2) The CARS-PLSR model, which is based on MP sample emissivity data training, is effective for MP sample prediction. In the U-PLSR model, the RPD values of the four minerals are 4.59, 4.06, 2.16, and 1.57, respectively, indicating that PLSR has a good prediction effect on MP. The calculation results of the U-CARS-PLSR model using the CARS method for sensitive wavelength selection show that the CARS method can effectively reduce the number of wavelengths, which is of great benefit to the practical application of TIH technology. With the CARS method, the number of selected wavelengths for the four minerals is reduced to 12, 11, 11, and 8, with RPD values of 4.24, Figure 12. RPD values of four minerals in U-PLSR and U-PCARS-PLSR models.

Conclusions
To solve the problem of rapid detection in the process of mineral production, this paper introduced the application of TIH technology in the detection of potassium salt samples as an example. The relationship between the potassium salt spectrum and mineral content was constructed, and the potassium salt content was successfully predicted under the PLSR method. The following conclusions are reached: (1) TIH imaging revealed highly prominent emissivity characteristics of the MP samples in the thermal infrared band. Furthermore, the temperature discrepancy between the inverted and actual temperatures of the samples was relatively small, with 71% of the samples exhibiting a temperature difference of less than 1 K. This suggests that the emissivity accuracy of the samples obtained in this experimental process is high and can be approximated as the emissivity spectrum corresponding to potassium salt samples. (2) The CARS-PLSR model, which is based on MP sample emissivity data training, is effective for MP sample prediction. In the U-PLSR model, the RPD values of the four minerals are 4.59, 4.06, 2.16, and 1.57, respectively, indicating that PLSR has a good prediction effect on MP. The calculation results of the U-CARS-PLSR model using the CARS method for sensitive wavelength selection show that the CARS method can effectively reduce the number of wavelengths, which is of great benefit to the practical application of TIH technology. With the CARS method, the number of selected wavelengths for the four minerals is reduced to 12, 11, 11, and 8, with RPD values of 4.24, 4.42, 1.73, and 1.37. The prediction accuracy of major minerals in MP is high (RPD > 4). The model has a good prediction effect on MP.
The study deduces that TIH technology has high accuracy for detecting potassium salts with limited mineral phases. It is reasonable to expect that TIH technology could have potential applications not only in detecting other minerals in production but also in other scenarios where mineral phases are limited, such as in the detection of allotropes.
TIH technology can accurately and rapidly detect mineral phases with limited mineral phases, compared with traditional detection methods. However, it is important to note that in this study, the TIH technology was used to collect data on dry potassium salt samples in a relatively open environment, which may be difficult to replicate under actual production conditions. Therefore, in subsequent studies, two approaches will be taken to address this issue: (1) Construct "shielding conditions" using low emissivity substances to simulate a more realistic environment; (2) Increase the sample size of potassium salt under aqueous conditions and quantitatively investigate the influence of water. These two approaches represent the focal points of follow-up research aimed at bringing the technology closer to practical application.

Data Availability Statement:
The data is unavailable due to privacy restrictions.