若转发本文章,请添加引用
前言
最近正好下载地层数据,刚开始以为需要手动矢量化,即下载tif后地理配准进行矢量化,研究区较小还好,若区域较大那将会是一个较大的工程,偶然有现成数据可以利用,因此整合网络资源教程给同学们节省时间。
一、地层、岩性
谈论到地层与岩性,需要知道这两个是不太一样的概念,虽然在一定程度上有关联。岩性一般诸如火成岩、沉积岩、变质岩、岩浆岩,地层的话一般按照时代划分,能够分为宇、界、系、统、阶,这其中就能包括岩石特性。
这里我们使用两个网站:
全国地质资料馆https://www.ngac.cn/125cms/c/qggnew/index.htm
地质科学数据出版系统:http://dcc.ngac.org.cn/cn/page/index
在全国地质资料馆可以查到不同比例尺的地质图,比如1:5W,1:20W,1:25W等,然后再地质科学数据出版系统中找到对应区幅进行下载(由于地质科学数据出版系统中公开了1:20W的资料,因此以1:20W地质图为例)
二、数据下载
1.确定数据地质图幅号
首先在全国地质资料馆—数据与服务—数字资料馆(地质图数据库)找到对应分比例尺地质图
然后打开网站:http://ngac.org.cn/Distribute/20wanQD.htm,根据研究区经纬度确定图幅及图幅号,比如汶川地区图幅号为H-48-02。然后对其进行检索。检索方式记得点击地址图库。
然后点击进入即可在线游览,利用WTMS服务也可以加载到arcgis中游览,不过速度会比较慢。
2.下载数据
在地质科学数据出版系统中找到全国1∶20W数字地质图(公开版)空间数据库http://dcc.ngac.org.cn/cn//geologicalData/details/doi/10.23650/data.A.2019.NGA120157.K1.1.1.V1,然后下滑找到对应地质图幅号,需要注意这里需要注册账号,且一个账号每天只能从该网站下载五个数据,若论文使用需要仔细阅读添加引用。下载解压后有三个文件夹。
3.数据转换
该数据格式为mapgis格式,下一步需要将其转化为我们常用的shp格式。
Mapgis转换方法:
(1)Mapgis软件转换(但是目前破解版Mapgis67较多,自己也尝试后发现转换后属性数据部分会有偏差,有些显示错误,改善方法是使用高版本Mapgis)
(2)使用FME软件转化。
这里用第二种方法进行转换,参考https://zhuanlan.zhihu.com/p/371786900,使用FME2020+MyFME插件进行转换(MyFME为必须,否则无法识别Mapgis格式)
按照文章及压缩包内教程后,打开FME Workbench:(记得修改时间!)
网盘链接
1.右上角文件—新建—空工作空间
2.工具栏选择读模块,数据集选择解压数据后MAPGIS文件夹下的所有数据。
3.工具栏选择写模块,格式选择esrishapeflie ,然后确定后勾选全选、排序。
4.连线,将对应的左右灰色箭头拉向红色箭头,最后点击运行即可,正常运行不会报错。
总结
最后打开arcgis查看效果,发现结果还不错,数据包括断层数据、地层数据、侵入岩数据等,较大程度上可以避免矢量化,但是转换时没有选择参考下,使用时需要定义投影CGS——WGS84。