李菲菲,高麗君

(四川省地震局測(cè)繪工程院，四川 雅安 625000)

1985年12月，四川省地震局地震測(cè)量隊(duì)周如品等同志根據(jù)搜集到的資料編纂了《川西大地形變測(cè)量網(wǎng)、點(diǎn)布設(shè)圖集》。圖集中記載了川西大地形變測(cè)量場(chǎng)地的基本情況及場(chǎng)地圖。其中，場(chǎng)地圖選用四川省測(cè)繪局1984年9月編制印刷的1∶150萬(wàn)《四川省地圖》作底圖。圖集中的場(chǎng)地圖多年來(lái)被廣泛應(yīng)用于會(huì)商報(bào)告、異常核實(shí)報(bào)告以及涉及跨斷層觀測(cè)資料的科研論文中。三十多年來(lái)，道孚、折多塘等場(chǎng)地的測(cè)網(wǎng)曾被重新布設(shè)，很多場(chǎng)地的觀測(cè)環(huán)境也發(fā)生了較大改變。2017年3月，由中國(guó)地震局監(jiān)測(cè)預(yù)報(bào)司統(tǒng)一組織管理，中國(guó)地震局第一監(jiān)測(cè)中心牽頭和協(xié)調(diào)管理，包括四川省地震局在內(nèi)的十多家單位分片具體實(shí)施的財(cái)政專(zhuān)項(xiàng)——全國(guó)跨斷層場(chǎng)地優(yōu)化改造項(xiàng)目啟動(dòng)。依托該項(xiàng)目，四川省地震局測(cè)繪工程院(以下簡(jiǎn)稱(chēng)“測(cè)繪院”)對(duì)28個(gè)跨斷層場(chǎng)地進(jìn)行了優(yōu)化改造和測(cè)點(diǎn)維護(hù)。至此，《川西大地形變測(cè)量網(wǎng)、點(diǎn)布設(shè)圖集》中的場(chǎng)地圖已幾乎不再適用，所以有必要繪制新的場(chǎng)地圖。一幅完整的場(chǎng)地圖包含了地形、地貌、地物等要素，實(shí)地測(cè)圖的工作量及工作成本都太大。高分辨率的地形數(shù)據(jù)可通過(guò)空間數(shù)據(jù)服務(wù)商訂購(gòu)，但費(fèi)用昂貴；即便在測(cè)繪地理信息局購(gòu)買(mǎi)，一幅1∶1萬(wàn)DEM數(shù)據(jù)的價(jià)格也需要2 000元。在不具備無(wú)人機(jī)這樣先進(jìn)的航測(cè)設(shè)備的情況下，測(cè)繪院利用免費(fèi)的DEM數(shù)據(jù)及Google earth、Global Mapper等軟件重新制作了28個(gè)跨斷層場(chǎng)地的矢量場(chǎng)地圖。

1 獲取DEM數(shù)據(jù)

目前國(guó)內(nèi)主流的免費(fèi)高程數(shù)據(jù)資源主要有SRTM、ASTER GDEM、DLR三種。公開(kāi)的SRTM地形數(shù)據(jù)的空間分辨率精度為90 m，可通過(guò)SRTM官方網(wǎng)站http://srtm.csi.cgiar.org/srtmdata/或地理空間數(shù)據(jù)云網(wǎng)站http://www.gscloud.cn/下載。ASTER GDEM V1原始數(shù)據(jù)局部地區(qū)存在異常；ASTER GDEM V2版則采用了一種先進(jìn)的算法對(duì)V1版GDEM影像進(jìn)行了改進(jìn)，提高了數(shù)據(jù)的空間分辨率精度和高程精度，其全球空間分辨率為30 m。DLR數(shù)據(jù)是2000年2月11日德國(guó)航天局在美國(guó)發(fā)射的“奮進(jìn)”號(hào)航天飛機(jī)上用X波段雷達(dá)進(jìn)行地形測(cè)繪所獲得的數(shù)據(jù)，數(shù)據(jù)呈網(wǎng)狀，有些地方覆蓋不全，但在能覆蓋到的地方精度比ASTER GDEM V2高，約為10 m。綜合考慮三種數(shù)據(jù)的優(yōu)缺點(diǎn)后，選擇ASTER GDEM V2版數(shù)據(jù)作為場(chǎng)地圖的DEM數(shù)據(jù)。以侏倭場(chǎng)地為例，其概略坐標(biāo)平均值為(100.3°E，31.6°N)。在地理空間數(shù)據(jù)云網(wǎng)站中，找到“GDEMV2 30M 分辨率數(shù)字高程數(shù)據(jù)”列表，檢索侏倭場(chǎng)地所在的數(shù)據(jù)文件并進(jìn)行下載，即可得到名為“ASTGTM2_N31E100.zip”的文件。

