陳鵬飛 徐新剛

1 中國科學(xué)院地理科學(xué)與資源研究所 / 資源與環(huán)境信息系統(tǒng)國家重點(diǎn)實(shí)驗(yàn)室, 北京100101; 2 江蘇省地理信息資源開發(fā)與利用協(xié)同創(chuàng)新中心, 江蘇南京 210023; 3 國家農(nóng)業(yè)信息化工程技術(shù)研究中心, 北京100097

遙感技術(shù)是反演作物生理、生化參數(shù)的重要工具[1-3], 其在作物長勢監(jiān)測、產(chǎn)量估測中有廣泛應(yīng)用。到目前為止, 人們最常用到的對(duì)地觀測平臺(tái)有衛(wèi)星、有人機(jī)和觀測塔。與其他平臺(tái)相比, 低空無人機(jī)遙感具有成本低、能根據(jù)天氣隨時(shí)起降等優(yōu)勢。它提供了一種新的、多樣化的對(duì)地觀測數(shù)據(jù)獲取手段。當(dāng)其攜帶高質(zhì)量、低成本傳感器開展對(duì)地觀測時(shí), 可實(shí)現(xiàn)對(duì)目標(biāo)區(qū)域高頻次、便捷快速監(jiān)測[4-6]。因此, 無人機(jī)遙感對(duì)及時(shí)獲取作物長勢狀況, 支撐精準(zhǔn)農(nóng)業(yè)管理具有重要意義。

目前, 基于無人機(jī)遙感已開展許多和農(nóng)業(yè)相關(guān)的研究工作[7-9]。但是在無人機(jī)的應(yīng)用過程中, 仍有許多問題需要解決, 其中之一就是無人機(jī)影像正射拼接問題[10-12]。由于無人機(jī)影像飛行高度比較低,它獲取的一景影像的覆蓋范圍很小[13]。要想獲取整個(gè)研究區(qū)的單幅影像, 需要拼接成千上百張的無人機(jī)影像。將海量的無人機(jī)數(shù)據(jù)拼接, 以代表地球表面的狀態(tài)信息是一件非常有挑戰(zhàn)性的工作。由于非常耗時(shí), 傳統(tǒng)航空測量中使用的基于自動(dòng)空中三角網(wǎng)和光束法平差的方法, 很難在無人機(jī)影像拼接中使用[14]。近年來, 研究者提出基于現(xiàn)代計(jì)算機(jī)機(jī)器視覺和運(yùn)動(dòng)結(jié)構(gòu)算法可以很好解決這一問題[15-17]。目前, 已有一系列軟件基于這一思想設(shè)計(jì), 比如Photoscan 軟件、Pix4D 軟件等。然而, 由于缺乏對(duì)這些新設(shè)計(jì)軟件系統(tǒng)比較和驗(yàn)證, 用戶在選擇時(shí)存在一定的盲目性。

幾何精度和光譜精度是衡量拼接效果的兩個(gè)重要指標(biāo)。對(duì)于幾何精度, 目前已有一些對(duì)比不同拼接軟件性能的研究[18]。Gross 和Heumann[14]比較了Pix4D 軟件、Photoscan 軟件和ICE (Microsoft Image Composite Editor)軟件拼接影像的幾何精度, 發(fā)現(xiàn)Pix4D 軟件的精度最高, 然后是 Photoscan 軟件;Turner 等[19]比較了Pix4D 軟件、Photoscan 軟件和一款同樣基于機(jī)器視覺和運(yùn)動(dòng)結(jié)構(gòu)算法開發(fā)的非商用軟件, 結(jié)果表明Photoscan 軟件的精度最高, 然后是Pix4D 軟件?？傮w而言, 已有研究表明在幾何精度方面, 所有比較的商用和開源軟件中 Pix4D 和Photoscan 的精度最高。對(duì)于光譜精度, 還較少有人詳細(xì)比較不同軟件的性能。僅見Gross 和Heumann[14]目視比較了Photoscan, Pix4D 和ICE 拼接影像的清晰度。不同拼接軟件在拼接過程中對(duì)無人機(jī)原始影像的影響如何？從光譜角度看, 哪個(gè)是最好的拼接軟件？這些問題對(duì)于無人機(jī)農(nóng)業(yè)應(yīng)用來說是非常重要的。

