劉 梅 馬啟良 原居林 倪 蒙 練青平 郭愛環(huán)

基于無人機(jī)高光譜遙感技術(shù)對(duì)內(nèi)陸?zhàn)B殖池塘水質(zhì)監(jiān)測(cè)的研究*

劉 梅1馬啟良2原居林1①倪 蒙1練青平1郭愛環(huán)1

(1. 農(nóng)業(yè)農(nóng)村部淡水漁業(yè)健康養(yǎng)殖重點(diǎn)實(shí)驗(yàn)室 浙江省魚類健康與營養(yǎng)重點(diǎn)實(shí)驗(yàn)室 浙江省淡水水產(chǎn)研究所 浙江湖州 313001; 2. 湖州師范學(xué)院信息技術(shù)中心 浙江湖州 313000)

無人機(jī); 高光譜; 養(yǎng)殖池塘; 機(jī)器算法; 水質(zhì)反演

隨著浙江省“五水共治”以及“剿滅劣V類水”戰(zhàn)略的深入實(shí)施, 養(yǎng)殖尾水亂排直排現(xiàn)象面臨巨大的環(huán)境壓力, 養(yǎng)殖尾水生態(tài)化處理迫在眉睫。截止2020年, 浙江省主要通過“沉淀池+過濾壩1+曝氣池+過濾壩2+生態(tài)池(簡稱“三池兩壩”)”技術(shù)完成了內(nèi)陸約5.33萬hm2養(yǎng)殖池塘尾水治理任務(wù)(劉梅等, 2021)。但是由于養(yǎng)殖池塘及尾水處理池小而散、涉及面積廣等問題, 且養(yǎng)殖池塘對(duì)水質(zhì)的監(jiān)測(cè)要求具有實(shí)時(shí)、動(dòng)態(tài)、快速的特點(diǎn), 仍缺乏有效的水質(zhì)監(jiān)測(cè)體系, 阻礙了現(xiàn)代養(yǎng)殖小區(qū)數(shù)字化養(yǎng)殖場的創(chuàng)建與發(fā)展(李道亮等, 2020)。

目前, 對(duì)養(yǎng)殖池塘或者養(yǎng)殖尾水處理效果采取的仍是常規(guī)的水質(zhì)監(jiān)測(cè)方法, 即對(duì)特定養(yǎng)殖池塘或尾水處理系統(tǒng)定期、定點(diǎn)進(jìn)行長年累月采樣及水質(zhì)監(jiān)測(cè), 這種方法受人力、物力、時(shí)間及天氣的限制, 采集的數(shù)據(jù)量不可能太多, 而且成本高、速度慢; 并且對(duì)于整個(gè)養(yǎng)殖小區(qū)而言, 這些采樣點(diǎn)數(shù)據(jù)只能代表局部水質(zhì)狀況, 難以獲得大范圍淡水養(yǎng)殖水域水質(zhì)參數(shù)的分布和變化情況, 不能滿足對(duì)養(yǎng)殖池塘水質(zhì)大尺度、快速、實(shí)時(shí)的監(jiān)測(cè)要求。因此, 迫切需求一種實(shí)時(shí)、快速監(jiān)測(cè)養(yǎng)殖池塘及養(yǎng)殖尾水處理池水質(zhì)動(dòng)態(tài)變化的有效手段。

鑒于此, 本文以淡水養(yǎng)殖池塘及配套尾水處理池為研究養(yǎng)殖小區(qū), 通過無人機(jī)搭載高光譜傳感器獲取光譜反射率數(shù)據(jù), 建立水質(zhì)指標(biāo)的反演模型并研究其濃度空間分布, 為養(yǎng)殖池塘水環(huán)境實(shí)時(shí)調(diào)控、養(yǎng)殖尾水池處理效果評(píng)價(jià)及構(gòu)建數(shù)字化養(yǎng)殖小區(qū)提供技術(shù)支撐。

1 材料與方法

1.1 研究區(qū)概況

