Performance Recognition for Sulphur Flotation Process Based on Froth Texture Unit Distribution

As an important indicator of flotation performance, froth texture is believed to be related to operational condition in sulphur flotation process. A novel fault detection method based on froth texture unit distribution (TUD) is proposed to recognize the fault condition of sulphur flotation in real time. The froth texture unit number is calculated based on texture spectrum, and the probability density function (PDF) of froth texture unit number is defined as texture unit distribution, which can describe the actual textual feature more accurately than the grey level dependence matrix approach. As the type of the froth TUD is unknown, a nonparametric kernel estimation method based on the fixed kernel basis is proposed, which can overcome the difficulty when comparing different TUDs under various conditions is impossible using the traditional varying kernel basis.Through transforming nonparametric description into dynamic kernel weight vectors, a principle component analysis (PCA) model is established to reduce the dimensionality of the vectors. Then a threshold criterion determined by the TQ statistic based on the PCA model is proposed to realize the performance recognition.The industrial application results show that the accurate performance recognition of froth flotation can be achieved by using the proposed method.


Introduction
Sulphur flotation is a complex physical process influenced by multiple operational variables such as inlet air flow, pulp level, and it is naturally hydrophobic to attach to the air bubbles.The objective of sulphur flotation is to separate valuable sulphur minerals from useless materials or other minerals so as to gain the upgraded sulphur minerals [1].Sulphur concentrate grade depends on flotation separation performance, and it is affected by the accuracy of the performance recognition.It is well recognized that froth visual appearance observed can characterize the combining effect of multiple process conditions on flotation [2], and it is also known as the indicator of flotation separation performance.Recent advances in image processing and computer vision based froth appearance monitoring systems contribute greatly to the feature extraction of visual descriptors [3][4][5].Computerbased vision technology is now moving out of the research laboratory and into the plant to become a useful means of monitoring and controlling flotation performance at the cell level [6][7][8].The development of base level process control (control of pulp level, air flow rate, etc.) has been significant progress, but automated advanced and optimization flotation control strategies based on computer vision have been more difficult to implement [9].The performance recognition is available for the optimal control of flotation [10], and flotation performance is closely related to the concentrate grade.Therefore, it is of great importance to improve the sulphur concentrate grade by developing an effective performance recognition method based on computer vision.
It is shown that the froth texture is a good indicator to the performance of the flotation cells [11], and texture information is believed to strongly associate with mineral grade [12].Numerous reported literatures are devoted to the extraction of froth image texture features.The grey-level cooccurrence matrix (GLCM) approach is one of the most popular statistical methods used in practice to measure the textural information of images.Most of the researchers calculated second-order statistics based on GLCM such as angular second moment, entropy, moment of inertia, and moments of deficit and relevance to recognize the flotation performance.GLCM was used as texture descriptor to classify different types of froths, and they also provide qualitative information on the changes in the visual appearance of the froth [13].However, there is only one angle of displacement used in the GLCM approach, with 0 ∘ , 45 ∘ , 90 ∘ , and 135 ∘ calculated to acquire average second-order statistics, neglecting the different variation in different directions and leading to large computing of the high-dimensional matrix.The actual flotation froth texture is more complex, so that a simple statistical property in the GLCM approach is difficult to accurately describe it, leading to the inaccuracy of performance recognition.Based on the proposed concept of texture unit, a new statistical approach to texture analysis, termed the texture spectrum approach, was proposed [14].The proposed method extracted the textural information of an image with a more complete respect of texture characteristics (simultaneously in all eight directions instead of only one displacement vector used in the GLCM approach).It is worth noticing that the PDF of texture unit number, defined as the texture unit distribution (TUD), is found to be nonnormal.Further exploring of the information indicated by froth structure has shown that the TUD is multipeaky and highly skewed, which neither belongs to any existing mathematical model based distribution.To depict the unknown continuous process of froth flotation, nonparametric estimation provides a credible solution.Commonly used nonparametric estimation techniques include histogram, frequency polygon, shift average histogram, kernel methods [15], wavelet method, and B-spline expansion models.Theoretical researches on tracking the output probability density distribution to a target distribution shape by using various control approaches [16] can also increase the possibility of froth visual features based process control.Nevertheless, the traditional kernel estimation cannot compare the various froth TUDs under different flotation conditions with the varying kernel basis.Therefore, the fixed kernel basis is proposed to describe the TUDs in various froth images.
To relate the flotation operation condition with flotation performance, Jampana et al. revealed that the increase in pulp level causes concentrate grade to decrease [17], as the variation of pulp level has great effect on the froth retention time in the flotation cell [18].The continuous decrease of froth retention time can lead to less collision time between mineral particles and bubbles with decreased gangue drop, which resulting in the deteriorating performance of mineral concentrate grade.Conventionally, industry process performance recognition heavily depends on the frequent inspection of froth views and manipulation of experienced human operators, which is often problematic with strong subjective and unable to regulate the fault performance timely, leading to the unstable flotation process and low concentrate grade.Along with the implementation of online monitoring system of froth visual appearance, quantitative performance recognition becomes highly desired and essential to maintain the operational variables at acceptable rates.Cilliers proposed a quantitative fault detection and diagnosis model which is successfully applied for hydrocyclones [19].In industrial case studies of aluminum flotation, Xu explored the froth structure by using kernel density estimation technique to approximate the output probability density of surface bubble size distribution rectified by the empirical formula and its application on process fault detection [20].
The froth texture characterizes the roughness of the froth surface, which indicates the mineral contents of froth.When the pulp level is too high with slurry overflow, froth texture is smooth; in this case, the middle value takes a large portion of the texture unit number in the whole image, which results in a high peak in the texture unit distribution curve.On the other hand, when the pulp level becomes too low, the froth cannot overflow, such that the mineral contents in the froth accumulate to a high level.Therefore, the texture becomes coarse and the middle value takes a small portion of the texture unit number in the whole image, resulting in a low peak in the texture unit distribution curve.By transforming the texture unit distribution into the weight vector using the fixed kernel estimator, a weight PCA model can be established to handle the variation in the texture unit distribution.The sulphur froth image contains a great deal of noise because of the acid fog in the sulphur flotation. 2 statistic based on the PCA model can reveal the major variation of the froth texture, and  statistic can reveal the noise contained in the image.Thus the proposed new statistical variable  is proposed to detect the sulphur flotation fault effectively by considering the influence of noise.
The main advantages of the proposed method in this paper are that (i) the texture unit distribution can describe the froth texture feature more completely, by considering eight directions of grey level variation information, compared to the GLCM method.(ii) The mathematical model of texture unit distribution is unknown, as it is nonnormal and multipeaky, so nonparametric estimation method is more suitable to approximate it.The fixed kernel basis can compare the different flotation performance reflected by the weight coefficients of texture unit distribution, compared with the traditional varying kernel basis.(iii) The new statistic  can reduce the influence of noise in the accuracy of performance recognition, compared to the traditional  2 and  statistics.
This work aims to explore the froth texture by using kernel density estimation technique to approximate the surface froth TUD and its application on sulphur flotation process performance recognition.A nonparametric kernel estimator by the fixed kernel basis is designed to approximate texture unit distribution, such that the output TUD is formulated in terms of dynamic weights, on which a principle component analysis (PCA) model is established.Then an effective performance recognition criterion is determined using the proposed  statistic based on PCA model.The fault condition is successfully detected on the industrial data of offline froth images.Next section introduces the froth texture unit number calculation oriented texture spectrum scheme.Section 3 presents the output TUD curve modeling by using designed kernel density estimators.The kernel weight vector based PCA model is established, and a threshold criterion determined by the  statistic based on the PCA model is proposed to realize the performance recognition in Section 4. Section 5 presents the experimental results and discussion.Conclusion is provided in the last section.

