權(quán)赫春， 房金鵬， 許 鎮(zhèn)， 林哲浩

（1.延邊大學(xué)地理與海洋科學(xué)學(xué)院，吉林延吉133002；2.黃河河務(wù)局，濟(jì)南250011）

0 引 言

地理信息系統(tǒng)（Geographical Information System，GIS）學(xué)科是一門實(shí)踐性很強(qiáng)的高新技術(shù)學(xué)科，因此，理論教學(xué)與實(shí)驗(yàn)教學(xué)的共同發(fā)展是GIS學(xué)科發(fā)展和人才培養(yǎng)的必由途徑［1-2］；是集計(jì)算機(jī)科學(xué)、地理科學(xué)、測繪學(xué)、遙感學(xué)、環(huán)境科學(xué)等多學(xué)科為一體的新興邊緣學(xué)科。它不但與全球定位系統(tǒng)（GPS）和遙感（RS）相結(jié)合，構(gòu)成3S集成系統(tǒng)，而且與CAD、多媒體、通信、因特網(wǎng)、辦公自動(dòng)化、虛擬現(xiàn)實(shí)等多種技術(shù)相結(jié)合，構(gòu)成了綜合的信息技術(shù)。地理信息科學(xué)專業(yè)跨自然、人文和技術(shù)三大領(lǐng)域，具有較強(qiáng)的綜合性和專業(yè)技術(shù)性。

隨著我國經(jīng)濟(jì)、社會(huì)和科學(xué)技術(shù)的發(fā)展，各企業(yè)對大學(xué)生的實(shí)踐和創(chuàng)新能力提出了新的要求。《國家中長期教育改革和發(fā)展規(guī)劃綱要（2010—2020年）》中指出：“優(yōu)化知識結(jié)構(gòu)，豐富社會(huì)實(shí)踐，強(qiáng)化能力培養(yǎng)。著力提高學(xué)生的學(xué)習(xí)能力、實(shí)踐能力、創(chuàng)新能力”?；趧?chuàng)新實(shí)驗(yàn)的創(chuàng)新實(shí)踐教學(xué)對于提高學(xué)生綜合素質(zhì)、培養(yǎng)學(xué)生的創(chuàng)新精神與實(shí)踐能力具有特殊的作用［3］。而高校對學(xué)生的培養(yǎng)歷來以課堂教育為主，對學(xué)生創(chuàng)新和實(shí)踐能力的重視不夠［4］。學(xué)生在校所學(xué)知識和創(chuàng)新實(shí)踐能力與企業(yè)、科研機(jī)構(gòu)等要求相差甚遠(yuǎn)，經(jīng)過一定的培訓(xùn)和再教育才能適應(yīng)企業(yè)要求。

為了有效培養(yǎng)本科生創(chuàng)新實(shí)踐能力，本文根據(jù)地理信息科學(xué)專業(yè)的特點(diǎn)，設(shè)計(jì)了基于無人機(jī)與3S技術(shù)的創(chuàng)新實(shí)驗(yàn)。試圖通過該項(xiàng)創(chuàng)新實(shí)驗(yàn)的操作與完成，培養(yǎng)地理信息科學(xué)專業(yè)學(xué)生的創(chuàng)新與實(shí)踐能力。

1 實(shí)驗(yàn)?zāi)康?/h2>

掌握攝影測量學(xué)理論知識、無人機(jī)的拍攝和專用軟件使用方法，了解相機(jī)校正、像點(diǎn)位移、控制點(diǎn)校正等處理及4D產(chǎn)品生成過程［5］；掌握GPS的使用方法、控制點(diǎn)和災(zāi)害點(diǎn)的定位方法；了解利用RS技術(shù)獲取災(zāi)害信息、監(jiān)督和無監(jiān)督分類方法及土地利用圖的制作過程，并掌握RS相關(guān)軟件使用方法；掌握GIS的緩沖、地形分析、重分類、疊加等分析功能，了解利用GIS技術(shù)制作災(zāi)害因子圖層和危險(xiǎn)性區(qū)劃圖的過程；提高學(xué)生創(chuàng)新實(shí)踐能力、動(dòng)手能力和解決問題能力。

2 實(shí)驗(yàn)方案

