馮 琰

（蘇州科技大學(xué) 環(huán)境科學(xué)與工程學(xué)院，江蘇 蘇州 215009）

隨著城市管理科學(xué)化、精細(xì)化水平的提高，城市各行各業(yè)對(duì)城市高精度DEM的應(yīng)用需求也有了更高的要求。在城市規(guī)劃行業(yè)，DEM數(shù)據(jù)是城市規(guī)劃的重要基礎(chǔ)數(shù)據(jù)，利用DEM數(shù)據(jù)可以進(jìn)行土地監(jiān)測(cè)[1]、水源保護(hù)區(qū)規(guī)劃、水系規(guī)劃、景觀規(guī)劃以及規(guī)劃方案[2]的三維效果顯示等；在市政管理部門(mén)，利用DEM數(shù)據(jù)可以提取等高線和坡度分布信息，進(jìn)行道路設(shè)計(jì)和施工、道路橫縱斷面分析、地下管線設(shè)計(jì)管理等；在水利部門(mén)，利用DEM數(shù)據(jù)可以進(jìn)行水流向分析、流域分析、城市積水分析等，對(duì)于城市洪水暴雨預(yù)警等有積極作用。

與全球和國(guó)家級(jí)別DEM相比，城市級(jí)別DEM在制作方法、表現(xiàn)方式、存儲(chǔ)方法和更新維護(hù)方法等方面具有以下特點(diǎn)：（1）由于城市中地形地貌破碎復(fù)雜，具有豐富的特征線，比如馬路牙、下沉式廣場(chǎng)、隧道出入口等都需要表現(xiàn)，因此，必須有高精度、高分辨率DEM表現(xiàn)城市級(jí)別DEM。（2）國(guó)內(nèi)各大中型城市都具有豐富的基礎(chǔ)地形數(shù)據(jù)，例如上海市擁有1∶500、1∶1000和1∶2000多種大比例尺地形數(shù)據(jù)，更新周期分別為0.5 a、1 a、1 a。大比例尺基礎(chǔ)地形數(shù)據(jù)具有豐富的特征點(diǎn)線數(shù)據(jù)，例如道路面、水域面、陡坎線和斜坡線等，這些特征點(diǎn)線可以直接用于城市DEM建模，從而可以代替?zhèn)鹘y(tǒng)的航空攝影測(cè)量方法，提高DEM的生產(chǎn)和更新效率。（3）DEM數(shù)據(jù)結(jié)構(gòu)包括規(guī)則格網(wǎng)（GRID）和不規(guī)則三角格網(wǎng)（TIN），針對(duì)城市地形地貌特點(diǎn)，TIN能較好的顧及城市地貌特征，逼真的表示復(fù)雜地形特征，并能克服地形起伏不大區(qū)域的數(shù)據(jù)冗余問(wèn)題。因此，對(duì)TIN組織管理DEM數(shù)據(jù)方法進(jìn)行了深入研究。（4）城市發(fā)展到一定程度，城市的地形地貌變化相對(duì)變緩，每年只是局部區(qū)域會(huì)有變化和調(diào)整，DEM數(shù)據(jù)的維護(hù)更新也須適應(yīng)這一特點(diǎn)。在原有城市TIN數(shù)據(jù)基礎(chǔ)上實(shí)現(xiàn)局部TIN的更新。不同于傳統(tǒng)DEM建模方法，對(duì)基于城市框架數(shù)據(jù)的大比例尺城市DEM建模方法進(jìn)行了研究，基于不規(guī)則三角格網(wǎng)（TIN）的城市DEM數(shù)據(jù)管理方法，建立城市DEM數(shù)據(jù)庫(kù)，并解決基于TIN的DEM局部增量動(dòng)態(tài)更新問(wèn)題。

1 大比例尺城市DEM建模方法

1.1 DEM建模方法概述

