劉芳，周毅，曾艷菲，陳峰，趙國慶

(浙江省水利水電勘測設(shè)計院，杭州 310002)

0 引言

水利水電工程移民實物調(diào)查工作主要是對工程涉及區(qū)域內(nèi)的各類實物指標(biāo)進(jìn)行數(shù)量和質(zhì)量上的測量、登記，是水利水電工程移民安置中的重要環(huán)節(jié)。目前，社會各階層逐漸意識到以往水利水電工程“重工程，輕移民；重搬遷，輕安置”的弊端[1]及其對移民生產(chǎn)生活帶來的不良影響，并且隨著人民生活水平和自身權(quán)益保護(hù)意識的提高，水利水電工程移民安置過程中涉及的實物補(bǔ)償成為移民關(guān)注的重中之重，稍有不妥便易引起上訪、聚集鬧事等群體性事件。此外，移民實物成果是水利水電建設(shè)項目論證規(guī)模與控制投資的決定性因素之一[2]。因此，作為移民安置基礎(chǔ)的移民實物調(diào)查便顯得尤為重要。

實物調(diào)查內(nèi)容涉及土地、人口、房屋、企事業(yè)單位、專業(yè)項目(交通、電力、通信等)等，目前通用的移民信息采集方式包括傳統(tǒng)紙質(zhì)調(diào)查方式和應(yīng)用調(diào)查信息系統(tǒng)2種[3]，但是2種方式均需要調(diào)查人員攜帶繁重設(shè)備深入項目所在地，長期駐扎進(jìn)行逐戶或逐項調(diào)查，工作內(nèi)容瑣碎復(fù)雜，工作方式機(jī)械單一，需要投入大量的人力、財力、物力[4]。對于關(guān)系到實物補(bǔ)償?shù)姆课菡{(diào)查，更需要調(diào)查人員對每一幢大小、新舊不一的房屋進(jìn)行現(xiàn)場測量，十分耗時費(fèi)力。因而，如何在水利水電移民工程前期階段快速并且準(zhǔn)確地調(diào)查移民實物信息，同時減少對項目地居民的影響，對于工程方案比選、移民安置概算、水庫水電站項目的導(dǎo)截流蓄水方案確定等項目決策具有重要意義。

近年來，隨著地理信息技術(shù)和低空航空攝影測量的快速發(fā)展，無人機(jī)航空攝影已形成專門的產(chǎn)業(yè)服務(wù)鏈，在大比例尺地形圖制作、建筑物三維建模、農(nóng)村地籍測量、城市規(guī)劃、交通選址等方面應(yīng)用廣泛[5-8]。無人機(jī)傾斜攝影是無人機(jī)航空攝影中的重要組成部分，通過在無人機(jī)上搭載多個傳感器，同時對地物的前后左右及頂部5個方向進(jìn)行攝影，獲取不同方向的具有一定傾斜角度的影像，同時記錄拍攝時間、拍攝點的位置信息、航速和航高等基礎(chǔ)信息，可以在短時間內(nèi)完整地獲取地物信息，建立三維模型。無人機(jī)傾斜攝影具有快速靈活、精度高、效率高、成本低等許多優(yōu)勢，對于城市建設(shè)、自然資源調(diào)查和社會發(fā)展具有重要推動作用。目前，無人機(jī)傾斜攝影技術(shù)在水利水電領(lǐng)域的應(yīng)用多集中在地形測繪、與BIM技術(shù)結(jié)合展示工程模型等方面，在移民安置過程中鮮有應(yīng)用[9-10]。本文以浙江省葛岙水庫建設(shè)為契機(jī)，探討無人機(jī)傾斜攝影技術(shù)在移民實物調(diào)查領(lǐng)域的應(yīng)用和精度，以期為水利水電工程移民實物調(diào)查工作的快速、準(zhǔn)確開展提供一種可行的新方法，為項目決策快速提供基礎(chǔ)資料。

1 工作流程

結(jié)合移民實物調(diào)查的特點和無人機(jī)傾斜攝影的主要環(huán)節(jié)制定了本試驗的主要工作流程，具體如圖1所示。

圖1 工作流程圖

1.1 無人機(jī)影像獲取

作業(yè)區(qū)確定是航空作業(yè)的第一步，一般根據(jù)研究對象或項目要求確定。飛行方案設(shè)計則是結(jié)合作業(yè)區(qū)的特點設(shè)計作業(yè)方案，通常先搜集作業(yè)區(qū)的地形圖和歷史衛(wèi)星影像資料，了解作業(yè)區(qū)的自然地理狀況，進(jìn)行飛行方案設(shè)計，初步確定航攝比例尺、航高、影像重疊度、航攝分區(qū)、航攝時間、航飛路線等技術(shù)參數(shù)，再進(jìn)行實地勘察，查看作業(yè)區(qū)內(nèi)是否有高大遮擋物、高壓電塔等，根據(jù)實際情況調(diào)整相關(guān)飛行技術(shù)參數(shù)。實際航飛時將事先規(guī)劃好的航線載入無人機(jī)飛行控制系統(tǒng)，飛控手按照航線設(shè)計控制無人機(jī)執(zhí)行航攝任務(wù)，航飛的同時，無人機(jī)搭載的多傳感器相機(jī)按照預(yù)先設(shè)定的拍攝時間間隔、傾斜角度等參數(shù)同時進(jìn)行多方位的地物攝影，無人機(jī)的定位定向系統(tǒng)同步記錄拍攝時的無人機(jī)位置信息，完成傾斜影像及同步位置數(shù)據(jù)的獲取。飛行過程中，需及時將各個架次獲取的影像及航飛記錄文件下載至電腦，檢查數(shù)據(jù)完整性及數(shù)據(jù)質(zhì)量，對未達(dá)要求的及時進(jìn)行補(bǔ)飛。

1.2 三維模型構(gòu)建

將無人機(jī)獲取的所有傾斜相片、相機(jī)傳感器及鏡頭參數(shù)、位置數(shù)據(jù)全部導(dǎo)入軟件，完成影像的檢查、合并等預(yù)處理步驟；然后導(dǎo)入外業(yè)獲取的作業(yè)區(qū)控制點，將控制點與影像人工進(jìn)行關(guān)聯(lián)，利用航攝相片、航攝地物與航攝相機(jī)三者的幾何位置關(guān)系進(jìn)行空中三角測量，獲取待求點的位置及像片外方位元素等，運(yùn)用多視角傾斜影像匹配技術(shù)生成三維不規(guī)則三角網(wǎng)(triangulated irregular network，TIN)，對整體數(shù)據(jù)進(jìn)行空三加密，獲取加密點云；最后生成數(shù)字表面模型(digital surface model，DSM)，對影像進(jìn)行糾正生成真正射影像(true digital ortho map，TDOM)，建立初始三維模型，并利用影像給初始三維模型添加紋理，形成三維實景模型。模型建立過程中關(guān)注檢查模型質(zhì)量，使之符合精度要求。

1.3 地物要素采集

將三維實景模型及真正射影像加載到專業(yè)的三維測圖及數(shù)據(jù)采集軟件中，利用軟件強(qiáng)大的圖形編輯及專業(yè)數(shù)據(jù)處理功能對作業(yè)區(qū)進(jìn)行圖形屬性一體化采集。借助比較成熟的無人機(jī)傾斜攝影測量地形圖技術(shù)，獲取居民點、地貌、交通、管線、水系、植被等主要地物要素的地理位置及屬性信息，并對各地物要素進(jìn)行分類分層存儲，同一類型要素均分為點、線、面、注記4個圖層。

1.4 移民信息提取

通過無人機(jī)傾斜攝影獲取的地物要素在分類上與現(xiàn)行的《水利水電工程建設(shè)征地移民實物調(diào)查規(guī)范》中的移民實物分類存在部分不一致的情況，而且地物的權(quán)屬關(guān)系無法確定，因此三維測圖完成后，需對采集的地物要素進(jìn)行移民信息提取。移民信息提取采用的方法為外業(yè)調(diào)繪與現(xiàn)場詢問相結(jié)合，具體流程如下：1)針對內(nèi)業(yè)采集的地物要素進(jìn)行修補(bǔ)測，確保屬于移民實物調(diào)查規(guī)范中的實物內(nèi)容清楚、明確、無遺漏；2)根據(jù)移民實物調(diào)查規(guī)范的分類要求對外業(yè)補(bǔ)調(diào)繪后的地物要素建立移民地理數(shù)據(jù)庫，實現(xiàn)圖屬關(guān)聯(lián)；3)通過公安戶籍部門將作業(yè)區(qū)的人口信息以行政村為單位逐戶輸出，通過工商行政主管部分將作業(yè)區(qū)的企事業(yè)單位的基本信息以行政村為單位逐一輸出；4)由作業(yè)區(qū)行政村的相關(guān)負(fù)責(zé)人協(xié)助，將房屋、附屬物、零星林果木信息與3)中的人口信息或企事業(yè)單位相對應(yīng)，明確房屋的權(quán)屬及邊界，現(xiàn)場確定土地的邊界及對應(yīng)的權(quán)屬人/單位；5)請涉及專業(yè)項目的相關(guān)主管部門對采集到的專業(yè)項目信息進(jìn)行核對、補(bǔ)充、確認(rèn)；6)通過作業(yè)區(qū)涉及的鄉(xiāng)鎮(zhèn)/街道獲取社會經(jīng)濟(jì)信息；7)根據(jù)上述成果進(jìn)行內(nèi)業(yè)整理，完善移民地理數(shù)據(jù)庫中的地物信息，關(guān)聯(lián)移民地理數(shù)據(jù)庫中各地物的權(quán)屬人或權(quán)屬單位，調(diào)整房屋和地塊的邊界；8)檢查核對形成最終的水利水電工程移民實物數(shù)據(jù)庫。