Surface Froth Texture Unit Number Calculation
Experimental setup consists of RGB camera with resolution of 1280×960 and lens of 35 mm, high frequency light source, cover hook protecting camera from dust, acid fog and ambient light, and optical fiber with length over 200 m for signal communication to industrial PC computer in operating room.The camera is mounted 96.5 cm vertically above the froth surface of the target cell, and froth images with window size 24.26 × 18.2 cm 2 are captured online at the rate of 15 frames/s.Meanwhile, the corresponding process operational and performance data are collected on industrial scale.
Froth images collected from industry field display that various froth texture feature leads to the different performance.The existing texture description method such as texture spectrum, spatial and neighboring grey-level cooccurrence matrix are derived from this fact.Froth image observed is a type of gradient images.Nevertheless, simple second-order statistical variables in the GLCM approach are difficult to accurately describe the froth texture, the texture unit (TU) oriented texture spectrum scheme proposed by [14] is used to describe the texture features.In a froth digital image, each pixel is surrounded by eight neighboring pixels.The local texture information for a pixel can be extracted from a neighborhood of 3 × 3 pixels called texture unit, which represents the smallest complete unit (in the sense of having eight directions surrounding the pixel).Given a neighborhood of 3 × 3 pixels (which will be denoted by a set containing nine elements:  = { 0 ,  1 , . . .,  8 }, where  0 represents the intensity value of the central pixel, and   { = 1, 2, . . ., 8} is the intensity value of the neighboring pixel ) defines the corresponding texture unit by a set containing eight elements, TU = { 1 ,  2 , . . .,  8 }, where   { = 1, 2, . . ., 8} is determined by the following formula: for  = 1, 2, . . ., 8, and the element   occupies the same position as the pixel .
As each element of TU has one of three possible values, the combination of all eight elements results in 3 8 = 6561 possible texture units in total.There is no unique way to label and order the 6561 texture units.In our study, the 6561 texture units are labeled by using the following formula: where  TU represents the texture unit number, and   is the ith element of texture unit set TU = { 1 ,  2 , . . .,  8 }.
In addition, the eight elements may be ordered differently.If the eight elements are ordered clockwise as shown in Figure 1, the first element may take eight possible positions from the top left (a) to the middle left (h), and then the 6561 texture units can be labeled by the above formula under eight different ordering ways (from a to h).TUD is defined as the occurrence frequency for every texture unit number, and it exhibits probability density function (PDF) distribution of froth texture unit number.The online acquired sulphur froth image in cleaner cell in normal condition is shown in Figure 3. Figure 4 shows the froth TUD.The froth texture unit probability density distribution is found to be nonnormal and multipeaky.

