任 順，張 雄，任 東，楊信廷,3，張 力

(1.三峽大學(xué) 計(jì)算機(jī)與信息學(xué)院，湖北 宜昌 443002；2.三峽大學(xué) 湖北省農(nóng)田環(huán)境監(jiān)測(cè)工程技術(shù)研究中心，湖北 宜昌 443002；3.農(nóng)產(chǎn)品質(zhì)量安全追溯技術(shù)及應(yīng)用國(guó)家工程實(shí)驗(yàn)室，北京 100097)

近年來(lái)，有毒重金屬對(duì)我國(guó)土壤環(huán)境質(zhì)量危害日益嚴(yán)重，耕地、農(nóng)產(chǎn)品等污染問題亟待解決,因此，重金屬監(jiān)測(cè)已成為農(nóng)產(chǎn)品保護(hù)和安全生產(chǎn)的重要任務(wù)[1]。目前，土壤重金屬污染檢測(cè)方法有傳統(tǒng)實(shí)驗(yàn)室檢測(cè)法與速測(cè)法，其中傳統(tǒng)實(shí)驗(yàn)室檢測(cè)方法主要有原子吸收與原子熒光光譜法、電感耦合等離子質(zhì)譜法等，該類方法雖然準(zhǔn)確度和精確度高，但儀器運(yùn)行條件要求高且不適用于現(xiàn)場(chǎng)快速判斷；速測(cè)法有X射線熒光光譜法(XRF)、生物傳感器技術(shù)和激光誘導(dǎo)擊穿光譜法[2-3]等。其中，XRF具有檢測(cè)速度快、成本低且可同時(shí)檢測(cè)多種元素等優(yōu)點(diǎn)，已被廣泛用于土壤重金屬污染和農(nóng)產(chǎn)品檢測(cè)。相比實(shí)驗(yàn)室測(cè)定方法，XRF法前處理簡(jiǎn)單，對(duì)主要的重金屬可有效監(jiān)測(cè)和快速篩查，可為農(nóng)田和農(nóng)作物重金屬污染的監(jiān)測(cè)和防控、制定合理的農(nóng)業(yè)發(fā)展規(guī)劃提供科學(xué)依據(jù)[4]。XRF法快速檢測(cè)重金屬時(shí)，針對(duì)不同類型和性質(zhì)的土壤，需建立重金屬濃度預(yù)測(cè)模型，再根據(jù)實(shí)際檢測(cè)結(jié)果對(duì)預(yù)測(cè)模型優(yōu)化。由于土壤成分復(fù)雜，采用XRF儀獲取的光譜數(shù)據(jù)具有很高的空間復(fù)雜程度。此外，所得到的光譜數(shù)據(jù)往往由數(shù)以千計(jì)的波長(zhǎng)點(diǎn)組成，而樣本數(shù)量常受制于實(shí)驗(yàn)條件和成本，在進(jìn)行化學(xué)建模算法時(shí)，常會(huì)出現(xiàn)統(tǒng)計(jì)學(xué)稱為維數(shù)災(zāi)禍的多項(xiàng)式復(fù)雜程度的非確定性(NP)難題。模型集群分析策略強(qiáng)調(diào)，要最大限度分析已有樣本集信息，通過隨機(jī)采樣利用大量子集模型信息，獲得數(shù)據(jù)集的內(nèi)在結(jié)構(gòu)，而不僅僅依靠單一模型信息，可更好地避免模型對(duì)樣本的依賴[5-7]。

本研究以土壤重金屬的X射線熒光光譜為研究對(duì)象，將0～26 keV范圍內(nèi)的光譜數(shù)據(jù)與濃度梯度法配制的銅(Cu)、鋅(Zn)、砷(As)、鉛(Pb)、鉻(Cr)5種重金屬污染土壤樣本含量值進(jìn)行關(guān)聯(lián)，將多特征串聯(lián)策略區(qū)間組合優(yōu)化-競(jìng)爭(zhēng)適應(yīng)性重加權(quán)采樣-連續(xù)投影(ICO-CARS-SPA)算法所提取出的特征變量作為輸入，建立偏最小二乘(PLS)模型預(yù)測(cè)土壤樣本中5種重金屬含量，并與一些常見變量選擇方法及其組合對(duì)比，結(jié)合光譜特征變量的優(yōu)選和模型的預(yù)測(cè)性能，實(shí)現(xiàn)了對(duì)農(nóng)田土壤重金屬含量高效、精準(zhǔn)的定量檢測(cè)，可為土壤重金屬治理和污染防控以及制定科學(xué)的農(nóng)業(yè)發(fā)展規(guī)劃提供依據(jù)。

