余思汗，劉 超，王建勇，單帥帥，楊 順

（1. 寧夏回族自治區(qū)地震局，寧夏 銀川 750001；2. 北方信息控制研究院集團有限公司，江蘇 南京 210000）

0 引言

發(fā)生破壞性地震后，快速、精準地開展地震災(zāi)害損失評估工作對初期應(yīng)急處置和應(yīng)急救援力量部署尤為重要，而詳細的建筑物數(shù)據(jù)是評估工作的基礎(chǔ)[1]，也是應(yīng)急基礎(chǔ)數(shù)據(jù)庫的重要組成[2]。建筑物數(shù)據(jù)的精細程度直接關(guān)系到地震災(zāi)害評估的準確性，繼而影響到現(xiàn)場救援的目標精準性[3]。然而建筑物數(shù)據(jù)由于數(shù)據(jù)量大、獲取難、更新難、精細度差等原因，導致地震應(yīng)急基礎(chǔ)數(shù)據(jù)庫中該數(shù)據(jù)的現(xiàn)實性、全面性等都和實際情況有差異，因此解決好建筑物空間分布數(shù)據(jù)是當前應(yīng)急基礎(chǔ)數(shù)據(jù)庫更新工作中的重中之重。

為了保證地震災(zāi)害評估結(jié)果的準確性，部分學者對建筑物空間分布格網(wǎng)化方法做了深入研究[4]，主要有以下四種方法：均勻分布法[5]、直接提取法[6]、現(xiàn)場調(diào)查法[7]和直接推算法[8]。雖然這些方法在一定程度上解決了建筑物空間分布的問題，但是每種方法都存在一些不足之處：均勻分布法雖容易實現(xiàn)，但其結(jié)果與實際情況相差甚遠，嚴重影響到了地震災(zāi)害評估的準確性；直接提取法需要大量的技術(shù)人員進行遙感目視解譯，工作繁重且成本較高，不利于大范圍開展；現(xiàn)場調(diào)查法費時費力，工作周期長且效率低，也難以推廣應(yīng)用；間接推算法過于依賴精細化的人口空間分布數(shù)據(jù)，而數(shù)據(jù)難獲取則成為了制約該方法的難題。為了高效解決建筑物空間分布數(shù)據(jù)問題的同時，又可以節(jié)約人力物力，本文以寧夏回族自治區(qū)石嘴山市平羅縣為例，借助于測繪部門提供的公眾服務(wù)數(shù)據(jù)，提出了基于多源數(shù)據(jù)的建筑物空間分布格網(wǎng)化研究方法。

1 研究區(qū)概況及數(shù)據(jù)準備

1.1 研究區(qū)概況

本文以寧夏回族自治區(qū)石嘴山市平羅縣為例，作為研究區(qū)域（圖1）。平羅縣地處國家絲綢之路經(jīng)濟帶、呼包銀能源金三角、寧夏沿黃城市帶等經(jīng)濟區(qū)前沿地帶，介于東經(jīng)105°57′42″～106°58′02″、北緯 38°36′18″～39°51′13″之間。截至2018 年年底，平羅縣常住總?cè)丝?1.6萬人，下轄13 個鄉(xiāng)鎮(zhèn)，縣域總面積2060 平方公里，城關(guān)鎮(zhèn)是平羅縣政府駐地，也是城區(qū)所在地。地形分為賀蘭山區(qū)、山前洪積扇區(qū)、西大灘碟形洼地、沖積平原區(qū)、靈鹽臺地和河灘區(qū)六大單元，主要活動斷裂有賀蘭山東麓斷裂[9]、銀川—平羅隱伏斷裂[10]、蘆花臺隱伏斷裂[11]等，其中賀蘭山東麓斷裂為1739 年銀川—平羅8 級地震的發(fā)震構(gòu)造[12]。

圖1 研究區(qū)域Fig.1 Study area

1.2 數(shù)據(jù)來源及處理

（1）遙感影像數(shù)據(jù)：遙感影像數(shù)據(jù)來源于天地圖網(wǎng)站，是我國國產(chǎn)的高分2 號衛(wèi)星遙感影像，影像的空間分辨率為1m，影像發(fā)布時間為2019 年第四季度。

（2）建筑物數(shù)據(jù)：建筑物數(shù)據(jù)來源于寧夏自然資源廳提供的《2017 年寧夏地理空間基礎(chǔ)信息共享庫（公眾版）》[13]。該建筑物數(shù)據(jù)基本上由數(shù)字城市建筑物數(shù)據(jù)、地理國情監(jiān)測建筑物數(shù)據(jù)組成，故在部分城市地區(qū)可能存在兩種數(shù)據(jù)重疊在一起的問題，這時就需要將地理國情建筑物數(shù)據(jù)刪除，保留精度更高的數(shù)字城市建筑物數(shù)據(jù)，這樣有利于提高建筑物的總體精度。由于數(shù)據(jù)是2017 年的，而近幾年城市化發(fā)展迅速，因此數(shù)據(jù)會存在現(xiàn)實性的問題，需對建筑物數(shù)據(jù)進行更新，以天地圖的影像服務(wù)作為底圖，采用目視解譯的方式進行更新[14]。具體數(shù)據(jù)處理流程如圖2 所示。

