简介：本资源为全国108288个小型矿产点位的高精度矢量SHP数据集，面向地质勘查、矿产资源管理、GIS空间分析及环境评估等领域的科研人员、高校师生与行业从业者，解决矿产空间分布建模、资源潜力评价与规划决策中基础地理数据缺失问题。压缩包共8个文件（6.8MB），含.shp（几何图形）、.dbf（属性表，含矿种、规模、成因类型、利用现状、地质工作程度等12项字段）、.prj（坐标系定义）、.shx（索引）及.cpg（编码声明）等核心组件，确保数据可直接导入ArcGIS、QGIS等平台完成叠加分析、统计制图与空间查询。已有74人学习下载，数据结构规范、字段完整、坐标统一，支持从点位定位、属性筛选到成矿规律挖掘的全流程应用，是开展区域矿产资源数字化评估与可持续开发研究的重要基础支撑。 干GIS这行，跟点数据打交道是家常便饭。前阵子帮一个做资源评估的团队处理全国小型矿产点位矢量数据，文件里有108288个点，每个点都是一处小型矿产地的空间位置。刚拿到手时我第一反应是：数据量不算大，但牵扯到的格式认知、坐标系校验、属性表清洗、后续三维化与格式转换这一整套流程，足够写一篇像样的实操笔记了。这篇文章就围绕“全国小型矿产点位矢量SHP数据”的获取、解析、清洗和转换展开，覆盖SHP打开、shp转3dtiles、shp转txt、批量shp转CAD这些高频需求。
适合看这篇内容的人有三类：一是刚接触GIS、手里只有SHP文件却不知道怎么下手的初学者；二是做矿产资源分析、国土空间规划、生态评价的从业者，需要把点位数据真正用起来；三是卡在“点位数据怎么转3dtiles、怎么转CAD”这种具体操作里的实操派。读完你可以照着流程走一遍，从裸数据做到可发布、可汇报的成果。
1. 看懂这批点位数据：108288个小型矿产点是什么
1.1 十万级点位数据的分量
108288不是随便凑出来的数字，它代表的是全国范围内已经登记、圈定的小型矿产地的数量级。这里“小型”是相对大型、中型矿产地而言的储量规模口径，小型矿产点通常出现在地方矿业权登记、矿产资源国情调查、储量库更新等数据源里。十万个点听起来不多，但当你把它画到全国地图上，再叠加行政区、流域、交通路网这些图层，它就是能支撑宏观决策和区域分析的基础数据。
这类数据的典型用法有几类：
- 按行政区统计小型矿产点的分布密度，判断区域资源禀赋。
- 与地质背景、断裂构造等图层叠加，做成矿规律分析与找矿预测。
- 为外业调查编制核查点位清单，按经纬度逐个跑现场。
- 在三维地球或数字孪生场景里做点位可视化展示。
我在实际项目里遇到最多的是前两类。尤其在做矿产资源规划时，十万个点位的空间分布直接决定了哪些区域需要重点管控、哪些区域适合集约开发。没有这份底图，规划方案基本就是拍脑袋。
1.2 “小型矿产”的数据口径与字段含义
在资源储量分类体系里，矿产地的规模一般按矿种不同有对应的储量或资源量阈值。比如煤矿、铁矿、铜矿的“小型”标准并不一样，所以点数据本身不直接告诉你“它有多小”，而是通过属性字段来区分。
拿到一份小型矿产点位SHP，属性表里大概率会出现这些字段：
| 字段名 | 含义 | 典型值示例 |
|---|---|---|
| NAME / MC | 矿产地名称 | 某某铜矿、某某铅锌矿 |
| MINERAL / KZ | 主要矿种 | 煤、铁、铜、铅锌、金等 |
| SCALE / GM | 规模等级 | 大型、中型、小型、矿点 |
| PROVINCE / CITY | 所在行政区 | 云南省、个旧市等 |
| LONGITUDE / X | 经度或投影X坐标 | 102.35（十进制度） |
| LATITUDE / Y | 纬度或投影Y坐标 | 24.36 |
| STATUS | 开采状态 | 开采、关闭、探矿权 |
| SOURCE | 数据来源 | 储量库、国情调查等 |
我见过的数据里，有的把经纬度直接给了十进制度，有的是度分秒字符串，还有的是高斯投影坐标。这一点在数据清洗时必须先确认，否则后面所有分析都是错的。
1.3 这批数据的典型应用场景
一个基于十万点位数据的项目，通常是怎么推进的？拿我最近经手的案例来说，团队要做的是“某区域小型矿产开发生态影响评估”。流程大概是：先把全国点位按省份字段筛选出目标区，再叠加县域轮廓、流域边界、保护区范围，做空间连接和缓冲区分析，最后按矿种分类导出专题图。
这个过程中，你需要的不只是点位本身，还经常要配合其他公开数据：县域轮廓SHP、流域边界SHP、气象站矢量点、道路交通线等等。很多人搜“云南省县域轮廓shp”“珠江流域shp”，本质上都是在做同一套空间分析工作流。点位数据只是起点，后续的叠加、筛选、可视化才是重头戏。
2. SHP文件的结构与坐标系：先搞懂底层逻辑，再动手
2.1 shapefile不是“一个文件”，而是一组文件
很多新手拿到SHP时，只看到一个“xxx.shp”图标，以为它是一个单文件，拷贝的时候只拷走了这个文件，结果到别的电脑上打开就报错。这是个非常经典的问题。shapefile本质上是ESRI定义的一种矢量数据存储格式，它由一组同名前缀的文件共同组成。
一个完整的shapefile至少包含：
| 扩展名 | 作用 | 缺失后果 |
|---|---|---|
| .shp | 几何信息：点、线、面的坐标 | 无法显示图形 |
| .shx | 几何索引，加速读取 | 数据打不开或读取极慢 |
| .dbf | 属性表，存字段值 | 只有图形，没有属性 |
| .prj | 坐标系定义文本 | 系统不知道点在哪 |
| .cpg | 字符编码描述 | 中文属性可能乱码 |
另外还可能有.sbn、.sbx、.xml、.lock等辅助文件，比如.sbn和.sbx是空间索引，在ArcGIS里编辑过之后会出现。拷贝数据时最好把整个文件夹打包带走，别只拖那个.shp。
为什么.cpg文件缺失会导致中文乱码？因为.dbf文件本身不带编码声明，系统只能靠猜。.cpg里写的是UTF-8还是GBK，决定了ArcGIS或QGIS用什么解码方式去读属性。很多从老系统导出的数据只有GBK编码的.dbf，没有.cpg，直接用新版ArcGIS Pro打开就会看到一堆乱码。处理方法后面第6节会讲。
2.2 坐标系问题：为什么你的点会跑到海里
点数据跑偏是矢量操作里最让人头疼的问题。全国范围的小型矿产点位，最合理的坐标系是WGS84或CGCS2000的地理坐标系，也就是经纬度。但也有不少历史数据用的是北京54、西安80，或者高斯-克吕格投影。一旦.prj文件缺失或写错，软件默认按WGS84经纬度解析，点位就会偏到几百上千公里外，甚至掉进太平洋。
我处理的一个典型场景：某团队拿到的SHP打开后点全在海上，排查发现.prj文件描述的是CGCS2000高斯投影，但属性表里的X、Y字段其实是经纬度。也就是说，数据本身没问题，是坐标系定义和实际几何坐标对不上。解决方式很简单，在ArcGIS里用Define Projection重新指定为WGS84地理坐标系即可。
所以拿到任何SHP的第一步，永远是看三个东西：
- .prj文件是否存在，里面写的什么坐标系。
- 图层范围是否合理，经纬度范围应该在73°E到135°E、18°N到54°N之间。
- 属性表里的X、Y字段和几何坐标是否一致。
2.3 属性表：打开数据后第一件要做的事
不要急着画图渲染，先打开属性表预览一遍。点开之后重点看几项：记录总数是否为108288，能否对上文件名里的数量级；字段名称是否清晰，有没有乱码；关键字段是否有大量空值；经纬度字段是数值型还是文本型。
如果字段是文本型，后面做坐标换算、空间计算时会非常麻烦。比如“102.3567”存成了文本，在ArcGIS里做数值运算之前得先转成双精度。QGIS里可以用Field Calculator转，ArcGIS里可以用Convert Field工具。遇到这种情况别嫌麻烦，第4节的清洗流程里我会详细说明。
3. 数据来源与获取：哪些渠道靠谱，哪些坑要避
3.1 官方资料库与公开数据平台
“全国108288个小型矿产点位矢量SHP数据下载”这类需求，我第一反应永远是：数据本身大概率来自国土或地质系统的公开成果。普通用户拿不到涉密坐标，但经过脱密处理的矿产地信息，在官方平台是可以申请或下载的。
靠谱渠道我实际用过这些：
- 全国地质资料馆与地质云平台：提供地质图、矿产地数据库等公开成果，部分数据需要注册申请。
- 各省自然资源厅地质资料馆：可以下载省级矿产地质相关的矢量数据。
- 全国矿业权人勘查开采信息公示系统：能查到每个矿业权的位置、矿种、状态，可结合坐标批量生成点位。
- 国家地球系统科学数据中心、地理空间数据云：地学领域开放数据的老牌平台，数据规范程度高。
这类官方渠道下载的数据一般带完整.prj、.cpg文件，坐标系规范、属性表结构清晰，后期处理省很多事。
3.2 学术共享与开放数据平台
除了官方渠道，学术社区和开放数据平台也是重要补充。比如OSGeo、OpenStreetMap能提供基础地理底图，但不是矿产点位数据的主力来源。科研项目发布的成果数据（比如某些国家重点研发计划共享数据）经常包含矿产地、灾害点、监测井等点位信息，可以用作参考或交叉验证。
需要提醒的是，第三方平台下载的数据质量参差不齐。有的数据直接从纸质图件数字化而来，坐标精度偏差可能达到几百米；有的属性表字段是拼音缩写，没有说明文档；还有的干脆缺少.prj文件。在学术研究里，这类数据只适合做趋势分析，不适合做高精度工程判断。
3.3 授权与合规：别在数据源上翻车
这是很多人在兴奋下载数据时忽略的一点。矿产地坐标数据，无论来自官方还是第三方，都可能涉及数据授权和版权问题。官方数据通常会有“使用者需注明来源”“不得用于商业用途”之类的限制条款；第三方爬取整理的数据更是处于灰色地带。
实操中我的建议有三条：
- 优先使用官方公开数据和有明确授权协议的共享数据。
- 如果项目涉及对外交付或商业用途，务必确认数据授权范围，必要时走正式申请流程。
- 不要把自己下载的数据二次打包传播，避免法律风险。
数据合规这个事，在朋友圈里可能没人提，但真出了问题就是大问题。你辛辛苦苦做了几个月的分析，如果数据来源不干净，最后成果是站不住脚的。
4. 实操：从打开到清洗，把十万点位整理成可用状态
4.1 在ArcGIS Pro / QGIS中正确打开SHP
ArcGIS Pro打开SHP很简单：打开软件后，在“目录”面板中找到文件夹位置，直接右键SHP文件选择“添加至当前地图”，或者用“添加数据”按钮浏览到文件位置。QGIS更直接，把.shp文件从文件夹拖到QGIS窗口里就能加载。
但有个细节，如果SHP文件所在路径包含中文，某些版本的ArcGIS Pro可能会识别异常。我习惯把所有输入数据统一放到英文路径下，比如D:\data\mineral_points\，文件命名也用英文或拼音。这一步能规避很多莫名奇妙的报错。
打开后第一件事不是制图，而是右键图层查看“属性”，找到“源”或“坐标系”选项卡，确认坐标系统。如果看到“未知坐标系”或者和预期不符，立即停下处理，别继续分析。
4.2 坐标校验：把异常点揪出来
即使坐标系正确，十万个点里也可能藏着“脏点”。比如坐标值写反了（经纬度互换）、坐标超出中国范围、重复点位、坐标精度明显异常（比如小数点后只剩一位）。我处理这批数据时先跑了一遍范围检查，发现其中有几十个点的经度在0到70之间，明显不对，很可能是属性表录入时经纬度填反，或者数据源本身存在异常值。
一个通用排查逻辑：
- 先看整体的几何范围，打开图层属性里的“范围”，确认大致在73°E到135°E、18°N到54°N。
- 用属性表筛选经纬度字段，找出超出合理范围的点。
- 使用QGIS的“按位置选择”或ArcGIS的“Select by Attributes”，把异常点单独导出，人工抽查。
- 确认异常点无法修复的，标记并剔除，不要把脏数据留到分析阶段。
4.3 点位去重与字段规整
十万点级别的数据，去重是绕不开的。重复可能来自数据源合并时未做处理，同一矿产地存在多条记录。去重思路有两种：一是按矿产地名称字段去重，二是按坐标位置去重。
按名称去重比较简单，直接在ArcGIS Pro里用“删除相同项”工具，选择NAME字段。按坐标去重更稳妥，因为同一个矿产地可能有多种叫法。可以在属性表里增加一个“X_Y”文本字段，把经度和纬度拼接起来，再用删除相同项按这个字段去重。
字段规整也是重要一步。常见问题包括：
- 经纬度字段是文本类型，需要转成双精度。
- 度分秒格式（如102°30′20″）需要转成十进制度。
- 矿种字段存在多种写法（如“铜矿”“Cu”“铜”），需要统一映射。
这些操作都可以在ArcGIS Pro的字段计算器或QGIS的字段计算器里完成。给一个简单的字段计算器示例，把度分秒转十进制度：
-- 假设原始字段"LON_DMS"格式为 "102°30'20\"" -- 转成十进制度 -- ArcGIS Python解析器 def dms2dd(dms): d, m, s = dms.replace('°',' ').replace('\'',' ').replace('"','').split() return float(d) + float(m)/60 + float(s)/3600在字段计算器里调用dms2dd(!LON_DMS!)即可。
5. 实操：shp转3dtiles、shp转txt、批量shp转CAD
5.1 shp转3dtiles，让十万个点飞上三维地球
“shp转3dtiles”是这两年特别高频的需求，尤其是配合Cesium、超图或者自研数字孪生平台做三维展示。点位数据转3dtiles，本质上就是把十万个经纬度坐标按四叉树规则切块组织成3D Tiles格式，浏览器端按需加载，避免一次性渲染十万个点卡死浏览器。
工具选择上，我常用两条路线：
- 路线一：商业软件，比如CesiumLab，导入SHP或GeoJSON，选择“点数据”，设置高度来源（比如从属性字段读高度，或者按固定高度抬升），输出3dtiles。操作路径短、出成果快，适合赶工项目。
- 路线二：开源免费方案，比如使用QGIS先把SHP导出为GeoJSON，再用lonlat2tiles或者基于Three.js的3dtiles工具转换。
具体到CesiumLab的操作，大致是：新建转换任务，选择SHP文件，数据分类选“点”；在空间参考设置里确认坐标系统；输出格式选“3D Tiles点云”；然后设置每个Tile的节点容量，一般点数据可以设5000到10000个点一个节点；最后生成。生成后在Cesium里加载tileset.json就能看到点位。
十万点级别的数据，在三维场景里默认样式是小白点，太小看不清，太大遮挡地形。我的经验是：先按矿种字段做分类着色，属性相同矿种用同一颜色，pointSize根据视点距离动态调整，或者用Cesium的pointPrimitive批量绘制。如果只是做“看个分布”的演示，可以把点转成小圆柱体或billboard图标，效果更直观。
5.2 shp转txt，外业设备和轻量应用都靠它
外业核查、无人机航线规划、或者给不装GIS软件的业务同事导数据，都需要把SHP转成最朴素的txt或CSV。这一步看似简单，但坑也不少。
方法一：ArcGIS属性表导出。打开属性表，选择所有记录，点击“导出”，可以选择文本文件，字段由自己勾选。这个办法最直接，但导出的经纬度默认精度可能不是特别高，需要设置字段格式。
方法二：脚本处理。如果是命令行玩家，用GDAL/OGR一行命令就能转：
ogr2ogr -f CSV points.csv mineral_points.shp -lco GEOMETRY=AS_XY这个命令会把几何坐标输出到X、Y两列。如果原始数据是投影坐标，想要经纬度，可以用-t_srs EPSG:4326参数先转换坐标系统再输出：
ogr2ogr -f CSV points_4326.csv mineral_points.shp -lco GEOMETRY=AS_XY -t_srs EPSG:4326如果你习惯用Python，直接用ogr读取并写出文件也很方便。脚本示例我放到下面：
from osgeo import ogr shp_path = r"D:/data/mineral_points.shp" txt_path = r"D:/data/mineral_points.txt" ds = ogr.Open(shp_path) lyr = ds.GetLayer(0) with open(txt_path, "w", encoding="utf-8") as f: f.write("fid,longitude,latitude,name,mineral_type\n") for feat in lyr: geom = feat.geometry() if geom is not None and geom.GetGeometryName() == "POINT": x = geom.GetX() y = geom.GetY() name = feat.GetField("NAME") mtype = feat.GetField("MINERAL") f.write(f"{feat.GetFID()},{x:.6f},{y:.6f},{name},{mtype}\n") ds = None print("done")注意：上述脚本读取的坐标是数据源自身的坐标系。如果数据源是CGCS2000投影坐标，想要WGS84经纬度，需要先做坐标转换，不能直接输出。
5.3 批量shp转CAD，给设计院和评审单位交底
做矿产资源规划或工程选址时，经常要把点位数据提交给设计院或评审专家，对方用的软件是AutoCAD而不是GIS。这时候需要把SHP转成DXF或DWG。
单个文件转换很简单，ArcGIS里有“导出为CAD”工具，选SHP输入，指定输出类型DXF/DWG就行。麻烦的是批量。举个例子，你需要按行政区分省导出，或者按矿种分类导出几十个SHP再转换。这时候用ArcGIS的“批量导出为CAD”工具最合适，它接受多个输入要素类，统一转换为一个或多个CAD文件。
我的批量做法是基于ArcGIS Pro的“地理处理-批量导出为CAD”，操作流程：
- 准备好所有需要转换的SHP文件，放在同一个文件夹。
- 打开“批量导出为CAD”工具，把文件夹下所有SHP拖入输入列表。
- 选择输出版本（DWG 2010还是DXF 2018），设置输出路径。
- 运行后得到单个或多个CAD文件，每个输入图层对应一个CAD图层名。
需要注意的坑是单位问题。CAD默认单位是毫米，而SHP里的坐标单位是米或度。如果直接把经纬度导入CAD，点位会显示在接近（0,0）的位置，因为经纬度数值太大，CAD视图无法自动缩放。我过去常用的处理方法是：先把SHP投影成高斯平面坐标（米制），再利用CAD的“地图图块定义”功能把坐标转换写进DWG文件，这样打开图纸时坐标位置才是对的。
如果是给外单位的设计图，建议先转成米制投影坐标，比如CGCS2000 3度带高斯投影，再导出DXF。这样至少在CAD里打开后能看到坐标落在正确区域。
6. 常见问题与避坑清单
6.1 高频问题速查表
这些年处理SHP点位数据，我遇到的高频问题基本可以汇总成一张表：
| 现象 | 可能原因 | 排查与解决 |
|---|---|---|
| SHP打开后只有点没有属性 | .dbf文件缺失或损坏 | 检查同目录下是否有.dbf文件，从源数据重新拷贝 |
| 属性表中文乱码 | 缺少.cpg或编码声明错误 | 在QGIS里设置图层编码为UTF-8或GBK，或者ArcGIS里设置读取编码 |
| 点位跑到海上或国外 | 坐标系定义错误 | 用Define Projection工具重新指定坐标系，检查经纬度字段是否反了 |
| 属性表记录数少了几百条 | 数据源本身存在坐标过滤或查询被误用 | 检查图层属性里的“定义查询”（Definition Query），看是否误加过滤条件 |
| shp转3dtiles后点看不见 | 高度设置错误或节点容量过大 | 检查高度字段是否为空，调小节点容量，调整点样式 |
| 导出CSV后坐标变成了一长串整数 | 数据源是投影坐标，没有转经纬度 | 用ogr2ogr加上-t_srs EPSG:4326参数重新导出 |
| dbf用Excel打开时中文乱码 | Excel默认按ANSI解析，不支持UTF-8编码的dbf | 不要直接用Excel编辑dbf，用DB Browser或QGIS编辑属性表 |
这张表背后对应的问题我几乎每个都踩过。特别是“投影坐标导出了整数坐标”这个问题，最容易出现在把CGCS2000高斯投影的SHP用ogc2ogr直接转CSV的时候。解决方式很简单，先转坐标系再导出即可。
6.2 几个实操里容易忽略的细节
第一个细节：dbf文件不要用Excel直接编辑保存。Excel打开dbf后如果另存为xlsx再转回dbf，字段类型、精度都可能被破坏。十万人里不一定有一个碰到这个问题，但碰到了就是灾难。我见过有人用Excel改了属性表，保存后图层打不开，最后靠备份才找回来。改属性表，请老老实实打开ArcGIS或QGIS的编辑模式。
第二个细节：检查数据时不要只看“点数”，要看“有效点数”。108288个点里可能有重复点、空几何点，也可能有坐标为零的异常点。统计“有效点数”通常比文件标题里的总数更重要。用QGIS的“矢量-几何工具-检查几何”能快速找出空几何和无效几何。
第三个细节：做全国尺度分析时，注意带地理坐标系还是投影坐标系。十个省份范围内的空间计算，用米制投影坐标可以保证距离准确；但如果只是看图，经纬度就够了。我默认的规范是：空间分析用投影坐标，发布服务用经纬度坐标，各取所长。
6.3 最后分享一个我个人的处理习惯
拿到任何一份SHP点位数据，我固定走四步：先看投影，再查属性，然后做范围检查，最后才进入可视化或转换环节。这个流程看起来慢，实际是效率最高的方式。很多新手拿到数据就急着配色出图，结果发现坐标系错了、属性少了，返工两三次，浪费的时间远比先检查多得多。
如果你手头正在处理类似“全国108288个小型矿产点位”这种体量的数据，我建议你第一遍先别做任何分析，单纯把属性表完整过一遍，把所有字段含义和空值情况搞清楚。这个动作花不了半小时，但能让你在后面所有分析里都底气十足。数据是别人的，但清洗完的数据是自己的。十万个点落在哪里、代表什么矿种、状态是什么，都心里有数之后，出图出表就只是时间问题了。
本文还有配套的精品资源，点击获取
