水體光譜測量與數據處理

白 璟，彭 寬

（1.武漢市測繪研究院，湖北 武漢 430000；2.武漢市勘察設計有限公司，湖北 武漢 430000）

1 概述

1.1 水體光譜特性研究的發展和現狀

近年來，遙感科學與技術有了突飛猛進的發展，上到衛星影像，下到航空攝影測量，遙感的發展正是對光譜數據利用能力的提高，水是生命之源，水體在自然界中重要地位也是毋庸置疑的，所以對水體的研究應該從多方面多層次考慮，以前我們僅僅利用遙感科學技術來發現水域，隨著技術的不斷革新，現在我們已經可以利用遙感技術來識別水域并判斷水質狀況。近年來，隨著對物質光譜特征研究的深入、算法的改進，遙感監測水質已從定性發展到定量，可通過遙感數據識別水質情況，并對水體的各種參數進行預測，對水體質量的判斷有了明顯的提高。

1.2 水體的反射光譜特性

水體的反射主要在藍綠光波段，其他波段吸收率很強。

2 光譜數據采集與預處理

2.1 實驗區概況

武漢市地處長江中游，位于長江、漢水的交匯處，素有“江城”之稱。市區江河縱橫，湖泊密布，僅長江、漢江、府環河三條水系縱貫市區就長達180.5 公里。武漢大多數湖泊呈富營養化狀態，一些湖泊已接近甚至達到極富營養化。為水體光譜研究提供了豐富的資源。

2.2 研究內容

2.2.1 光譜數據的采集

光譜數據主要來源于可見光和近紅外波段的地面光譜儀測量獲得的。在光譜的采集過程中，影響數據質量的因素很多，有儀器因素，人為因素，環境因素，這些相互影響較為復雜，不但會產生誤差甚至會出現錯誤數據，有些誤差可以通過數據處理的方法解決，有些則無法彌補，所以光譜的采集至關重要，一定要按儀器的操作規程和對環境的要求來進行。

2.2.2 光譜數據的預處理

利用光譜儀測量得到的數據由于各種原因的相互影響客觀上會存在較大的噪聲，這樣的數據在實際研究時意義不大，為地物反演研究的需要，需要對數據進行預處理。

2.2.3 光譜數據的分析

對采集來的水體光譜數據進行分析，找出水體光譜特征點，以及光譜特征與哪些因素有關。

2.3 研究方案

光譜數據主要來源于2020年5月間通過可見光和近紅外波段的地面光譜儀測量獲得的。

選取自然界中最具有代表性的水體作為光譜采樣目標，實驗中我們主要測量了有代表性的水域東湖，東湖水體富營養化程度比較高，在做水質反演方面有比較實際的作用。

3 水體光譜數據處理

我們利用光譜儀測得的光譜數據由于每個點上環境因素的影響、人為操作的不同、測量次數的差異，造成了數據的數量和質量都有差別，因此需要對原始數據進行處理，才能為后續的特征分析所使用。目前數據處理主要可以用求平均值、波段合并，歸一化和光譜微分等方法來處理，以提高數據的信噪比，本文主要介紹求平均值和光譜微分的處理方法。

3.1 光譜數據求平均值

通過大量的光譜測量得到較多的光譜曲線，有些數據能很好的反映地物的基本特征但是有些數據在某些波段效果比較差，需要選擇吻合較好的曲線進行數據處理，在選擇吻合較好的光譜曲線的方法上可以選擇最小二乘法，他能比較有效的剔除誤差太大的光譜數據。最小二乘法比觀察法的精度要高一些但實現起來難度要大一些。

以下是2020年5月測得的東湖水體光譜數據：

圖1 吻合較好的7 組光譜數據

圖2 7 組數據平均后得到平均光譜

由以上效果圖可以看出，對同一點測量得到的光譜取平均值后能很好的反映該點的光譜特征的同時還可以減弱由于個別測量數據不準引起的較大誤差。光譜測量中，在條件允許的情況下建議對每一點進行多次測量。

3.2 光譜曲線一階微分

對光譜的一階微分處理可以去除部分線性或接近線性的背景、噪聲光譜對目標光譜的影響，R（λi）′光譜數據的一階微分可以用以下公式近似計算：

其中，λi-1，λi，λi+1為相鄰波長，R（λi）′為波長λi的一階微分反射光譜。計算得到的一階微分光譜如圖3 所示。一階光譜為反演水質參數，植被參數，土壤含水量等工作提供了新的數據。

圖3 水體光譜反射率一階微分圖

4 水體光譜特征位置選取與影響因素分析

純凈的水體在可見光波段的反射率曲線是接近線性的，且隨著波長向紅外波段的逐漸增大，反射率呈現逐漸減小的趨勢，直線化特征明顯。自然狀態下的水體則會由于其中含有各種物質對光輻射的吸收和散射性質的不同而顯示出不同的光譜特性。在本次光譜測量測得的水體中，決定其光譜特征的物質主要有：

（1）浮游植物（主要是藻類）和由浮游植物死亡而產生的有機碎屑。

（2）水體底泥經再懸浮而產生的無機懸浮顆粒，簡稱懸浮物。

（3）由黃腐酸、腐殖酸組成的有機物，通常稱黃色物質。

可以用以下的簡單模型近似描述各種物質對水體光譜反射率RW的影響：

其中 RW為水面反射率，bW、bS和 bP分別為水、無機懸浮物和藻類物質的后向散射系數，aW、aS、aP和 aY分別為水、無機懸浮物、藻類物質和黃色物質的吸收系數。

本次測得的水體的光譜特征主要以浮游植物引起的富養化特征（如東湖）為主。為了進一步研究這些特征參數，通過對水體光譜曲線的分析，得出了東湖水體光譜的幾個特征點。

由于各種物質的影響，水體光譜曲線具有特征明顯且穩定的特征峰谷。我們在400-800nm 波段范圍選擇了4 個峰谷表征水體的光譜特征。

圖4 典型水體光譜（2006年5月東湖）

（1）藍紫波段吸收峰B：為光譜在藍紫波段的反射率最小值，λ∈420-525nm。

（2）綠反射峰（綠峰位置）G：為光譜在綠波段的反射率最大值，λ∈520-560nm。

（3）紅波段吸收峰R：為光譜在紅波段的反射率最小值，λ∈660-690nm。

（4）紅波段反射峰N：為光譜在紅波段向近紅外波段過度區間的波段的反射率最大值，λ∈690-720nm。

以上是水體的4 個最基本的特征位置，本次實驗測得的數據都能很好的反映水體的以上波譜特征。

5 研究進展

5.1 研究進展

通過對水體光譜數據的分析，提出了水體的4 個相關特征位置。并對其產生的原因作了定性的分析。由于受條件限制，對地物采樣次數和采樣點數量較少，進一步的工作需要加大采樣點的數量和范圍，對多個季節的數據進行分析比較。

5.2 研究展望

（1）改進和加強對地物光譜數學分析方法。以數學分析方法為基礎，結合物理方法，使對地物特征的定性分析向定量分析發展。

（2）對水體光譜的研究可以從可見光波段范圍向近紅外和熱紅外波段范圍發展。