圖2 數(shù)據(jù)處理流程圖Fig.2 Data processing flow chart

1.3 城市建筑物調(diào)查

由于城區(qū)內(nèi)建筑物數(shù)據(jù)缺乏層數(shù)屬性，不能夠精確表達建筑物面積，而城區(qū)內(nèi)又為多層建筑物密集的區(qū)域，因此就需要對城區(qū)內(nèi)建筑物進行調(diào)查。采用了人工外業(yè)和無人機相結(jié)合的調(diào)查手段，主要是利用小型無人機成本低、數(shù)據(jù)采集效率高、機動靈活等優(yōu)點[15]，對城區(qū)內(nèi)建筑物進行航拍調(diào)查，從而判斷建筑物的層數(shù)。將航拍的照片時間加小數(shù)點后3 位的經(jīng)緯度來命名，在GIS 軟件中對建筑物數(shù)據(jù)屬性表新增舉證照片字段，使用GIS 軟件中的超級鏈接功能，把相對應(yīng)的無人機航拍照片的存儲位置及名稱填充進該字段，便可在GIS 軟件中使用超級鏈接功能進行查看建筑物的舉證照片，為建筑物的層數(shù)提供照片舉證依據(jù)。利用GIS 軟件中的字段計算器功能，將調(diào)查的結(jié)果賦值給同一建筑物群的所有單體建筑物，這樣就完成了城區(qū)內(nèi)建筑物層數(shù)的調(diào)查工作。圖3 為調(diào)查后的城區(qū)內(nèi)建筑物分布圖。

圖3 城區(qū)內(nèi)建筑物分布圖Fig.3 Distribution of buildings in urban area

1.4 農(nóng)村建筑物面積估算

對于農(nóng)村的建筑物數(shù)據(jù)，絕大多數(shù)都是采用了地理國情監(jiān)測的建筑物數(shù)據(jù)，該數(shù)據(jù)在農(nóng)村地區(qū)的采集標準較為粗略，往往將相連的房屋包括院落也一并采集為建筑物用地，這將大大影響到了建筑物的精度，因此需隨機選取15個自然村落，對其建筑物進行采集，同樣采用以天地圖為底圖的目視解譯的方法，采集的標準是以建筑物的幾何輪廓為邊界，盡可能采集的精確。樣本點建筑物數(shù)據(jù)采集完成之后，將新采集的農(nóng)村建筑物數(shù)據(jù)與該村落的地理國情檢測數(shù)據(jù)進行對比，得到農(nóng)村建筑物的實際比值，對15 個樣本點的實際比值求平均數(shù)，該平均值則為農(nóng)村建筑物的實際比值，最后獲得的平均實際比值為0.46。利用GIS 軟件的字段計算器功能，對所有的農(nóng)村建筑物面積進行重新賦值，把原來的建筑物圖斑幾何面積乘以實際比值，就是估算的農(nóng)村建筑物的實際面積。

2 研究方法及過程

2.1 研究思路

本文利用數(shù)字城市建筑物數(shù)據(jù)、地理國情監(jiān)測建筑物數(shù)據(jù)，通過GIS 軟件的空間分析功能找出重疊部分，并刪除掉重疊部分中較粗略的監(jiān)測數(shù)據(jù)，并對數(shù)據(jù)進行拓撲檢查，徹底排除掉面重疊的數(shù)據(jù)錯誤，再利用天地圖影像服務(wù)作為底圖，對建筑物進行目視解譯更新，從而得到現(xiàn)實性較強的建筑物數(shù)據(jù)。根據(jù)調(diào)查的城區(qū)內(nèi)建筑物的層數(shù)和農(nóng)村建筑物實際比值，建立建筑物面積的數(shù)學模型，從而確定建筑物的實際面積，最終通過矢量數(shù)據(jù)空間運算的方式得到研究區(qū)建筑物空間分布格網(wǎng)。

2.2 建筑物空間格網(wǎng)化的面積計算方法