TUD Curve Modeling
The surface sulphur froth TUD is nonnormal.Unlike traditional method applying singular feature such as mean or variance with the assumption that the distribution is normal, probability density distribution is suggested to accurately describe statistical feature of froth texture.The fact that the mathematical model of TUD is unknown makes nonparametric estimation method fitting to depict the unknown continuous process of froth flotation.

Nonparametric Kernel Estimation.
Consider a probability density function () describing the probability distribution of  in (, ) as follows: Density estimation accomplishes the fitting of ().Though classic nonparametric histogram estimator is good for data presentation, its discontinuity causes difficulty if derivatives of the estimates are required.A continuous version of the histogram is the frequency polygon formed by interpolating the midpoints of a histogram.Histogram based methods seek the balance between estimation accuracy and feature dimensionality, which can be very expensive for large samples.Apart from the histogram, the kernel estimator is most commonly used [15], which is given by where (( −   )/ℎ) is the ith kernel function, and   is the center of the ith kernel function along the horizontal axis.

Output TUD Kernel Estimation.
Supposing there is a dynamic stochastic system with input () ∈   and output () ∈ [, ], the probability of output () lying in [, ) is defined as where the  ker (, ) represents the output TUD after froth texture unit number calculation.() is control input such as the input amount of pulp level which is a dominant operational condition in the sulphur flotation system.The  ker (, ) can be approximated by kernel estimators designed in formula ( 5) and the corresponding weights   .
However, the traditional kernel estimation cannot compare the various froth TUDs under different flotation conditions with the varying kernel basis.Therefore, the fixed kernel basis is proposed to describe the TUDs in various froth images, such that the TUD curves can be transformed into dynamic kernel weight vectors, based on which the fault condition can be detected in sulphur flotation.Meanwhile, the computational complexity is also reduced using the designed fixed kernel basis.
Adjusting to the range of froth texture unit number, a number of kernel bases are selected to depict the TUD  in Figure 5.Its window width ℎ is fixed across the entire sample.As the TUD for sulphur froth image is multipeaky and complicated, 25 kernel bases with fixed window width are used to approximate TUD, which is plotted in dashed line in Figure 5.One dotted curve presents the first and second kernel basis multiplying the corresponding weight coefficients.And the estimation result of Figure 3 froth TUD is plotted in solid line.Figure 6 presents the kernel density estimation methods to approximate actual texture unit distribution of sulphur froth image in Figure 3.The results have shown the kernel estimation can accomplish the description of froth texture unit probability density distribution with general low feature dimensionality and high accuracy.