考慮到Photoscan 和Pix4D 這兩款拼接軟件是目前幾何精度最高和最常被人們使用的軟件, 本研究通過設(shè)計(jì)實(shí)驗(yàn), 對(duì)比分析了它們在拼接過程中對(duì)原始光譜的影響。

1 材料與方法

1.1 田間試驗(yàn)

本研究數(shù)據(jù)來自北京市昌平區(qū)小湯山精準(zhǔn)農(nóng)業(yè)示范基地開展的冬小麥氮肥梯度試驗(yàn)(116°26′35″E,40°10′46″N)。研究區(qū)處于半濕潤大陸性季風(fēng)氣候區(qū),陽光充足, 溫度溫和, 降雨適中。冬小麥和夏玉米輪作是該地區(qū)典型的耕作模式。田間試驗(yàn)面積4800 m2。選用2 個(gè)小麥品種, 4 個(gè)氮肥處理和4 個(gè)重復(fù)。一共32 個(gè)田間試驗(yàn)小區(qū)。小區(qū)大小9 m × 15 m。小麥品種是Lunxuan 167 和Jingdong 18。氮肥處理分別為0、90、180 和270 kg hm-2。除氮肥外, 其他田間管理措施相同。

1.2 數(shù)據(jù)獲取

1.2.1 無人機(jī)數(shù)據(jù)獲取 在冬小麥生長的挺身期進(jìn)行無人機(jī)飛行。飛行時(shí), 無人機(jī)搭載RedEdge-M(MicaSense, 美國華盛頓州)相機(jī)來獲取多光譜影像。所用無人機(jī)為3DR Solo (3DR, 美國加利福尼亞州)四旋翼無人機(jī)。RedEdge-M 相機(jī)包含5 個(gè)波段,中心波長分別位于475、560、668、717 和840 nm, 波段寬度為20、20、10、10 和40 nm。鏡頭焦距為5.5 mm,水平視場角為47.2°, 垂直視場角為35.4°, 采用全局快門模式成像, 分辨率為1280×960 像素。無人機(jī)在中午12：00 左右, 晴空無云時(shí)飛行, 飛行高度為30 m,對(duì)應(yīng)影像空間分辨率約為1.99 cm。無人機(jī)飛行時(shí),所拍影像的旁向和航向重疊率均設(shè)定為75%。無人機(jī)飛行前, 采用手動(dòng)模式拍攝白板。白板信息用來后期將影像DN 值轉(zhuǎn)換為反射率。本次飛行共獲得329 景影像, 其中無人機(jī)飛行平臺(tái)及一景獲取的影像如圖1 所示。

圖1 無人機(jī)飛行平臺(tái)及獲取的一景影像Fig. 1 UAV platform and an acquired image

1.2.2 無人機(jī)數(shù)據(jù)預(yù)處理 Pix4D 軟件(Pix4D,洛桑, 瑞士)和Photoscan 軟件(Agisoft, 圣彼得堡,俄羅斯)被分別用來拼接獲取的無人機(jī)影像。本研究使用了Pix4D AG 的4.2.25 版和Photoscan 的1.4.2版。此外, 為了檢驗(yàn)軟件在拼接處理中對(duì)原始影像光譜信息的保留程度, 基于相機(jī)制造商提供的開源程序, 本研究制作了每次相機(jī)拍攝獲取5 個(gè)波段影像所合成的多光譜影像。相關(guān)Python 程序可以從網(wǎng)站https：//github.com/micasense/imageprocessing 下載,包括各波段影像的畸變校正、反射率計(jì)算和不同波段影像的對(duì)齊處理。基于以上處理, 共有3 種類型的數(shù)據(jù)源, 包含2 個(gè)用不同拼接軟件拼接的影像和沒有拼接的單張多光譜影像。為了保證它們有相同的空間基準(zhǔn), 用Trimble Geo 7X (Trimble, 加利福尼亞州, 美國)差分GPS 采集厘米級(jí)誤差GPS 點(diǎn)數(shù)據(jù),來對(duì)Pix4D 拼接的影像和Photoscan 拼接的影像進(jìn)行幾何校正。另外, 在每個(gè)小區(qū)隨機(jī)布設(shè)樣點(diǎn), 選擇對(duì)應(yīng)樣點(diǎn)天頂?shù)膯螐垷o人機(jī)影像, 參照校正好的拼接影像進(jìn)行幾何校正。其中, 各小區(qū)布設(shè)的樣點(diǎn)信息如圖2 所示。