圖1 ICO算法流程示意圖Fig.1 Flowchart of ICO algorithm

1 原理與算法

1.1 試驗(yàn)區(qū)間組合優(yōu)化算法

區(qū)間組合優(yōu)化(ICO)是基于模型集群分析策略(MPA)框架下的化學(xué)建模算法[8]，MPA強(qiáng)調(diào)從數(shù)據(jù)集中統(tǒng)計(jì)分析時(shí)，應(yīng)通過隨機(jī)采樣，分析大量隨機(jī)產(chǎn)生的子集模型信息，從不同角度考察數(shù)據(jù)集的內(nèi)在性質(zhì)，獲得數(shù)據(jù)集的內(nèi)在結(jié)構(gòu)。ICO算法流程如圖1所示。

1.2 競(jìng)爭(zhēng)性自適應(yīng)重加權(quán)采樣算法

競(jìng)爭(zhēng)性自適應(yīng)重加權(quán)采樣(CARS)算法選擇波數(shù)的方法是基于PLS模型回歸系數(shù)，衡量波長(zhǎng)選擇的重要依據(jù)是PLS回歸系數(shù)絕對(duì)值的大小，通過挑選出交叉驗(yàn)證均方根誤差(RMSECV)值最小的子集來(lái)對(duì)應(yīng)獲得與目標(biāo)值最優(yōu)的波長(zhǎng)組合。計(jì)算則采用指數(shù)衰減函數(shù)，逐次保留絕對(duì)值大的回歸系數(shù)對(duì)應(yīng)的波數(shù)點(diǎn)[9-10]。

1.3 連續(xù)投影算法

連續(xù)投影算法(SPA)是一種前向變量循環(huán)特征波長(zhǎng)選取方法，該算法利用向量的投影分析，能找到含有最低限度冗余信息的變量組。SPA一般能夠較為有效地消除波長(zhǎng)之間共線性，同時(shí)提高模型精確度[11-12]。

1.4 多特征串聯(lián)算法

X射線熒光光譜的高維度無(wú)法規(guī)避地含有較多噪聲和冗余。ICO采用波長(zhǎng)區(qū)間進(jìn)行變量?jī)?yōu)選，能排除大量無(wú)用信息和噪聲波長(zhǎng)點(diǎn)，且在每次優(yōu)化后的聯(lián)合區(qū)間上進(jìn)行局部搜索策略，更好地體現(xiàn)了柔性收縮策略優(yōu)勢(shì)，減少偶然誤差的發(fā)生，最終可得到一組位置、組合、寬度等都經(jīng)過優(yōu)化的有效波長(zhǎng)區(qū)間。但I(xiàn)CO算法選擇的波長(zhǎng)數(shù)量一般較多，區(qū)間內(nèi)部不可避免地存在共線性和冗余問題，最終選中的變量區(qū)間可采用單一變量選擇算法進(jìn)行精簡(jiǎn)[13]。

CARS算法作為一種高效的波長(zhǎng)選擇算法，通過選擇少量波長(zhǎng)變量子集，可以給出更佳或相當(dāng)?shù)念A(yù)測(cè)效果。由于該算法在迭代過程中同時(shí)引入了蒙特卡羅采樣(MCS)和自適應(yīng)加權(quán)采樣(ARS)兩個(gè)隨機(jī)因素，其單獨(dú)使用的穩(wěn)定性難以令人滿意[14-15]。

SPA算法生成的每一個(gè)波長(zhǎng)組合中，新入選的波長(zhǎng)點(diǎn)與上一個(gè)波長(zhǎng)點(diǎn)的相關(guān)性均為最低，這也體現(xiàn)了SPA算法能夠有效消除波長(zhǎng)之間的共線性[16]。但對(duì)于X射線熒光光譜而言，有效變量間的投影距離并不一定最大，通過SPA算法挑選出變量子集很可能存在部分無(wú)關(guān)信息變量甚至是干擾變量。