首先是創(chuàng)建格網(wǎng)，利用GIS 軟件在研究區(qū)行政區(qū)劃矢量數(shù)據(jù)上建立最小邊界幾何，再使用GIS 軟件中的漁網(wǎng)工具生成500m×500m 的格網(wǎng)。建筑物空間格網(wǎng)的面積是格網(wǎng)內(nèi)所有建筑物面積的累加，因此建筑物空間格網(wǎng)的面積計算要基于每個單體建筑物的面積。假設(shè)格網(wǎng)單元內(nèi)的建筑物面積之和為S，第i 個單體建筑物的占地面積為Si，單元區(qū)域內(nèi)共有n 個單體建筑物，第i 個單體建筑物的層數(shù)為m，則格網(wǎng)內(nèi)的面積S 為：

對于研究區(qū)的農(nóng)村地區(qū)而言，幾乎全部為單層建筑物，只有少數(shù)城鄉(xiāng)結(jié)合部有兩層或三層的建筑物，但是相對于整體的農(nóng)村建筑物數(shù)據(jù)來說，幾乎可以忽略不計，因此全部按單層建筑物來處理。根據(jù)建立的模型，即農(nóng)村地區(qū)m 值恒為1。結(jié)合該方法建立的模型，在GIS軟件中運用空間鏈接功能，將格網(wǎng)單元內(nèi)的建筑物面積全部累加賦值給格網(wǎng)。在運算過程中，可能存在一個建筑物分別處于兩個或者兩個以上的格網(wǎng)單元內(nèi)，解決辦法是以建筑物的幾何中心為目標進行賦值。即根據(jù)建筑物的幾何中心所處的格網(wǎng)單元，將建筑物面積賦值給該格網(wǎng)單元作為一個累加單體建筑物。將建筑物的面積全部賦值給格網(wǎng)單元之后，就完成了建筑物空間分布格網(wǎng)化的數(shù)據(jù)。圖4 為研究區(qū)500m×500m 尺度的建筑物格網(wǎng)面積粒度。

3 結(jié)果驗證

為了驗證建筑物空間分布格網(wǎng)化方法的可行性，先是抽樣了單個格網(wǎng)進行驗證，隨機抽取了25 個格網(wǎng)，在GIS 軟件中使用裁剪功能對建筑物進行裁剪，重新統(tǒng)計了樣本格網(wǎng)的建筑物面積，并與原格網(wǎng)面積進行對比，其相對誤差率最大為4.22%，說明使用建筑物幾何中心的方法是可行的。對于格網(wǎng)數(shù)據(jù)整體的驗證，則是采用了2018 年寧夏統(tǒng)計局發(fā)布的統(tǒng)計年鑒人均居住面積和研究區(qū)內(nèi)常住人口數(shù)據(jù)來進行驗證。截至2018 年年底，平羅縣的常住人口為31.6 萬人，統(tǒng)計部門發(fā)布的人均住房面積為35平方米[16]，根據(jù)常住人口和人均面積得出研究區(qū)內(nèi)的居住建筑物面積為1106 萬平方米，而格網(wǎng)化后的建筑物面積約為1507.68 萬平方米，居住建筑物面積占格網(wǎng)化后面積的73.36%。考慮到除了居住建筑物之外，還有企業(yè)廠房、商業(yè)綜合體、學校、醫(yī)院、辦公樓、體育場館等建筑物，因此驗證建筑物空間分布格網(wǎng)化的結(jié)果比較合理，建立的模型和研究方法也是可行的。

4 結(jié)論與討論

本文采用基于多源數(shù)據(jù)建筑物面積空間分布累加的方法，獲得了寧夏平羅縣的空間分布格網(wǎng)建筑物數(shù)據(jù)。利用GIS 軟件的空間分析能力，確定了唯一的單體建筑物，將高分辨率遙感影像作為底圖，采用了目視解譯的方法進行了建筑物的更新，通過現(xiàn)場人工調(diào)繪和無人機調(diào)查的方式，確定了城區(qū)內(nèi)建筑物群的層數(shù)，從而計算出了單體建筑物面積，農(nóng)村地區(qū)則采用了抽樣核查的方式，給出了農(nóng)村建筑物的實際比值系數(shù)，最終通過空間運算得到了研究區(qū)500m×500m 尺度的空間分布格網(wǎng)建筑物數(shù)據(jù)。該數(shù)據(jù)為震后災(zāi)害評估提供了可靠的空間分布建筑物數(shù)據(jù)，使得評估結(jié)果的準確性得到了提高，也為特大災(zāi)害救援力量部署和應(yīng)急處置決策提供依據(jù)。

由于城區(qū)內(nèi)建筑物比較復雜，雖然利用無人機進行調(diào)查提高了效率，但是整體的工作量依然較大，不利于單體建筑物面積的快速提取。在今后的研究工作中，應(yīng)該注重于城市建筑物面積模型的精確估算，減少建筑物調(diào)查的工作量，從而提高基于多源數(shù)據(jù)的建筑物空間分布格網(wǎng)方法應(yīng)用價值。