本研究選擇湖州市典型主要養(yǎng)殖品種之一的青蝦養(yǎng)殖小區(qū)尾水處理區(qū)及周邊池塘為研究區(qū)域, 位于湖州市菱湖鎮(zhèn)楊港現(xiàn)代生態(tài)漁業(yè)園內(nèi)(120°7′8.7″N, 30°46′32.7″E)。2018年該養(yǎng)殖小區(qū)采用“三池兩壩”技術(shù)建立完成尾水處理區(qū), 截止2019年湖州近4萬hm2的養(yǎng)殖池塘均采用該技術(shù)實(shí)現(xiàn)尾水的生態(tài)化處理, 減輕了對(duì)周邊河網(wǎng)水體及太湖的富營養(yǎng)化影響。但是由于在后期的出水水質(zhì)監(jiān)測(cè)均采用傳統(tǒng)方法, 費(fèi)時(shí)費(fèi)力, 也不利于養(yǎng)殖小區(qū)數(shù)字化建設(shè)及管控。

1.2 采樣點(diǎn)布設(shè)和水質(zhì)檢測(cè)

圖1 研究區(qū)地理位置及采樣點(diǎn)空間分布

1.3 高光譜影像的獲取與處理

同時(shí)利用大疆無人機(jī)M600 Pro搭載高光譜成像儀GaiaSky-mini獲取尾水處理區(qū)及周邊養(yǎng)殖池塘的高光譜影像。無人機(jī)飛行高度為100 m, 無人機(jī)高光譜相機(jī)的掃描范圍為400～1 000 nm, 光譜通道為176, 光譜分辨率為3.5 nm。所選實(shí)驗(yàn)時(shí)間均晴朗微風(fēng)(1～2級(jí)), 無云層覆蓋。

無人機(jī)高光譜影像的預(yù)處理主要包括鏡像變換、黑白幀校正、場地校正等, 均在高光譜系統(tǒng)自帶軟件完成。然后對(duì)拍攝的大量數(shù)據(jù)進(jìn)行研究區(qū)圖像拼接, 將拼接后的影像導(dǎo)入ENVI5.3軟件, 根據(jù)經(jīng)緯度坐標(biāo), 定位49個(gè)水面采樣點(diǎn)的具體位置, 以采樣點(diǎn)為中心的6×6像素密度作為計(jì)算區(qū)域, 再將該區(qū)域內(nèi)光譜反射率求平均值, 該平均值即可作為該采樣點(diǎn)的反射率數(shù)據(jù), 共獲得49組與地面水質(zhì)數(shù)據(jù)相對(duì)應(yīng)的光譜反射率, 隨機(jī)選取39組數(shù)據(jù)用于模型訓(xùn)練樣本, 其余10組用于模型檢驗(yàn)樣本。

1.4 高光譜指數(shù)分析

將水樣中各水質(zhì)指標(biāo)與采樣點(diǎn)在遙感影像中的反射率進(jìn)行相關(guān)性分析, 獲得研究區(qū)不同水質(zhì)參數(shù)的敏感單波段。諸多學(xué)者研究結(jié)果顯示, 將反射率進(jìn)行適當(dāng)?shù)臄?shù)學(xué)變換, 構(gòu)建光譜任意兩波段組合而成的差值光譜指數(shù)(difference spectral indices, DSI)、歸一化光譜指數(shù)(normalized spectral indices, NSI)和比值光譜指數(shù)(ratio spectral indices, RSI), 可以降低背景信息的干擾, 將更有效提取光譜信息(黃宇等, 2020; 黃昕晰, 2020), 計(jì)算如式(1)、式(2)和式(3)所示。

DSI =R1–R2, (1)

NSI = (R1–R2) / (R1+R2), (2)

RSI =R1/R2, (3)

式中,R1和R2分別代表400～1 000 nm波段范圍內(nèi)任意2個(gè)波段的反射率。

1.5 水質(zhì)反演模型構(gòu)建流程