目前國(guó)內(nèi)外常用的DEM制作方法有：（1）全野外測(cè)量數(shù)據(jù)構(gòu)建，這種方法雖然精度很高，但耗時(shí)長(zhǎng)，成本高，更新困難，應(yīng)用范圍也比較?。唬?）將現(xiàn)有的地形圖數(shù)字化，生成數(shù)字網(wǎng)絡(luò)模型，該方法雖然節(jié)約成本，但對(duì)于城市中的細(xì)節(jié)無(wú)法準(zhǔn)確表達(dá)；（3）攝影測(cè)量和遙感的方法，即利用兩相鄰像對(duì)，在全自動(dòng)攝影測(cè)量系統(tǒng)的支持下，通過(guò)輸入地面控制點(diǎn)、相對(duì)定向和絕對(duì)定向以及影像自動(dòng)匹配，生成DEM數(shù)據(jù)；（4）激光雷達(dá)法，這種方法采集到的數(shù)據(jù)為DSM點(diǎn)云數(shù)據(jù)[3]，利用DSM數(shù)據(jù)制作DEM數(shù)據(jù)也是近來(lái)研究的熱點(diǎn)，但是此方法成本高，技術(shù)難度大，難以完全濾掉建筑、樹(shù)木等噪聲的影響。以上幾種建模方法各有優(yōu)缺點(diǎn)，目前比較常用的是攝影測(cè)量與遙感方法。這種建模方法可以得到高精度、高分辨率的DEM數(shù)據(jù)，在有高分辨率影像前提下，可以準(zhǔn)確表達(dá)城市地貌細(xì)節(jié)特征。但是這種方法需要大量的人工干預(yù)，建模效率較低。

以上幾種DEM建模方法中，攝影測(cè)量和遙感方法是目前比較常用的方法，較廣泛地應(yīng)用于國(guó)家和省級(jí)DEM的建模中；近幾年由于激光雷達(dá)技術(shù)的發(fā)展，利用該技術(shù)進(jìn)行DEM建模也在一些地區(qū)得以開(kāi)展研究，但由于其成本高，一些關(guān)鍵技術(shù)一直未得到較好解決，該種建模方法也未得到廣泛推廣。充分利用城市中已有的數(shù)據(jù)資源，大比例尺城市框架數(shù)據(jù)，結(jié)合城市地形地貌特點(diǎn)，開(kāi)展城市DEM建模方法研究。

1.2 建模方法

城市框架數(shù)據(jù)包括高程散點(diǎn)、道路面、道路附屬設(shè)施、鐵路中心線、水域面、街坊面、等高線、陡坎線和斜坡線。利用以上數(shù)據(jù)，在城市DEM中可以準(zhǔn)確描述道路、水域、陡坎、下沉式廣場(chǎng)、隧道出入口等復(fù)雜城市地形地貌?；诖蟊壤叱鞘锌蚣軘?shù)據(jù)的DEM建模方法包括以下幾步：數(shù)據(jù)準(zhǔn)備、數(shù)據(jù)篩查、特征線采集、三維特征線插值和TIN網(wǎng)構(gòu)建。

1.2.1 數(shù)據(jù)準(zhǔn)備

該方法所需數(shù)據(jù)源包括1∶500、1∶1000、1∶2000大比例尺基礎(chǔ)地形數(shù)據(jù)及航空數(shù)字正射影像數(shù)據(jù)。由于上海市1∶500、1∶1000和1∶2000數(shù)據(jù)覆蓋不同區(qū)域范圍，因此，在數(shù)據(jù)準(zhǔn)備階段需要將三種比例尺數(shù)據(jù)進(jìn)行合并，形成覆蓋上海全域的、無(wú)縫拼接的大比例尺城市框架數(shù)據(jù)。此外，為了保證數(shù)據(jù)現(xiàn)勢(shì)性，需要利用最新航空正射影像數(shù)據(jù)對(duì)城市框架數(shù)據(jù)進(jìn)行更新，主要包括：（1）對(duì)與影像上不符的已竣工通車(chē)的市政道路及其以上等級(jí)道路，進(jìn)行現(xiàn)勢(shì)性更新；（2）對(duì)與影像上不符的通航河道，進(jìn)行現(xiàn)勢(shì)性更新；（3）對(duì)城市重點(diǎn)區(qū)域，進(jìn)行現(xiàn)勢(shì)性更新。