2 實驗案例

本實驗結(jié)合葛岙水庫建設(shè)需求，選定浙江省寧波市葛岙水庫的建設(shè)征地范圍為作業(yè)區(qū)進(jìn)行無人機(jī)傾斜攝影測量。

2.1 作業(yè)區(qū)概況

作業(yè)區(qū)位于浙江省寧波地區(qū)，屬于低山丘陵地帶，為集中式居民地。作業(yè)區(qū)范圍根據(jù)葛岙水庫建設(shè)征地范圍確定，面積約3.12 km2。

2.2 無人機(jī)航攝

本實驗采用自研的ZN1000+五鏡頭無人機(jī)航拍系統(tǒng)，該系統(tǒng)包括六旋翼無人機(jī)飛行器、傾斜攝影相機(jī)、飛行控制系統(tǒng)和地面控制站，是一套靈活性高、穩(wěn)定性好的專業(yè)化航測系統(tǒng)，其中傾斜攝影相機(jī)為睿鉑RIY-D2五鏡頭相機(jī)，正攝21.8 mm定焦鏡頭，傾斜36.2 mm定焦鏡頭，傾斜角45°視場角，垂直鏡頭為航向65.4°、旁向56.2°，傾斜鏡頭為航向43.8°、旁向37.2°，可實現(xiàn)多方向影像獲取，具有體積小、重量輕、穩(wěn)定性高和操作便捷的特點。ZN1000+五鏡頭無人機(jī)航拍系統(tǒng)可滿足測繪、電力、水利、國土等行業(yè)的基礎(chǔ)應(yīng)用。

