DEM栅格文件提取带高程属性的点文件——ARCGIS

格式转化的目的

VISUAL MODFLOW FLEX（VMF）中构建模型时，需要地表高程文件（surface），文件导入格式可为 .DEM，.GRD，.TXT，.ASC等；也可以通过point文件，直接在VMF中创建surface文件，主要需要点的X，Y，Z数据。

但本学渣目前下载到的DEM栅格文件多为.img或.tif格式的文件，这里只会通过ARCGIS提取DEM栅格数据中点的XYZ信息。【如果大家有另外的好办法，请一定告诉我！】

DME数据下载地址

1. 根据省、市下载：http://gaohr.win/site/special/2019/2019-09-14-dem-30m-cn.html 

2. 根据经纬度下载（地理空间数据云）：http://www.gscloud.cn/sources/accessdata/aeab8000652a45b38afbb7ff023ddabb?pid=302

打开ARCGIS中的3D Analyst、Spatial Analyst扩展工具

1. 若未打开扩展工具，在使用这两个功能时会出现如下提示。

2. 解决方法：自定义→扩展模块→选中3D Analyst、Spatial Analyst，关闭窗口即可。

DEM栅格文件提取带高程属性的点文件——ARCGIS

1. 在ARCGIS导入DEM文件：.img或.tif。（可能也可以不导入🤣）

2. 栅格转为多点文件：ArcToolBox→3D Analyst工具→转换→由栅格转出→栅格转多点。（此过程需要一段时间，转化成功后，会在内容列表中显示该文件，格式为.SHP，文件大，点无法编辑）

3. 要素转点：ArcToolBox→数据管理工具→要素→要素转点。（此过程同样需要一段时间，转化成功后，会在内容列表中显示该点文件，格式为.SHP，文件小，点可编辑）

4. 点属性增加高程项：ArcToolBox→Spatial Analyst工具→提取分析→值提取至点。（输入点要素中选择上一步生成的shp文件，输入栅格中选择原始DEM文件，完成操作后，右键该文件→打开属性表，此时已有高程属性RASTERVALUE）

5. 点属性增加XY（经纬度）：ArcToolBox→数据管理工具→要素→添加XY坐标。（右键该文件→打开属性表，此时已增加XY坐标值）

6. 将点属性复制出来，创建TXT文件粘贴进去保存即可拥有XYZ的point文件。

注：可能会出现点数量过少、密度过低的情况，可以在“生成多点文件”步骤之前，通过“重采样”功能（ArcToolBox→数据管理工具→栅格→栅格处理→重采样），降低输出像元大小XY，重采样技术选择CUBIC，确认即可。（另外，建议在此之前先将DEM文件裁剪到所需大小，否则后续操作会导致文件占用内存飙升，甚至软件出现闪退的情况。）