1.2.2 數(shù)據(jù)篩查

數(shù)據(jù)篩查是對(duì)城市框架數(shù)據(jù)中不合理部分、缺失部分、邏輯關(guān)系不正確部分進(jìn)行處理，使其滿足DEM建模需求。主要包括以下工作：

（1）數(shù)據(jù)完整性排查：將所有道路面、水域面、街坊面及池塘面按作業(yè)區(qū)塊合并分割，排查出作業(yè)范圍內(nèi)存在的空白區(qū)域，并根據(jù)正射影像賦予空白區(qū)域相應(yīng)的屬性。

（2）水域面篩查：水域面包括河流、湖泊和水庫(kù)等，通過(guò)水域面篩查，保留大面積的，貫通的水域面，刪除其中零星的水域面。

（3）池塘面篩查：通過(guò)池塘面篩查，保留大面積的，連片的池塘面，刪除其中小面積的、零星的池塘面。

（4）道路標(biāo)高點(diǎn)篩查：刪除高架、橋梁及引橋上的標(biāo)高點(diǎn)。

（5）地形地貌高程點(diǎn)篩查：刪除高程特殊的高程散點(diǎn)（如-900 m或9 999 m），篩選每一個(gè)街坊面內(nèi)部的高程點(diǎn)，保留具有代表性的點(diǎn)，刪除一些微地貌上的點(diǎn)。

1.2.3 特征線采集

這里特征線采集是指框架數(shù)據(jù)中不包含的，對(duì)表現(xiàn)城市形態(tài)非常重要的地形地貌特征線，比如下沉式廣場(chǎng)、大堤大壩、隧道出入口等。這些特征線需要在航空影像立體像對(duì)環(huán)境下采集，利用城市框架數(shù)據(jù)中的陡坎線和斜坡線，分割原街坊面，使得建筑物附近的地表盡可能平緩，構(gòu)成新街坊面，并保留該街坊面內(nèi)所有地形地貌高程點(diǎn)。對(duì)沿海地區(qū)的街坊面，利用城市框架的陡坎線和斜坡線，采集大堤邊線及大堤旁的斜坡線。

1.2.4 三維特征線插值

三維特征線插值是指利用城市框架數(shù)據(jù)中的高程散點(diǎn)數(shù)據(jù)對(duì)城市框架數(shù)據(jù)特征線插值賦予高程信息，從而得到三維特征線數(shù)據(jù)。其中等高線的處理方法不同于其他特征線，由于等高線屬性字段中包含高程值，因此，根據(jù)屬性字段可以直接生成三維等高線。

1.2.5 TIN網(wǎng)構(gòu)建

圖1 TIN網(wǎng)

通過(guò)以上四步，可以得到高程散點(diǎn)、三維道路線、三維街坊線、三維等高線、三維水域線等三維特征點(diǎn)線，在第五步中，將利用這些三維特征點(diǎn)線，根據(jù)一定規(guī)則連接成覆蓋整個(gè)區(qū)域且互不重疊的但相鄰接的一系列三角網(wǎng)，以構(gòu)成一個(gè)不規(guī)則三角網(wǎng)TIN。采用三角網(wǎng)增長(zhǎng)算法構(gòu)建TIN網(wǎng)，最終得到以TIN為數(shù)據(jù)結(jié)構(gòu)的DEM數(shù)據(jù)。如圖1所示，為某區(qū)域的TIN三角網(wǎng)數(shù)據(jù)。

2 大比例尺城市DEM空間數(shù)據(jù)庫(kù)

對(duì)于一個(gè)城市DEM數(shù)據(jù)而言，DEM是連續(xù)的、無(wú)縫的。大比例尺城市DEM空間數(shù)據(jù)庫(kù)建庫(kù)目的就是將所有相關(guān)的特征點(diǎn)線數(shù)據(jù)和不規(guī)則格網(wǎng)TIN數(shù)據(jù)按照一定方法有效組織起來(lái)，并根據(jù)其地理分布特點(diǎn)，建立統(tǒng)一空間索引，進(jìn)而可以快速調(diào)度空間數(shù)據(jù)庫(kù)中的任意范圍數(shù)據(jù)，達(dá)到對(duì)全市DEM數(shù)據(jù)的快速無(wú)縫瀏覽和查詢。

