简介中国一级河流分布数据集基于1:100万河流数据提取涵盖长江、黄河、黑龙江等主要一级河流的空间走向与边界信息适合GIS用户、科研人员用于水资源管理、生态评估、流域规划及防洪减灾等专题制图。包体共70个文件压缩包大小156.15MB核心文件包括SHP标准矢量数据、MXD可编辑制图模板、TIF标准成图及配套的dbf、shx、prj等辅助文件。其中SHP可直接在ArcGIS、QGIS中进行空间分析与可视化MXD预置河流样式、流向箭头、图例、比例尺等要素便于快速编辑出图TIF为600dpi高清地图可直接用于报告与演示。资源还附带中国行政区、九段线、世界底图等基础地理数据方便叠加制图。目前已有169人学习适合需要规范河流分布底图及一体化成图方案的研究与规划人员。1. 中国一级河流分布shape数据到手不代表能用mxd与成图才是完整的可交付成果多数人拿到一份中国一级河流分布的shape文件第一反应是拖进ArcMap看能不能显示。结果往往有两种要么河流显示出来但错位严重要么线型和标注完全没法直接出图。原因很简单——shape文件只是几何和属性的载体坐标系是否正确、符号系统是否配置、图面是否按制图规范排版这些都不包含在.shp里。真正能交付给同事、领导或甲方直接用的是「标准shape文件 可编辑mxd工程 标准成图」三件套。这套东西解决的不只是有没有数据的问题而是拿到数据后能不能少走弯路、直接进入生产流程的问题。适合谁用水文规划、国土空间规划、GIS数据生产、地图制图岗位的从业者以及准备把河流基础数据做成项目底图的技术人员。2. 从原始河流数据到标准shape文件坐标系、拓扑与口径一次对齐2.1 一级河流怎么定义先定口径再做字段筛选中国一级河流这个说法在公开数据里比较含糊。水文专业里常说的一级支流是指直接汇入干流的支流而国家一级河流在不同部委的清单里又有差异。做数据生产时不能拿着概念就开工第一步是先确定口径否则后面筛选出来的河流数量、名称、分级都会对不上。我一般以两个口径交叉确认一是水利行业标准《中国河流代码》SL 249里对主要江河及其支流的编码体系二是全国水资源调查评价中采用的河流分级表。实操层面先在原始数据里找到分级字段常见的有LEVEL、GRADE、CLASS几种命名值为1的通常是干流或一级河流。如果原始数据里没有分级字段就按河流名称匹配已知清单——长江、黄河、珠江、松花江、淮河、海河、辽河、塔里木河、怒江、澜沧江等入海或入内陆盆地的骨干河流单独挑出来。这个环节最容易翻车的是把一级支流当成一级河流。比如长江的支流岷江、嘉陵江名称上是一级支流但在国家河流分级里属于二级甚至三级。所以字段筛选前一定要把数据生产说明或元数据里的分级定义读一遍别凭直觉按支流层级筛选。2.2 统一到CGCS2000定义投影与重投影脚本拿到手的公开河流数据坐标系千奇百怪。常见的有WGS84地理坐标系、北京54、西安80投影坐标系甚至有些数据没有.prj文件。中国境内做标准成果现在统一用CGCS2000国家大地坐标系2000它是国内测绘和自然资源行业的标准基准。地理坐标系用CGCS2000EPSG:4490投影坐标系按所在中央经线选高斯-克吕格3度带比如EPSG:4521到4530之间。用ArcPy做批处理比较省事。下面这段脚本是把一份无坐标系或坐标系错误的原始线层先定义为CGCS2000地理坐标系再投影到3度带import arcpy arcpy.env.workspace rD:\river_work arcpy.env.overwriteOutput True # 原始数据假设它实际是CGCS2000但缺prj文件 raw_fc river_raw.shp defined_fc river_defined.shp # 1. CGCS2000地理坐标系EPSG 4490 sr_geo arcpy.SpatialReference(4490) arcpy.DefineProjection_management(raw_fc, sr_geo) # 2. 投影到CGCS2000 3度带中央经线117°E sr_proj arcpy.SpatialReference(4528) arcpy.Project_management(raw_fc, defined_fc, sr_proj)逻辑说明第一步DefineProjection_management是给没有坐标系信息的要素类贴标签前提是你确认数据本身的几何坐标值在正确范围内——比如经纬度值大约在73°E到135°E之间。第二步Project_management做真正的坐标转换会把经纬度换算成平面坐标。4528是CGCS2000_3_Degree_GK_Zone_16不对要核对中央经线117°E对应的3度带带号是39EPSG代码是4528没问题。如果你拿捏不准用arcpy.ListSpatialReferences查一下名称再填代码。参数上最需要注意的是中央经线。全国范围出图时用地理坐标系就行出分幅图或入库才按区域转投影。跨带河流比如长江跨了3个3度带强行用一个投影带会导致边缘变形这种情况下我一般保留地理坐标系做全图分带投影只用于局部成果。2.3 拓扑检查与线的简化shape文件能被ArcMap直接读取的前提坐标统一后还有个隐藏问题原始河流线层经常有自相交、重复线段、断线悬挂点。这些拓扑错误在屏幕上看不出来但会导致长度计算偏大、标注位置异常、出图时河流线有毛刺。检查拓扑可以用ArcGIS的拓扑工具也可以直接用QGIS的v.clean一步搞定。ArcGIS的操作路径是新建要素数据集 → 导入河流线层 → 新建拓扑规则不能自相交、不能有伪节点 → 验证 → 按错误批量修复。如果是用代码QGIS的Processing框架里v.clean很高效import processing from qgis.core import QgsVectorLayer layer QgsVectorLayer(rD:\river_work\river_defined.shp, river, ogr) result processing.run(grass:v.clean.polylines, { input: layer, threshold: 1.0, # 容差1米 output: rD:\river_work\river_clean.shp })threshold参数是关键意思是把距离小于1米的节点合并。设太小没用设太大会把河流的弯曲细节吃掉。我一般取数据精度的2倍1:100万数据用2米容差比较稳妥。几何清理后再检查一下要素数量如果数量变化超过1%回头看一下是不是容差过大误删了短小支流。3. 属性表设计落在shape文件上字段、代码、长度计算的完整做法3.1 字段设计名称、编码、长度、流域面积一个都不能少标准shape文件的生产属性表设计决定了后续制图和统计能不能自动化。只带几何和名称的河流数据到出图阶段会发现图例没法做、注记没法统一、筛选条件只能靠手选。我一般给一级河流线层设计这么几个字段RIVER_NAME河流中文全名、RIV_CODE按SL 249编码、LEVEL短整型1表示一级河流、ST_LEN_KM双精度长度单位公里、BASIN_AREA双精度流域面积平方公里按需填写、SOURCE_INFO文本填写数据来源和年份。字段名的设计有讲究。shape文件的DBF属性表字段名有长度限制而且不能用中文字段名——大多数GIS软件对中文字段名的支持是历史遗留灾难。所以我全部用英文短名中文含义放到字段别名里展示。RIV_CODE这一项大家容易漏但它是河流图斑关联统计年鉴、水文站数据的唯一标识没有编码后期做查询就要靠文字匹配非常痛苦。3.2 用ArcPy补字段、算长度、改别名一条命令一个坑字段设计好后用ArcPy把结构和计算一步到位import arcpy fc rD:\river_work\shape\River_Level1.shp # 1. 添加字段 arcpy.AddField_management(fc, ST_LEN_KM, DOUBLE, , , 16, 河流长度(km)) arcpy.AddField_management(fc, RIVER_NAME, TEXT, , , 50, 河流名称) arcpy.AddField_management(fc, RIV_CODE, TEXT, , , 20, 河流编码) # 2. 计算长度shape.lengthKILOMETERS 返回公里数 arcpy.CalculateField_management(fc, ST_LEN_KM, !shape.lengthKILOMETERS!, PYTHON3) # 3. 修改字段别名 arcpy.AlterField_management(fc, ST_LEN_KM, ST_LEN_KM, 河流长度(km)) arcpy.AlterField_management(fc, RIVER_NAME, RIVER_NAME, 河流名称)逻辑说明第二步用!shape.lengthKILOMETERS!求线长度返回的是投影坐标系下的平面长度。所以字段添加一定要放在坐标系转换和拓扑清理之后否则数据实际是经纬度时shape.length返回的是十进制度数算出来的长度毫无意义。PYTHON3指定表达式解析器保证!字段名!语法有效。第三步AlterField_management改的是别名不影响字段名本身。这一步容易被忽略但mxd里图例、标注表达式、属性表表头都读别名不改别名图例上就会显示ST_LEN_KM这种硬编码显得很不专业。要注意的是在ArcGIS Pro里AlterField支持更好ArcMap 10.6以上的版本才稳定支持这个函数。3.3 DBF编码与.cpg河流名称乱码的根源和解决办法shape文件中属性表存储在.dbf文件里没有内置编码信息全靠同目录下的.cpg文件声明编码。ArcMap读取.dbf时默认按系统语言代码页处理Windows中文系统下经常把UTF-8的数据按GBK读结果就是长江变成闀挎睙且无法直接改。解决办法有两个看你的数据从哪来。如果源数据是QGIS另存出来的UTF-8 shape文件就在同目录放一个River_Level1.cpg内容写UTF-8ArcMap再打开就不乱码。如果数据是DWG或MapGIS转出来的GBK编码需要在QGIS里做一次另存为转成UTF-8再入库。实际生产里我强烈建议把属性表从DBF里抽出来入库到FileGDB或PostGIS做编辑最后导出成shape时才转DBF格式这样能避开很多历史问题。4. mxd可编辑工程文件搭建模板化的制图环境4.1 mxd里到底存了什么为什么没有shape文件也能打开工程mxd是ArcMap的工程文件很多人误以为数据被存在mxd里其实它只存储引用和配置数据框坐标系、图层顺序、符号化设置、标注引擎参数、图面整饰元素的位置以及数据源的相对/绝对路径。没有mxd你每次开ArcMap都要重新拖数据、调符号、加图例有了mxd打开就是一套排好的图。关键理解点是可编辑。交付的mxd应当是数据源指向明确的工程文件双击能打开、图层能增删、符号能改、标注能调。要做到这一点必须保证shape文件和mxd的相对路径关系稳定并且符号系统不是靠内存里的临时配置而是保存在每个图层的.lyr文件或mxd固有的图层设置里。4.2 手动搭建一次mxd数据框、标注、图廓一整套我自己搭建mxd有固定的顺序步骤不对后面返工很麻烦第一步新建空白mxd在数据框属性里把坐标系设为CGCS2000地理坐标系这样无论图层是什么坐标系数据框都会按统一基准做动态投影。第二步把清理、编码好的河流shape文件拖进地图右键图层选择属性→符号系统用类别按LEVEL字段分级渲染一级河流线宽1.5mm其余0.6mm颜色用标准水系蓝。第三步打开标注选项卡勾选标注此图层中的要素标注字段选RIVER_NAME字体用黑体或宋体字号设8pt按Place label along line沿河流方向放置。第四步切换到布局视图插入图例、比例尺、指北针、图名和资料来源说明。图例只显示RIVER_NAME要素的符号比例尺按出图比例设为1:200万。第五步最关键的一步文件→地图文档属性→勾选存储数据源的相对路径名然后保存mxd。这一步不做整个mxd换个目录就废了。4.3 用arcpy批量生成mxd适合批量出图的工程化做法如果只是做一两个省的一级河流手动搭mxd很快。但要做全国分幅出图比如按1:200万标准分幅切成几十个图幅手动搭建会累到怀疑人生。用arcpy.mapping模块可以批量生成mxdimport arcpy # 模板mxd先手动搭好符号、图廓和图例存为template.mxd template arcpy.mapping.MapDocument(rD:\river_work\template.mxd) df arcpy.mapping.ListDataFrames(template)[0] # 要出图的份数按分幅范围循环 frames [ (H-48, rD:\river_work\tiles\H-48.shp), (J-50, rD:\river_work\tiles\J-50.shp), ] for sheet_name, layer_path in frames: mxd_copy template # 清除旧图层添加当前分幅河流数据 for lyr in arcpy.mapping.ListLayers(mxd_copy): arcpy.mapping.RemoveLayer(df, lyr) new_layer arcpy.mapping.Layer(layer_path) arcpy.mapping.AddLayer(df, new_layer, TOP) new_layer.symbology arcpy.mapping.Layer(rD:\river_work\river_symbology.lyr) out_mxd rD:\river_work\mxd\River_{0}.mxd.format(sheet_name) arcpy.mapping.SaveToACopy(mxd_copy, out_mxd)这段脚本的坑不少。第一arcpy.mapping.MapDocument是ArcMap的环境如果你用的是ArcGIS Pro必须换成arcpy.mp.ArcGISProjectAPI完全不同。第二每轮循环要先把上一轮的图层移除干净否则会叠加。第三符号系统单独存成.lyr文件的好处是改线宽、改颜色不用重新跑脚本只改lyr文件再跑一遍循环即可。saveToACopy生成的是莫副本不会污染模板。5. 河流制图生产中的5个常见问题与避坑记录5.1 河流整体错位、跑到省界外坐标系混用现象一级河流图层在ArcMap里和底图比如省界对不上河流整体偏移几百米到几公里。原因通常是原始数据是WGS84经纬度但.prj文件错误地标成CGCS2000或者数据框坐标系与图层坐标系不一致。动态投影在屏幕上会强行显示但导出成图、量算长度时错误就暴露了。解决不依赖动态投影用Project_management把数据重投影到目标坐标系并检查.prj文件内容确保和元数据一致。这一步偷懒后面全白干。5.2 移动工程目录后mxd图层红色感叹号相对路径没开现象把整个工作目录从D盘拷到移动硬盘打开mxd所有图层名前出现红色感叹号数据全部丢失显示。原因ArcMap默认存储绝对路径mxd里的数据源还指向D盘原路径。解决保存前在地图文档属性中勾选存储相对路径名同时保证数据目录结构一致——比如mxd在project\mxd\下shape在project\shape\下保持相对关系拷贝即可。如果你是接手别人的mxd遇到红色感叹号也可以右键图层→属性→源→设置数据源重新指向。5.3 属性表河流名称乱码、字段丢失DBF字符集现象打开属性表RIVER_NAME字段显示或乱码个别字段完全不存在。原因DBF的编码声明依赖.cpg文件源数据是GBK编码但没带.cpgWindows按UTF-8去读或者字段名超过10个字符被截断。解决用QGIS打开乱码的shape在另存为里选择UTF-8编码自动生成.cpg文件字段名控制在10个字符以内超过的用别名表达。5.4 标注挤成一团压住河面标注引擎参数没调现象用默认标注方式长江、黄河这种长河道的名称标注要么重叠在支流标注上要么跑到图面之外。原因标注放置策略没有按线要素优化没有设置标注优先级和字符间距。解决标注引擎里选择沿几何放置设置标注层次——一级河流标注权重最高支流权重降低再不行就把标注转成注记Convert labels to annotation手动拖动位置。对标准成图来说转注记才能保证最终输出的图面完全可控。5.5 线宽跟着缩放变蜘蛛网参考比例尺缺失现象在布局视图里缩小或放大显示河流线宽要么细成头发丝要么粗成一条黑带。原因符号单位默认是像素屏幕缩放时线宽不随地图比例尺变化但出图DPI和比例尺一变效果和预览完全不同。解决在布局视图设置参考比例尺Reference Scale我习惯设为1:200万然后符号线宽用毫米制mm设定。这样无论缩放多少倍最终成图线宽都是物理尺寸打印出来和屏幕上效果一致。6. 标准成图验证与一个可复用的出图习惯输出前查一遍这5项标准成图的验证不是看起来没问题而是用数据检查说话。我每次交付前必查的5项是要素数量是否和源数据一致用arcpy.GetCount对比属性表是否有空值按RIVER_NAME IS NULL做选择集坐标系名是否和文件名声明一致mxd用相对路径再次打开是否正常导出PDF后放大到200%检查线型、注记压盖、图廓线是否完整。一个值得长期保留的习惯是三个目录加自解释文件名。原始数据放00_source只读目录中间过程放10_process_gdbFileGDB工程库最终交付放20_delivery。文件名写成River_Level1_CGCS2000_3DGK_117E.shp坐标系、投影带、中央经线全部写进文件名看到文件名就知道数据基准避免别人用的时候还要猜坐标系。mxd同理命名带版本日期。这个习惯帮我避免了很多次这数据到底是WGS84还是CGCS2000的血泪追问。希望帮到你。本文还有配套的精品资源点击获取