简介全国3300多个自然保护区空间分布数据包面向GIS分析人员、生态科研与环保管理者可用于保护区分布格局研究、生态规划与制图。数据截至2024年覆盖全国3300余个保护区属性字段含级别、名称、始建年份、所属省市县、面积、编号、类型、保护对象等。压缩包共10个文件约154MB包含标准Shapefile空间数据shp/dbf/shx/prj/xml、可编辑mxd工程文件、标准成图TIF以及Excel属性表另附2个zip压缩包含中国行政区等辅助数据。其中shp适用于ArcGIS等平台做叠加与缓冲区分析mxd已按中国行政区精细配准TIF可直接用于成果展示。已有167人浏览学习。通过该数据包可快速获取全国保护区空间位置与属性信息免去逐一手动收集整理的麻烦mxd便于个性化编辑Excel表格便于统计汇总适合撰写研究报告、制作专题地图或开展生态保护决策支持。1. 全国3300多个自然保护区空间数据一套能直接开工的标准底图包做生态评价、环评或者林地一张图的同行应该都有这种体会画图前最耗时间的不是模型本身而是找底图、理边界、统一坐标系。网上能下的全国尺度数据不少真正能直接开工的很少——要么是教程扣出来的局部Demo要么坐标系东一块西一块要么属性表字段残缺。全国3300多个自然保护区的空间分布数据是把边界、名称、类型、出图一次配齐的三件套MXD工程文件管排版、标准Shapefile管数据和调用、标准TIF管汇报展示。适合搞规划、写环评报告、做GIS分析和生态评价的从业者拿到手就能拖进ArcMap开工不用再从在线地图描边凑底图。2. MXD、SHP、TIF三种文件怎么选工程文件管排版、矢量文件当本体、成图当底图这个资源包解压后你能看到三种核心格式.mxd、.shp和.tif。不少新手第一反应是双击MXD然后发现图层红叹号、图形飞走就开始怀疑数据有问题。其实不是数据坏了而是没搞懂这三种格式的分工。MXD是ArcMap的工程文件存的是图层路径、符号样式、版面布局不存几何SHP是真正的边界和属性本体TIF是已经渲染好的标准底图。三者作用不同使用顺序也不同。下面按这三类格式分别说清楚定位、打开方式和适用场景。2.1 MXD不是数据本体先分清“打开什么”和“改什么”MXD在ArcMap里是“现场排版文件”它记录了图层怎么叠加、符号怎么配、比例尺放在哪、图例用什么字体。拿到数据包先打开MXD看到的是数据作者预设好的全国保护区一张图效果。但这不意味着边界就存在MXD里——MXD只是指路牌真正的几何在SHP里。我一般建议工作流反过来先新建一个空白地图把SHP直接拖进去确认数据完整、坐标正确之后再打开MXD看排版。如果一上来就依赖MXD换一台电脑后很容易遇到数据源丢失的问题。MXD里记录的绝对路径指向某位数据制作人员的C盘目录到你机器上路径不存在图层全部变红叹号。遇到这种情况别急着重建先右键红叹号图层打开Properties→Source看看“Set Data Source”能不能重新指到解压目录下的SHP。能指就指指不上就去改MXD的文档属性。在ArcMap里点File→Map Document Properties勾选“Store relative pathnames to data sources”以后整包拷贝MXD和SHP到其他盘符路径变了也能自动相对查找。这个操作每次交付前值得强制走一遍。MXD本身还有一个容易翻车的点ArcGIS 10.x创建的MXD拿到ArcGIS 10.2以下版本打开标注位置大概率会错位。如果你办公室用的是老版本先把MXD另存为较低版本兼容格式再继续编辑。2.2 Shapefile才是干活的核心属性表里有边界、名称和级别SHP是这套数据里最值得细看的本体。一个完整的Shapefile由至少三个文件组成.shp存几何、.dbf存属性、.shx存索引再加上可选的.prj存坐标系。解压后如果发现某个图层缺了.dbf属性表打开就是空的缺了.prjArcMap会猜坐标系边界偏到海里都查不出来。这套数据的SHP属性表里一般会包含保护区名称、级别国家级/省级/县级、类型森林生态/湿地/荒漠/野生动物等、所在行政区域和面积字段。做筛选时可以直接按“级别”字段做符号化也可以按“类型”字段把全国森林类保护区单独提出来分析。实际动手之前建议先用预览窗口看一遍坐标范围。全国陆地经纬度范围大致在73°E到135°E、18°N到54°N之间如果预览时看到坐标是8位数的米制坐标例如4540000之类的数值说明SHP已经是投影坐标系不是经纬度。两种情况都能用但叠加时别混用。遇到混用要么用投影工具统一到CGCS2000/Albers等距圆锥投影要么统一回WGS84经纬度再分析各图层保持一致就行。2.3 TIF标准成图定位是“出图”能看不能改TIF是栅格成图通常已经配好颜色、边界和注记供汇报、插图或者快速预览使用。它是GeoTIFF格式文件里带地理参考信息拖进ArcMap/QGIS能直接对到正确坐标位置。但它和矢量SHP的本质差别是TIF的边界是像素不能选中、不能编辑、不能单独提某个保护区属性。使用场景一般是两种一是放到底图上当背景快速查看全国保护区的整体分布格局二是导出报告插图时直接用TIF保证配色统一不用再花时间重调符号。你要是想在TIF上改一个边界线位置那就得回到SHP编辑完成后重新导出成图而不是直接拿PS去涂TIF。TIF只能做“结果”不能当“源数据”。另外注意检查TIF是否带透明通道。部分标准成图为了美观会留白边或者画框放到PPT深色背景下会出现一个大白块。后面第四部分会讲怎么处理这类背景问题。3. 从解压到第一张出图坐标系检查、字段复核和图幅调整资源包下载后不是双击MXD就万事大吉第一步永远是“验数据”。我习惯把流程拆成三步先看目录结构、再查坐标系和要素基数、最后才动排版。文本步骤看起来多但做顺了五分钟内能完成。3.1 解压后第一步检查目录、坐标参考和要素基数先把整个包解压到一个不带空格和中文的路径下比如D:\gis\data\reserve。路径带中文在ArcMap某些旧版本里会导致数据源连不上、SHP字段编码错乱。解压后先看一眼根目录确认MXD指向的SHP文件是否都在同名文件夹内。如果包比较大、图层多逐个手点太慢。可以用一个简单的ArcPy脚本批量读取每个SHP的信息输出坐标系名称、要素数量和几何类型几秒钟就把整个包摸清底细。import arcpy import os workspace rD:\gis\data\reserve arcpy.env.workspace workspace for shp in arcpy.ListFeatureClasses(): desc arcpy.Describe(shp) count int(arcpy.management.GetCount(shp).getOutput(0)) sr desc.spatialReference print(文件: {}, 类型: {}, 要素数: {}, 坐标系: {}.format( shp, desc.shapeType, count, sr.name))这段ArcPy脚本的逻辑是ListFeatureClasses()遍历工作空间下所有矢量要素类拿Describe读取几何类型拿GetCount取要素总数最后把坐标系名称打印出来。参数上要注意arcpy.env.workspace必须指向SHP所在目录路径里不能有中文字符。输出里如果看到多个SHP的坐标系名称不一致先用投影工具统一再到下一步否则后面叠加会出缝。3.2 用ArcPy批量合并与计算面积一点脚本代码做全国尺度的统计时经常要把3300多个面合并成按类型或按级别汇总的结果。ArcMap里的Dissolve工具能实现但每次手动选字段容易漏。用脚本参数化之后换一套SHP还能复用。import arcpy arcpy.env.workspace rD:\gis\data\reserve input_fc 保护区边界.shp output_fc 保护区按类型合并.shp arcpy.management.Dissolve( in_featuresinput_fc, out_feature_classoutput_fc, dissolve_field[类型], multi_partMULTI_PART )这个小节里多数读者并不需要严格复制但如果你恰好要用务必注意dissolve_field参数列表里必须写属性表里真实存在的字段名写成中文或类别名都会报错。合并之后检查一下输出要素数量全国森林生态类、湿地类、荒漠类保护区通常不超过十类如果输出要素数和输入差不多多半是“类型”字段没有值得先处理空白字段再合并。计算面积时我推荐用测地面积不要用默认的平面面积。在ArcGIS里添加几何属性时选择AREA_GEODESIC单位设为SQUARE_KILOMETERSarcpy.management.AddGeometryAttributes( in_features保护区边界.shp, geometry_propertyAREA_GEODESIC, length_unit, area_unitSQUARE_KILOMETERS )用测地面积的原因很简单全国范围跨多个投影带平面面积会把高纬度地区的面积算大测地面积基于椭球体计算跨带不偏。参数area_unitSQUARE_KILOMETERS直接控制输出面积字段的单位需要平方米就改成SQUARE_METERS。3.3 改图例、比例尺和标注在MXD里完成“排版”数据本身没问题之后才轮到动MXD。打开MXD后找到左下角的“Layout View”切换按钮进入版面视图。此时页面上的指北针、比例尺和图例都是占用版面空间的独立元素可以拖动位置和调整大小。改指北针时右键选择Properties常见制图规范要求指北针方向为正北不要选斜角样式。比例尺我这里给个经验值出A3横向图幅比例尺放在1:600万到1:1200万左右比较合适再大就显得空再小字挤在一起看不清。如果保护区分布图上要同时体现国家级和省级图例按“级别”字段分四到五类配色不要用连续色带否则图例范围超出类别含义。字体系列统一用黑体或思源黑体。ArcMap默认宋体导成PDF后在没装中文字体的电脑上会临时替换标注位置可能偏移。检查方式很简单只改图例标题字体整体预览一遍即可。注意MXD里导出的图片默认不压缩导出为PDF时会自动嵌入字体最好在导出设置里勾选“Convert text to curves”或“Export text as graphics”这样交付给别人的PDF不会因缺字体而乱版。4. 常见问题排查坐标系、记录数与TIF背景的五个实战坑以下五个问题是用这套数据时最容易踩的。每一条都是我实际遇到或复现过的按“现象→原因→解决”梳理排查时可以直接对着看。4.1 打开MXD提示“未知的空间参考”图形飞到非洲现象双击MXD后ArcMap正常打开但数据显示在地图窗口右下角或黑乎乎一片缩放范围一看竟在几内亚湾附近。原因MXD里的Data Frame坐标系与SHP的.prj坐标参考不一致或者SHP缺少.prj文件ArcMap按默认WGS84去猜坐标结果自然对不上。解决先查看SHP的Properties→Source选项卡确认坐标系名称。如果SHP没有.prj文件用“Define Projection”工具手动设定坐标系不要直接使用“Project”工具因为“Project”是在已有坐标系基础上做转换不是赋予坐标系。设定后再回到Data Frame属性里把坐标系切换到与SHP一致图形位置就会恢复正常。4.2 属性表记录数对不上3300多个多半是重名或缺值现象统计要素数量时结果比说明里写的“全国3300多个”少一截或者字段里出现大量空值。原因部分保护区分片管理同一名称的面要素被拆成多个部分属性表里占多行也有部分SHP在制作时漏了属性录入导致按名称筛选时统计不上。解决先用“频数统计”工具以“名称”字段为统计字段输出记录重复次数。重复值不是错误要区分“同一保护区的多个图斑”和“不同保护区恰巧重名”看“级别所在省份”组合字段是否一致。缺值字段用ArcMap的字段计算器按行政区域或保护区类型做条件赋值不要直接在Excel里改dbf因为Excel对dbf编码支持很差保存后容易出现中文乱码。4.3 保护区边界叠到“一张图”上有缝投影统一再跑拓扑现象将保护区SHP叠加到林地一张图或国土空间规划底图上边界错位几十米到几百米看起来像两张图纸套不上。原因保护区SHP的坐标参考是WGS84经纬度底图用的是CGCS2000或者某个地方坐标系坐标基准不同必然错位。解决统一用“投影”工具把两套数据转换到同一参考系。常见做法是转成CGCS2000 / 3-degree Gauss-Kruger zone或Albers等积圆锥投影具体选哪种看底图本身的定义。转换完成后如果还有局部边界差异打开编辑器选中落单的边界顶点做小幅吸附。3300多个面不可能挨个手改一般只处理你分析区域范围内的相交部分即可全国层叠加时更看重整体格局而非厘米级重合。4.4 TIF插到汇报图里带黑边/白边背景不透明很突兀现象把TIF标准成图直接拉到PowerPoint或Word里发现图幅外围一圈黑底或白底盖住了下面要展示的区域地形。原因TIF是栅格图没有矢量那样的透明背景图片边缘的无效值像素默认按纯黑或纯白渲染。解决在ArcMap里右键TIF图层→Layer Properties→Symbology选项卡打开“Display Background Value”选项把背景值设为0或255并勾选“Display as RGB”下的透明背景。如果TIF自带NoData值还要在“Display NoData”选项里勾选透明。设置完成后导出图片记得选PNG格式PNG保留透明通道JPG会强行压掉透明度再次出现白底。这里有个小坑不同TIF的无效值不一样有的用0有的用255全图一片白时改成0试试全黑就改255。4.5 用Excel直接改dbf字段导致中文乱码现象在Excel里打开SHP对应的.dbf表把保护区名称批量修改后保存回到ArcMap打开属性表字段值变得难以辨认。原因dbf文件使用常规编码Excel保存为默认格式后改变了编码结构非中文环境的ArcGIS再读就乱了套。解决不在Excel里直接改dbf改用ArcGIS属性表编辑器逐行修改或者用Python脚本批量处理import arcpy fc rD:\gis\data\reserve\保护区边界.shp with arcpy.da.UpdateCursor(fc, [名称]) as cursor: for row in cursor: if row[0] is not None: row[0] row[0].strip().replace( , ) cursor.updateRow(row)这段脚本处理的是去掉名称字段首尾空格和全角空格。UpdateCursor按行读取字段修改后调用updateRow写回。注意replace( , )第一参数是全角空格第二个是半角空格在文本编辑器里两个字符肉眼可能看起来一样实际编码不同。如果不小心把全角空格替换成空字符串字段内容会粘在一起影响后续表关联。5. 对接新名录做增量更新字段清洗、关联键与拓扑保持这套数据是基础底图但你手头很可能还有一份更新的保护区名录要把新批复的保护区加进来或者把某些边界调整后的结果替换进去。增量更新的核心不是“画一笔存一次”而是字段清洗、关联键构建和更新检查三件事。处理不好新增的几十个保护区会和你原有的3300多个打架。5.1 先清洗名称字段空格、全半角和括号统一两份数据的行政区写法经常不一致一份写“内蒙古”另一份写“内蒙古自治区”括号有的是全角“”有的是半角“()”名称尾部还可能带不可见空格。直接拿名称做关联十有八九对不上。先把外部名录读进来用pandas做标准化import pandas as pd df pd.read_csv(名录.csv, encodingutf-8) df[名称] df[名称].str.strip() df[名称] df[名称].str.replace(, () df[名称] df[名称].str.replace(, )) df[名称] df[名称].str.replace(r\s, , regexTrue) df.to_csv(名录_清洗后.csv, indexFalse)这段代码是按顺序做的strip()去掉首尾空白replace(, ()和replace(, ))把全角括号统一成半角方便后续字符串匹配最后\s用来吃掉名称内部所有空白字符。regexTrue参数必须显式写否则pandas会把正则表达式当成普通字符串处理清洗结果不如预期。清洗完以后再看一眼“名称”和“级别”两列有没有空值这个动作虽小但能省去后面关联结果里一大堆NaN排查。5.2 用“名称类型”作关联键做增量更新只用名称关联最容易翻车因为不同保护区可能重名同一保护区又可能因为类型字段不同被分到不同行。安全的做法是把“名称级别”或者“名称类型”拼成组合键。import geopandas as gpd import pandas as pd old gpd.read_file(保护区边界.shp) new pd.read_csv(名录_清洗后.csv, encodingutf-8) old[key] old[名称].astype(str) | old[类型].astype(str) new[key] new[名称].astype(str) | new[类型].astype(str) merged old.merge( new.drop_duplicates(subsetkey), onkey, howleft ) merged.to_file(保护区边界_增量.shp, encodingutf-8, driverESRI Shapefile)这段脚本把旧SHP里的“名称类型”列和新名录里的“名称类型”列拼接成单字段之后执行左连接。旧数据的几何和属性都保留新名录中匹配上的字段会加新列写入。参数上注意sep|的选择避免字段内容本身包含“|”符号造成误分裂实际操作里用短横线-或下划线_更稳妥。连接完成后把merged[merged[新字段].isna()]打印出来看这批行就是新旧名录对不上的记录多半原因是名称措辞不同或类型划分不同需要人工确认。匹配上之后再把新增图斑补充进SHP。新增图斑不能只贴坐标数据还要同时补上边界几何和属性字段否则面积统计时会把新图斑算成无属性记录。常见做法是把新增图斑单独建成一个SHP用Append工具或Merge工具合并到原文件然后重新计算面积字段。5.3 更新后重新出图的四个检查项更新完成后别急着关电脑先跑一遍出图检查避免交付出了问题再返工。第一项检查要素基数。合并前后对比总数如果新增了50个图斑总数应该精确增加50多出来的说明关联时产生了重复图斑。第二项检查几何有效性。在ArcToolbox里跑“Check Geometry”工具重点看“Null geometry”和“Self intersections”更新操作中最容易写坏的就是几何顶点顺序。第三项检查字段完整性。新增图斑的“级别”“类型”“面积KM2”三个字段不能为空面积字段必须重新计算否则汇总统计时新增图斑会拖低平均值。第四项检查符号化配色。图例如果按“级别”字段分类级别值为空的图斑会落到图层最底层不显示要把等级字段缺值的情况提前处理掉。四个检查都过再沿用第三章的方法改字体、调比例尺出一版更新后的TIF成图。6. 用脚本把属性表变成交付清单坐标、面积一次性导出数据更新完最终交付往往不是SHP而是一张坐标表、面积表或图斑清单。每次手工复制属性表导出Excel格式不一样不说坐标还是6位小数领导看不懂。日常我会用一个批处理脚本把SHP的关键字段连同质心坐标一次性导出成CSV再用它做汇报表。import arcpy import csv arcpy.env.workspace rD:\gis\data\reserve shp 保护区边界_增量.shp out_csv 保护区清单.csv sr_wgs84 arcpy.SpatialReference(4326) with arcpy.da.SearchCursor( shp, [名称, 级别, 类型, 面积KM2, SHAPEXY], spatial_referencesr_wgs84 ) as cursor: with open(out_csv, w, newline, encodingutf-8-sig) as f: writer csv.writer(f) writer.writerow([名称, 级别, 类型, 面积KM2, 经度, 纬度]) for name, level, typ, area, xy in cursor: writer.writerow([ name, level, typ, round(area, 2), round(xy[0], 6), round(xy[1], 6) ])这里两个细节是经验积累出来的。第一spatial_referencesr_wgs84参数强制输出坐标为WGS84经纬度不管SHP内部是什么投影导出结果都是统一的经纬度格式后续丢进Excel或数据库都好处理。第二文件写入时用encodingutf-8-sig而不是utf-8为此我翻过一次车直接用utf-8写出的CSVExcel双击打开中文全是乱码改成utf-8-sig后Excel能识别文件头并自动切到UTF-8编码。这两行参数看着不起眼实际交付时能省掉大量来回沟通。脚本跑完生成的CSV通常还要追加两列保护区所在地级市名称和批复年份这个信息属性表里不一定有可以用名录表再关联一次。关联之后继续导出CSV最终作为年度任务台账的附件或者配合ArcGIS Online地图服务发布时使用。数值单位再确认一次面积字段如果来自AddGeometryAttributes的AREA_GEODESIC测地面积单位是平方公里导出后不需要缩放。如果是平面面积计算工具算出来的很多默认单位是平方米导出前先除以1,000,000避免汇总时面积偏大三个数量级。我早期接过一个数据整理任务当时刚理解即可想当然认为SHP源文件里的面积字段已经正确结果没验证单位就导出了汇表汇报时发现某省保护区总面积比官方发布的数据翻了几十倍场面一度很尴尬。后来养成的习惯是任何第三方SHP提供方的面积字段先抽查五个图斑的量算值确认单位后再脚本化导出。从那以后我每次交付数据前都强制走一遍这套流程验坐标、查要素基数、复算面积、导出CSV后再手工抽核两行。前期多花十分钟后面少解释两小时希望这套流程和踩坑记录能帮你少走一段弯路。本文还有配套的精品资源点击获取