無人機(jī)飛行前對葛岙水庫建設(shè)征地范圍進(jìn)行了實地踏勘，了解作業(yè)區(qū)環(huán)境和地形起伏情況，并收集了帶征地紅線的地形圖和現(xiàn)有控制點坐標(biāo)。無人機(jī)航線是利用Zero-tech軟件結(jié)合地形及建設(shè)征地范圍形狀設(shè)計的最佳航線，本次實驗共飛行15個架次，圖2所示為第1架次的設(shè)計航線。本次航飛相對航高為168 m，航向重疊度為80%，旁向重疊度為75%。像元尺寸大小3.9 μm，像幅為6 000像素×4 000像素。航攝期間結(jié)合航線、地形情況、成圖精度及像控點布設(shè)規(guī)范在作業(yè)區(qū)采用區(qū)域網(wǎng)法布設(shè)像控點，采用選取作業(yè)區(qū)內(nèi)的道路交叉口、地物拐角點等易于辨認(rèn)的特征點及現(xiàn)場人工布設(shè)標(biāo)志2種方式布設(shè)控制點，人工布標(biāo)為每間隔350 m布設(shè)一個十字標(biāo)(白色油漆或石灰)，并在水庫淹沒紅線的凸角轉(zhuǎn)折處予以加密。采用GPS-RTK測量技術(shù)施測，即在作業(yè)區(qū)內(nèi)已知點架設(shè)基準(zhǔn)站和移動站測量計算參數(shù)并校點合格后進(jìn)行像控點測量，每個控制點測量3次取平均值作為觀測結(jié)果，3次測量值的較差要求平面小于2 cm，高程小于3 cm。控制采用CGCS2000平面坐標(biāo)系及1985國家高程基準(zhǔn)。