將ICO選擇的優(yōu)化組合區(qū)間結(jié)合與競(jìng)爭(zhēng)適應(yīng)性重加權(quán)采樣法(CARS)和連續(xù)投影算法(SPA)串聯(lián)(ICO-CARS-SPA)進(jìn)行特征波長(zhǎng)提取，既能有效地鎖定特征波段區(qū)間，同時(shí)優(yōu)選出的最優(yōu)波長(zhǎng)子集又能夠充分減少最優(yōu)區(qū)間內(nèi)的共線性和冗余，形成優(yōu)勢(shì)互補(bǔ)。

2 實(shí)驗(yàn)部分

2.1 材 料

在周邊1 km范圍內(nèi)無(wú)污染源的農(nóng)田采集制樣土壤，為防止實(shí)驗(yàn)器具影響和干擾樣本重金屬的濃度，樣本制作過程采用陶瓷用具。將土壤烘干研磨過0.45 mm孔篩后，參考《土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)》[17]中對(duì)Ⅰ、Ⅱ、Ⅲ類土壤中各元素含量的規(guī)定并結(jié)合農(nóng)田土壤實(shí)際情況制作實(shí)驗(yàn)樣本。每個(gè)樣本制作30 g土壤，共91個(gè)樣本。經(jīng)初步測(cè)量，所采集農(nóng)田土壤中As、Cu、Cr、Pb、Zn的含量依次為3.0、15、30、18、39 mg/kg。綜合國(guó)家土壤等級(jí)劃分標(biāo)準(zhǔn)和所采集土壤的具體情況，設(shè)I為各個(gè)元素初始濃度梯度，其中1～40號(hào)樣本的濃度梯度為I；41～50號(hào)樣本的濃度梯度為2*I；51～60號(hào)樣本濃度梯度為3*I；61～70號(hào)樣本濃度梯度為4*I；71～91號(hào)樣本濃度梯度為J*I；其中J=1,2,3,4,……,n，具體實(shí)施方案參照表1。按土壤需求濃度，將重金屬標(biāo)準(zhǔn)溶液用丙酮稀釋后混入土壤中，混合均勻后置于通風(fēng)櫥中自然揮發(fā)，待土壤完全風(fēng)干后，再研磨混勻，按指定標(biāo)號(hào)放入樣品盒中，密封保存。

表1 土壤樣本配制方案Table 1 Soil sample allocation scheme

2.2 光譜數(shù)據(jù)采集

每次取一個(gè)樣本置于便攜式X射線土壤重金屬檢測(cè)儀(由三峽大學(xué)協(xié)同創(chuàng)新中心與北京農(nóng)業(yè)質(zhì)量標(biāo)準(zhǔn)與檢測(cè)技術(shù)研究中心聯(lián)合研制)上，獲取樣本在0～26 keV范圍內(nèi)共4 096個(gè)通道內(nèi)的光譜信息，每個(gè)樣本經(jīng)過轉(zhuǎn)動(dòng)方向測(cè)量3次并求取其平均光譜。由于36和37號(hào)樣本的光譜數(shù)據(jù)遺失，同時(shí)根據(jù)配制的樣品濃度繪制91個(gè)樣本的實(shí)際濃度曲線，發(fā)現(xiàn)60號(hào)樣品中的Zn濃度和79號(hào)樣品中Pb濃度分布與配制濃度差異較大，分析原因可能是溶液配制時(shí)出現(xiàn)了問題，故剔除這4個(gè)樣本。剩余87個(gè)樣本的平均光譜圖見圖2，其中配置Ag靶微型X光管、電子冷卻Si-PIN探測(cè)器，設(shè)置電壓為30 kV，電流為40 μA，積分時(shí)間為300 ms。

圖2 87個(gè)樣本的平均光譜圖Fig.2 Average spectrogram of 87 samples

2.3 建模方法與評(píng)價(jià)

2.3.1 光譜數(shù)據(jù)集的劃分該方法以樣本集被測(cè)指標(biāo)的理化參照值作為劃分標(biāo)準(zhǔn)，將剔除4個(gè)異常樣本后的87個(gè)不同濃度重金屬的土壤樣本通過濃度梯度法進(jìn)行校正集和驗(yàn)證集的劃分，將樣本按照2∶1分成2組，得校正集58個(gè)，驗(yàn)證集29個(gè)，樣本分布如圖3。

圖3 校正集(A)和驗(yàn)證集(B)樣本分布圖Fig.3 Distribution figure of calibration set(A) and prediction set(B)

3 結(jié)果與討論

3.1 ICO算法選擇特征波長(zhǎng)子區(qū)間

