根據(jù)縣級(jí)國土調(diào)查數(shù)據(jù)庫特點(diǎn)及第三次全國國土調(diào)查現(xiàn)狀分類,該文提出了以土地利用數(shù)據(jù)為基礎(chǔ),生產(chǎn)1 km格網(wǎng)路網(wǎng)密度的總體設(shè)計(jì)流程,如圖1所示。生產(chǎn)1 km格網(wǎng)路網(wǎng)密度的總體流程可以概括為數(shù)據(jù)預(yù)處理、骨架線抽取、路網(wǎng)密度計(jì)算和產(chǎn)品質(zhì)檢4部分。
本數(shù)據(jù)集具體生產(chǎn)步驟如下:
(1)數(shù)據(jù)預(yù)處理
數(shù)據(jù)預(yù)處理包括分縣抽取、投影轉(zhuǎn)換、數(shù)據(jù)拼接、屬性整合4個(gè)部分。首先從2800多個(gè)縣區(qū)的地類圖斑層中分別抽取鐵路用地、公路用地、城鎮(zhèn)村道路用地及軌道交通用地圖斑。由于各縣區(qū)為高斯克呂格投影3度分帶,不同經(jīng)度劃分帶號(hào)不同,為便于數(shù)據(jù)整合將分縣抽取后的圖層統(tǒng)一轉(zhuǎn)換為Albers坐標(biāo)系投影。再利用矢量轉(zhuǎn)柵格系統(tǒng)提供的數(shù)據(jù)拼接功能,將分縣區(qū)的數(shù)據(jù)以省為單位拼接,最后用融合工具消除坐落權(quán)屬的影響,避免軌道交通用地及道路圖斑破碎化。
圖1
?
數(shù)據(jù)集生產(chǎn)技術(shù)流程 (2)骨架線抽取
將數(shù)據(jù)預(yù)處理后的鐵路用地、公路用地、城鎮(zhèn)村道路用地及道路軌道交通用地圖斑分省組織目錄,借助WJ-III地圖工作站軟件,通過修改分省對應(yīng)的索引文件讀取綜合知識(shí)庫,分別抽取鐵路用地骨架線、公路用地骨架線、城鎮(zhèn)村道路用地骨架線及軌道交通用地骨架線。線狀網(wǎng)絡(luò)數(shù)據(jù)拓?fù)潢P(guān)系的優(yōu)劣是決定線狀網(wǎng)絡(luò)數(shù)據(jù)質(zhì)量的重要因素。WJ-III地圖工作站相較于其他軟件,基于拓?fù)湎嗨菩缘牡缆肪W(wǎng)綜合算法,為線狀網(wǎng)絡(luò)數(shù)據(jù)的連通性及正確性提供保障,在可視化效果上更貼合線狀網(wǎng)絡(luò)數(shù)據(jù)的應(yīng)用,且運(yùn)行效率快
[17]。
(3)路網(wǎng)密度計(jì)算
①算法及模型
由定義可知,路網(wǎng)密度在數(shù)量上表示為單位區(qū)域面積的道路里程,即:
\(Md=\sum L/F\) (1)
式中,\(Md\)為路網(wǎng)密度(km/km2),\(\sum L\)為道路總長度(km),F為1 km格網(wǎng)面積。
②公里格網(wǎng)匯總計(jì)算
使用GDAL(地理空間抽象數(shù)據(jù)庫)作為讀寫空間數(shù)據(jù)的底層構(gòu)建庫,調(diào)用Java接口開發(fā)路網(wǎng)密度格網(wǎng)計(jì)算工具。首先生成全國統(tǒng)一的標(biāo)準(zhǔn)公里格網(wǎng)框架的1 km柵格底圖,再將柵格轉(zhuǎn)為面要素分別與鐵路用地、公路用地、城鎮(zhèn)村道路用地、軌道交通用地4個(gè)圖層的骨架線相交獲取長度屬性,對各圖層骨架線長度按格網(wǎng)id分組求和,計(jì)算網(wǎng)格內(nèi)骨架線總長度。生成鐵路網(wǎng)密度、公路網(wǎng)密度、城鎮(zhèn)村道路網(wǎng)密度及軌道交通網(wǎng)密度柵格圖后,通過柵格計(jì)算器對以上4幅柵格圖相加求和,再按照第三次全國國土調(diào)查2019年行政區(qū)邊界裁剪柵格,生成最終的全國1 km\(×\)1 km格網(wǎng)路網(wǎng)密度數(shù)據(jù)集。
(4)產(chǎn)品質(zhì)檢
對路網(wǎng)密度產(chǎn)品進(jìn)行質(zhì)量檢查,檢查產(chǎn)品空間完整性、數(shù)據(jù)值域、影像位數(shù),投影參數(shù),以及與全國標(biāo)準(zhǔn)格網(wǎng)空間吻合程度。質(zhì)檢無誤后,利用相關(guān)軟件進(jìn)行產(chǎn)品制圖,編制產(chǎn)品數(shù)據(jù)集文檔。