圖2 樣點(diǎn)空間分布圖Fig. 2 Distribution map for sampling points

1.3 數(shù)據(jù)分析

首先, 提取每個(gè)樣點(diǎn)在Pix4D 拼接影像、Photoscan拼接影像和單張影像上的光譜信息。提取時(shí), 在每個(gè)采樣點(diǎn)提取半徑為0.75 m 圓形區(qū)域內(nèi)所有像元的光譜值。其次, 對(duì)于每一來源的影像數(shù)據(jù), 以各樣點(diǎn)提取的所有像元光譜值的平均值代表樣點(diǎn)光譜值,用以分析不同來源影像光譜反射率之間的相關(guān)性與差異。其中, 差異分析采用了散點(diǎn)圖分析和成對(duì)數(shù)據(jù)t檢驗(yàn)。最后, 對(duì)于每個(gè)來源的影像數(shù)據(jù), 基于提取的各樣點(diǎn)處所有像元光譜反射率計(jì)算變異系數(shù)(CV), 對(duì)比不同來源影像變異系數(shù)的差異。變異系數(shù)反映了影像的對(duì)比度。同樣, 用散點(diǎn)圖分析和成對(duì)數(shù)據(jù)t檢驗(yàn)來分析它們之間的差異。流程圖如圖3所示。

2 結(jié)果與分析

2.1 不同數(shù)據(jù)源光譜反射率的比較

從圖4 可以看出, 用P i x 4 D 拼接的影像和Photoscan拼接的影像均沒有出現(xiàn)模糊和偽影區(qū)域?？紤]到相關(guān)性分析結(jié)果, 樣點(diǎn)處兩種拼接影像反射率與單張影像反射率間均存在較好的相關(guān)性, 相關(guān)系數(shù)0.68以上, 達(dá)到0.001顯著水平(圖5)。但從它們的散點(diǎn)圖看, 代表Pix4D拼接影像的樣點(diǎn)基本都在1∶1線以下, 說明Pix4D拼接影像的反射率在各個(gè)波段均小于單張影像反射率, 所用方法傾向于將原始影像各波段反射率變小。在藍(lán)光波段、綠光波段、紅光波段, 代表Photoscan拼接影像的樣點(diǎn)基本都在1∶1線周圍, 而在紅邊波段、近紅外波段, 代表Photoscan拼接影像的樣點(diǎn)大部分在1∶1線的上方,說明Photoscan拼接軟件所用方法在反射率值較低波段, 能較好地保留原始影像光譜信息, 但在反射率值較高波段, 會(huì)增大波段反射率值。總體說來, 與Pix4D相比,Photoscan能保留更多的原始影像光譜反射率信息。

圖3 數(shù)據(jù)分析流程圖Fig. 3 Flow chart for data analysis

圖4 Pix4D 拼接影像(a)、Photoscan 拼接影像(b)、一些單張影像(c)的假彩色合成影像和小區(qū)施氮量空間分布圖(d)Fig. 4 False color images for Pix4D mosaicked image (a), Photoscan mosaicked image (b), some single images (c), and map of applied nitrogen in each zone (d)

圖5 單張影像與拼接影像在各波段的反射率散點(diǎn)圖Fig. 5 Reflectance scatter plots of single image and mosaicked image in each band

從成對(duì)數(shù)據(jù)t檢驗(yàn)的結(jié)果看, Pix4D拼接影像的反射率與單張影像反射率在各波段均存在顯著性差異, 達(dá)0.001顯著水平(表1)。對(duì)于Photoscan拼接影像來說其在紅邊波段和近紅外波段的反射率與單張影像對(duì)應(yīng)波段的反射率存在顯著性差異(P＜0.001)。但在藍(lán)光波段、綠光波段和紅光波段, Photoscan拼接影像的波段反射率與單張影像的波段反射率不存在顯著性差異(P＜0.001)。以上結(jié)果表明, 在藍(lán)、綠、紅波段, Photoscan軟件在拼接影像時(shí)能比Pix4D軟件在拼接影像時(shí)保留更多的原始光譜反射率信息。除此之外, 當(dāng)比較Pix4D拼接影像和Photoscan拼接影像之間光譜反射率的差異時(shí), 它們分布于散點(diǎn)圖的不同區(qū)域, 且各個(gè)波段均表現(xiàn)出顯著性差異(P＜0.001)。說明Pix4D和Photoscan在拼接單張影像時(shí)采用了不同的數(shù)據(jù)處理策略。