圖2 一架次飛行航線

2.3 三維建模與地物要素采集

本實驗采用Smart 3D商業(yè)軟件進(jìn)行數(shù)據(jù)處理與三維建模[11]。主要過程為：建立模型工程；對數(shù)據(jù)進(jìn)行預(yù)處理，預(yù)處理包括檢查影像及控制點質(zhì)量，影像勻光勻色處理等；將獲取的影像數(shù)據(jù)導(dǎo)入軟件并添加相機(jī)參數(shù)；選擇空間坐標(biāo)系并設(shè)置其他有關(guān)參數(shù)；添加外業(yè)獲取的控制點，進(jìn)行影像自動匹配并將影像與控制點進(jìn)行人工關(guān)聯(lián)；利用光束平差法進(jìn)行空中三角測量以加密數(shù)據(jù)；生成高密度點云，在其基礎(chǔ)上生成三維TIN網(wǎng)格，再生成白模，最后對白模添加紋理，建立真三維模型[12]，完整并逼真地呈現(xiàn)作業(yè)區(qū)現(xiàn)狀。生成三維實景模型(圖3)的過程中同步獲取真正射影像(圖4)、數(shù)字高程模型、數(shù)字表面模型(圖5)等產(chǎn)品。本實驗生成的真正射影像精度為0.02 m，三維傾斜模型精度為0.06 m。

圖3 實景三維模型(局部)

圖5 數(shù)字表面模型

地物要素采集在商業(yè)軟件EPS中由專業(yè)三維測圖人員進(jìn)行矢量化操作來獲取。首先在軟件中同時加載三維傾斜模型和超大真正射影像，在三維傾斜模型中利用繪點、繪線、繪面、線劃法、五點房等功能快速繪制不規(guī)則地物及標(biāo)準(zhǔn)房屋。繪制房屋的過程同時鍵入房屋屬性信息，包括結(jié)構(gòu)類型、層高、層數(shù)、房屋外立面裝修情況等[13-14]，繪制附屬物、道路、管線、植被、河流湖泊等其他地物時也同步采集相關(guān)屬性信息(圖6)，無需房檐改正的房屋可直接在正射影像上快速繪制。操作過程中充分借助軟件提供的快捷鍵來提高數(shù)據(jù)采集效率，數(shù)據(jù)采集過程中嚴(yán)格執(zhí)行《數(shù)字航空攝影測量測圖規(guī)范》相關(guān)測圖要求。