采用ICO算法進(jìn)行特征變量選擇，并建立PLS模型預(yù)測(cè)重金屬含量。不同區(qū)間劃分下試驗(yàn)結(jié)果如表2所示，當(dāng)子區(qū)間劃分?jǐn)?shù)為20時(shí)，5種重金屬均獲得最小的均方根誤差值。

表2 不同子區(qū)間數(shù)PLS建模對(duì)比Table 2 Comparison of PLS modeling for different subinterval numbers

(續(xù)表2)

以生物毒性顯著的重金屬元素Cr為特征提取對(duì)象的結(jié)果如下：采用ICO算法，進(jìn)行最優(yōu)聯(lián)合區(qū)間的選取，利用5折交叉驗(yàn)證法建立PLS模型選擇特征變量。選擇最佳波長(zhǎng)區(qū)間劃分?jǐn)?shù)量20，PLS模型中最大潛變量數(shù)為10，所選子模型的比例為0.05，加權(quán)自舉采樣次數(shù)為1 000。在ICO算法迭代過程中，隨著迭代次數(shù)的增加，每個(gè)波長(zhǎng)區(qū)間的采樣權(quán)重變化情況見圖4A，采樣權(quán)重值越接近1時(shí)，顏色越接近深紅色；采樣權(quán)重值越接近0時(shí)，顏色越接近深藍(lán)色；顏色介于深藍(lán)色和深紅色之間表明采樣權(quán)重值處于0和1之間。圖中還可見，第6個(gè)波長(zhǎng)區(qū)間的采樣權(quán)重值在第2次迭代過程中已經(jīng)約為1，其采樣權(quán)重值在第3次迭代過程中依然有機(jī)會(huì)變得小于1，最終該區(qū)間由于權(quán)重過小被剔除；第9個(gè)波長(zhǎng)區(qū)間在第1～5次迭代過程中顏色偏綠藍(lán)，權(quán)重系數(shù)約為0.3～0.5，但在第6次迭代過程中權(quán)重系數(shù)仍有機(jī)會(huì)變大，最終被選中。由ICO每次迭代中提取的子模型的RMSECV圖可見，經(jīng)過11次迭代后均方根誤差趨于穩(wěn)定(圖4B)，此時(shí)其值為22.624 5，最終挑選的波長(zhǎng)區(qū)間為[4，9，13，20](圖4C)。在ICO最終選中的聯(lián)合區(qū)間內(nèi)引入了局部搜索策略，進(jìn)行寬度的自動(dòng)優(yōu)化，最終挑選了805個(gè)特征波長(zhǎng)。

圖4 重金屬Cr進(jìn)行ICO變量選擇(區(qū)間數(shù)20)Fig.4 ICO variable selection for heavy metal Cr(interval number 20)A.sampling weight change graph of each interval during the iteration process(迭代過程中各區(qū)間的采樣權(quán)重變化圖),B.RMSECV value of the sub-model extracted in each iteration of the ICO(ICO 每次迭代中提取的子模型的RMSECV值),C.ICO algorithm for wavelength range selection(ICO算法進(jìn)行波長(zhǎng)區(qū)間選擇)

圖5 CARS運(yùn)算提取變量原理圖Fig.5 Schematic diagram of CARS operation extraction variablesA.variation trend of the wavelength variable(波長(zhǎng)變量的變化趨勢(shì)),B.variation trend of the RMSEcv(RMSEcv的變化趨勢(shì)),C.trend of wavelength regression coefficient(波長(zhǎng)回歸系數(shù)的變化趨勢(shì))

3.2 ICO-CARS算法選擇特征波長(zhǎng)

將ICO算法挑選出的805個(gè)特征波長(zhǎng)進(jìn)一步使用CARS剔除區(qū)間波長(zhǎng)中的無(wú)關(guān)變量。由于CARS方法在迭代過程引入蒙特卡羅采樣(MCS)和自適應(yīng)加權(quán)采樣(ARS)兩個(gè)隨機(jī)因素，造成每次運(yùn)行挑選的波長(zhǎng)結(jié)果不盡相同，為增加波長(zhǎng)選擇的穩(wěn)定性和可靠性，對(duì)CARS進(jìn)行100次重復(fù)計(jì)算，得到其RMSECV平均值為8.642 0，標(biāo)準(zhǔn)偏差(STD)為0.686 3。圖5為基于CARS算法的波長(zhǎng)變量篩選過程圖，其MCS抽樣運(yùn)行次數(shù)N為50，PLS主成分?jǐn)?shù)為10，采用5折交叉驗(yàn)證。