以吉林省延邊州為研究區(qū)域（見圖1），基于無人機(jī)與3S技術(shù)，利用BP神經(jīng)網(wǎng)絡(luò)、信息量法等數(shù)學(xué)模型，進(jìn)行泥石流災(zāi)害危險(xiǎn)性評價(jià)。該實(shí)驗(yàn)主要包括：①利用無人機(jī)制作DOM（數(shù)字正射影像）、DEM（數(shù)字高程模型）等4D產(chǎn)品；②選擇、量化及歸一化處理泥石流災(zāi)害因子；③基于模型的危險(xiǎn)性評價(jià)；④ 精度評價(jià)及繪制危險(xiǎn)性區(qū)劃圖等4個(gè)部分，實(shí)驗(yàn)流程如圖2所示。

圖1 研究區(qū)域及泥石流災(zāi)害點(diǎn)分布圖

圖2 實(shí)驗(yàn)流程圖

3 實(shí)驗(yàn)儀器和軟件

儀器：雷凌SL700高精度航測無人機(jī)（見圖3）、GPS、GIS數(shù)據(jù)采集儀、全站儀、計(jì)算機(jī)等。

圖3 雷凌SL700高精度航測無人機(jī)

軟件：Arcgis10.2（GIS數(shù)據(jù)處理軟件）、PGSTEAMER（遙感數(shù)據(jù)處理軟件）、Skyphoto-Super（無人機(jī)數(shù)據(jù)處理軟件）、SkyPhoto-Map3D（三維測圖軟件）等。

4 實(shí)驗(yàn)步驟

4.1 利用無人機(jī)制作DOM和DEM

首先對局部地區(qū)進(jìn)行地面像控點(diǎn)的布設(shè)，根據(jù)參考文獻(xiàn)及實(shí)踐經(jīng)驗(yàn)，一般需要布設(shè)5～12個(gè)控制點(diǎn)。并利用全站儀或GPS對像控點(diǎn)進(jìn)行聯(lián)測。布設(shè)像控點(diǎn)之后，利用高精度航測無人機(jī)對研究區(qū)域進(jìn)行拍攝。拍攝時(shí)主要考慮精度、拍攝比例尺、航攝分區(qū)的劃分、航線敷設(shè)、航攝時(shí)間等問題［6］。其次，利用Skyphoto-Super軟件對無人機(jī)相片進(jìn)行處理，制作DEM和DOM。處理過程主要包括：相機(jī)、像控點(diǎn)、創(chuàng)建航帶參數(shù)設(shè)置；控三匹配、評差解算、更新元數(shù)據(jù)；模型定向、影像自動(dòng)匹配、匹配編輯、生成及拼接DEM；控三編輯、點(diǎn)云加密、DOM生成、DOM糾正及編輯［7］。

4.2 基于RS和GIS制作泥石流災(zāi)害因子圖層

通過文獻(xiàn)考察和災(zāi)害分析，選擇高程、坡度、坡向、降水量、河流、地層巖性、土地利用類型、土壤類型、斷層9個(gè)因子作為泥石流災(zāi)害的影響因子。其中高程、坡度和坡向圖利用GIS空間分析功能和DEM數(shù)據(jù)生成。土地利用圖利用RS技術(shù)中的監(jiān)督分類法進(jìn)行分類，分成森林、居民地、農(nóng)業(yè)用地及河流4個(gè)種類［8-10］，并利用無人機(jī)影像生成的DEM和DOM數(shù)據(jù)進(jìn)行編輯和糾正。其他數(shù)據(jù)從延邊州數(shù)據(jù)庫中獲取。

通過ArcGIS10.2軟件對9種影響因子進(jìn)行狀態(tài)等級劃分，利用GIS分析功能統(tǒng)計(jì)各狀態(tài)等級的泥石流災(zāi)害點(diǎn)比例，依據(jù)泥石流災(zāi)害點(diǎn)比例對影響因子的各狀態(tài)等級賦值。

4.3 基于模型的泥石流災(zāi)害危險(xiǎn)性計(jì)算