圖6 ArcMap中主要地物要素圖(局部)

2.4 移民信息提取

采集完各類地物要素后，對部分內(nèi)業(yè)判斷有疑問的地物進(jìn)行外業(yè)修補(bǔ)測，以確保實物內(nèi)容齊全。修補(bǔ)測程序完成后得到完整可靠的地形圖，在地形圖的基礎(chǔ)上，根據(jù)《移民實物調(diào)查規(guī)范》對各地物要素進(jìn)行分類存儲及掛接圖形的屬性數(shù)據(jù)，轉(zhuǎn)換為地理數(shù)據(jù)庫。通過相關(guān)部門及單位獲取人口、企事業(yè)單位等信息，形成以戶或家為單元的數(shù)據(jù)表單，并在作業(yè)區(qū)村鎮(zhèn)負(fù)責(zé)人的協(xié)助下與三維模型獲取的房屋、附屬物、零星林果木、地類地塊等信息關(guān)聯(lián)對應(yīng)，同時明確地物要素采集過程中無法區(qū)分的房屋邊

界、地塊邊界等界址信息，確保各實物信息均確權(quán)到權(quán)屬人或權(quán)屬單位。此外，由專業(yè)項目主管部門確認(rèn)專業(yè)項目權(quán)屬邊界，并補(bǔ)充相關(guān)屬性信息。通過上述過程獲取的實物數(shù)據(jù)及社會經(jīng)濟(jì)信息組成完整的矢量型水利水電工程移民實物數(shù)據(jù)庫。

2.5 移民實物成果精度分析

《移民實物調(diào)查規(guī)范》中規(guī)定了不同階段的人口、房屋、耕園地及林草地等主要實物的調(diào)查精度，其中項目建議書階段的調(diào)查精度要求分別是人口(±10%)、房屋(±10%)、耕園地(±10%)、林草地(±15%)。由于三維傾斜模型獲取的移民信息主要為房屋、土地、附屬物和線狀專業(yè)項目，因此本實驗綜合選取房屋及地類地塊進(jìn)行精度比對。結(jié)合項目關(guān)注實物尺寸、規(guī)模等數(shù)量信息的特點，本實驗在整個作業(yè)區(qū)范圍內(nèi)均勻選取了102幢房屋與喜利得激光測距儀測量的房屋面積進(jìn)行比對，選取24個地塊與采用GPS-RTK技術(shù)獲取的土地勘測定界成果進(jìn)行對比，用于檢查通過三維傾斜模型測量的地物精度。三維傾斜模型測量的房屋面積與現(xiàn)場實際測量獲取的房屋面積對比分析情況如表1和圖7所示。三維傾斜模型測量的地塊面積與土地勘測定界的地塊面積對比如圖8所示。

表1 房屋面積差值范圍表

圖7 房屋誤差絕對差值散點圖

圖8 地塊面積相對誤差

根據(jù)表1中的平均值及中誤差，可以看出本次實驗房屋面積的精度較高，房屋面積相對誤差最大5.74%，最小的僅0.07%，各幢房屋面積的平均相對誤差比例為1.72%，達(dá)到實物調(diào)查可行性研究報告階段的±3%精度要求。圖7顯示三維傾斜模型測量的房屋面積與實際獲取的房屋面積差值集中在±3 m2之間，占64.71%。