2.2 不同數(shù)據(jù)源影像對(duì)比度的比較

表1 Pix4D 拼接影像、Photoscan 拼接影像和單張影像之間反射率成對(duì)數(shù)據(jù)t 檢驗(yàn)結(jié)果Table 1 Paired t-test results of reflectance data among the Pix4D mosaicked image, the Photoscan mosaicked image, and the single image

Pix4D拼接影像、Photoscan拼接影像和單張影像樣點(diǎn)處光譜反射率平均變異系數(shù)和散點(diǎn)圖, 如表2、圖6 所示。與P i x 4 D 拼接影像的變異系數(shù)相比,Photoscan 拼接影像每個(gè)波段的變異系數(shù)與單張影像相應(yīng)波段的變異系數(shù)更加接近(表2)。從散點(diǎn)圖上看, 對(duì)于Pix4D 拼接影像, 其在藍(lán)光波段和紅光波段的變異系數(shù)與對(duì)應(yīng)單張影像相應(yīng)波段的變異數(shù)據(jù)有較好的對(duì)應(yīng)關(guān)系, 樣點(diǎn)均分布在1∶1 線周圍, 但其在綠光波段、紅邊波段和近紅外波段的變異系數(shù)明顯低于單張影像相應(yīng)波段的變異系數(shù), 絕大部分樣點(diǎn)分布在1∶1 線的下部; 對(duì)于Photoscan 拼接影像,除在紅邊波段, 樣點(diǎn)略偏向1∶1 線的下部, 其在其他各個(gè)波段的變異系數(shù)與對(duì)應(yīng)單張影像相應(yīng)波段的變異系數(shù)均有較好的對(duì)應(yīng)關(guān)系, 樣點(diǎn)均分布在1∶1 線周圍。從成對(duì)數(shù)據(jù)t檢驗(yàn)結(jié)果看, 除紅光波段和藍(lán)光波段外,Pix4D 拼接影像在其他波段的變異系數(shù)與單張影像顯著不同, 顯著性大于0.001 水平(表3); Photoscan 拼接影像除在紅邊波段與單張影像對(duì)應(yīng)波段的變異系數(shù)存在顯著性差異外, 在其他波段的變異系數(shù)均與單張影像不存在顯著性差異(P＜0.001)。以上結(jié)果說明, 與Pix4D拼接影像相比, Photoscan 拼接影像能保留更多的單張影像的對(duì)比度信息。此外, 除藍(lán)光波段外, Photoscan 拼接影像各波段的變異系數(shù)與Pix4D 拼接影像相應(yīng)波段的變異系數(shù)在散點(diǎn)圖中處于不同位置, 它們之間存在顯著性差異(P＜ 0.001)。

3 討論

表2 Pix4D 拼接影像、Photoscan 拼接影像和單張影像樣點(diǎn)處光譜反射率平均變異系數(shù)(CV)Table 2 Averaged coefficient of variation (CV) data for Pix4D, Photoscan, and the single image

圖6 單張影像與拼接影像在各波段的反射率變異系數(shù)散點(diǎn)圖Fig. 6 Reflectance coefficient of variation (CV) scatter plots of single image and mosaicked image in each band

Gross 和Heumann[14]通過目視觀測發(fā)現(xiàn)Pix4D拼接影像和Photoscan 拼接影像均有模糊和偽影區(qū)域, 從側(cè)面說明兩個(gè)軟件在拼接無人機(jī)數(shù)據(jù)時(shí)均改變了原始影像光譜信息。如前所述, 本研究中用Pix4D 拼接的影像和Photoscan 拼接的影像均沒有出現(xiàn)模糊和偽影區(qū)域。但是, 當(dāng)進(jìn)行定量比較時(shí),Pix4D 拼接影像、Photoscan 拼接影像和單張影像間在光譜反射率和其變異系數(shù)上均存在差異。這與Gross 和Heumann 的結(jié)果一致。此外, 3 種不同來源數(shù)據(jù)的光譜反射率和其變異系數(shù)的比較表明, Photoscan 在拼接無人機(jī)影像過程中, 比Pix4D 保留了更多的原始影像信息。