利用BP神經(jīng)網(wǎng)絡(luò)和信息量模型計(jì)算各影響因子間的相對權(quán)重。其中BP神經(jīng)網(wǎng)絡(luò)利用70%的隨機(jī)訓(xùn)練數(shù)據(jù)進(jìn)行訓(xùn)練，確定各因子的最佳權(quán)重［11-13］。本文以9種災(zāi)害相關(guān)因子作為輸入變量，輸入層節(jié)點(diǎn)個(gè)數(shù)為9。是否發(fā)生泥石流災(zāi)害為輸出層，輸出層節(jié)點(diǎn)個(gè)數(shù)為1，若發(fā)生泥石流災(zāi)害，輸出因子為1，未發(fā)生泥石流災(zāi)害，輸出因子為2。神經(jīng)元之間的傳遞函數(shù)采用的是正弦S型傳遞函數(shù)，選取16組不同數(shù)據(jù)作為本研究輸入層、隱藏層和輸出層間的連接強(qiáng)度。利用RMSE分析各預(yù)測模型，尋找誤差最小的的連接強(qiáng)度及隱藏層節(jié)點(diǎn)數(shù)。由誤差分析圖表明：誤差最小的隱藏層節(jié)點(diǎn)數(shù)為18（見圖4），誤差最小的連接強(qiáng)度為0.3/0.7，確定最優(yōu)模型為9-18-1結(jié)構(gòu)。

圖4 連接強(qiáng)度誤差分析圖

信息量模型不需要對影響因子的各狀態(tài)等級賦值，直接計(jì)算各災(zāi)害相關(guān)因子不同狀態(tài)等級的信息量值［14-15］。單獨(dú)計(jì)算各因素對泥石流災(zāi)害發(fā)生事件提供的信息量Ixi（y，xi），在實(shí)際運(yùn)算過程中可采取頻率進(jìn)行條件概率的估算，

式中：S為研究區(qū)內(nèi)柵格單元總數(shù)；N為研究區(qū)內(nèi)分布有泥石流災(zāi)害點(diǎn)的柵格單元數(shù)量；Si為研究區(qū)內(nèi)含有評價(jià)因素xi的柵格單元數(shù)量；Ni為分布在因素xi內(nèi)含有泥石流災(zāi)害點(diǎn)的柵格單元數(shù)量。通過上述公式計(jì)算各單元在各影響因子分別作用下的信息量值后，求和得到各柵格單元的信息量值。

4.4 基于GIS技術(shù)繪制泥石流災(zāi)害危險(xiǎn)性區(qū)劃圖

為了泥石流災(zāi)害危險(xiǎn)性區(qū)劃，本研究根據(jù)BP神經(jīng)網(wǎng)絡(luò)計(jì)算的權(quán)重和信息量模型計(jì)算的信息量值，采用柵格數(shù)據(jù)疊置分析法進(jìn)行了計(jì)算和分析。在柵格數(shù)據(jù)內(nèi)部，疊加運(yùn)算是通過像元之間的各種運(yùn)算來實(shí)現(xiàn)的。設(shè)x1，x2，…，xn分別表示第1～n層上同一坐標(biāo)屬性值，f函數(shù)表示各層上屬性與用戶需求之間的關(guān)系，E為疊置后屬性輸出層的屬性值，則

疊加操作的輸出結(jié)果是各層屬性數(shù)據(jù)的算術(shù)運(yùn)算結(jié)果。在本研究中x1，x2，…，x9分別是高程、坡度、坡向、降水量、河流、地層巖性、土地利用類型、土壤類型、斷層等災(zāi)害因子圖層中柵格的屬性值，f的運(yùn)算方式是權(quán)重疊加運(yùn)算法，結(jié)果E是泥石流災(zāi)害危險(xiǎn)性值。結(jié)合信息量法和BP神經(jīng)網(wǎng)絡(luò)模型與GIS技術(shù)繪制的危險(xiǎn)性區(qū)劃圖如圖5和6所示。利用自然斷點(diǎn)法將泥石流災(zāi)害危險(xiǎn)性分為Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ5個(gè)等級，級別越高表明危險(xiǎn)性越高。

圖5 泥石流災(zāi)害危險(xiǎn)性區(qū)劃圖（信息量模型）

圖6 泥石流災(zāi)害危險(xiǎn)性區(qū)劃圖（BP神經(jīng)網(wǎng)絡(luò)）

5 結(jié)果與討論

5.1 信息量模型評價(jià)結(jié)果

由信息量值表明：在延邊州當(dāng)高程＜600 m，坡度為10°～45°，地層巖性為軟弱巖層，年平均降水量在665～991 mm，距離斷裂帶距離小于4 000 m，距離河流距離小于500 m時(shí)，這部分影響因子的狀態(tài)等級對泥石流災(zāi)害的發(fā)生貢獻(xiàn)較大?？傮w而言，高降水量、軟弱巖層以及臨近斷層、河流的區(qū)域泥石流災(zāi)害發(fā)生頻率較高。