Weight PCA Model Based Performance Recognition
A fault performance is defined as the departure from an acceptable range of an observed output or operating variable.Timeous detection of fault can determine whether the abnormal condition occurs [21].The information indicated by froth characteristics is a combining effect of multiple operational variables such as pulp level and inlet air flow in sulphur flotation.By retaining the variance of inlet air flow during a short period of time, froth texture in cleaner cell is closely related to concentrate grade which is determined by the regulation of pulp level.Human operators are in capable of performing timeous monitoring of various process variables, and the process manipulation mostly relies on heuristics of their froth view observation.Hodouin used PCA to analyze and interpret the behavior of mineral flotation and grinding circuits in a large mineral processing plant [22].Kourti summarized the latest developments in multivariate statistical process control (MSPC) and its application for fault detection and isolation (FDI) in industrial processes [23].

Weight PCA Model.
PCA is a multivariate statistical technique used in MSPC and FDI perspectives [23].PCA uses latent variables instead of every measured variable in the process because they can better explain the behavior of the process.By monitoring the sulphur froth appearance such as froth texture, the process fault performance can be inferred and identified based on established PCA model.The output TUD for sulphur flotation froth can be transformed to dynamic kernel weight vectors () through formula (8).PCA reduces the dimensionality of the original weights () by projecting it onto a lower dimensionality space.It obtains the principal causes of variability in the sulphur flotation process.If some of these causes change, it can be due to a fault in the process.
Consider the weight matrix W ∈  ×(−1) , containing  samples of  − 1 dynamic kernel weight coefficients collected under normal operation in sulphur flotation.This matrix must be normalized to zero mean and unit variance with the scale parameter vectors  and  as the mean and variance vectors, respectively.Next step to calculate PCA is to construct the covariance matrix S: and performing the SVD decomposition on S: where Λ is a diagonal matrix that contains in its diagonal eigenvalues of S sorted in decreasing order ( 1 ≥  2 ≥ ⋅ ⋅ ⋅ ≥  −1 ≥ 0).Columns of matrix Λ are the eigenvectors of S.
The transformation matrix P ∈  × is generated choosing  eigenvectors or columns of V corresponding to  principal eigenvalues.Matrix P transforms the space of the measured variables into the reduced dimension space as follows: Columns of matrix P are called loadings, and elements of T are called scores.Scores are the values of the original measured variables that have been transformed into the reduced dimension space.
Operating in (11), the scores can be transformed into the original space as follows: The residual matrix E is calculated as Finally the original data space can be calculated as It is very important to choose the number of principal components , because TP  represents the principal sources of variability in the process, and E represents the variability corresponding to process noise.There is Cumulative Percent Variance (CPV) approach for determining the number of components to be retained in a PCA model as [24].The measure of the percent variance (CPV() ≥ 95%) captured by the first  principal components is adopted as follows:

A New Statistical Variable Based Fault Performance
Recognition.Having established a PCA model based on historical data collected when only common cause variation are present, multivariate control charts based on Hotelling's  2 and square prediction error (SPE) or  can be plotted.The fault performance recognition can be reduced to this two traditional variables ( 2 and ) characterizing two orthogonal subsets of the original space.However, some of sulphur froth images contain a great deal of noise because of the acid fog in the sulphur flotation.The traditional  2 statistic can only describe the variation in the texture information, therefore, the normal image of containing noise caused by acid fog may be detected as a fault image for its disability to handle the noise influence.As the  statistic can represent the random noise in the froth texture, by combining  2 statistic and  statistic, the new statistic  is proposed to detect sulphur flotation fault performance more accurately: where  is the regulation factor controlling the value range of . takes values between 99% and 100%. 2 can be calculated as the sum of squares of a new process weight vector  as follows: where Λ  is a squared matrix formed by the first  rows and columns of Λ.
The sulphur flotation process is considered normal for a given significance level  if where   (,  − ) is the critical value of the Fisher-Snedecor distribution with  and  −  degrees of freedom and  the level of significance. takes values between 90% and 95%. 2 is based on the first  principal components, so that it provides a test for derivations in the latent variables that are of greatest importance to the variance of the sulphur flotation process.This statistic will only detect an event if the variation in the latent variables is greater than the variation explained by common causes.
New events can be detected by calculating the SPE or  of the residuals of a new observation. statistic is calculated as the sum of squares of the residuals.The scalar value  is a measurement of goodness of fit of the sample to the model and is directly associated with the noise as follows: with e = (I − PP  ) .
The upper limit of this statistic can be computed as the next form: with where   is the value of the normal distribution with , and the level of significance and   are the eigenvalues of the PCA residual covariance matrix E  E/( − 1).When an unusual event occurs and it produces a change in the covariance structure of the model, it will be detected by a high value of .
According to the formulae ( 18) and ( 21), the critical value of the new statistical variable can be calculation as Through using the output TUD weight based PCA model, a criterion can be designed to detect the fault.The new statistical variable  is calculated on the weight PCA model.Then the critical value   of statistical variable is set as the threshold value.When the value of  for the new sample is larger than the threshold value   evaluated by formula (24), the fault can be detected.8) to approximate the bubble texture unit distributions of froth images in Figure 7, the 3D mesh plot of the output TUD when fault occurs is presented.