2 制作草圖

以場(chǎng)地坐標(biāo)為制圖中心，在MapInfo軟件中繪制600 m×600 m的矩形框，記錄矩形框的邊界坐標(biāo)(用于控制生成等高線的區(qū)域邊界)，同時(shí)在矩形框中繪制比例尺、棄測(cè)測(cè)點(diǎn)、在用測(cè)點(diǎn)、測(cè)線等要素(如圖1)。

3 生成等高線

解壓下載到的數(shù)據(jù)文件，用Global Mapper打開(kāi)文件名中帶“dem”的tif格式文件。通過(guò)Analysis-Generate Contours (from Terrain Grid)，生成目標(biāo)區(qū)域的等高線。通過(guò)File-Export-Export Elevation Grid Format，輸出Surfer Grid格式的高程網(wǎng)格數(shù)據(jù)。用Surfer軟件打開(kāi)高程網(wǎng)格數(shù)據(jù)，發(fā)現(xiàn)等高線中存在一些小的彎曲，可用軟件自帶的低通濾波器剔除這些小的彎曲(如圖2)。

圖3 濾波前后等高線對(duì)比

4 數(shù)字化地物

從MapInfo軟件中將矩形框?qū)С龅紾oogle earth軟件中(由于將矩形框的海拔高度設(shè)置為“貼近地面”，導(dǎo)致顯示出的矩形框有變形)。在Google earth軟件中將矩形框內(nèi)的道路、河流、房屋的輪廓線及斷層線、植被分界線繪制出來(lái)(如圖3)，并分別保存為kml格式文件。在Global Mapper中，將所有用“線”表示的要素逐一打開(kāi)，并分別輸出為Surfer BLN矢量格式的文件。

5 編輯顏色、標(biāo)注、圖例等

在Surfer中，將等高線、道路、河流、房屋、測(cè)點(diǎn)、測(cè)線、斷層線等各要素分圖層打開(kāi)，并“按坐標(biāo)疊加圖形”。編輯各圖層的線條屬性、填充屬性，添加標(biāo)注文字。選中所有對(duì)象，導(dǎo)出為emf格式文件。用Adobe Illustrator軟件打開(kāi)并編輯emf格式文件。需要編輯的內(nèi)容包括房屋的填充圖案、小路的虛線表示、圖框的加粗、標(biāo)注文字的移動(dòng)等，編輯完成后的成果圖如圖4(b)所示。

圖4 侏倭場(chǎng)地新舊場(chǎng)地圖對(duì)比(圖4(a)來(lái)自《川西大地形變測(cè)量網(wǎng)、點(diǎn)布設(shè)圖集》)

6 討論

重新繪制后的場(chǎng)地圖，測(cè)點(diǎn)、測(cè)線及各個(gè)地物的位置均由坐標(biāo)控制，保證了各要素之間相對(duì)位置的準(zhǔn)確性。場(chǎng)地圖中各要素均以矢量格式保存，方便編輯及不同要求的圖形輸出。由于Google earth影像圖并非實(shí)時(shí)更新，依靠Google earth影像圖提取的部分地物信息也具有一定的滯后性，滯后時(shí)間短則幾個(gè)月，多則幾年。盡管如此，用上述方法繪制出的場(chǎng)地圖已能滿足各類(lèi)報(bào)告中對(duì)觀測(cè)場(chǎng)地的展示需求。每年更新一次場(chǎng)地圖可改善地物信息的滯后性，但要獲取實(shí)時(shí)的、精細(xì)的場(chǎng)地圖，建議配備專(zhuān)業(yè)航測(cè)設(shè)備及配套的內(nèi)業(yè)處理軟件。