統(tǒng)計(jì)延邊州泥石流災(zāi)害危險(xiǎn)性各等級面積占延邊州面積比例分別為18%、23%、20%、23%、16%。各危險(xiǎn)等級泥石流災(zāi)害點(diǎn)比例分別為3%、5%、9%、27%、56%，各等級面積及災(zāi)害點(diǎn)比例圖見圖7。此結(jié)果與延邊州歷史泥石流災(zāi)害發(fā)生區(qū)域較為吻合，說明本次評價(jià)結(jié)果精確度較高。

圖7 各等級面積及災(zāi)害點(diǎn)比例

5.2 BP神經(jīng)網(wǎng)絡(luò)模型評價(jià)結(jié)果

由人工BP神經(jīng)網(wǎng)絡(luò)模型分析得出的相對權(quán)重值表明：9種影響因子對泥石流災(zāi)害的影響為河流＞高程＞坡度＞坡向＞降水量＞土地利用類型＞斷層＞土壤類型＞地層巖性。

各等級面積及泥石流比例如圖8所示。危險(xiǎn)性等級較高的區(qū)域主要分布在安圖縣北部、汪清縣西部、敦化市北部、敦化市南部、和龍市東部、圖們市北部、琿春市中部。

圖8 各等級面積及泥石流比例

延邊州泥石流災(zāi)害危險(xiǎn)性等級圖與泥石流災(zāi)害點(diǎn)相疊加并分析其吻合程度，延邊州Ⅳ級危險(xiǎn)區(qū)、Ⅴ級危險(xiǎn)區(qū)與延邊州泥石流災(zāi)害點(diǎn)吻合程度達(dá)到81%，這說明人工BP神經(jīng)網(wǎng)絡(luò)模型對于延邊州地區(qū)的泥石流災(zāi)害危險(xiǎn)性評價(jià)較為精確。

5.3 模型對比分析

兩種模型的各等級災(zāi)害點(diǎn)與面積的比例如圖9所示。通過兩種模型預(yù)測的危險(xiǎn)區(qū)域大致相同，但也存在著細(xì)微的區(qū)別。兩種模型得到的高危險(xiǎn)區(qū)域都為安圖縣北部、汪清縣西部、敦化市北部及南部地區(qū)、和龍市東部和圖們市北部?，q春市中部地區(qū)在人工BP神經(jīng)網(wǎng)絡(luò)模型的預(yù)測中為Ⅴ級危險(xiǎn)區(qū)，在信息量模型的預(yù)測中主要為Ⅳ級危險(xiǎn)區(qū)。分析其原因主要是，人工BP神經(jīng)網(wǎng)絡(luò)模型側(cè)重于分析不同影響因子之間的關(guān)聯(lián)，而信息量模型則主要考察同一個(gè)影響因子各分級狀態(tài)之間的關(guān)聯(lián)，這也是兩種模型的內(nèi)在區(qū)別。

圖9 各等級災(zāi)害點(diǎn)與分級面積比例

6 結(jié) 語

通過無人機(jī)與3S技術(shù)設(shè)計(jì)的創(chuàng)新實(shí)驗(yàn)，學(xué)生們在完成整個(gè)實(shí)驗(yàn)過程中，了解到無人機(jī)的拍攝及4D產(chǎn)品的生成過程、控制點(diǎn)布設(shè)及GPS在控制點(diǎn)測量中的應(yīng)用、基于遙感技術(shù)繪制災(zāi)害相關(guān)圖層的過程、利用數(shù)學(xué)模型和GIS的空間分析功能相結(jié)合評價(jià)災(zāi)害危險(xiǎn)性的整個(gè)流程。在實(shí)際操作過程中學(xué)生們還可以自主選擇不同的數(shù)學(xué)模型、遙感分類方法以及選擇多種類型的災(zāi)害相關(guān)圖層，試圖提高災(zāi)害危險(xiǎn)性評價(jià)精度。該項(xiàng)試驗(yàn)使學(xué)生們對GIS、RS、GPS及攝影測量等所學(xué)專業(yè)知識的綜合應(yīng)用有了更深層次的了解，有效提高了地信科學(xué)專業(yè)學(xué)生的創(chuàng)新能力、動(dòng)手能力和解決問題的能力，同時(shí)也引起了學(xué)生們對科研創(chuàng)新和專業(yè)知識的濃厚興趣。