由CARS算法篩選有效波長(zhǎng)變量的變化趨勢(shì)可看到變量數(shù)呈指數(shù)函數(shù)下降(圖5A)，選擇變量的個(gè)數(shù)也從急劇減少到緩慢遞減，最終趨于穩(wěn)定。波長(zhǎng)選擇過程采用5折交叉驗(yàn)證得到的RMSECV變化趨勢(shì)，若RMSECV值減小，表明剔除了無(wú)關(guān)信息變量；若RMSECV值增大則表明剔除了有效信息變量。由圖5B可見，在1～20次的采樣過程中，RMSECV呈遞減趨勢(shì)，第20次時(shí)的RMSECV值最低，此后開始增加，表明采樣20次時(shí)有效剔除了ICO提取光譜區(qū)間中的無(wú)關(guān)信息。圖5C表示各波長(zhǎng)回歸系數(shù)隨著采樣次數(shù)增加而變化，圖中“*”為RMSECV值最低點(diǎn)。由圖可見，當(dāng)運(yùn)行20次時(shí)，RMSECV值最低，此時(shí)保留的變量為79個(gè)。

3.3 ICO-CARS-SPA算法選擇特征波長(zhǎng)

經(jīng)過ICO-CARS變量篩選后，再通過SPA算法進(jìn)行波長(zhǎng)精簡(jiǎn)(圖6)。通過最小誤差均方根值(RMSE)來(lái)確定最終所選特征波長(zhǎng)個(gè)數(shù)，RMSE值越小，表明模型穩(wěn)定性越好、精度越高。

圖6 ICO-CARS-SPA變量選擇Fig.6 ICO-CARS-SPA variable selectionA.RMSE of SPA model varies the number of variables(SPA模型RMSE隨變量個(gè)數(shù)的變化),B.optimal characteristic wavelength selected by SPA(SPA模型選擇的最優(yōu)特征波長(zhǎng))

結(jié)果顯示，RMSE隨著波長(zhǎng)個(gè)數(shù)的增加呈逐漸減小的趨勢(shì)，當(dāng)波長(zhǎng)個(gè)數(shù)大于33時(shí)，RMSE值變化不再顯著，此時(shí)RMSE值為15.439 9(圖6A)。說明最終優(yōu)選的33個(gè)特征波長(zhǎng)保留了更多有效信息，作為最優(yōu)的波長(zhǎng)點(diǎn)個(gè)數(shù)，選取的波長(zhǎng)點(diǎn)在實(shí)驗(yàn)光譜中的索引如圖6B所示，所選波段占原始光譜信息的0.81%。

采用ICO-CARS-SPA最終優(yōu)選出的33個(gè)敏感波長(zhǎng)為輸入變量，利用PLS對(duì)重金屬元素Cr含量建模，并按照相同特征提取步驟對(duì)Cu、Zn、As、Pb進(jìn)行特征優(yōu)選和PLS建模。結(jié)果表明，采用多特征串聯(lián)策略ICO-CARS-SPA-PLS建模后5種重金屬的R2、RMSE和MRB均達(dá)到了滿意的效果(表3)。

3.4 不同特征波長(zhǎng)選擇方法比較

本文建立的ICO-CARS-SPA方法與ICO、最小絕對(duì)值收斂和選擇算法(LASSO)、CARS、ICO-CARS和ICO-SPA等[18-19]其他波長(zhǎng)選擇算法的對(duì)比結(jié)果見表3。其中LASSO是一種系數(shù)壓縮變量選擇方法，即通過引入懲罰項(xiàng)使影響較小或者無(wú)影響的自變量系數(shù)趨近于零，進(jìn)而可實(shí)現(xiàn)只保留與響應(yīng)變量最相關(guān)的解釋變量，將大量無(wú)關(guān)變量剔除，但高維數(shù)據(jù)中常會(huì)面臨統(tǒng)計(jì)學(xué)稱為維數(shù)災(zāi)禍的NP難題，因此不可避免地造成LASSO特征波長(zhǎng)選擇個(gè)數(shù)將小于樣本個(gè)數(shù)，無(wú)法保證數(shù)據(jù)集中所有重要信息被完全保留。由于土壤中重金屬的X射線熒光光譜數(shù)據(jù)樣本數(shù)據(jù)集的光譜波長(zhǎng)數(shù)遠(yuǎn)大于樣本個(gè)數(shù)，LASSO算法不適合直接應(yīng)用于波長(zhǎng)選擇問題。