偏最小二乘法(partial least square, PLS)通過將自變量和因變量的高維數(shù)據(jù)空間投影到相應(yīng)的低維空間, 得到相互正交的特征向量, 再建立線性回歸模型, 不僅克服多重共線性問題, 即在回歸的同時(shí)強(qiáng)調(diào)自變量對(duì)因變量的解釋和預(yù)測(cè)作用, 而且能較好地辨識(shí)系統(tǒng)信息和噪聲。

徑向基神經(jīng)網(wǎng)絡(luò)(radial basis function network, RBF)屬于多層前饋神經(jīng)網(wǎng)絡(luò), 具有自組織、自學(xué)習(xí)、知識(shí)推理的特點(diǎn), 對(duì)大量非結(jié)構(gòu)性、非精確性規(guī)律具有自適應(yīng)、信息記憶和非線性逼近能力以及優(yōu)化計(jì)算等優(yōu)點(diǎn)。RBF在Matlab中實(shí)現(xiàn), 功能函數(shù)采用Sigmoid函數(shù)。

支持向量機(jī)(support vector machine, SVM)是以統(tǒng)計(jì)學(xué)習(xí)理論為基礎(chǔ)、訓(xùn)練誤差為約束條件、以結(jié)構(gòu)風(fēng)險(xiǎn)最小化為優(yōu)化目標(biāo)的一種專門用來研究有限樣本預(yù)測(cè)的高效機(jī)器學(xué)習(xí)方法。具有較強(qiáng)的泛化能力、較好的魯棒性以及沒有局部最小點(diǎn)的優(yōu)勢(shì)。本研究SVM采用R語言中的e1017包構(gòu)建預(yù)測(cè)模型, 使用radial (徑向基函數(shù))為核函數(shù), 成本函數(shù)cost和gamma選擇模型中最佳參數(shù)。

最后根據(jù)模型的相關(guān)評(píng)價(jià)標(biāo)準(zhǔn)篩選各水質(zhì)要素最佳反演模型, 再將最優(yōu)的監(jiān)測(cè)模型反演到無人機(jī)高光譜影像上, 分別制作各水質(zhì)參數(shù)的空間分布圖, 具體計(jì)算路線如圖2所示。

1.6 模型篩選評(píng)價(jià)指標(biāo)

水質(zhì)參數(shù)預(yù)測(cè)模型采用決定系數(shù)(coefficient of determination,2)、均方根誤差(root mean square error, RMSE)、相對(duì)分析誤差(relative prediction deviation, RPD)以及平均絕對(duì)誤差(relative absolute error, MAE)對(duì)比模型穩(wěn)定性及精度。其中2越大, 說明模型越穩(wěn)定, 一般大于0.8即說明模型穩(wěn)定。RMSE越小, 則模型精度越高。RPD>2.0時(shí)表示模型可靠性高, 具有極好的預(yù)測(cè)能力, RPD值介于1.5～2.0之間, 表明模型可對(duì)樣品做粗略估測(cè), 模型有待改進(jìn), RPD<1.5時(shí), 表明模型預(yù)測(cè)能力差(葉勤等, 2017)。MAE<20%可用于養(yǎng)殖水質(zhì)的預(yù)測(cè), 滿足養(yǎng)殖水體水質(zhì)反演誤差要求。其計(jì)算公式如下:

圖2 無人機(jī)高光譜水質(zhì)反演模型構(gòu)建的技術(shù)路線

2 結(jié)果與分析

2.1 養(yǎng)殖池塘和尾水處理區(qū)光譜特征