2.1 大比例尺城市DEM數(shù)據(jù)組織與管理

對(duì)于TIN模型的存儲(chǔ)結(jié)構(gòu)，國(guó)內(nèi)外的不少學(xué)者已提出了許多種方法，其中最早的是由Peucker及Flower等人提出的按網(wǎng)點(diǎn)坐標(biāo)和高程表及地形點(diǎn)鄰接指針鏈存儲(chǔ)，網(wǎng)點(diǎn)鄰接指針鏈用每個(gè)點(diǎn)所有鄰接點(diǎn)的編號(hào)按順時(shí)針（或逆時(shí)針）方向存儲(chǔ)。這種結(jié)構(gòu)的存儲(chǔ)量小，編輯方便，但計(jì)算量大，不利于TIN的快速顯示。Gold，McCullagh，Tarvelas提出了用網(wǎng)點(diǎn)坐標(biāo)和高程、三角形表、鄰接三角形的存儲(chǔ)方法，該方法具有拓?fù)湫矢叩忍攸c(diǎn)，但存儲(chǔ)量大，編輯不方便。McKenna提出了混合表示網(wǎng)點(diǎn)及三角形鄰接關(guān)系結(jié)構(gòu)。國(guó)內(nèi)也有學(xué)者提出了相應(yīng)的存儲(chǔ)結(jié)構(gòu)，陳曉勇提出了將TIN轉(zhuǎn)化為規(guī)則三角網(wǎng)來(lái)存儲(chǔ)。劉學(xué)軍提出了用點(diǎn)、邊、三角形結(jié)構(gòu)存儲(chǔ)，相對(duì)于前幾種存儲(chǔ)方法而言，采用這種結(jié)構(gòu)占用的存儲(chǔ)量更大，但隨著計(jì)算機(jī)硬件技術(shù)的不斷發(fā)展，計(jì)算速度相對(duì)而言更為重要，由于這種結(jié)構(gòu)拓?fù)潢P(guān)系明顯，便于管理，因而被廣泛的應(yīng)用在數(shù)字高程模型的數(shù)據(jù)組織中。因此，文中也采用這種數(shù)據(jù)組織方法。

圖2 TIN數(shù)據(jù)組織結(jié)構(gòu)

圖2的三個(gè)表，分別用來(lái)記錄組成三角形的頂點(diǎn)號(hào)、邊號(hào)和三角形號(hào)。這三個(gè)表之間，通過(guò)索引指針相互發(fā)生聯(lián)系，即由某一三角形（面），可以檢索出構(gòu)成該三角形的三條邊（線），從而又可檢索出該三角形的三個(gè)頂點(diǎn)（點(diǎn)），由某條邊又可以很方便地檢索出共用該邊的兩個(gè)三角形。這種結(jié)構(gòu)，不僅便于TIN數(shù)據(jù)存儲(chǔ)，建立拓?fù)涔芾恚曳奖銓?shí)現(xiàn)TIN數(shù)據(jù)編輯修改。

2.2 城市DEM空間數(shù)據(jù)庫(kù)管理系統(tǒng)

