Soil Salinization Level Monitoring and Classifying by Mixed Chaotic Systems

: Soil salinization process is a complex non-linear dynamic evolution. To classify a system with this type of non-linear characteristic, this study proposed a mixed master/slave chaotic system based on Chua’s circuit and a fractional-order Chen-Lee chaotic system to classify soil salinization level. The subject is the soil in Xinjiang with different levels of human interference. A fractional-order Chen-Lee chaotic system was constructed, and the spectral signal processed by the Chua’s non-linear circuit was substituted into the master/slave chaotic system. The chaotic dynamic errors with different fractional orders were calculated. The comparative analysis showed that 0.1-order has the largest chaotic dynamic error change, which produced two distinct and divergent results. Thus, this study converted the chaotic dynamic errors of fractional 0.1-order into chaotic attractors to build an extension matter-element model. Finally, we compared the soil salt contents (SSC) from the laboratory chemical analysis with the results of the extension theory classiﬁcation. The comparison showed that the combination of fractional order mixed master/slave chaotic system and extension theory has high classiﬁcation accuracy for soil salinization level. The results of this system match the result of the chemical analysis. The classiﬁcation accuracy of the calibration set data was 100%, and the classiﬁcation accuracy of the validation set data was 90%. This method is the ﬁrst use of the mixed master/slave chaotic system in this ﬁeld and can satisfy certain soil salinization monitoring needs as well as promote the application of the chaotic system in soil salinization monitoring.