研究養(yǎng)殖小區(qū)共49個(gè)采樣點(diǎn)的光譜反射率曲線如圖3所示, 其中圖3a列出了40個(gè)養(yǎng)殖池塘采樣點(diǎn)的反射率曲線, 圖3b列出了9個(gè)尾水區(qū)采樣點(diǎn)的反射率曲線?？梢钥闯龀靥了w的光譜特征在400～ 560 nm范圍內(nèi), 光譜反射率呈上升趨勢(shì), 至560～ 580 nm附近形成一個(gè)反射高峰, 這主要是是由水體中葉綠素和胡蘿卜素在此波段具有較弱的吸收特點(diǎn), 同時(shí)加上藻類和懸浮物的散射作用而引起反射高峰(Gitelson, 1993)。而到580 nm后, 反射率曲線呈下降趨勢(shì), 直至680～710 nm處附近則形成一個(gè)峰谷, 主要是由于水體中的葉綠素在此波段具有強(qiáng)吸收的特點(diǎn)(唐軍武等, 2004); 而緊接著在790～810 nm范圍內(nèi)形成一個(gè)峰值, 主要是由于水中懸浮物的散射作用(黃宇業(yè)等, 2019)。而尾水處理區(qū)由于經(jīng)過過濾、沉淀、曝氣等凈化過程, 水體中懸浮物和藻類去除率可達(dá)到70%以上, 因此, 尾水處理區(qū)在680～710 nm和790～810 nm范圍內(nèi)均未形成明星的峰值(圖3b), 而其他波段的變化規(guī)律和養(yǎng)殖池塘具有一致性。綜合分析光譜曲線可知, 不同池塘水體的光譜曲線變化趨勢(shì)總體一致, 但是由于不同采樣點(diǎn)具有不同的水質(zhì)指標(biāo)濃度, 導(dǎo)致其峰谷值及曲線高低變換存在一定差異。

2.2 水質(zhì)參數(shù)與單波段反射率的相關(guān)性

圖3 養(yǎng)殖池塘(a)和尾水處理池(b)采樣點(diǎn)光譜反射率

圖4 養(yǎng)殖小區(qū)池塘和尾水處理池水質(zhì)指標(biāo)與單波段光譜反射率的相關(guān)系數(shù)

2.3 水質(zhì)參數(shù)與雙波段反射率的相關(guān)性

表1 各水質(zhì)指標(biāo)與單波段及雙波段反射率指數(shù)最高相關(guān)系數(shù)

*注: 相關(guān)系數(shù)(對(duì)應(yīng)波段)

圖5 各水質(zhì)參數(shù)與最佳雙波段反射率指數(shù)相關(guān)系數(shù)分布圖

2.4 最佳水質(zhì)反演模型篩選與檢驗(yàn)

表2 各水質(zhì)參數(shù)反演模型評(píng)價(jià)

圖6 各水質(zhì)參數(shù)實(shí)測(cè)值與最佳反演模型預(yù)測(cè)值對(duì)比圖

2.5 水質(zhì)參數(shù)的高光譜圖像反演

圖7 養(yǎng)殖小區(qū)池塘各水質(zhì)指標(biāo)濃度空間反演效果

3 討論

無人機(jī)高光譜技術(shù)在定量反演水質(zhì)指標(biāo)時(shí), 由于具有較高的光譜分辨率和較強(qiáng)的波段連續(xù)性等特點(diǎn), 可以獲得更為全面的光譜特征信息, 但同時(shí)由于光譜測(cè)量中的某些人為和自然因素的干擾, 光譜數(shù)據(jù)需要進(jìn)行不同的數(shù)學(xué)變換以減少系統(tǒng)誤差和背景噪聲的影響, 同時(shí)降低信息重疊度和冗余度, 進(jìn)一步增強(qiáng)信噪比, 提高光譜數(shù)據(jù)與水質(zhì)參數(shù)的相關(guān)性, 有利于提高線性模型的預(yù)測(cè)精度(趙慶展等, 2016)。本研究中3種數(shù)值變換方法明顯提高了TSS和CODMn與光譜數(shù)據(jù)的相關(guān)性, 均達(dá)到0.81以上, 差值指數(shù)對(duì)TN與光譜數(shù)據(jù)的相關(guān)性提高最大, 比值指數(shù)則對(duì)TP與光譜數(shù)據(jù)的相關(guān)性提高最大, 而氨氮是例外, 數(shù)值變換方法甚至降低了其與光譜數(shù)據(jù)的相關(guān)性, 這可能是因?yàn)樵谠鰪?qiáng)敏感波段信息的同時(shí)也增強(qiáng)了噪聲信息, 給整體數(shù)據(jù)帶來了較大的反演誤差, 降低了模型預(yù)測(cè)精度(譚玉敏等, 2014)。結(jié)果說明不同波段數(shù)值變換方法的處理效果不同, 應(yīng)合理選擇處理方法以進(jìn)行后續(xù)建模, 這與相關(guān)研究結(jié)論也較為一致(董廣香, 2007; 馬娜等, 2010)。