傳統(tǒng)的DEM等空間數(shù)據(jù)的管理一般是使用文件系統(tǒng)的方式進(jìn)行的，在這種方式下，文件按照生產(chǎn)時(shí)的圖幅進(jìn)行管理，每個(gè)圖幅對(duì)應(yīng)一個(gè)文件。和關(guān)系數(shù)據(jù)庫(kù)相比，文件系統(tǒng)在數(shù)據(jù)安全控制、數(shù)據(jù)共享、并發(fā)操作等方面具有難以克服的缺陷，由于網(wǎng)絡(luò)的發(fā)展，這些問(wèn)題顯得尤為突出。DEM空間數(shù)據(jù)庫(kù)管理系統(tǒng)的建設(shè)，可以控制整個(gè)數(shù)據(jù)庫(kù)系統(tǒng)的運(yùn)行，控制用戶的并發(fā)性訪問(wèn)，執(zhí)行對(duì)數(shù)據(jù)庫(kù)的安全、保密、完整性檢驗(yàn)，實(shí)施對(duì)數(shù)據(jù)的檢索、插入、刪除、修改等操作。（1）數(shù)據(jù)顯示與瀏覽：通過(guò)導(dǎo)航窗口能瀏覽數(shù)據(jù)庫(kù)中數(shù)據(jù)情況，能夠任意開(kāi)窗、放大、縮小和漫游；（2）查詢：可通過(guò)圖幅編號(hào)、區(qū)域邊界（如行政區(qū)域）、名稱(chēng)（如地名）、平面坐標(biāo)等在任意范圍內(nèi)進(jìn)行DEM數(shù)據(jù)的查詢；（3）數(shù)據(jù)入庫(kù)和更新：可將數(shù)據(jù)導(dǎo)入數(shù)據(jù)庫(kù)中，經(jīng)入庫(kù)工具檢查無(wú)問(wèn)題后，對(duì)無(wú)數(shù)據(jù)的區(qū)域數(shù)據(jù)直接入庫(kù)，對(duì)已有數(shù)據(jù)區(qū)域，更新現(xiàn)有數(shù)據(jù)，并與周邊數(shù)據(jù)進(jìn)行融合；（4）數(shù)據(jù)分發(fā)與提?。嚎梢园凑?qǐng)D幅編號(hào)、行政區(qū)域名、經(jīng)緯度范圍、大地坐標(biāo)范圍進(jìn)行分幅，矩形分塊提取數(shù)據(jù)，也可以通過(guò)多幅圖編號(hào)進(jìn)行數(shù)據(jù)分發(fā)和提取的批處理；（5）歷史版本管理：對(duì)DEM更新情況進(jìn)行記錄，保存DEM歷史更新情況和數(shù)據(jù)，實(shí)現(xiàn)不同歷史版本的DEM數(shù)據(jù)管理；（6）日志管理：數(shù)據(jù)庫(kù)變更信息、人員登錄情況等相關(guān)活動(dòng)記錄到日志中進(jìn)行管理；（7）用戶管理：數(shù)據(jù)庫(kù)用戶管理。

2.3 城市DEM空間數(shù)據(jù)庫(kù)維護(hù)和更新方法

隨著城市不斷建設(shè)發(fā)展，城市的DEM也在緩慢發(fā)生變化，這就需要定期對(duì)城市DEM數(shù)據(jù)進(jìn)行更新。在DEM數(shù)據(jù)的生產(chǎn)和更新的作業(yè)過(guò)程中，由于DEM數(shù)據(jù)的復(fù)雜性，該課題將研究相關(guān)的DEM數(shù)據(jù)檢查和維護(hù)更新工具，輔助作業(yè)人員進(jìn)行作業(yè)，從而提高數(shù)據(jù)生產(chǎn)質(zhì)量，保證數(shù)據(jù)庫(kù)數(shù)據(jù)的正確性。具體的功能包括：

（1）局部數(shù)據(jù)的分幅、提取和導(dǎo)出等操作；

（2）數(shù)據(jù)拓?fù)潢P(guān)系等的檢查；

（3）DEM數(shù)據(jù)的完整性、正確性、一致性等檢查；

（4）局部數(shù)據(jù)更新和周邊數(shù)據(jù)進(jìn)行接邊融合；

（5）數(shù)據(jù)生產(chǎn)作業(yè)信息記錄。

3 大比例尺城市DEM空間數(shù)據(jù)庫(kù)應(yīng)用實(shí)例

