【行業應用】高光譜檢測技術助力礦業領域研究
背景
experimental design試驗設計 本研究所用的鈦鐵礦原料來自中***四川（由攀枝花學院提供）。在這項工作中，所有樣品均在實驗室球磨機(南大QM-3SP2行星式)中研磨，之后以200目粒徑進行篩分。使用Innov-X Delta DS6000手持式XRF裝置測定兩份原料的鈦品位分別約為3.6780%和31.4530%。為了校準模型的開發提供廣泛的鈦濃度變化范圍，用電子分析天平（Sartorius BSA224S）將兩份原料以1%(Ti)的梯度制備5份，在球磨機中分別以300轉/分鐘的速度充分混合30分鐘，再次使用XRF測定制備的鈦鐵礦粉末樣本的多組分品位。礦石樣本制備完成后，通過GaiaField-N17E-HR近紅外高光譜相機（江蘇科技有限公司）采集礦石樣本高光譜圖像。 提取光譜數據 圖1b（藍色線條）為HSI捕獲到的原始平均光譜曲線，可以明顯看到在整個光譜區域沒有特異的波峰波谷。為了找出峰谷突出的區域，分析了光譜波段之間的相關性。對于維度為（x, y, λ）的高光譜數據，轉換為（z = x × y, λ），并對z進行相關性分析，結果如a所示。圖中顯示紅色區域的波段在整個波段中占比94.34%，其余顏色區域總共占比5.66%。為了***大限度地保留 “純凈數據”，防止其它顏色區域數據的干擾，通過查詢顏色欄數據并結合波段得知紅色與其余顏色區域的相關性分界點為0.8（波段：483）。隨后，將相關性小于0.8的30個波段從數據分析中剔除，從而消除了1685.47-1735.34 nm的光譜區域。得到的光譜曲線如圖1b（紅線）所示。該過程獲得的光譜數據（150×483）用于后續的模型建立和分析。 （a）波段相關性分析 （b）原始平均光譜 Conclusion結論 通過XRF檢測出了鈦鐵礦中品位較高的9種金屬元素，根據其品位高低對樣本進行分組。其中，1%及以上為主要品位（Ti, Fe），0.02-1%為次要品位（V, Mn, Co, Cu），0.02%及以下為微量品位（Zn, Zr, Pb）。鈦鐵礦中各組分品位的統計分布如表1所示，可以清晰看出，不同組分的品位有明顯的差異，這對回歸模型的建立至關重要。 表1 鈦鐵礦樣品品位測定的統計分析 ▼ 鈦鐵礦樣本的光譜反射率值在0.03 - 0.27之間，證實粉狀鈦鐵礦的整體反射率較低，并隨品位的增加而降低。圖1b顯示在NIR區域，該鈦鐵礦的光譜特征分別在940nm、1020nm和1300~1650nm處出現了波谷。具體地，光譜反射率值在940 nm處達到谷值，隨后急劇下降，這主要歸因于O-Ti-O的拉伸和變形。該鈦鐵礦較寬的吸收波段主要位于1300 ~ 1650 nm，Izawa指出這是由八面體配位的亞鐵在此附近產生晶體場躍遷分裂而引起的。 采用SavitZky***Golay卷積平滑法(簡稱SG平滑算法)、多元散射校正算法、包絡線去除算法、二階導數算法、多元散射校正+SG平滑算法等5種預處理方法對原始光譜數據進行處理。基于樹突網絡模型對比原始光譜數據，探索出HSI光譜數據的***佳預處理方法，結果如表2所示。預測集中，與原始數據相比，MSC+SG方法下的各指標R2P提高了5.88%，RMSEP降低了34.39%，RPD和RPIQ分別提高了28.86%和38.93%。 表2 DD模型下不同預處理方法的結果 ▼ 為直觀地看到各預處理方法對數據的改變，繪制了HSI采集的鈦鐵礦樣本三維原始光譜如圖2a所示，圖2b-d所示三維光譜為CR、SD、MSC+SG預處理的結果。SG預處理后的光譜反射率曲線趨勢與鈦鐵礦原始光譜反射率曲線相似。MSC和MSC+SG預處理的光譜變化不大，但有效地消除了光譜散射，突出了光譜峰谷區域。相比之下，CR和SD預處理后的光譜變化明顯，兩種方法的光譜反射率范圍差異過大，導致原始數據的***些特征丟失。然而，***佳MSC+SG預處理圖譜顯示，光譜特征區域顯著增強，反射率數據得以明顯集中，有效增強了光譜特征信息。因此，MSC+SG可以作為***種有效的光譜預處理方法來預測鈦鐵礦多組分的品位信息。 ▲ 圖2 HSI采集鈦鐵礦樣品的三維光譜：(a) RAW光譜；(b) CR預處理后的光譜；(c) SD預處理后的光譜；(d) MSC + SG預處理后的光譜。 采用t分布式隨機鄰居嵌入、區間偏***小二乘法-VCPA-IRIV和區間偏***小二乘法-VIP-IRIV三種不同的降維算法對近紅外光譜數據進行不同程度地約簡。圖3為不同降維算法的數據分布情況。t-SNE改變了原始數據集，經該方法輸出的高維（483維）原始光譜數據特征映射到了三維平面之上。特征選擇方法iPLS-VCPA-IRIV和iPLS-VIP-IRIV三步混合策略算法提取的波長數分別為28和38，比原始維度相應減少了94.19 %和92.12 %。總的來說，三種降維算法有效地消除了原始光譜矩陣中的冗余信息，為減少計算時間和確保精度奠定了基礎。 表3為鈦鐵礦樣本中各組分的VIF檢驗分數。在相關程度上，L代表相關度低，H代表相關度高。對于每個變量，VIF得分均不超過2，說明各變量之間的共線性現象不突出，相關程度符合研究要求。因此，可以采用數據預處理，降維等方法對各組分品位數據進行處理，并參與模型預測。 表3 鈦鐵礦樣品中各組分的VIF指數 ▼ 提取礦石樣本的顏色特征是否可以用于構建品位預測模型，還有待進***步驗證。基于鈦鐵礦多組分品位（平均品位）與顏色特征參數進行Pearson相關性分析（圖4）。從圖中可以看出，除S顏色通道的特征參數與多組分品位存在弱相關（R < 0.3）外，其余顏色特征參數對于品位預測模型都是合理的。造成這種現象的原因是，鈦鐵礦本身是***種不透明礦物，而HSI采集到的圖像過于單***，整體顏色偏暗。 ▲ 圖4 品位和顏色特征參數的皮爾遜相關性分析 多元回歸模型DD、PLSR和SSA-KELM可能對兩個維度中某***特定數據集的表現***好，但基于多組分品位模型穩定性的需求，BPNN才是***佳的模型選擇。圖5為基于BPNN兩個效果***好的光譜和圖像數據結果，可以看出，兩個數據的指標變化差異不大，證實了BPNN模型即使在小樣本量下也具有很強的泛化性能和魯棒性。綜合評價表明，BPNN對兩組數據集的預測都是成功的，但利用圖像數據預測鈦鐵礦多組分品位的可靠性略低于利用特征選擇的光譜數據。 (a) iPLS-VCPA-IRIV選擇的光譜數據 (b)全波長圖像數據 從iPLS-VCPA-IRIV特征選擇的擬合光譜數據圖（圖6）可以看出，在校正集中，鈦鐵礦多組分品位都均勻準確地分布在理想曲線上。然而，對于預測集來說，明顯可以看出來的是，除擬合良好的主要品位（Ti、Fe）和次要品位（V、Mn、Co、Cu）外，微量品位（Zn、Zr、Pb）的擬合似乎并不令人滿意。因此，在下***步工作中，可以對少樣本量下礦石內部微量品位的變化進行更深入的研究。 (a) 校準集表現 (b)預測集表現 本研究從光譜和空間維度出發，探討了HSI與化學計量學相結合預測鈦鐵礦多組分品位的能力。 PRODUCT 產品推薦 ▲ GaiaField 便攜式成像光譜系統
主要功能： ● ***鍵實現自動曝光、自動調焦、自動掃描速度匹配、自動采集并保存數據 ● 輔助取景攝像頭實現對拍攝區域的監控 ● 內置電池，可連續工作4小時以上 ● 數據預覽及校正功能：輻射度校正、反射率校正、區域校正、鏡頭校準、均勻性校準 ● 鏡頭可更換 ● 只需***根USB線實現連接通信 ● 數據格式完美兼容Evince、Envi、SpecSight等數據分析軟件 ● 支持Win7~10-32位或64位系統 ● 高配版可內置微型處理器，利用手機、Ipad通過Wifi熱點遠程控制
(文章來源于儀器網)