根據(jù)圖8結(jié)果，地塊面積相對誤差最小為0.09%，最大為8.37%(溝渠)，平均相對誤差為2.30%，相對誤差超過3%的地塊共有7個(均為林地)，占選取的檢查地塊個數(shù)的29.17%。從圖中可以看出相對誤差主要取決于地塊面積，地塊面積越小，誤差越大。經(jīng)過檢查，高大樹木遮擋及房屋密集分布是造成部分房屋及地塊誤差較大的主要原因。此外，信息采集與數(shù)據(jù)精度分析過程中，發(fā)現(xiàn)如下特點：1)三維傾斜攝影測量對于有弧形墻、缺角、非矩形、上下層墻體不對應(yīng)等不規(guī)則的房屋特別適用，快捷方便遠(yuǎn)超現(xiàn)場測量；2)樓間距過窄會導(dǎo)致三維實景模型變形，從而增大誤差；3)房屋墻體被綠色植物覆蓋時誤差增加；4)聯(lián)排房屋不易區(qū)分，測量時容易將多個權(quán)屬人/單位的房屋合并測量；5)祖?zhèn)骼吓f房屋一般單幢房屋分屬不同權(quán)屬人，基于三維傾斜模型測量時無法區(qū)分，必須依賴現(xiàn)場問詢；6)溝渠、道路等線狀地類受樹木遮擋，測量時需仔細(xì)；7)各地民居各具特點，測量前制定針對性方案極為必要。

綜上，本實驗基于無人機(jī)傾斜攝影開展的水利水電工程移民實物調(diào)查工作完全可以滿足《移民實物調(diào)查規(guī)范》中對項目建議書階段移民實物調(diào)查的精度要求，基本滿足可行性研究報告階段移民實物調(diào)查的精度要求(人口、房屋、耕園地精度要求均為±3%，林草地精度要求為±5%)，但在房屋權(quán)屬劃分到戶、地塊劃分到組等權(quán)屬界定方面限制性較多，需綜合現(xiàn)場問詢及修補(bǔ)測。

3 結(jié)束語

本文詳細(xì)闡述了應(yīng)用無人機(jī)傾斜攝影技術(shù)開展移民實物調(diào)查的方法及工作流程，并以葛岙水庫移民實物調(diào)查為例，測試了無人機(jī)傾斜攝影技術(shù)在水利水電工程建設(shè)征地范圍三維實景建模及移民實物信息提取方面的可行性。無人機(jī)傾斜攝影測量技術(shù)在應(yīng)對任務(wù)繁重、時間緊迫且入戶調(diào)查難度大的移民實物調(diào)查項目(尤其是未開展大比例尺地形圖測繪的項目)時具有明顯優(yōu)勢，影像獲取簡單靈活，數(shù)據(jù)采集方便快捷，可大量節(jié)省外業(yè)調(diào)查的人力物力，提高移民信息采集的效率和準(zhǔn)確率，為移民投資的準(zhǔn)確估算及移民安置工作的穩(wěn)步開展奠定良好基礎(chǔ)。而且在應(yīng)用無人機(jī)傾斜攝影技術(shù)開展移民實物調(diào)查的過程中間接實現(xiàn)了水利水電工程建設(shè)征地范圍歷史風(fēng)貌的高保真留存，工程建設(shè)后為移民保留了一份永久的影像記憶檔案。需要指出的是，本文僅測試了低山丘陵地區(qū)，并未對高山地區(qū)及高大建筑物密集的城區(qū)進(jìn)行驗證，這2種區(qū)域的可行性需要進(jìn)一步測試驗證；本實驗中影像地面分辨率為0.2 m，與精度要求較高的工程可行性研究階段的移民實物調(diào)查尚有差距，將來可以進(jìn)一步驗證影像地面分辨率高達(dá)5 cm甚至3 cm時該方法的適用性。本文僅探討了無人機(jī)傾斜攝影技術(shù)在移民實物調(diào)查中的應(yīng)用可行性，未來還可以探究無人機(jī)傾斜攝影技術(shù)在水庫導(dǎo)截流方案比選、與BIM技術(shù)結(jié)合輔助庫周交通改復(fù)建及移民安置區(qū)建設(shè)等移民安置其他方面的應(yīng)用。