簡介滹沱河流域 shp 面文件是一份可直接用于 ArcGIS 等地理信息軟件的矢量數(shù)據(jù)包主要反映流域范圍與邊界特征。資源面向水文學、地理學、城鄉(xiāng)規(guī)劃以及環(huán)境保護工作者解決這類使用者缺少基礎流域底圖、需要手工勾畫或采集矢量邊界的問題。壓縮包共 8 個文件完整覆蓋 Shapefile 格式的主要組成部分.shp 存儲面要素幾何坐標.dbf 保存屬性信息.prj 記錄投影坐標系.sbn、.sbx、.shx 等承擔空間索引與幾何加速讀取功能還附有 .xml 元數(shù)據(jù)和 .cpg 編碼聲明包體僅 36KB數(shù)據(jù)輕便規(guī)范。已有 396 人學習瀏覽說明該資源對入門級 GIS 用戶有一定參考價值。加載后可快速查看滹沱河流域邊界也可疊加土地利用、氣象、水文觀測等數(shù)據(jù)進行空間分析或作為底圖制作專題地圖同時可直接用于流域面要素的裁剪、疊加、提取等常見操作適合課堂教學、區(qū)域研究和項目演示等多種用途。1. 滹沱河流域shp格式面文件想做流域分析先從這張“面”說起滹沱河流域shp格式面文件是水利、環(huán)保和規(guī)劃項目里非常常見的一份基礎底圖。它的核心不是那幾條河道線而是一張完整的、閉合的流域邊界多邊形你可以在ArcGIS里直接疊加到DEM、行政區(qū)、土地利用或者遙感影像上。拿到這張面文件后最直接的用處就是圈定研究范圍算流域面積、統(tǒng)計高程分布、跟縣界做相交分析、裁出子流域這些操作都建立在“面SHP能被正確處理”的前提下。這篇筆記我會按自己拆數(shù)據(jù)的習慣來寫先說怎么在ArcGIS里把它加載成一份“能信任”的圖層再講裁剪、疊置、格式轉(zhuǎn)換、修復以及真正容易翻車的幾個地方。適合剛接觸SHP的水利從業(yè)者也適合被各種流域邊界數(shù)據(jù)折磨過的老手看看有沒有哪條坑你沒踩過。2. 把SHP面文件導入ArcGIS加載前的三件套檢查和坐標識別2.1 打開前先確認是“三件套”不是單個文件很多新手拿到SHP文件后習慣只看后綴帶.shp的那個文件雙擊或者拖進ArcMap結(jié)果提示“不能打開”第一反應是數(shù)據(jù)壞了。這里要先補一個基礎認知Shapefile不是一個文件而是一組文件。最少要有.shp幾何信息、.shx索引、.dbf屬性表這三個少任何一個都打不開或丟屬性。你在微信、QQ或者網(wǎng)盤里接收這種文件時如果對方只發(fā)了一個.shp那基本等于沒發(fā)。我一般會把SHP文件單獨放在一個文件夾里比如D:\hutuo_data\然后通過ArcCatalog連接到這個文件夾先在“預覽”頁簽里看一眼圖形確認它是不是面文件。為什么不推薦直接在ArcMap里拖因為ArcMap的圖層加載失敗時報錯信息很籠統(tǒng)你在里面分辨不出是文件缺失、坐標系統(tǒng)未知還是幾何本身有問題。而ArcCatalog的預覽視圖可以快速區(qū)分“點、線、面”三種ShapeType。這一點對流域數(shù)據(jù)尤其關鍵因為有些來源會把流域邊界做成線文件把流域畫成閉合線框看起來是面積實際上只是線要素后面做面積統(tǒng)計就徹底沒法用。在ArcGIS里也可以用一小段Python來快速判斷SHP的基本情況import arcpy arcpy.env.workspace rD:\hutuo_data shp hutuo_river_basin.shp desc arcpy.Describe(shp) print(要素類型, desc.shapeType) # 輸出 Polygon 才是面文件 print(空間參考, desc.spatialReference.name) print(是否含投影, desc.spatialReference.type) # Projected 或 Geographic這段代碼用arcpy.Describe讀取SHP文件的元數(shù)據(jù)。shapeType返回Polygon說明這是面圖層后面空間分析工具才認spatialReference.name會告訴你坐標系名稱比如“WGS_1984_UTM_Zone_50N”或者“CGCS2000_3_Degree_GK_CM_114E”。如果這里顯示Unknown說明.prj文件缺失這種情況下做面積計算和投影轉(zhuǎn)換都會出問題。2.2 坐標系怎么選滹沱河流域適合用什么投影面文件的坐標系決定了你算面積、量距離、做Buffer時結(jié)果對不對。滹沱河流域大部分位于山西省東部到河北省中西部東西跨度大約在112.5°E到115.5°E之間南北約38°N到39.5°N。這個范圍用經(jīng)緯度地理坐標顯示沒問題但如果直接用SHAPE.AREA字段計算面積或者用“測量工具”去量距離你得到的數(shù)值會非常奇怪因為經(jīng)緯度不是等距投影。最穩(wěn)妥的做法是在分析前把SHP轉(zhuǎn)換到合適的高斯-克呂格投影。這個流域的經(jīng)度范圍適合用CGCS2000坐標系下的3度分帶中央經(jīng)線114°E。在ArcGIS里打開“投影”工具輸入圖層選擇“CGCS2000_3_Degree_GK_CM_114E”輸出即可。轉(zhuǎn)換完成后如果單位變成米坐標值會出現(xiàn)“3開頭”或者“4開頭”的大數(shù)那就說明投影生效了。有些場合別人會發(fā)給你“北京54”或者“西安80”坐標系的數(shù)據(jù)。遇到這種情況不要直接拿來做疊加先通過“定義投影”弄清楚原始數(shù)據(jù)到底是哪一種再借助ArcGIS自帶的轉(zhuǎn)換參數(shù)轉(zhuǎn)到CGCS2000。這里特別提醒如果SHP文件沒有.prj文件ArcGIS無法自動識別坐標系你需要先手工把正確的坐標系定義進去否則后續(xù)做動態(tài)投影會遇到圖層明明在附近卻對不齊的情況。2.3 屬性表里能挖到什么先看字段再談分析面SHP的價值不止在邊界還在屬性表。打開屬性表后常見字段有OBJECTID、Shape_Leng、Shape_Area以及河流名稱、流域編碼、所屬水系這類業(yè)務字段。不同的制圖單位給出的字段命名差異很大有的叫NAME有的叫XZQMC。你在做字段關聯(lián)或制圖標注前必須先確認字段名的含義不要拿OBJECTID去關聯(lián)外部Excel表那幾乎一定會錯位。下面是一個典型的流域?qū)傩员斫Y(jié)構(gòu)示例字段名類型說明OBJECTID長整型系統(tǒng)唯一IDShape_Area雙精度投影平面面積NAME文本河流或流域名稱BASIN_CODE文本流域編碼SOURCE文本數(shù)據(jù)來源單位注意Shape_Area字段的值會隨當前圖層坐標系變化而變。如果圖層是地理坐標這個字段的值不是平方米只有投影坐標下它才接近真實的平方米數(shù)值。所以我的習慣是在開始任何分析前給SHP增加一個自己的面積字段用Geometry Calculator計算一次正式面積后面不管坐標系怎么折騰都以這個字段為準。3. 把流域SHP變成能分析的圖層裁剪、分區(qū)統(tǒng)計和關聯(lián)屬性3.1 用Clip工具提取子流域或局部范圍拿到完整的滹沱河流域面SHP后你可能只需要滹沱河石家莊段的匯水范圍或者只想保留某個水文站以上的區(qū)域。最常用的操作就是用另一個面圖層去做裁剪即ArcToolbox里的“分析工具→提取→裁剪”。裁剪相當于拿一個“餅干模子”去摳原始面輸出結(jié)果只保留模子內(nèi)的部分幾何邊界自動更新。操作上需要給“裁剪要素”提供一個閉合面這個面可以是你手工繪制的矩形也可以是按字段篩選后的子集。如果你已經(jīng)有子流域邊界直接用那個子流域面去裁剪如果沒有可以先用“按位置選擇”把主流域SHP里跟某條河流相交的面選中然后導出新圖層。這里要注意很多人喜歡先“合并”再裁剪但如果數(shù)據(jù)量大合并會帶來不必要的性能損耗和幾何重疊問題慎用。用ArcPy也能完成同樣的裁剪操作方便批量跑多個流域import arcpy arcpy.env.workspace rD:\hutuo_data arcpy.analysis.Clip( in_featureshutuo_river_basin.shp, clip_featuressjz_section.shp, out_feature_classhutuo_sjz_clip.shp ) print(裁剪完成)這里的in_features是被裁剪的流域面clip_features是裁剪范圍out_feature_class是輸出文件名。ArcGIS的Clip對輸入和裁剪要素必須是同類幾何面裁面沒有問題但如果你拿一個線圖層去裁剪面結(jié)果會什么也不剩。報錯時先檢查這兩個參數(shù)的幾何類型是否匹配這是最常見的翻車原因。3.2 讓流域SHP和DEM說上話分區(qū)統(tǒng)計面SHP放入ArcGIS后單獨看邊界意義有限真正能出成果的是把流域范圍作為掩膜去統(tǒng)計它覆蓋的柵格數(shù)據(jù)。比如你有滹沱河流域的30米DEM想知道流域內(nèi)平均高程、最高點、最低點不必手動描范圍直接用空間分析工具里的“分區(qū)統(tǒng)計”即可。分區(qū)統(tǒng)計的核心邏輯是把面SHP的每個要素當作一個“區(qū)”對落入該區(qū)內(nèi)的所有柵格像元做統(tǒng)計輸出一個dbf表格。這個工具的輸入?yún)?shù)有兩個關鍵點柵格數(shù)據(jù)或要素區(qū)域數(shù)據(jù)選擇SHP面文件區(qū)域字段選一個唯一ID字段。下面這段腳本演示批量場景import arcpy from arcpy.sa import * arcpy.env.workspace rD:\hutuo_data arcpy.env.extent hutuo_river_basin.shp arcpy.env.snapRaster dem30.tif outTable ZonalStatisticsAsTable( hutuo_river_basin.shp, FID, dem30.tif, basin_zonal.dbf, DATA, ALL ) print(分區(qū)統(tǒng)計結(jié)果, outTable)FID是原始SHP里的唯一編號用它做區(qū)域字段最保險因為不會出現(xiàn)同名要素導致統(tǒng)計結(jié)果被合并的情況。DATA參數(shù)表示忽略柵格中的NoData值ALL表示輸出平均值、最大值、最小值、標準差等全部統(tǒng)計項。統(tǒng)計完成后把輸出的basin_zonal.dbf通過公共字段連接回原SHP屬性表就能在圖層上按平均高程做符號化。3.3 給面文件補屬性計算面積和手工新建SHP有時候你拿到的SHP屬性表非常簡陋沒有面積字段甚至連坐標信息都不完整。我的習慣是先給原文件做一個備份然后添加一個雙精度字段專門用來記錄正式面積。這里有一個容易踩的坑如果你直接在SHP上添加字段并計算可能會把原始數(shù)據(jù)改壞所以更規(guī)范的做法是先復制一份再在副本上操作。復制SHP并計算面積可以用Python腳本import arcpy src rD:\hutuo_data\hutuo_river_basin.shp dst rD:\hutuo_data\hutuo_area.dbf arcpy.Copy_management(src, rD:\hutuo_data\hutuo_copy.shp, ) arcpy.AddField_management(dst, AREA_KM2, DOUBLE, 10, 2) arcpy.CalculateField_management( dst, AREA_KM2, !SHAPE.GEODESICAREA! / 1000000, PYTHON3 )SHAPE.GEODESICAREA是ArcGIS的Python表達式表示要素在橢球體上計算的測地面積單位是平方米除以1000000后得到平方千米。這里之所以不用!Shape_Area!是因為Shape_Area可能受當前坐標系影響當圖層處于經(jīng)緯度坐標時Shape_Area的數(shù)值會失真而GEODESICAREA只跟要素的真實幾何位置相關不受地圖投影影響。如果你工作區(qū)域跨度不大兩者差異不明顯如果做跨帶的大流域分析這個細節(jié)就可能決定數(shù)據(jù)是否可信。如果你沒有現(xiàn)成的SHP想從零開始新建一個流域面文件在ArcGIS里用“創(chuàng)建要素類”可以完成選擇面幾何類型指定坐標系然后打開編輯器繪制邊界。但手工繪制流域邊界精度很低我一般不推薦除非你只是做一個示意草圖。真實項目里還是以水文部門提取的流域SHP為準手工繪制僅用于快速演示。4. 格式轉(zhuǎn)換與修復DWG轉(zhuǎn)SHP、SHP轉(zhuǎn)KML和ShapeCHK修復實操4.1 DWG轉(zhuǎn)SHPCAD線稿要變成面先解決閉合問題很多歷史資料里滹沱河的河道和流域邊界都是DWG格式的CAD圖紙。要把DWG轉(zhuǎn)成SHP面文件不是簡單地“另存為”就能完成。CAD里的流域邊界通常是多段線如果這條多段線沒有完全閉合ArcGIS轉(zhuǎn)換后只會生成線要素不會生成面要素。這個坑在我過去處理河道邊界時幾乎每次都遇到。常見做法分兩步第一步把DWG導入ArcGIS用“CAD到地理數(shù)據(jù)庫”工具轉(zhuǎn)換成要素類第二步打開屬性表查看要素的幾何類型。如果導入后是Polyline需要用“要素轉(zhuǎn)面”工具把閉合線轉(zhuǎn)換成面。但“要素轉(zhuǎn)面”有個前提線要素不能有懸掛點和自相交否則生成的面會破碎。在CAD里先處理好數(shù)據(jù)比ArcGIS里補救高效得多。我會在CAD中用BOUNDARY命令對閉合區(qū)域重新生成面域或者用PEDIT命令檢查多段線是否閉合。確保每條邊界終點和起點重合后再導入ArcGIS。這個經(jīng)驗幾乎是處理“dwg轉(zhuǎn)shp”唯一可靠的路徑。下面是一個參考操作流程1. CAD中打開原圖凍結(jié)無關圖層只保留流域邊界線。 2. 全選邊界線執(zhí)行 PEDIT - 輸入 J 合并 - 閉合檢查。 3. CAD另存為 R2010 以下版本避免高版本不兼容。 4. ArcToolbox - 轉(zhuǎn)換工具 - CAD轉(zhuǎn)地理數(shù)據(jù)庫。 5. 對生成的線要素執(zhí)行要素轉(zhuǎn)面。 6. 用第七章提到的 ShapeCHK 檢查面文件。4.2 SHP轉(zhuǎn)KML把流域邊界拿到Google Earth里看很多非GIS專業(yè)的同事習慣把SHP轉(zhuǎn)成KML放到Google Earth里查看。ArcGIS里最常用的工具是“圖層轉(zhuǎn)KML”但SHP文件不能直接作為輸入它需要先有一個圖層文件或者地圖文檔。所以我一般先給SHP做一個簡單的符號化然后在內(nèi)容列表里右鍵圖層選擇“屬性→圖層屬性→圖層轉(zhuǎn)KML”。如果寫成Python邏輯是這樣import arcpy arcpy.MakeFeatureLayer_management( rD:\hutuo_data\hutuo_river_basin.shp, hutuo_lyr ) arcpy.LayerToKML_conversion( hutuo_lyr, rD:\hutuo_data\hutuo_basin.kmz, MAP_BOUNDARY, clampedToGround, COMPOSITE, 1024, 96, PNG8 )MAP_BOUNDARY表示只導出圖層可見范圍clampedToGround表示把要素貼在地面上而不是懸浮在某個高度。輸出文件后綴可以用.kmz壓縮包Google Earth能直接讀取。需要注意圖層里如果有中文標注導出前字體要選擇支持中文的字體否則KML在手機上打開會出現(xiàn)方塊字。反過來KML轉(zhuǎn)SHP也很常見。用“KML轉(zhuǎn)圖層”工具導入后輸出的是一個要素圖層但KML里的面可能被拆成了多個片段原因是Google Earth里的多邊形如果中間有鏤空或者復雜邊界轉(zhuǎn)成SHP后會生成多部件面。多部件面雖然也是合法SHP但做面積統(tǒng)計和拓撲檢查時容易忽略內(nèi)部孔洞。處理這類數(shù)據(jù)時我一般先用“多部件轉(zhuǎn)單部件”把幾何拆開再逐塊檢查。4.3 修復損壞的SHP用ShapeCHK檢查幾何和索引Shapefile是上世紀九十年代的老格式結(jié)構(gòu)簡單但也脆弱。文件斷電、非法關機、U盤拷貝中斷都可能讓SHP產(chǎn)生幾何錯誤、.shx索引和.shp內(nèi)容不匹配。ArcGIS打開時報“創(chuàng)建要素類失敗”或者“無法識別文件”大概率是SHP本身損壞了。這里強烈建議使用ShapeCHK這個老牌命令行工具。它是Esri圈子里流傳很廣的SHP診斷修復工具專門檢查.shp、.shx、.dbf三者一致性、幾何記錄長度、記錄頭是否完整。使用方法也不復雜打開命令行窗口把ShapeCHK和要修復的SHP放在同一目錄然后執(zhí)行shapechk hutuo_river_basin.shp工具會逐條掃描文件中的幾何記錄發(fā)現(xiàn)錯誤時輸出錯誤類型和大致位置。部分版本支持參數(shù)自動修復shapechk hutuo_river_basin.shp /fix/fix參數(shù)會嘗試重建索引文件并修正頭部信息但如果幾何本身已經(jīng)破損修復程度有限。所以我通常把ShapeCHK當作“體檢工具”而不是“萬能后悔藥”跑一遍沒報錯說明文件基礎結(jié)構(gòu)沒問題報了錯趕緊找原始備份不要抱著修復后的文件繼續(xù)做分析。順帶一提如果你想把SHP屬性表轉(zhuǎn)成txt用于核對可以在ArcGIS里右鍵屬性表選擇“導出”把dbf內(nèi)容另存為CSV或文本文件。CSV文件用Excel打開時字段編碼經(jīng)常亂我會選擇用VSCode或Notepad打開并切換到UTF-8編碼查看。5. 實操避坑滹沱河流域SHP處理時的五個高頻問題5.1 計算面積偏差大數(shù)值小得不像流域現(xiàn)象面SHP明明是完整流域用面積字段統(tǒng)計結(jié)果只有幾千平方米跟真實流域面積差了至少三個數(shù)量級。原因圖層坐標系處于經(jīng)緯度地理坐標狀態(tài)SHAPE.AREA按度計算不是平方米。在沒有投影的情況下Shape_Area字段不具備實際物理意義。解決先確認.prj內(nèi)容通過“投影”工具把SHP轉(zhuǎn)到適合滹沱河流域的高斯-克呂格投影。滹沱河流域我常用CGCS2000_3_Degree_GK_CM_114E。轉(zhuǎn)換完成后重新計算SHAPE.AREA得到的數(shù)值才接近真實平方米。如果不想動原數(shù)據(jù)計算字段時用!SHAPE.GEODESICAREA!也能繞開投影問題。5.2 面重疊和自相交做裁剪時報錯或結(jié)果為空現(xiàn)象用流域SHP去裁別的數(shù)據(jù)系統(tǒng)提示“幾何未正確初始化”或者輸出空圖層打開屬性表發(fā)現(xiàn)要素數(shù)為0。原因SHP里的多邊形存在自相交、重復節(jié)點、懸空邊界等拓撲問題。面文件在Web端下載或者從CAD轉(zhuǎn)換時沒有經(jīng)過拓撲檢查這樣的幾何在ArcGIS的許多分析工具中會被判定為不可信。解決用“數(shù)據(jù)管理工具→要素→修復幾何”工具跑一遍。修復幾何時勾選“刪除重復節(jié)點”和“修復自相交”。修復前先備份因為這類工具對數(shù)據(jù)有侵入性。如果修復后還有問題用ArcGIS的“拓撲”工具集建一個拓撲規(guī)則指定“不能重疊”再逐個檢查。5.3 屬性連接總是匹配不上中文名稱對不上號現(xiàn)象把Excel統(tǒng)計表按鄉(xiāng)鎮(zhèn)名稱連接到流域SHP屬性表結(jié)果大量記錄顯示“空”明明兩個表里都有同一個鄉(xiāng)鎮(zhèn)名。原因最常見的是兩個表里鄉(xiāng)鎮(zhèn)名存在不可見字符Excel從其他系統(tǒng)導出時字符串字段前面或后面帶著空格。另一個原因是SHP的.dbf文件編碼是GBK而Excel另存為UTF-8連接時字符不一致。解決連接前先用“表轉(zhuǎn)Excel”導出SHP屬性表然后跟外部Excel對比特別注意字段長度和空格。Excel里用TRIM()函數(shù)清理兩邊的空格再用VLOOKUP人工核對一遍。對編碼問題檢查SHP目錄下是否有.cpg文件如果沒有用QGIS打開原SHP并另存為UTF-8編碼的新SHP再回到ArcGIS里做連接。5.4 拷貝到別人電腦上打不開提示文件缺失現(xiàn)象把SHP從自己電腦傳到同事電腦對方說圖層能加載但屬性表空白或者干脆報“打開SHP文件失敗”。原因只拷貝了.shp文件.shx和.dbf沒跟上。更隱蔽的是三個文件名被改了其中一部分比如.shp叫“hutuo.shp”.shx卻叫“hutuo_new.shx”引擎按文件名匹配不到索引文件。解決傳輸SHP時把文件夾整體壓縮打包不要單獨拖一個文件。如果已經(jīng)有散亂的三件套改用“ArcCatalog→復制粘貼”方式重新生成一套確保命名統(tǒng)一。我的習慣是每次發(fā)送前用ShapeCHK檢查一遍確認三個文件齊全再發(fā)。5.5 省界線文件省1和省2的SHP有什么區(qū)別疊加時對不齊現(xiàn)象拿到國家基礎地理數(shù)據(jù)里的省1和省2兩個SHP加載到同一地圖中發(fā)現(xiàn)滹沱河流域附近的省界線條并不完全重合做疊加分析后出現(xiàn)很窄的縫隙或重疊。原因省1和省2是不同生產(chǎn)精度的界線數(shù)據(jù)。省1一般指經(jīng)過概化、用于中小比例尺制圖的省級行政區(qū)界線省2是更精細的邊界包含更復雜的細節(jié)和國界要素。兩者比例尺和拓撲規(guī)則不同幾何自然不完全一致。解決做流域空間分析時只選一套省界作為底圖不要同時疊加兩套并直接做Union。要以省2即精度更細的版本作為參考邊界省1可以用來做示意但不要參與拓撲計算。如果你發(fā)現(xiàn)流域SHP邊界和省界有沖突優(yōu)先以流域?qū)嶋H地貌范圍為準行政界線只是參考圖層。6. 再用深一點漁網(wǎng)分割、邊界簡化與三維場景過渡6.1 用“創(chuàng)建漁網(wǎng)”把大流域切成分塊有時你不需要整個滹沱河流域作為一個面做分析而需要把它劃分成2公里×2公里的網(wǎng)格逐格統(tǒng)計土地利用類型或者做水環(huán)境評價。手動畫方格不現(xiàn)實正確做法是用ArcGIS的“創(chuàng)建漁網(wǎng)”工具。漁網(wǎng)工具的本質(zhì)是生成一個規(guī)則的格網(wǎng)覆蓋指定范圍。關鍵參數(shù)有三個漁網(wǎng)原點、像元寬高和模板范圍。模板范圍直接選擇SHP面文件時漁網(wǎng)輸出邊界會跟流域的范圍框一致。下圖場景里我創(chuàng)建2公里格網(wǎng)然后與流域SHP疊加保留與流域相交的網(wǎng)格import arcpy arcpy.env.workspace rD:\hutuo_data arcpy.management.CreateFishnet( out_featuresgrid_2km.shp, origin_coord490000 4250000, y_axis_coord490000 4250010, cell_width2000, cell_height2000, number_rowsNone, number_columnsNone, labelsNO_LABELS, templatehutuo_river_basin.shp, geometry_typePOLYGON )參數(shù)origin_coord是漁網(wǎng)左下角的投影坐標y_axis_coord決定網(wǎng)格旋轉(zhuǎn)方向通常讓它與正北方向一致。template填流域SHP生成的漁網(wǎng)會自適應流域的最大外接矩形。生成網(wǎng)格后還要用“相交”工具把漁網(wǎng)和流域SHP相交這樣網(wǎng)格邊緣就貼合流域邊界而不是簡單的矩形外框。這一步如果不做你統(tǒng)計的格網(wǎng)會有一大半在流域外結(jié)果毫無意義。6.2 面邊界簡化與三維場景過渡流域SHP從專業(yè)水文機構(gòu)下載時邊界往往包含海量節(jié)點一個流域面可能帶著幾萬個點。這種精細數(shù)據(jù)在小比例尺出圖時不僅讓電腦卡頓還讓ArcMap的符號化變得遲鈍。我的經(jīng)驗是保存兩套數(shù)據(jù)一套原始精度用于分析一套簡化版用于制圖和發(fā)布。ArcGIS有條“簡化面”工具位于制圖工具箱中import arcpy arcpy.cartography.SimplifyPolygon( in_featureshutuo_river_basin.shp, out_feature_classbasin_simplify.shp, algorithmPOINT_REMOVE, tolerance50 Meters, minimum_area0 SquareMeters )POINT_REMOVE算法通過刪除與相鄰點夾角接近直線的點來降低節(jié)點密度50 Meters的容差表示小于50米的邊界彎曲會被忽略。簡化閾值不是越大越好如果設成500米流域邊界的細小河道特征會全部丟失。出圖給領導看容差可以放寬做水文分析保持原始精度不要隨意簡化。如果你準備把流域SHP推送到三維場景比如生成3DTiles做Web端可視化水汽效果會清楚反映出模型對流域面的支撐力?；驹硎前衙鍿HP先轉(zhuǎn)成GeoJSON再疊加上高程數(shù)據(jù)最后切片成b3dm。ArcGIS里可以用“轉(zhuǎn)出至GeoJSON”工具后續(xù)步驟交給開源工具完成。注意SHP屬性表里的中文名稱在GeoJSON里可能編碼錯亂轉(zhuǎn)之前先在屬性表里增加一個英文字段。我從實際項目里積累的習慣是拿到任何SHP第一件事復制一個原版?zhèn)浞萑缓笥肧hapeCHK檢查文件結(jié)構(gòu)用arcpy.Describe確認坐標系和幾何類型再檢查屬性表最后才開始裁剪和轉(zhuǎn)換。這套流程看起來繁瑣但確實幫我避免了很多在交付前一刻才發(fā)現(xiàn)幾何錯誤的窘境。尤其是滹沱河流域這類邊界跨度較大的面文件坐標系和拓撲問題幾乎是必然存在的早檢查早安心。希望這篇筆記里的“三件套檢查投影確認幾何修復”能幫到你至少下次再遇到“面文件算不出面積”的怪事時你知道該往哪幾個方向排查。本文還有配套的精品資源點擊獲取 滹沱河流域SHP面文件處理全指南：從ArcGIS加載到修復避坑
簡介滹沱河流域 shp 面文件是一份可直接用于 ArcGIS 等地理信息軟件的矢量數(shù)據(jù)包主要反映流域范圍與邊界特征。資源面向水文學、地理學、城鄉(xiāng)規(guī)劃以及環(huán)境保護工作者解決這類使用者缺少基礎流域底圖、需要手工勾畫或采集矢量邊界的問題。壓縮包共 8 個文件完整覆蓋 Shapefile 格式的主要組成部分.shp 存儲面要素幾何坐標.dbf 保存屬性信息.prj 記錄投影坐標系.sbn、.sbx、.shx 等承擔空間索引與幾何加速讀取功能還附有 .xml 元數(shù)據(jù)和 .cpg 編碼聲明包體僅 36KB數(shù)據(jù)輕便規(guī)范。已有 396 人學習瀏覽說明該資源對入門級 GIS 用戶有一定參考價值。加載后可快速查看滹沱河流域邊界也可疊加土地利用、氣象、水文觀測等數(shù)據(jù)進行空間分析或作為底圖制作專題地圖同時可直接用于流域面要素的裁剪、疊加、提取等常見操作適合課堂教學、區(qū)域研究和項目演示等多種用途。1. 滹沱河流域shp格式面文件想做流域分析先從這張“面”說起滹沱河流域shp格式面文件是水利、環(huán)保和規(guī)劃項目里非常常見的一份基礎底圖。它的核心不是那幾條河道線而是一張完整的、閉合的流域邊界多邊形你可以在ArcGIS里直接疊加到DEM、行政區(qū)、土地利用或者遙感影像上。拿到這張面文件后最直接的用處就是圈定研究范圍算流域面積、統(tǒng)計高程分布、跟縣界做相交分析、裁出子流域這些操作都建立在“面SHP能被正確處理”的前提下。這篇筆記我會按自己拆數(shù)據(jù)的習慣來寫先說怎么在ArcGIS里把它加載成一份“能信任”的圖層再講裁剪、疊置、格式轉(zhuǎn)換、修復以及真正容易翻車的幾個地方。適合剛接觸SHP的水利從業(yè)者也適合被各種流域邊界數(shù)據(jù)折磨過的老手看看有沒有哪條坑你沒踩過。2. 把SHP面文件導入ArcGIS加載前的三件套檢查和坐標識別2.1 打開前先確認是“三件套”不是單個文件很多新手拿到SHP文件后習慣只看后綴帶.shp的那個文件雙擊或者拖進ArcMap結(jié)果提示“不能打開”第一反應是數(shù)據(jù)壞了。這里要先補一個基礎認知Shapefile不是一個文件而是一組文件。最少要有.shp幾何信息、.shx索引、.dbf屬性表這三個少任何一個都打不開或丟屬性。你在微信、QQ或者網(wǎng)盤里接收這種文件時如果對方只發(fā)了一個.shp那基本等于沒發(fā)。我一般會把SHP文件單獨放在一個文件夾里比如D:\hutuo_data\然后通過ArcCatalog連接到這個文件夾先在“預覽”頁簽里看一眼圖形確認它是不是面文件。為什么不推薦直接在ArcMap里拖因為ArcMap的圖層加載失敗時報錯信息很籠統(tǒng)你在里面分辨不出是文件缺失、坐標系統(tǒng)未知還是幾何本身有問題。而ArcCatalog的預覽視圖可以快速區(qū)分“點、線、面”三種ShapeType。這一點對流域數(shù)據(jù)尤其關鍵因為有些來源會把流域邊界做成線文件把流域畫成閉合線框看起來是面積實際上只是線要素后面做面積統(tǒng)計就徹底沒法用。在ArcGIS里也可以用一小段Python來快速判斷SHP的基本情況import arcpy arcpy.env.workspace rD:\hutuo_data shp hutuo_river_basin.shp desc arcpy.Describe(shp) print(要素類型, desc.shapeType) # 輸出 Polygon 才是面文件 print(空間參考, desc.spatialReference.name) print(是否含投影, desc.spatialReference.type) # Projected 或 Geographic這段代碼用arcpy.Describe讀取SHP文件的元數(shù)據(jù)。shapeType返回Polygon說明這是面圖層后面空間分析工具才認spatialReference.name會告訴你坐標系名稱比如“WGS_1984_UTM_Zone_50N”或者“CGCS2000_3_Degree_GK_CM_114E”。如果這里顯示Unknown說明.prj文件缺失這種情況下做面積計算和投影轉(zhuǎn)換都會出問題。2.2 坐標系怎么選滹沱河流域適合用什么投影面文件的坐標系決定了你算面積、量距離、做Buffer時結(jié)果對不對。滹沱河流域大部分位于山西省東部到河北省中西部東西跨度大約在112.5°E到115.5°E之間南北約38°N到39.5°N。這個范圍用經(jīng)緯度地理坐標顯示沒問題但如果直接用SHAPE.AREA字段計算面積或者用“測量工具”去量距離你得到的數(shù)值會非常奇怪因為經(jīng)緯度不是等距投影。最穩(wěn)妥的做法是在分析前把SHP轉(zhuǎn)換到合適的高斯-克呂格投影。這個流域的經(jīng)度范圍適合用CGCS2000坐標系下的3度分帶中央經(jīng)線114°E。在ArcGIS里打開“投影”工具輸入圖層選擇“CGCS2000_3_Degree_GK_CM_114E”輸出即可。轉(zhuǎn)換完成后如果單位變成米坐標值會出現(xiàn)“3開頭”或者“4開頭”的大數(shù)那就說明投影生效了。有些場合別人會發(fā)給你“北京54”或者“西安80”坐標系的數(shù)據(jù)。遇到這種情況不要直接拿來做疊加先通過“定義投影”弄清楚原始數(shù)據(jù)到底是哪一種再借助ArcGIS自帶的轉(zhuǎn)換參數(shù)轉(zhuǎn)到CGCS2000。這里特別提醒如果SHP文件沒有.prj文件ArcGIS無法自動識別坐標系你需要先手工把正確的坐標系定義進去否則后續(xù)做動態(tài)投影會遇到圖層明明在附近卻對不齊的情況。2.3 屬性表里能挖到什么先看字段再談分析面SHP的價值不止在邊界還在屬性表。打開屬性表后常見字段有OBJECTID、Shape_Leng、Shape_Area以及河流名稱、流域編碼、所屬水系這類業(yè)務字段。不同的制圖單位給出的字段命名差異很大有的叫NAME有的叫XZQMC。你在做字段關聯(lián)或制圖標注前必須先確認字段名的含義不要拿OBJECTID去關聯(lián)外部Excel表那幾乎一定會錯位。下面是一個典型的流域?qū)傩员斫Y(jié)構(gòu)示例字段名類型說明OBJECTID長整型系統(tǒng)唯一IDShape_Area雙精度投影平面面積NAME文本河流或流域名稱BASIN_CODE文本流域編碼SOURCE文本數(shù)據(jù)來源單位注意Shape_Area字段的值會隨當前圖層坐標系變化而變。如果圖層是地理坐標這個字段的值不是平方米只有投影坐標下它才接近真實的平方米數(shù)值。所以我的習慣是在開始任何分析前給SHP增加一個自己的面積字段用Geometry Calculator計算一次正式面積后面不管坐標系怎么折騰都以這個字段為準。3. 把流域SHP變成能分析的圖層裁剪、分區(qū)統(tǒng)計和關聯(lián)屬性3.1 用Clip工具提取子流域或局部范圍拿到完整的滹沱河流域面SHP后你可能只需要滹沱河石家莊段的匯水范圍或者只想保留某個水文站以上的區(qū)域。最常用的操作就是用另一個面圖層去做裁剪即ArcToolbox里的“分析工具→提取→裁剪”。裁剪相當于拿一個“餅干模子”去摳原始面輸出結(jié)果只保留模子內(nèi)的部分幾何邊界自動更新。操作上需要給“裁剪要素”提供一個閉合面這個面可以是你手工繪制的矩形也可以是按字段篩選后的子集。如果你已經(jīng)有子流域邊界直接用那個子流域面去裁剪如果沒有可以先用“按位置選擇”把主流域SHP里跟某條河流相交的面選中然后導出新圖層。這里要注意很多人喜歡先“合并”再裁剪但如果數(shù)據(jù)量大合并會帶來不必要的性能損耗和幾何重疊問題慎用。用ArcPy也能完成同樣的裁剪操作方便批量跑多個流域import arcpy arcpy.env.workspace rD:\hutuo_data arcpy.analysis.Clip( in_featureshutuo_river_basin.shp, clip_featuressjz_section.shp, out_feature_classhutuo_sjz_clip.shp ) print(裁剪完成)這里的in_features是被裁剪的流域面clip_features是裁剪范圍out_feature_class是輸出文件名。ArcGIS的Clip對輸入和裁剪要素必須是同類幾何面裁面沒有問題但如果你拿一個線圖層去裁剪面結(jié)果會什么也不剩。報錯時先檢查這兩個參數(shù)的幾何類型是否匹配這是最常見的翻車原因。3.2 讓流域SHP和DEM說上話分區(qū)統(tǒng)計面SHP放入ArcGIS后單獨看邊界意義有限真正能出成果的是把流域范圍作為掩膜去統(tǒng)計它覆蓋的柵格數(shù)據(jù)。比如你有滹沱河流域的30米DEM想知道流域內(nèi)平均高程、最高點、最低點不必手動描范圍直接用空間分析工具里的“分區(qū)統(tǒng)計”即可。分區(qū)統(tǒng)計的核心邏輯是把面SHP的每個要素當作一個“區(qū)”對落入該區(qū)內(nèi)的所有柵格像元做統(tǒng)計輸出一個dbf表格。這個工具的輸入?yún)?shù)有兩個關鍵點柵格數(shù)據(jù)或要素區(qū)域數(shù)據(jù)選擇SHP面文件區(qū)域字段選一個唯一ID字段。下面這段腳本演示批量場景import arcpy from arcpy.sa import * arcpy.env.workspace rD:\hutuo_data arcpy.env.extent hutuo_river_basin.shp arcpy.env.snapRaster dem30.tif outTable ZonalStatisticsAsTable( hutuo_river_basin.shp, FID, dem30.tif, basin_zonal.dbf, DATA, ALL ) print(分區(qū)統(tǒng)計結(jié)果, outTable)FID是原始SHP里的唯一編號用它做區(qū)域字段最保險因為不會出現(xiàn)同名要素導致統(tǒng)計結(jié)果被合并的情況。DATA參數(shù)表示忽略柵格中的NoData值ALL表示輸出平均值、最大值、最小值、標準差等全部統(tǒng)計項。統(tǒng)計完成后把輸出的basin_zonal.dbf通過公共字段連接回原SHP屬性表就能在圖層上按平均高程做符號化。3.3 給面文件補屬性計算面積和手工新建SHP有時候你拿到的SHP屬性表非常簡陋沒有面積字段甚至連坐標信息都不完整。我的習慣是先給原文件做一個備份然后添加一個雙精度字段專門用來記錄正式面積。這里有一個容易踩的坑如果你直接在SHP上添加字段并計算可能會把原始數(shù)據(jù)改壞所以更規(guī)范的做法是先復制一份再在副本上操作。復制SHP并計算面積可以用Python腳本import arcpy src rD:\hutuo_data\hutuo_river_basin.shp dst rD:\hutuo_data\hutuo_area.dbf arcpy.Copy_management(src, rD:\hutuo_data\hutuo_copy.shp, ) arcpy.AddField_management(dst, AREA_KM2, DOUBLE, 10, 2) arcpy.CalculateField_management( dst, AREA_KM2, !SHAPE.GEODESICAREA! / 1000000, PYTHON3 )SHAPE.GEODESICAREA是ArcGIS的Python表達式表示要素在橢球體上計算的測地面積單位是平方米除以1000000后得到平方千米。這里之所以不用!Shape_Area!是因為Shape_Area可能受當前坐標系影響當圖層處于經(jīng)緯度坐標時Shape_Area的數(shù)值會失真而GEODESICAREA只跟要素的真實幾何位置相關不受地圖投影影響。如果你工作區(qū)域跨度不大兩者差異不明顯如果做跨帶的大流域分析這個細節(jié)就可能決定數(shù)據(jù)是否可信。如果你沒有現(xiàn)成的SHP想從零開始新建一個流域面文件在ArcGIS里用“創(chuàng)建要素類”可以完成選擇面幾何類型指定坐標系然后打開編輯器繪制邊界。但手工繪制流域邊界精度很低我一般不推薦除非你只是做一個示意草圖。真實項目里還是以水文部門提取的流域SHP為準手工繪制僅用于快速演示。4. 格式轉(zhuǎn)換與修復DWG轉(zhuǎn)SHP、SHP轉(zhuǎn)KML和ShapeCHK修復實操4.1 DWG轉(zhuǎn)SHPCAD線稿要變成面先解決閉合問題很多歷史資料里滹沱河的河道和流域邊界都是DWG格式的CAD圖紙。要把DWG轉(zhuǎn)成SHP面文件不是簡單地“另存為”就能完成。CAD里的流域邊界通常是多段線如果這條多段線沒有完全閉合ArcGIS轉(zhuǎn)換后只會生成線要素不會生成面要素。這個坑在我過去處理河道邊界時幾乎每次都遇到。常見做法分兩步第一步把DWG導入ArcGIS用“CAD到地理數(shù)據(jù)庫”工具轉(zhuǎn)換成要素類第二步打開屬性表查看要素的幾何類型。如果導入后是Polyline需要用“要素轉(zhuǎn)面”工具把閉合線轉(zhuǎn)換成面。但“要素轉(zhuǎn)面”有個前提線要素不能有懸掛點和自相交否則生成的面會破碎。在CAD里先處理好數(shù)據(jù)比ArcGIS里補救高效得多。我會在CAD中用BOUNDARY命令對閉合區(qū)域重新生成面域或者用PEDIT命令檢查多段線是否閉合。確保每條邊界終點和起點重合后再導入ArcGIS。這個經(jīng)驗幾乎是處理“dwg轉(zhuǎn)shp”唯一可靠的路徑。下面是一個參考操作流程1. CAD中打開原圖凍結(jié)無關圖層只保留流域邊界線。 2. 全選邊界線執(zhí)行 PEDIT - 輸入 J 合并 - 閉合檢查。 3. CAD另存為 R2010 以下版本避免高版本不兼容。 4. ArcToolbox - 轉(zhuǎn)換工具 - CAD轉(zhuǎn)地理數(shù)據(jù)庫。 5. 對生成的線要素執(zhí)行要素轉(zhuǎn)面。 6. 用第七章提到的 ShapeCHK 檢查面文件。4.2 SHP轉(zhuǎn)KML把流域邊界拿到Google Earth里看很多非GIS專業(yè)的同事習慣把SHP轉(zhuǎn)成KML放到Google Earth里查看。ArcGIS里最常用的工具是“圖層轉(zhuǎn)KML”但SHP文件不能直接作為輸入它需要先有一個圖層文件或者地圖文檔。所以我一般先給SHP做一個簡單的符號化然后在內(nèi)容列表里右鍵圖層選擇“屬性→圖層屬性→圖層轉(zhuǎn)KML”。如果寫成Python邏輯是這樣import arcpy arcpy.MakeFeatureLayer_management( rD:\hutuo_data\hutuo_river_basin.shp, hutuo_lyr ) arcpy.LayerToKML_conversion( hutuo_lyr, rD:\hutuo_data\hutuo_basin.kmz, MAP_BOUNDARY, clampedToGround, COMPOSITE, 1024, 96, PNG8 )MAP_BOUNDARY表示只導出圖層可見范圍clampedToGround表示把要素貼在地面上而不是懸浮在某個高度。輸出文件后綴可以用.kmz壓縮包Google Earth能直接讀取。需要注意圖層里如果有中文標注導出前字體要選擇支持中文的字體否則KML在手機上打開會出現(xiàn)方塊字。反過來KML轉(zhuǎn)SHP也很常見。用“KML轉(zhuǎn)圖層”工具導入后輸出的是一個要素圖層但KML里的面可能被拆成了多個片段原因是Google Earth里的多邊形如果中間有鏤空或者復雜邊界轉(zhuǎn)成SHP后會生成多部件面。多部件面雖然也是合法SHP但做面積統(tǒng)計和拓撲檢查時容易忽略內(nèi)部孔洞。處理這類數(shù)據(jù)時我一般先用“多部件轉(zhuǎn)單部件”把幾何拆開再逐塊檢查。4.3 修復損壞的SHP用ShapeCHK檢查幾何和索引Shapefile是上世紀九十年代的老格式結(jié)構(gòu)簡單但也脆弱。文件斷電、非法關機、U盤拷貝中斷都可能讓SHP產(chǎn)生幾何錯誤、.shx索引和.shp內(nèi)容不匹配。ArcGIS打開時報“創(chuàng)建要素類失敗”或者“無法識別文件”大概率是SHP本身損壞了。這里強烈建議使用ShapeCHK這個老牌命令行工具。它是Esri圈子里流傳很廣的SHP診斷修復工具專門檢查.shp、.shx、.dbf三者一致性、幾何記錄長度、記錄頭是否完整。使用方法也不復雜打開命令行窗口把ShapeCHK和要修復的SHP放在同一目錄然后執(zhí)行shapechk hutuo_river_basin.shp工具會逐條掃描文件中的幾何記錄發(fā)現(xiàn)錯誤時輸出錯誤類型和大致位置。部分版本支持參數(shù)自動修復shapechk hutuo_river_basin.shp /fix/fix參數(shù)會嘗試重建索引文件并修正頭部信息但如果幾何本身已經(jīng)破損修復程度有限。所以我通常把ShapeCHK當作“體檢工具”而不是“萬能后悔藥”跑一遍沒報錯說明文件基礎結(jié)構(gòu)沒問題報了錯趕緊找原始備份不要抱著修復后的文件繼續(xù)做分析。順帶一提如果你想把SHP屬性表轉(zhuǎn)成txt用于核對可以在ArcGIS里右鍵屬性表選擇“導出”把dbf內(nèi)容另存為CSV或文本文件。CSV文件用Excel打開時字段編碼經(jīng)常亂我會選擇用VSCode或Notepad打開并切換到UTF-8編碼查看。5. 實操避坑滹沱河流域SHP處理時的五個高頻問題5.1 計算面積偏差大數(shù)值小得不像流域現(xiàn)象面SHP明明是完整流域用面積字段統(tǒng)計結(jié)果只有幾千平方米跟真實流域面積差了至少三個數(shù)量級。原因圖層坐標系處于經(jīng)緯度地理坐標狀態(tài)SHAPE.AREA按度計算不是平方米。在沒有投影的情況下Shape_Area字段不具備實際物理意義。解決先確認.prj內(nèi)容通過“投影”工具把SHP轉(zhuǎn)到適合滹沱河流域的高斯-克呂格投影。滹沱河流域我常用CGCS2000_3_Degree_GK_CM_114E。轉(zhuǎn)換完成后重新計算SHAPE.AREA得到的數(shù)值才接近真實平方米。如果不想動原數(shù)據(jù)計算字段時用!SHAPE.GEODESICAREA!也能繞開投影問題。5.2 面重疊和自相交做裁剪時報錯或結(jié)果為空現(xiàn)象用流域SHP去裁別的數(shù)據(jù)系統(tǒng)提示“幾何未正確初始化”或者輸出空圖層打開屬性表發(fā)現(xiàn)要素數(shù)為0。原因SHP里的多邊形存在自相交、重復節(jié)點、懸空邊界等拓撲問題。面文件在Web端下載或者從CAD轉(zhuǎn)換時沒有經(jīng)過拓撲檢查這樣的幾何在ArcGIS的許多分析工具中會被判定為不可信。解決用“數(shù)據(jù)管理工具→要素→修復幾何”工具跑一遍。修復幾何時勾選“刪除重復節(jié)點”和“修復自相交”。修復前先備份因為這類工具對數(shù)據(jù)有侵入性。如果修復后還有問題用ArcGIS的“拓撲”工具集建一個拓撲規(guī)則指定“不能重疊”再逐個檢查。5.3 屬性連接總是匹配不上中文名稱對不上號現(xiàn)象把Excel統(tǒng)計表按鄉(xiāng)鎮(zhèn)名稱連接到流域SHP屬性表結(jié)果大量記錄顯示“空”明明兩個表里都有同一個鄉(xiāng)鎮(zhèn)名。原因最常見的是兩個表里鄉(xiāng)鎮(zhèn)名存在不可見字符Excel從其他系統(tǒng)導出時字符串字段前面或后面帶著空格。另一個原因是SHP的.dbf文件編碼是GBK而Excel另存為UTF-8連接時字符不一致。解決連接前先用“表轉(zhuǎn)Excel”導出SHP屬性表然后跟外部Excel對比特別注意字段長度和空格。Excel里用TRIM()函數(shù)清理兩邊的空格再用VLOOKUP人工核對一遍。對編碼問題檢查SHP目錄下是否有.cpg文件如果沒有用QGIS打開原SHP并另存為UTF-8編碼的新SHP再回到ArcGIS里做連接。5.4 拷貝到別人電腦上打不開提示文件缺失現(xiàn)象把SHP從自己電腦傳到同事電腦對方說圖層能加載但屬性表空白或者干脆報“打開SHP文件失敗”。原因只拷貝了.shp文件.shx和.dbf沒跟上。更隱蔽的是三個文件名被改了其中一部分比如.shp叫“hutuo.shp”.shx卻叫“hutuo_new.shx”引擎按文件名匹配不到索引文件。解決傳輸SHP時把文件夾整體壓縮打包不要單獨拖一個文件。如果已經(jīng)有散亂的三件套改用“ArcCatalog→復制粘貼”方式重新生成一套確保命名統(tǒng)一。我的習慣是每次發(fā)送前用ShapeCHK檢查一遍確認三個文件齊全再發(fā)。5.5 省界線文件省1和省2的SHP有什么區(qū)別疊加時對不齊現(xiàn)象拿到國家基礎地理數(shù)據(jù)里的省1和省2兩個SHP加載到同一地圖中發(fā)現(xiàn)滹沱河流域附近的省界線條并不完全重合做疊加分析后出現(xiàn)很窄的縫隙或重疊。原因省1和省2是不同生產(chǎn)精度的界線數(shù)據(jù)。省1一般指經(jīng)過概化、用于中小比例尺制圖的省級行政區(qū)界線省2是更精細的邊界包含更復雜的細節(jié)和國界要素。兩者比例尺和拓撲規(guī)則不同幾何自然不完全一致。解決做流域空間分析時只選一套省界作為底圖不要同時疊加兩套并直接做Union。要以省2即精度更細的版本作為參考邊界省1可以用來做示意但不要參與拓撲計算。如果你發(fā)現(xiàn)流域SHP邊界和省界有沖突優(yōu)先以流域?qū)嶋H地貌范圍為準行政界線只是參考圖層。6. 再用深一點漁網(wǎng)分割、邊界簡化與三維場景過渡6.1 用“創(chuàng)建漁網(wǎng)”把大流域切成分塊有時你不需要整個滹沱河流域作為一個面做分析而需要把它劃分成2公里×2公里的網(wǎng)格逐格統(tǒng)計土地利用類型或者做水環(huán)境評價。手動畫方格不現(xiàn)實正確做法是用ArcGIS的“創(chuàng)建漁網(wǎng)”工具。漁網(wǎng)工具的本質(zhì)是生成一個規(guī)則的格網(wǎng)覆蓋指定范圍。關鍵參數(shù)有三個漁網(wǎng)原點、像元寬高和模板范圍。模板范圍直接選擇SHP面文件時漁網(wǎng)輸出邊界會跟流域的范圍框一致。下圖場景里我創(chuàng)建2公里格網(wǎng)然后與流域SHP疊加保留與流域相交的網(wǎng)格import arcpy arcpy.env.workspace rD:\hutuo_data arcpy.management.CreateFishnet( out_featuresgrid_2km.shp, origin_coord490000 4250000, y_axis_coord490000 4250010, cell_width2000, cell_height2000, number_rowsNone, number_columnsNone, labelsNO_LABELS, templatehutuo_river_basin.shp, geometry_typePOLYGON )參數(shù)origin_coord是漁網(wǎng)左下角的投影坐標y_axis_coord決定網(wǎng)格旋轉(zhuǎn)方向通常讓它與正北方向一致。template填流域SHP生成的漁網(wǎng)會自適應流域的最大外接矩形。生成網(wǎng)格后還要用“相交”工具把漁網(wǎng)和流域SHP相交這樣網(wǎng)格邊緣就貼合流域邊界而不是簡單的矩形外框。這一步如果不做你統(tǒng)計的格網(wǎng)會有一大半在流域外結(jié)果毫無意義。6.2 面邊界簡化與三維場景過渡流域SHP從專業(yè)水文機構(gòu)下載時邊界往往包含海量節(jié)點一個流域面可能帶著幾萬個點。這種精細數(shù)據(jù)在小比例尺出圖時不僅讓電腦卡頓還讓ArcMap的符號化變得遲鈍。我的經(jīng)驗是保存兩套數(shù)據(jù)一套原始精度用于分析一套簡化版用于制圖和發(fā)布。ArcGIS有條“簡化面”工具位于制圖工具箱中import arcpy arcpy.cartography.SimplifyPolygon( in_featureshutuo_river_basin.shp, out_feature_classbasin_simplify.shp, algorithmPOINT_REMOVE, tolerance50 Meters, minimum_area0 SquareMeters )POINT_REMOVE算法通過刪除與相鄰點夾角接近直線的點來降低節(jié)點密度50 Meters的容差表示小于50米的邊界彎曲會被忽略。簡化閾值不是越大越好如果設成500米流域邊界的細小河道特征會全部丟失。出圖給領導看容差可以放寬做水文分析保持原始精度不要隨意簡化。如果你準備把流域SHP推送到三維場景比如生成3DTiles做Web端可視化水汽效果會清楚反映出模型對流域面的支撐力?；驹硎前衙鍿HP先轉(zhuǎn)成GeoJSON再疊加上高程數(shù)據(jù)最后切片成b3dm。ArcGIS里可以用“轉(zhuǎn)出至GeoJSON”工具后續(xù)步驟交給開源工具完成。注意SHP屬性表里的中文名稱在GeoJSON里可能編碼錯亂轉(zhuǎn)之前先在屬性表里增加一個英文字段。我從實際項目里積累的習慣是拿到任何SHP第一件事復制一個原版?zhèn)浞萑缓笥肧hapeCHK檢查文件結(jié)構(gòu)用arcpy.Describe確認坐標系和幾何類型再檢查屬性表最后才開始裁剪和轉(zhuǎn)換。這套流程看起來繁瑣但確實幫我避免了很多在交付前一刻才發(fā)現(xiàn)幾何錯誤的窘境。尤其是滹沱河流域這類邊界跨度較大的面文件坐標系和拓撲問題幾乎是必然存在的早檢查早安心。希望這篇筆記里的“三件套檢查投影確認幾何修復”能幫到你至少下次再遇到“面文件算不出面積”的怪事時你知道該往哪幾個方向排查。本文還有配套的精品資源點擊獲取