Introduction
Soil is made up of a complex mix of minerals, organic matter, living organisms, moisture, and air [1,2]. One of the global soil problems is soil salinization, especially in arid and semi-arid regions where evaporation is strong, rainfall is scarce, and soil salinization is particularly serious. Soil salinization leads to serious consequences such as soil fertility decline, soil compaction, and reduced crop productivity. Therefore, a timely, fast, and accurate grasp of the soil salinization status is essential to regional ecological stability [3][4][5][6][7][8]. At the same time, the salt content of saline soil is related to multiple components of the soil sample, and there is a complicated nonlinear relationship between the salt content and the hyperspectral reflectance curve. Therefore, it is difficult to reveal the nonlinear dynamic evolution process of soil salinization.
Xinjiang is located in an arid/semi-arid area rich in land and resources that is currently an important agricultural backup resource area for China [9][10][11]. However, the soil salinization in Xinjiang mainly refers to the salt in underground water rising to the soil surface through capillaries. After the water evaporates, the salt accumulates on the soil surface. Saline soil is also called saline-alkaline soil and mainly occurs in soil with high evaporation levels in arid/semi-arid areas and areas with high dissolved salt content in underground water [12,13]. Salinization not only damages resources and causes huge losses in agricultural production [14][15][16] but also threatens the environment and ecology. As a result, salinization has become a global problem [17,18]. Changes in soil salt content (SSC) are affected by the natural environment and human activities. A better soil salinization assessment requires spatial complexity and law of change that sufficiently describe salt content.
Soil salinization development is a complex non-linear dynamic evolutionary process that is closely related to factors such as geology, geomorphology, hydrology, climate, and human activities [19,20]. These dynamic change features and laws have always been popular research topics. Some studies have concentrated on predicting the salt content or nutritional content of saline soils. Farifteh et al. [21] prepared six types of soil with different MgCl 2 , NaCl, KCl, K 2 SO 4 , MgSO 4 , and Na 2 SO 4 concentrations. The salts are irrigated into three types of soil samples, silty clay loam, sandy loam, and sand texture soil. The relationship between the soil spectrum features and salt content was then analyzed. The result shows that as salt concentration decreases, the number and clarity of the diagnostic absorption waveband also decrease. The spectrum has the highest prediction accuracy towards MgCl 2 and MgSO 4 but has the lowest prediction accuracy towards KCl and Na 2 SO 4 . The absorption feature of wavebands above 1300 nm became wider. Srivastava et al. [22] used visible-near infrared (VNIR) spectrum technology to evaluate the soil salt content of saline soil in the Indo-Gangetic plains (IGP). It found that the area's soil salt content is sensitive to wavebands in the 1390-2400 nm range. The model proposed in this study can explain 80% or more of the electrical conductivity (EC), Na + , Cl − , Mg 2+ , sodium adsorption ratio (SAR), and Ca 2+ changes in the validation set. Wang et al. [23] collected saline soil spectrums in the Yellow River Delta region of China and found that the correlation coefficient between salinity spectral index (SSI) and SSC reached 0.91. Estimation with SSI indicated that the SSC single variable regression model's coefficient of determination was 0.873, and the root-mean-square deviation was 0.986. Mouazen et al. [24] studied soil samples in Belgium and Northern France and used three types of model building methods to predict phosphorus, sodium, magnesium, and potassium content. The results showed that by using the partial least squares regression's (PLSR) latent variables as a back-propagation neural network (BPNN), the predictive capability of the BPNN-LV model was the best. The RPD (ratio of the performance to deviation) value for organic carbon (OC(g/kg)) and magnesium was 2.54, and the RPD values for phosphorus, sodium, and potassium were between 1.77 to 1.94. Wang et al. [25] used Xinjiang's saline soil as their subject and analyzed the relationship between main salt ions and spectrums of fresh, air dried, and dried soil samples. The results show that the smaller the diameter of soil particles, the better the spectrum predictive result towards soil ion content. The inversed ion content of the fitted model built for dried soil was the most precise, and the RPD for K + , Na + , Ca 2+ , and SO 4 2− reached above 2. Small changes in soil parameters or the environment can cause differences in spectral reflectance, and such subtle differences are difficult to distinguish. At present, few studies have been conducted on the classification of soil salinization degree, and there is currently no report of chaos theory being used for salinization study. However, chaos theory has pseudo-randomness, ergodicity, regularity, and sensitivity towards initial conditions; thus, it is often applied to non-linear systems [26][27][28][29][30]. This paper is the first to use mixed chaos theory to describe the spectrum dynamic change law of soil with different salt content. The objective is to explore the feasibility of using chaos in the spectrum field. In this study, we analyzed saline soil from areas with different levels of human interference in Fukang Remote Sens. 2021, 13, 3819 3 of 21 City (Xinjiang) with a mixed fractional-order master/slave chaotic system (combination of Chua's chaotic circuit and fractional order Chen-Lee master/slave chaotic system). We used chaotic dynamic error, chaotic attractor, and extension theory to classify spectrum data from the soil with different salinization levels, which allowed us to determine the study area's soil salinization level. This research is divided into five sections, the first content is the Introduction, the second is Materials and Methods, the third is Simulation results, the fourth is Discussion, and the fifth is Conclusions. The aims of this research are: (1) Explain the feasibility of the application for the Chen-lee chaotic system in the field of soil hyperspectral sensing. (2) Propose a new intelligent classification system architecture to classify the degree of soil salinization. (3) Provide a scientific basis for soil salinization management in the study area, as well as provide a time and labor-saving soil salinization data extraction method.

Soil Sample Collection
There is a severe soil salinization problem in the agricultural development of Xinjiang. The total saline-alkaline soil area is approximately 8.476 million hm 2 , including cultivated land. Overall, 31.1% of the area is jeopardized by saline-alkaline conditions, making it the largest saline-alkaline soil area in China. Salinization and development have long been important factors that limit sustainable development of agriculture, resources, environment, and social economics in Xinjiang. The soil samples used in this study were collected in the desert area of Fukang in May 2017. The Fukang desert area is located in Fukang City, Xinjiang Provence, its geographic coordinates are 87 • 44 -88 • 46 E, 43 • 29 -45 • 45 N, which is a typical temperate continental desert climate [31], where the winter is cold, summer is hot, sunlight is sufficient, annual average temperature is 6.6 • C, annual and daily temperature differences are large, precipitation is sparse, annual average precipitation is 145 mm, evaporation is strong, and the average annual evaporation potential is 2292 mm. The soil types include gray desert soil, cracked soil, and sandy soil. The soil texture in this area is sandy soil, and the soil salinization is serious, which is represented in the arid area of northwestern China.
In order to arrange the sampling points reasonably, it is necessary to conduct a field investigation to determine the scope of the research area; in this paper, the research is divided into Areas A and B. Area A is an area without human activities, which is far away from human settlements, and the soil basically maintains the original ecological features. Area B is an area affected by human activities, which is located near Xinjiang's 102 Production and Construction Army Corp and has had heavier human activity in recent years. Most of the land has been developed into elm forests and growing areas.
Five sampling lines were laid out in Area A, with the spacing of the sampling lines is about 500 to 700 m, and five representative sampling points were set on each sampling line. Six sampling lines were arranged in Area B. The interval between the sampling lines was about 400 to 600 m, and five representative sampling points were set on each sampling line. The soil is mainly collected from the top 0-10 cm of soil, 25 sampling points are taken from Area A, and 30 sampling points are taken from Area B. Approximately 1 kg of soil was taken from each sampling point. In this way, a grid that can cover different areas is formed in the entire study area. The sampling points can not only better reflect the soil properties of each area but can also better control the uniformity of the sampling points in each area.
The samples were taken back to the laboratory and air-dried naturally. Pebbles, animal, and vegetation remains were removed. The remaining soil was grounded, put through a 1 mm sieve, and sent to Xinjiang Institute of Ecology and Geography. Chemical analysis was used to determine the SSC; the salt content of the soil sample was measured in the soil clear liquid using the Multi 3420 SET B portable multi-parameter analyzer, and the salt content was converted from the salinity data measured by this instrument. The location of the study area and the soil sampling points are shown in Figure 1. the salt content was converted from the salinity data measured by this instrument. The location of the study area and the soil sampling points are shown in Figure 1.

Soil Hyperspectral Measurement
The field hyperspectral of the soil samples was measured by FieldSpec®3Hi-Res spectrometer, manufactured by an ASD (Analytical spectral device) company, which is affiliated to PANalytical, the Netherlands. The measurement range of all bands was 350-2500 nm, and the number of output bands was 2151. The wavelength range of 350-1000 nm is in the visible light with a spectral resolution of 3 nm, and the sampling interval is 1.4 nm. The range of 1001-2500 nm is near-infrared with a spectral resolution of 8 nm and a sampling interval of 1.1 nm. The re-sampling interval was 1 nm. The spectrometer is widely used in remote sensing, geology, agriculture, environment, surveying, mapping, and other fields. The measurement objects can be soil, ice and snow, cultural relics, minerals, plants, rocks, metals, and so on. The relevant technical parameters are shown in Table 1.
The requirement for field spectrum measurement is a clear day with no cloud and wind. The probe is placed 15 cm vertically above the soil sample to collect the spectrum. Before each measurement, a white board is used for calibration to remove dark current effects. To make the spectrum data representative, five locations near each sampling point were used to collect the spectrums. Each location was measured 10 times to obtain 50 spectrum reflectance curves. The curves are averaged to obtain the final sample spectrum reflectance data.

Soil Hyperspectral Measurement
The field hyperspectral of the soil samples was measured by FieldSpec ® 3Hi-Res spectrometer, manufactured by an ASD (Analytical spectral device) company, which is affiliated to PANalytical, the Netherlands. The measurement range of all bands was 350-2500 nm, and the number of output bands was 2151. The wavelength range of 350-1000 nm is in the visible light with a spectral resolution of 3 nm, and the sampling interval is 1.4 nm. The range of 1001-2500 nm is near-infrared with a spectral resolution of 8 nm and a sampling interval of 1.1 nm. The re-sampling interval was 1 nm. The spectrometer is widely used in remote sensing, geology, agriculture, environment, surveying, mapping, and other fields. The measurement objects can be soil, ice and snow, cultural relics, minerals, plants, rocks, metals, and so on. The relevant technical parameters are shown in Table 1. The requirement for field spectrum measurement is a clear day with no cloud and wind. The probe is placed 15 cm vertically above the soil sample to collect the spectrum. Before each measurement, a white board is used for calibration to remove dark current effects. To Remote Sens. 2021, 13, 3819 5 of 21 make the spectrum data representative, five locations near each sampling point were used to collect the spectrums. Each location was measured 10 times to obtain 50 spectrum reflectance curves. The curves are averaged to obtain the final sample spectrum reflectance data.

Spectrum Data Pre-Processing
Because the spectrometer itself has a margin of error and there is noise in the measuring environment, the spectral reflectance can experience baseline drift and multiple scattering effects. These factors can affect the reliability of the final model. Thus, before building the model, it is necessary to pre-process the spectrum data. First, the Savitzky-Golay smoothing is used to treat the spectrum and remove errors that may be caused by spectrum curve noise. However, the 350-399 nm and 2401-2500 nm waveband has a lower spectrum. Thus, the spectral reflectance data we chose to use in this study is the 400-2400 nm wave band with higher signal noise and does not include the wavebands on the two ends. Because of the effects from the moisture absorption belt, the 1400 nm, 1800 nm, and 2200 nm spectrum reflectance curves show clear fluctuation. Therefore, the wavebands near the moisture absorption belt, 1355-1410 nm, and 1820-1942 nm were also removed. Figure 2 shows the spectral reflectance curves collected in Area A and B sampling points. Near the 1430, 1950, and 2180 nm bands, there is a deeper absorption valley, and at the 700, 900, 1100, 1750, and 2300 nm, there is a shallower absorption valley. As the sampling points change, the 600-1700 nm and 2200-2400 nm spectrum reflectance curves have a higher soil spectrum curve overlap, which makes it difficult to directly identify the objective. Thus, the soil salinization level is divided into five major categories based on the soil salinization level determination standards [32], as shown in Table 2. This study screened the spectral reflectance of sampling points with different salinization levels in Area A and B, as shown in Figure 3. Figure 3 shows that the change in spectral reflectance curve for different salinization levels is basically consistent. However, differences in soil salt content cause changes in spectrum reflectance. As the soil salt content increases, the spectrum reflectance also increases. This phenomenon is consistent with the result of previous studies [33]. Table 2. Salinization classification standards.

Degree of Soil Salinization Soil Salt Content (g/kg)
Non

Chua's Circuit
Chua's circuit is a type of simple non-linear electronic circuit design [34,35], as shown in Figure 4. It is formed by four linear components, the inductor L, resistor R, capacitor

Chua's Circuit
Chua's circuit is a type of simple non-linear electronic circuit design [34,35], as shown in Figure 4. It is formed by four linear components, the inductor L, resistor R, capacitor C 1 , C 2 , and non-linear Chua's diode N r .

Chua's Circuit
Chua's circuit is a type of simple non-linear electronic circuit design [34,35], as shown in Figure 4. It is formed by four linear components, the inductor L, resistor R, capacitor  According to Kirchhoff's Law, Chua's circuit status equation is shown as Equation (1).
where, V 1 , V 2 , and I 3 are the voltage on the two ends of the capacitor and the current that passes the inductor, respectively; f (V 1 ) is the non-linear resistor; N r is the non-linear Chua's diode. The characteristic curve of the non-linear Chua's diode N r is as shown in Figure 5. The polynomial reported in Figure 2 is, as shown in Equation (2): where, G a and G b are the slopes of the inside and outside line turn segment, E is the turning point voltage, I is current, and V is voltage. Where, 1 V , 2 V , and 3 I are the voltage on the two ends of the capacitor and the current that passes the inductor, respectively; 1 () fV is the non-linear resistor; N r is the nonlinear Chua's diode.
The characteristic curve of the non-linear Chua's diode N r is as shown in Figure 5. Where, G a and G b are the slopes of the inside and outside line turn segment, E is the turning point voltage, I is current, and V is voltage.
The polynomial reported in Figure 2 is, as shown in Equation (2):

Chen-Lee Chaotic System
Chaos theory is a type of non-linear system theory. Every chaotic system has chaotic attractors, and the signal outputted by the attractors displays an orderly but non-periodic motion trajectory. Every tiny change inputted can significantly affect the outputted motion track and chaotic attractor pattern. This is the butterfly effect. The chaotic self-synchronization system is composed of the master system (MS) and slave system (SS) and can be separately expressed as Equations (3) and (4).

Chen-Lee Chaotic System
Chaos theory is a type of non-linear system theory. Every chaotic system has chaotic attractors, and the signal outputted by the attractors displays an orderly but non-periodic motion trajectory. Every tiny change inputted can significantly affect the outputted motion track and chaotic attractor pattern. This is the butterfly effect. The chaotic selfsynchronization system is composed of the master system (MS) and slave system (SS) and can be separately expressed as Equations (3) and (4).
where, X ∈ R N and Y ∈ R N are the state vectors, X represents the master system states, Y represents slave system states. A is an N × M matrix, f (X) and f (Y) are non-linear vectors, and U is the designed non-linear control component. By making the master/slave mutually pursue each other and produce dynamic errors, we can determine system feature changes. When the master/slave receives a signal, subtraction of the master/slave can obtain the dynamic error status. This study uses a fractional-order master/slave Chen-Lee chaotic system [36][37][38][39], and its dynamic equation can be expressed as Equation (5): where, α, β, and γ are the system parameters. When they are 5, −10, and −3.8, the system is guaranteed to have a chaotic attractor [40]. x, y, and z represent state variables of the Chen-Lee chaotic system. The Chen-Lee chaotic system track is as shown in Figure 6. According to the master/slave chaotic dynamic error produced by the Chen-Lee chaotic system, if T denotes the transpose of a matrix, then X is the MS, and Y is the SS. Equations (3) and (4) can be revised as Equations (6) and (7).
Remote Sens. 2021, 13, 3819 Where,  ,  , and  are the system parameters. When they are 5, -10, and -3.8, the system is guaranteed to have a chaotic attractor [40]. X, Y, and Z represent state variables of the Chen-Lee chaotic system. The Chen-Lee chaotic system track is as shown in Figure 6. According to the master/slave chaotic dynamic error produced by the Chen-Lee cha- u , T denotes the transpose of a matrix, then X is the MS, and Y is the SS. Equations (3) and (4) can be revised as Equations (6) and (7). To make the master/slave pursue each other and produce dynamic errors, which we can use to determine system feature changes, we defined the errors e 1 = x 1 − y 1 , e 2 = x 2 − y 2 , e 3 = x 3 − y 3 . Thus, the dynamic error can be expressed as Equation (8).
In order to guarantee that the error dynamics also demonstrate a chaotic behavior, we set the control component U = [u 1 , u 2 , u 3 ] T as Equation (9), and the controlled error dynamic equation is shown in Equation (10). .
It shows that the dynamic error Equation (10) has the same structure as Equation (5). Therefore, we can guarantee the dynamic error Equation (10) will also show a stranger attractor with the same system parameter. We can reform Equation (10) to a matrix form as in Equation (11).
To make the dynamic error produced by the Chen-Lee chaotic system be clearer and conform more to non-linearity, this study used the fractional order differential theory proposed by Grünwald-Letnikov [41,42]. Its relative dynamic error e can be expressed as Equation (13).
To be able to clearly express spectrum signal features, Equation (13) is re-written as Equation (14).
To realize the dynamic errors with a digital system, its dynamic errors can be set as . In view of this, the dynamic self-synchronization error system mapping can be revised as in Equation (19), below: In this study, we will use the Self-Synchronization Error System Mapping to extract the feature differences in spectral reflectance of soil salinization level.

Extension Theory
Extension theory studies the possibility of extending events and pioneer innovative laws and methods; it is also used to solve contradictions [43,44]. Extension theory can be divided into two parts, the first part is the matter-element theory, and the second part is the extension set. Matter-element theory mainly studies the extendibility of matter-element and the nature of matter-element transformations. The extension set quantifies during the contradiction solving process. Finally, we build a corresponding quantifying tool and plan according to the correlation level. The fuzzy set range is extended from <0, 1> to <−∞, ∞>. The extension set is as shown in Figure 7.
Extension theory uses the matter-element model to describe an event. The extension theory builds a matter-element model for specific items/events. The matter-element model can be expressed as Equation (20).
where R represents an event, N is the name of the event, c is the characteristic of the event, and v is the event measurement.
The equation of extension theory can be seen in a previous study [45]. Extension theory uses the matter-element model to describe an event. The extension theory builds a matter-element model for specific items/events. The matter-element model can be expressed as Equation (20).
Where R represents an event, N is the name of the event, C is the characteristic of the event, and v is the event measurement.
The equation of extension theory can be seen in a previous study [45].

Intelligent Classification System Architecture
The spectrum of saline soil is an extremely complex evolution process, which is affected by many inter-influencing factors, and the process exhibits strong nonlinear characteristics. Small changes in soil parameters or the environment will cause differences in spectral reflectance, making such subtle spectral differences difficult to identify. The Chen-Lee chaotic system is used as a signal analysis tool. The intelligent classification system architecture can be shown in Figure 8. The specific eight steps can be described as follows: One, we use an ASD spectrometer to collect the spectral signal of the soil in the study area and observe the hyperspectral signal of the non-saline soil, which is used as a part of the input signal of the X variable.
Two, the hyperspectral signal of the saline soil is decomposed by Chua's circuit, and the Chua's circuit status equation is shown as Equation (1). The spectral sequence processed by Chua's circuit is used as a part of the input signal of the Y variable.
Three, it generates the main system signal and the slave system signal with the selfsynchronization of the chaotic system. The hyperspectral signal is substituted into the master-slave chaotic system, and the slave system can be synchronized with the master system under nonlinear control; the chaotic self-synchronization system can be separately expressed as Equation (3) and (4).
Four, the hyperspectral signal and salinity information of the saline soil are substituted into the constructed fractional-order Chen-Lee chaotic system, further subtraction is performed, and the dynamic equation of fractional order master/slave Chen-Lee chaotic system can be expressed as Equation (5).
Five, we further obtained the dynamic error of the master-servant chaotic system, and the dynamic self-synchronization error system mapping can be revised as in Equation (19). Spectrum signals from the soil with different salinization levels are introduced into different fractional-order mixed master/slave chaotic systems to produce chaotic dynamic error distribution.
Six, the hyperspectral signal and salinity information of the saline soil are substituted into the constructed fractional-order chaotic system, its dynamic behavior characteristics

Intelligent Classification System Architecture
The spectrum of saline soil is an extremely complex evolution process, which is affected by many inter-influencing factors, and the process exhibits strong nonlinear characteristics. Small changes in soil parameters or the environment will cause differences in spectral reflectance, making such subtle spectral differences difficult to identify. The Chen-Lee chaotic system is used as a signal analysis tool. The intelligent classification system architecture can be shown in Figure 8. The specific eight steps can be described as follows: One, we use an ASD spectrometer to collect the spectral signal of the soil in the study area and observe the hyperspectral signal of the non-saline soil, which is used as a part of the input signal of the X variable.
Two, the hyperspectral signal of the saline soil is decomposed by Chua's circuit, and the Chua's circuit status equation is shown as Equation (1). The spectral sequence processed by Chua's circuit is used as a part of the input signal of the Y variable.
Three, it generates the main system signal and the slave system signal with the selfsynchronization of the chaotic system. The hyperspectral signal is substituted into the master-slave chaotic system, and the slave system can be synchronized with the master system under nonlinear control; the chaotic self-synchronization system can be separately expressed as Equation (3) and (4).
Four, the hyperspectral signal and salinity information of the saline soil are substituted into the constructed fractional-order Chen-Lee chaotic system, further subtraction is performed, and the dynamic equation of fractional order master/slave Chen-Lee chaotic system can be expressed as Equation (5).
Five, we further obtained the dynamic error of the master-servant chaotic system, and the dynamic self-synchronization error system mapping can be revised as in Equation (19). Spectrum signals from the soil with different salinization levels are introduced into different fractional-order mixed master/slave chaotic systems to produce chaotic dynamic error distribution.
Six, the hyperspectral signal and salinity information of the saline soil are substituted into the constructed fractional-order chaotic system, its dynamic behavior characteristics are analyzed, and the chaotic dynamic error and attractor difference of different sampling points under different fractional orders are discussed.
Seven, the fractional-order chaotic system with the largest dynamic error is extracted, and the extension matter-element model is established based on its chaotic attractor features. The extension theory is used to intelligently classify the system. The matter-element model can be expressed as Equation (20).
Eight, we can draw a phase-plane distribution map based on the state phasors to achieve the purpose of classifying salinization.

Mixed Fractional Order Master/Slave Chaotic System Dynamic Error Distribution
This study combined the Chua's circuit with the fractional order Chen-Lee master/slave chaotic system to analyze saline soil, severe saline soil, moderate saline soil, mild

Mixed Fractional Order Master/Slave Chaotic System Dynamic Error Distribution
This study combined the Chua's circuit with the fractional order Chen-Lee master/slave chaotic system to analyze saline soil, severe saline soil, moderate saline soil, mild saline soil, and non-saline soil. Distribution change trend in a different order number is used to search for status change features. Each fractional-order difference is compared, and order numbers with greater dynamic error changes were chosen. Saline soil, severe saline soil, and mild saline soil were chosen for Area A. The chaotic dynamic error distributions for non-integer orders 0.1, 0.3, 0.6, 0.9, and integer-order 1 are shown in Figure 9. Saline soil, severe saline soil, and moderate saline soil were chosen for Area B. The chaotic dynamic error distribution for non-integer orders 0.1, 0.3, 0.6, 0.9, and integer-order 1 are shown in Figure 10. Figures 9 and 10 show that the chaotic dynamic error distribution change trend of the integer and non-integer orders generally form a line. Soil samples with lower salinization had greater Xand Y-axis component changes in chaotic dynamic error distributions and vice versa. Observations also show that of all the fractional orders, the chaotic dynamic error of orders 0.1 and 0.3 shows two divergent lines. However, the change in order 0.1 is more obvious, so we chose non-integer order 0.1 for follow-up analysis.

Mixed Fractional Order Master/Slave Chaotic Attractor
To make the building of an extension matter-element model easier, this study used the extension theory as a smart classification tool. The chaotic dynamic error distribution is converted to the chaotic attractor pattern for easier extension matter-element model building. Simulations in Section 3.1 show that the dynamic error change in the non-integer order 0.1 mixed fractional-order chaotic system is larger. Thus, we chose the chaotic attractor produced by order 0.1 mixed fractional-order chaotic system for analysis. The chaotic attractor distribution for five sampling points in Area A and B is shown in Figure 11. Figure 11 shows that the chaotic attractor coordinate distribution exhibits a linear trend because of the different salt contents. The bottom left corners of Areas A and B have sampling points with the highest salt content. The top right corners have sampling points with the lowest salt content. The salt content going from the bottom left to the top right corner exhibits a decrement status. By using this chaotic attractor coordinate distribution map, we can easily identify soil with different salinization levels.

Building an Extension Matter-Element Model and Extension Classic Domain
To build a smart classification system, this study used the extension theory as the backend classification system. Chaotic attractors are introduced into the extension theory. After building a matter-element model, this study established the extension classic domain range and calculated the correlation level to determine the accuracy of salinization level classification.
Based on the order 0.1 chaotic attractor coordinates shown in Figure 11, the matterelement model vectors from Area A and B are expressed as Equations (21)

Mixed Fractional Order Master/Slave Chaotic Attractor
To make the building of an extension matter-element model easier, this study used the extension theory as a smart classification tool. The chaotic dynamic error distribution is converted to the chaotic attractor pattern for easier extension matter-element model building. Simulations in Section 3.1 show that the dynamic error change in the non-integer otic attractor distribution for five sampling points in Area A and B is shown in Figure 11. Figure 11 shows that the chaotic attractor coordinate distribution exhibits a linear trend because of the different salt contents. The bottom left corners of Areas A and B have sampling points with the highest salt content. The top right corners have sampling points with the lowest salt content. The salt content going from the bottom left to the top right corner exhibits a decrement status. By using this chaotic attractor coordinate distribution map, we can easily identify soil with different salinization levels.

Building an Extension Matter-Element Model and Extension Classic Domain
To build a smart classification system, this study used the extension theory as the backend classification system. Chaotic attractors are introduced into the extension theory. After building a matter-element model, this study established the extension classic domain range and calculated the correlation level to determine the accuracy of salinization level classification.
Based on the order 0.1 chaotic attractor coordinates shown in Figure 11, the matterelement model vectors from Area A and B are expressed as Equations (21) After building the matter-element model, the extension classic domain and section domain can be established based on the model. The interval range defined by the extension classic domain and section domain is used as the basis for determining soil salinization. The extension classic domain range is shown in Figure 12. The yellow area in Figure  12 is the saline soil classic domain, the orange area is the severe saline soil classic domain, the green area is the moderate salinization classic domain, the red area is the mild saline soil classic domain, and the blue is the section domain of all the saline soil.  Figure 12. The yellow area in Figure 12 is the saline soil classic domain, the orange area is the severe saline soil classic domain, the green area is the moderate salinization classic domain, the red area is the mild saline soil classic domain, and the blue is the section domain of all the saline soil.

Calculating Correlation Level
After the classic domain is established, we can calculate the correlation level. Tables  3 and 4 show the soil sample extension correlation levels for Areas A and B. Extension correlation calculation results in Tables 3 and 4 show that A1, B1, and B2 is saline soil, A2, A3, and B3 are severe saline soil, B4 and B5 are moderate saline soil, and A4 and A5 are mild saline soil. If the correlation level is 1, it is the indicated status; if the reverse, then it is not the indicated status.

Calculating Correlation Level
After the classic domain is established, we can calculate the correlation level. Tables 3 and 4 show the soil sample extension correlation levels for Areas A and B. Extension correlation calculation results in Tables 3 and 4 show that A1, B1, and B2 is saline soil, A2, A3, and B3 are severe saline soil, B4 and B5 are moderate saline soil, and A4 and A5 are mild saline soil. If the correlation level is 1, it is the indicated status; if the reverse, then it is not the indicated status.

Extension Classification Results
To verify the accuracy of the extension theory classification, we compared the laboratory chemical analysis results to the computer-calculated extension correlation level results. Results for Areas A and B are shown in Tables 5 and 6, respectively. Simulations show that the extension correlation classification results are consistent with the laboratory chemical analysis results; the classification accuracy of calibration set data is 100%. This study used a combination of Chua's chaotic circuit and a fractional-order Chen-Lee master/slave chaotic system to analyze soil sample spectrums. The extension theory is used to make classification smarter. The results obtained in this study matched the soil salinization standard proposed by Qifei Li [42]. This shows that the chaotic system can be applied to spectrum analysis. The mixed chaotic system matched with the extension theory proposed in this study can rapidly and intuitively identify soil salinization levels.

Classification Result of Validation Set
In order to verify the robustness of the system, this paper randomly selected 10 samples from the research areas as the validation set data, substituted the spectral reflectance data into the master/slave chaotic system, and extracted the chaotic attractor into the extension classification system established in this paper. The extension classic domain range of the validation set is shown in Figure 13. The yellow area in Figure 13 is the saline soil extension classic domain; the orange area is the severe saline soil extension classic domain; the blue area is the moderate salinization extension classic domain; the green area is the mild saline soil extension classic domain.
ples from the research areas as the validation set data, substituted the spectral refl data into the master/slave chaotic system, and extracted the chaotic attractor into tension classification system established in this paper. The extension classic domain of the validation set is shown in Figure 13. The yellow area in Figure 13 is the sal extension classic domain; the orange area is the severe saline soil extension classic d the blue area is the moderate salinization extension classic domain; the green are mild saline soil extension classic domain.
. Figure 13. Extension classic domain range of the validation set.
The extension correlation degree value and salt content value of 10 soil samp shown in Table 7. In Table 7, the correlation degree value is one as the correspondin It can be noted that the correlation degrees of T6, T8, T9, and T10 are 0.985, 0.973 and 0.964, respectively. Thus, indicating that the attractor phase coordinates of t samples do not fall within the classical domain. However, the extension distance close to the classical domain of soil. This article can consider these four samples saline soil state. The attractor phase coordinates of sample T4 fall within the class main of moderate saline, but the actual salt content is 20.4 g/kg, which belongs to th soil state, so the sample T4 is a misclassified chaotic attractor. The extension correlation degree value and salt content value of 10 soil samples are shown in Table 7. In Table 7, the correlation degree value is one as the corresponding state. It can be noted that the correlation degrees of T6, T8, T9, and T10 are 0.985, 0.973, 0.987, and 0.964, respectively. Thus, indicating that the attractor phase coordinates of the four samples do not fall within the classical domain. However, the extension distance is very close to the classical domain of soil. This article can consider these four samples as the saline soil state. The attractor phase coordinates of sample T4 fall within the classical domain of moderate saline, but the actual salt content is 20.4 g/kg, which belongs to the saline soil state, so the sample T4 is a misclassified chaotic attractor. According to the above validation results, the system randomly extracts the spectral reflectance of 10 soil samples for chaotic analysis and extension classification, and the classification accuracy is 90%, which proves that the system is very robust. In the future, the state phasor of the master/servant chaotic dynamic error can be extracted in the way of machine learning for learning classification, or the convolutional neural network can be used to directly identify the coordinate map of the master/servant chaotic dynamic error to achieve a more accurate and robust classification.

Discussion
Some scholars use Chen-Lee chaotic systems to study wind power systems. Ying-Che Kuo et al. [46] used the Chen-Lee chaotic system to load the normal and fault signals of the bearings into the chaos synchronization error dynamics system. Results showed that the proposed chaos synchronization method could separate the characteristics completely, which had a better fault diagnosis rate than the methods. Chih-Jer Lin et al. adopted [47] the Chen-Lee chaotic system for analysis and took its chaotic attractor as the feature of the ball bearing malfunction recognition, the method combing Chen-Lee chaotic system, and BPNN (Back Propagation Neural Network) had a high identification rate and overall efficiency for the ball bearing malfunction. Ying-Che Kuo et al. used [48] the unified chaotic system-based design to study the wind turbine blade fault diagnostic system, Lorenz system, the Lu system, and the Chen system was used to identify which one responds best to this system. The experimental results showed that the unified chaotic system had a significant effect on the vibration signal analysis. However, all the researchers did not use a fractional-order Chen-Lee chaotic system. At the same time, there have not yet been reports on the Chen-Lee chaotic system in the field of soil hyperspectral. However, in our paper, we proposed a new method of classifying the degree of soil salinization to achieve a higher-precision classification using the fractional Chen-Lee chaotic system and coupled it with Chua's circuit and control theory.
Fractional chaotic systems can not only accurately describe the physical model of the system due to the use of fractional orders but also accurately distinguish small differences in signals due to their high sensitivity to initial values. In this paper, the fractional-order chaos system is used to improve the resolution between the peaks of hyperspectral signals, and the effect of the fractional-order chaos system to improve the pretreatment of soil hyperspectral signals is discussed, which provides a new idea for the pretreatment of hyperspectral soil data. Based on the spectrum and salt content information of different sampling points, this paper analyzes the mathematical model and dynamic behavior characteristics of the fractional chaotic system and compares the corresponding chaotic motion trajectories under different fractional orders and identifies the characteristics of the hyperspectral signal. Thereby, the resolution between the peaks of the hyperspectral signal is enhanced, and the best fractional chaotic system is found.
Typical fractional-order chaotic systems include Chen, Lorenz, Sprott, etc. The dynamic behavior characteristics of different combination systems varies. This research proposes a mixed master/slave chaotic system based on the Chua's circuit and fractional order Chen-Lee chaotic system for classifying saline soil salinization level. It is found that the 0.1-order chaotic dynamic error change is the most obvious, and the fractional order was further analyzed. Based on the characteristics of the 0.1-order chaotic attractor, an extension matter-element model was established, the correlation degree was calculated, and the calculated extension classification result was compared with the salt content data, and the accuracy rate was high. The simulation results show that the method proposed in this research can accurately classify the degree of salinization. However, the classification of saline soil based on other combinations of chaotic systems requires further research and discussion.

Conclusions
This study proposes a soil spectrum data analysis method based on fractional order mixed master/slave chaotic system. This is the first time this new method is used to analyze large amount of spectrum data. Soil salinization is a very complex dynamic evolution process that exhibits a very strong non-linear characteristic. Chaos theory can completely and accurately describe this type of very complex non-linear process. The salinization classification method proposed by this study is simple, fast, and highly expandability. In the future, this research method can combine UAV (Unmanned Aerial Vehicle/Drones) hyperspectral images and hyperspectral satellite remote sensing images to further analyze the classification accuracy of soil salinization in a large area. We can extract the fractionalorder with the largest dynamic error to improve the classification accuracy by analyzing the chaotic attractor characteristics of the soil hyperspectral at different fractional orders. The results of this paper can be transplanted to embedded devices and used to develop intelligent system processing software for predicting soil elements, so that the entire set of smart classification systems can be portable and faster. It can integrate hyperspectral data collection, storage, analysis, modeling, and other modules, and automatically realize the feature extraction and model prediction of hyperspectral data in batches. This has important application prospects for the development of resource and environmental remote sensing, and creates certain economic benefits.