Application Results and Discussion
To evaluate the proposed weight PCA model based fault detection approach, a series of industrial experiments are carried out in a Chinese sulphur froth flotation plant.In the test runs, froth image videos are captured through the previously introduced monitoring system in the last cleaner flotation cell.Subsequently, the froth videos are processed by the developed image analysis software which is capable of extracting froth features such as TUD online.Figure 7 presents the three types of froth images in different performance, which are collected and analyzed under the same condition in terms of resolution, angle, light condition, position, view scale, and so forth.
In practical sulphur flotation process experiments, the air flow rate and feed-in conditions are kept at a steady state so as to stabilize the production process.The adjustment of pulp level (or froth depth) becomes the major manipulating parameter, which directly determines flotation performance.As an indication of flotation performance, the froth texture feature is one of determinants of mineral separation efficiency.Bubbles with relative complex texture generally carry more valuable mineral particles, whose corresponding pulp level value is to be maintained to an acceptable bounded range.When one of the dominant operating variable pulp level is fluctuated, in this case the regulation of slurry underflow, froth surface visual features such as froth texture and color spectral information are reacting to the change of pulp level value.An increase in pulp level was considered, such that its simultaneous effect on froth texture unit distribution can be identified.As is shown in Figure 7, the froth images evolved as pulp level value varied gradually during a period, and the corresponding operational conditions were measured at the same time.
As for texture unit number calculation, normal kernel with following basis functions is selected according to formula (5).The window width ℎ is set to be 200 as a smoothing parameter, and centered points of each kernel   = 250 * , ( = 1, . . ., 25).Since the froth texture unit involved ranges from 0 to 6560, the kernel functions with fixed window width are supposed to cover the entire texture unit value range.Thus, the froth TUD can be approximated by (8), where  = 25.The weights () of normal kernel expansion have dimension of 25, and only 24 of which are independent.
By applying the kernel estimation on the TUDs of froth images in Figure 7, the 3D mesh plot of the output TUD is shown in Figure 8.At an hourly interval, the froth image video is captured at the point since it is reasonable to consider that the froth TUD is representative during a short-time period in this study case.Meanwhile, the process operational conditions are measured correspondingly.As can be seen, the froth TUD tends to shift dramatically with occurring a low peak when the slurry underflow increased at 9:00,  which resulted in froth depth value increased from 190 mm to 350 mm in response.Then the excessive decrease of froth depth to 30 mm produces a corresponding upward change of the peak of TUD curve.Accordingly, the separation performance mineral grade deteriorated from 81% at 8:00 to 50% at 11:00.
The weight PCA model applied in this case is established as where  = 9.
According to the formula (24), threshold value can be calculated as   = 18.11 by setting  = 0.95,  = 0.99.Setting Figure 7(a) as the normal TUD, Figure 9 shows the threshold with solid line and the  statistics for froth images in Figure 7.As can be seen, Figures 7(b)-7(c) are clearly identified as fault status because the statistic  >   , which are consistent with the observation results from human operators.
Attempts have been made in calculating false alarm rate on a testing database.The testing data consist of 243 offline froth videos captured from the sulphur flotation industry during August of 2011.The fault detection is accomplished by a threshold criterion calculated from formula (24), according to which the  statistic above the threshold value indicates that a fault occurs.Table 1 gives the detection performance of the testing database.As can be seen in Figure 9, solid line represents the threshold for fault detection, and asterisks are  statistics of each normal video sample.Triangles and diamonds are the samples with fault A and fault B. The total fault detection accuracy on the database is 93.83%.It is possible that the false detection alarm ascribes to the texture unit number calculation malfunction of the captured froth images.

Conclusion
In this paper the description of texture unit number probability density distribution and its relationship to pulp level operational status are investigated.Unlike traditional discussion of froth texture feature focusing mostly on second-order statistics based on GLCM including angular second moment, entropy, moment of inertia, and moments of deficit and relevance, a nonparametric estimation method is proposed to describe the TUD more accurately based on the fixed normal kernel basis, and the fault performance is detected through the proposed  statistic.Desired fault detection for pulp level regulation in froth flotation industry is achieved using the proposed method.

Figure 1 :
Figure 1: Eight clockwise, successive ordering ways of the eight elements of the texture unit.The first element  1 may take eight possible positions from a to h.

Figure 2
Figure 2  gives an example of transforming a neighborhood in sulphur flotation froth image to a texture unit with the texture unit number under the ordering way a.TUD is defined as the occurrence frequency for every texture unit number, and it exhibits probability density function (PDF) distribution of froth texture unit number.The online acquired sulphur froth image in cleaner cell in normal condition is shown in Figure3.Figure4shows the froth TUD.The froth texture unit probability density distribution is found to be nonnormal and multipeaky.

Figure 2 :
Figure 2: Example of transforming a neighborhood to a texture unit with the texture unit number.

Figure 3 :
Figure 3: The online acquired sulphur froth images in cleaner cell.

Figure 5 :
Figure 5: Normal kernel estimation and the weight coefficients.

Figure 6 :
Figure 6: Normal kernel based method to approximate the actual texture unit density estimation for Figure 3.

Figure 7 :
Figure 7: Industrial froth images captured in a continuous time: (a) normal froth image, (b) fault A froth image, and (c) fault B froth image.

Figure 8 :
Figure8: The 3D mesh plot of the output TUD.By applying the kernel estimation proposed in formula(8) to approximate the bubble texture unit distributions of froth images in Figure7, the 3D mesh plot of the output TUD when fault occurs is presented.

Figure 9 :
Figure 9: The  statistics of the testing froth video database.

Table 1 :
The detection performance of the testing database.