4 結(jié)論

馬娜, 胡云鋒, 莊大方, 等, 2010. 基于最佳波段指數(shù)和J-M距離可分性的高光譜數(shù)據(jù)最佳波段組合選取研究——以環(huán)境小衛(wèi)星高光譜數(shù)據(jù)在東莞市的應(yīng)用為例[J]. 遙感技術(shù)與應(yīng)用, 25(3): 358-365.

葉勤, 姜雪芹, 李西燦, 等, 2017. 基于高光譜數(shù)據(jù)的土壤有機(jī)質(zhì)含量反演模型比較[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 48(3): 164-172.

史舟, 梁宗正, 楊媛媛, 等, 2015. 農(nóng)業(yè)遙感研究現(xiàn)狀與展望[J].農(nóng)業(yè)機(jī)械學(xué)報(bào), 46(2): 247-260.

呂航, 馬蔚純, 周立國, 等, 2013. 淀山湖微量水質(zhì)參數(shù)衛(wèi)星高光譜遙感估算[J]. 復(fù)旦學(xué)報(bào)(自然科學(xué)版), 52(2): 238-246.

劉彥君, 夏凱, 馮海林, 等, 2019. 基于無人機(jī)多光譜影像的小微水域水質(zhì)要素反演[J]. 環(huán)境科學(xué)學(xué)報(bào), 39(4): 1241-1249.

劉梅, 原居林, 倪蒙, 等, 2021. “三池兩壩”多級(jí)組合工藝對(duì)內(nèi)陸池塘養(yǎng)殖尾水的處理[J]. 環(huán)境工程技術(shù)學(xué)報(bào), 11(1): 97-106.

李俊生, 2007. 高光譜遙感反演內(nèi)陸水質(zhì)參數(shù)分析方法研究——以太湖為例[D]. 北京: 中國科學(xué)院遙感應(yīng)用研究所: 6-10.

李道亮, 劉暢, 2020. 人工智能在水產(chǎn)養(yǎng)殖中研究應(yīng)用分析與未來展望[J]. 智慧農(nóng)業(yè)(中英文), 2(3): 1-20.

楊振, 盧小平, 武永斌, 等, 2020. 無人機(jī)高光譜遙感的水質(zhì)參數(shù)反演與模型構(gòu)建[J]. 測(cè)繪科學(xué), 45(9): 60-64, 95.

張兵, 申茜, 李俊生, 等, 2009. 太湖水體3種典型水質(zhì)參數(shù)的高光譜遙感反演[J]. 湖泊科學(xué), 21(2): 182-192.

張武肖, 孫盛明, 戈賢平, 等, 2015. 急性氨氮脅迫及毒后恢復(fù)對(duì)團(tuán)頭魴幼魚鰓、肝和腎組織結(jié)構(gòu)的影響[J]. 水產(chǎn)學(xué)報(bào), 39(2): 233-244.

陳俊英, 邢正, 張智韜, 等, 2019. 基于高光譜定量反演模型的污水綜合水質(zhì)評(píng)價(jià)[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 50(11): 200-209.

趙慶展, 劉偉, 尹小君, 等, 2016. 基于無人機(jī)多光譜影像特征的最佳波段組合研究[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 47(3): 242-248.

胡震天, 周源, 2020. 基于低空多光譜遙感的城市水質(zhì)監(jiān)測(cè)方法研究[J]. 地理空間信息, 18(7): 4-8.

段洪濤, 張柏, 宋開山, 等, 2006. 長春南湖水體透明度高光譜定量模型研究[J]. 中國科學(xué)院研究生院學(xué)報(bào), 23(5): 633-639.

徐福祥, 高志強(qiáng), 尚偉濤, 等, 2018. 基于無人機(jī)的海濱景區(qū)綠潮監(jiān)測(cè)[J]. 海洋與湖沼, 49(5): 1061-1067.

唐軍武, 田國良, 汪小勇, 等, 2004. 水體光譜測(cè)量與分析Ⅰ: 水面以上測(cè)量法[J]. 遙感學(xué)報(bào), 8(1): 37-44.

黃宇, 陳興海, 劉業(yè)林, 等, 2020. 基于無人機(jī)高光譜成像技術(shù)的河湖水質(zhì)參數(shù)反演[J]. 人民長江, 51(3): 205-212.

黃宇業(yè), 付東洋, 劉大召, 等, 2019. 珠江口水體表觀光譜特性與類型分析[J]. 海洋環(huán)境科學(xué), 38(6): 891-897.

黃昕晰, 應(yīng)晗婷, 夏凱, 等, 2020. 基于無人機(jī)多光譜影像和OPT-MPP算法的水質(zhì)參數(shù)反演[J]. 環(huán)境科學(xué), 41(8): 3591-3600.

董廣香, 2007. 三峽庫區(qū)高光譜水質(zhì)參數(shù)反演數(shù)據(jù)處理及分析研究[D]. 北京: 中國測(cè)繪科學(xué)研究院: 10-17.

譚玉敏, 夏瑋, 2014. 基于最佳波段組合的高光譜遙感影像分類[J]. 測(cè)繪與空間地理信息, 37(4): 19-22.

BANSOD B, SINGH R, THAKUR R, 2018. Analysis of water quality parameters by hyperspectral imaging in Ganges River [J]. Spatial Information Research, 26(2): 203-211.

GITELSON A, GARBUZOV G, SZILAGYI F,, 1993Quantitative remote sensing methods for real-time monitoring of inland waters quality [J]. International Journal of Remote Sensing, 14(7): 1269-1295.

GUIMARAES T T, VERONEZ M R, KOSTE E C,, 2019. Evaluation of regression analysis and neural networks to predict total suspended solids in water bodies from unmanned aerial vehicle images [J]. Sustainability, 11(9): 2580.

KOPONEN S, PULLIAINEN J, KALLIO K,, 2002. Lake water quality classification with airborne hyperspectral spectrometer and simulated MERIS data [J]. Remote Sensing of Environment, 79(1): 51-59.

ROMERO-TRIGUEROS C, NORTES P A, ALARCóN J J,, 2017. Effects of saline reclaimed waters and deficit irrigation onphysiology assessed by UAV remote sensing [J]. Agricultural Water Management, 183: 60-69.

YAO X, HUANG Y, SHANG G Y,, 2015. Evaluation of six algorithms to monitor wheat leaf nitrogen concentration [J]. Remote Sensing, 7(11): 14939-14966.

WATER QUALITY MONITORING OF INLAND AQUACULTURE PONDS BASED ON UAV HYPERSPECTRAL REMOTE SENSING TECHNOLOGY

LIU Mei1, MA Qi-Liang2, YUAN Ju-Lin1, NI Meng1, LIAN Qing-Ping1, GUO Ai-Huan1

(1. Agriculture Ministry Key Laboratory of Healthy Freshwater Aquaculture, Key Laboratory of Fish Health and Nutrition of Zhejiang Province, Zhejiang Institute of Freshwater Fisheries, Huzhou 313001, China; 2. Center of Information Technology, Huzhou University, Huzhou 313000, China)

UAV; hyperspectral technology; inland ponds; machine algorithm; water quality inversion

S951.2; S966; X87

10.11693/hyhz20210600133

*浙江省公益技術(shù)研究計(jì)劃項(xiàng)目, LGN20C190004號(hào); 國家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目, 2019YFD0900302號(hào)。劉 梅, 博士, 助理研究員, E-mail: liumei@zju.edu.cn

原居林, 高級(jí)工程師, E-mail: yuanjulin1982@126.com

2021-06-08,

2021-07-27