表3 Pix4D 拼接影像、Photoscan 拼接影像和單張影像之間反射率變異系數(shù)的成對(duì)數(shù)據(jù)t 檢驗(yàn)結(jié)果Table 3 Paired t test results of reflectance coefficients of variation (CV) among the Pix4D mosaicked image, the Photoscan mosaicked image, and the single image

幾何精度和光譜精度是兩個(gè)衡量拼接軟件性能的關(guān)鍵技術(shù)指標(biāo)。在幾何精度方面, Gross 和Heumann 的研究表明與其他軟件相比, Pix4D 軟件拼接影像的幾何精度最高[14], 但是Sona 等[18]和Turner等[19]的研究則表明Photoscan 軟件拼接影像的幾何精度要好于Pix4D。顯然, 這兩款拼接軟件拼接影像的幾何精度取決于具體的應(yīng)用場景。不過, 基于這些研究的結(jié)果看, Pix4D 和Photoscan 拼接影像的幾何精度差異并不大。此外, 即使拼接影像存在細(xì)微的幾何校正誤差, 后期還可以用野外采集的高精度的GPS 點(diǎn)對(duì)它們進(jìn)行幾何精校正。本研究僅僅比較了兩款拼接軟件拼接影像光譜信息的差異?？紤]到光譜信息, 盡管兩款拼接軟件拼接影像的光譜反射率都與單張影像光譜反射率有很高的相關(guān)性, 但與Pix4D 軟件拼接的影像相比, Photoscan 軟件拼接影像的反射率與單張未拼接影像的反射率具有更好的對(duì)應(yīng)關(guān)系, 且反射率的變異系數(shù)與原始影像也更加接近。說明Photoscan 在拼接影像時(shí)保留了更多的原始影像信息, 特別是對(duì)比度信息。好的對(duì)比度, 會(huì)使基于Photoscan 軟件拼接的影像在復(fù)雜區(qū)域進(jìn)行目標(biāo)分類時(shí)比基于Pix4D 軟件拼接的影像表現(xiàn)更好。除了幾何精度和光譜精度這兩個(gè)標(biāo)準(zhǔn), 價(jià)格也是用戶在使用拼接軟件時(shí)應(yīng)該考慮的因素。Photoscan 專業(yè)版的售價(jià)為3499 美元, Pix4D Mapper 專業(yè)版的售價(jià)為8700 美元。Pix4D AG 為農(nóng)業(yè)版, 它的售價(jià)為5185 美元。顯然, Photoscan 更加便宜。

本研究只是比較了兩款軟件在平原區(qū)農(nóng)業(yè)地塊應(yīng)用的情況。今后, 需要在不同地形條件, 對(duì)結(jié)果進(jìn)行更廣泛的驗(yàn)證。目前來說, 本研究提供了一個(gè)比較軟件拼接影像光譜精度的實(shí)例, 其結(jié)論適合于類似區(qū)域。

4 結(jié)論

Pix4D 軟件和Photoscan 軟件在拼接無人機(jī)影像時(shí), 盡管拼接影像的光譜反射率均與單張影像光譜反射率有較高的相關(guān)性, 但在拼接處理影像的過程中均改變了原始單張影像的反射率。其中, Pix4D 軟件會(huì)導(dǎo)致各波段反射率變低, 而Photoscan 軟件能較好地處理具有低反射率值的光譜波段, 但在具有高反射率值的光譜波段會(huì)導(dǎo)致處理后的反射率變大。與Pix4D 相比, 基于Photoscan 拼接的影像, 其反射率及影像對(duì)比度與原始單張影像更加接近。對(duì)于與本研究類似區(qū)域的應(yīng)用來說, 推薦Photoscan 為性價(jià)比較好的軟件。

致謝 感謝孫剛、劉金然、徐志濤在田間實(shí)驗(yàn)過程中給予的支持與幫助。