利用大比例尺城市DEM數(shù)據(jù)進(jìn)行智能化災(zāi)害預(yù)警發(fā)布，下雨時(shí)，通過(guò)對(duì)雨量的預(yù)估，利用DEM數(shù)據(jù)判斷全市易積水路段和易積水街坊，并對(duì)街坊積水程度進(jìn)行定量評(píng)估，并將這些分析結(jié)果實(shí)時(shí)發(fā)布服務(wù)政府決策和方便市民生活。如圖3所示，當(dāng)下雨量為70 mm/h時(shí)，不同街區(qū)的積澇風(fēng)險(xiǎn)分析圖，分別用不同顏色表示積澇風(fēng)險(xiǎn)程度，分別為很危險(xiǎn)、危險(xiǎn)、比較危險(xiǎn)和一般危險(xiǎn)。如圖4所示，為城市中易積水路段分析圖，圖中黑色加粗標(biāo)識(shí)路段為易積水路段。

圖3 積澇風(fēng)險(xiǎn)分析圖

圖4 易積水道路分析圖

4 結(jié)語(yǔ)

大比例尺城市DEM數(shù)據(jù)是“智慧城市”建設(shè)的基礎(chǔ)，是建設(shè)三維城市模型，進(jìn)行城市建設(shè)、規(guī)劃和土地管理、市政工程、氣象和水務(wù)管理等都必不可缺的資源，大比例尺城市DEM空間數(shù)據(jù)庫(kù)的建設(shè)，能為相關(guān)領(lǐng)域的智能管理、智能分析、預(yù)警預(yù)報(bào)、科學(xué)決策和分析提供有力的支持，能在相關(guān)領(lǐng)域產(chǎn)生良好社會(huì)效益和經(jīng)濟(jì)效益。

[1]馮琰，潘琛.結(jié)合地表模型和遙感影像的地表變化監(jiān)測(cè)[J].測(cè)繪科學(xué)，2016，41（3）：44-48.

[2]汪旻琦，馮琰，顧星曄，等.上海市3維城市規(guī)劃輔助審批系統(tǒng)建設(shè)與應(yīng)用研究[J].測(cè)繪與空間地理信息，2013，36（7）：97-100.

[3]顧星曄，馮琰，潘琛，等.DSM數(shù)據(jù)存儲(chǔ)與變化檢測(cè)算法研究[J].測(cè)繪與空間地理信息，2015，38（1）：169-171.

[4]邱文雙，白海麗，趙飛.一種基于數(shù)字線劃圖的大比例尺城市DEM的建立方法[J].礦山測(cè)量，2008（3）：40-41.

[5]諸敏秋，竇煒，趙偉.南京市DEM構(gòu)建及其在城市規(guī)劃中的應(yīng)用[J].城市勘測(cè)，2009（1）：76-79.

[6]盧小平，田繼輝，張劍偉，等.多尺度DEM數(shù)據(jù)的組織與管理方法研究[J].測(cè)繪通報(bào)，2010（6）：12-15.

[7]楊剛，李朝奎，王芳，等.基于ArcSDE的DEM空間數(shù)據(jù)庫(kù)的設(shè)計(jì)與實(shí)現(xiàn)[J].中南林業(yè)科技大學(xué)學(xué)報(bào)，2011，31（4）：169-173.

[8]但尚銘，安海鋒，但玻，等.基于AVHRR和DEM的重慶城市熱島效應(yīng)分析[J].長(zhǎng)江流域資源與環(huán)境，2009，18（7）：680-685.

[9]彭亮，鄒崢嶸，徐濤，等.AutoCAD數(shù)據(jù)在基于Imagis建模中的預(yù)處理與質(zhì)量控制研究[J].測(cè)繪科學(xué)，2008，21（3）：38-42.

[10]吳煥萍，潘懋，胡金星，等.基于空間索引的規(guī)則格網(wǎng)DTM內(nèi)插算法研究[J].地理與地理信息科學(xué)，2004，20（1）：43-46.

[11]季建華，楊軍生.大比例尺DEM生產(chǎn)過(guò)程中的粗差檢測(cè)與質(zhì)量控制的研究[J].現(xiàn)代測(cè)繪，2004，27（3）：31-33.