表3 不同波長(zhǎng)選擇算法PLS建模性能比較Table 3 Comparison of PLS modeling performance of different wavelength selection

由表3可知，盡管ICO算法可將局部搜索策略應(yīng)用于最終挑選出的最優(yōu)聯(lián)合區(qū)間進(jìn)行寬度的自動(dòng)優(yōu)化，但區(qū)間挑選算法選擇的波長(zhǎng)數(shù)量一般較多，不利于快速預(yù)測(cè)，且忽視了聯(lián)合區(qū)間內(nèi)部光譜間存在的相關(guān)性和冗余問題，因此可對(duì)優(yōu)選出的聯(lián)合區(qū)間進(jìn)行進(jìn)一步波長(zhǎng)篩選。而ICO-CARS建模效果優(yōu)于ICO-SPA是因?yàn)閷?duì)于X射線熒光光譜而言，有效變量間的投影距離并不一定最大，在初選算法波長(zhǎng)數(shù)量較多的情況下，SPA篩選出的變量子集中可能包含一些無(wú)關(guān)信息甚至是干擾變量。

通過對(duì)ICO、LASSO、CARS、ICO-CARS、ICO-SPA和ICO-CARS-SPA波長(zhǎng)優(yōu)選對(duì)比分析可知，采用波長(zhǎng)優(yōu)選算法建模后5種重金屬的R2、RMSE和MRB相對(duì)于單一PLS建模大多有不同程度地提升，其中以ICO-CARS-SPA-PLS提升最明顯。表明ICO-CARS-SPA-PLS的多特征串聯(lián)提升策略從整體上比其他模型優(yōu)化方法建模效果更好，更能兼顧到其他元素，能夠增強(qiáng)上一輪變量選擇的過程與下一輪變量選擇的關(guān)聯(lián)性，根據(jù)每輪變量選擇的好壞進(jìn)行動(dòng)態(tài)調(diào)整，在一定程度上解決了變量區(qū)間選擇的“筑巢效應(yīng)”。

3.5 模型廣適性比較研究

為驗(yàn)證ICO-CARS-SPA變量篩選模型具有較好的廣適性，加入100次隨機(jī)分組進(jìn)行實(shí)驗(yàn)，5種重金屬元素的評(píng)價(jià)參數(shù)及標(biāo)準(zhǔn)偏差結(jié)果見表4。

表4 5種重金屬100次隨機(jī)分組ICO-CARS-SPA-PLS定量模型預(yù)測(cè)結(jié)果Table 4 Prediction results of ICO-CARS-SPA-PLS quantitative model of 100 random groups of five heavy metals

由表3～4可見，由于Zn元素在土壤本底中的濃度較大，導(dǎo)致在制備Zn樣本時(shí)，其配制濃度跨度較大，建模結(jié)果顯示驗(yàn)證集上Zn的決定系數(shù)提高，但RMSEp效果較差，這可能是由于出現(xiàn)了過擬合現(xiàn)象。因此在后續(xù)實(shí)驗(yàn)中通過擴(kuò)大樣本數(shù)量，對(duì)樣本中的重金屬濃度劃分更細(xì)致，對(duì)Zn元素做進(jìn)一步改進(jìn)，可提高其預(yù)測(cè)效果。由表4可知，ICO-CARS-SPA-PLS模型具有較好的廣適性，能夠?qū)崿F(xiàn)對(duì)樣本中5種重金屬的定量檢測(cè)。

4 結(jié) 論

通過多特征串聯(lián)策略的ICO-CARS-SPA算法，建立了X射線熒光光譜重金屬含量的定量檢測(cè)模型。實(shí)驗(yàn)制備了91個(gè)土壤樣本，剔除4個(gè)問題樣本后，選取了87個(gè)樣本，并對(duì)比分析了幾種常見的光譜特征變量選擇方法及其組合。結(jié)果表明，以ICO-CARS-SPA算法提取的特征波長(zhǎng)建立的PLS重金屬含量預(yù)測(cè)模型精度高、誤差小、模型推廣性能好?；赬射線熒光光譜結(jié)合多特征串聯(lián)策略算法ICO-CARS-SPA是有效的特征提取方法，降低了模型復(fù)雜度，有利于推動(dòng)X射線熒光光譜技術(shù)在農(nóng)業(yè)信息領(lǐng)域上的進(jìn)一步應(yīng)用。