lisflood地形DEM数据处理

1.DEM下载

分别尝试用12.5m和30m的DEM转化成ASCII文件，发现12.5米生成的网格过多可能导致运行速度太慢，因此使用30mDEM进行研究，DEM数据下载链接：

https://pan.baidu.com/s/17VZrzprLUI0hmG__W4Hpdw?pwd=ompt
提取码：ompt
哥白尼全国DEM

下载得到的为.tif文件，需要进行以下两步转换才能得到所需的ASCII文件。

2.修改坐标系（将经纬度坐标转换为投影坐标即可，对使用哪一种投影坐标没有限制）

转换工具：ArcToolbox-数据管理工具-投影和变换，若为矢量图层则可以直接使用投影工具，若为栅格图层则需要进一步点开栅格，使用投影栅格进行坐标系转换。（适用于已知坐标系正确的情况下，转换为另一个坐标系）。

注：定义投影适用于没有坐标系或者坐标系错误的情况下，该功能对矢量图层和栅格图层均适用。

以栅格图层为例，将WGS84转为WGS_1984_World_Mercator投影，点击倒三角选择输入栅格文件，再选择输出坐标系，这里没有选择CGCS2000，因为会提示无地理（坐标）变换（不同椭球体之前的转换），目前暂未找到解决方案，故先转换为WGS_1984_World_Mercator投影，不需要选择地理（坐标）变换。注意命名的时候要在文件后面加.tif，否则可能生成arcmap 打不开的文

 点击确定即可在右下角看到以m为单位的xy坐标。（点击右上角切换成html编辑器即可缩放图片大小）：

将WGS84转为Gauss Kruger CGCS2000投影：

• 先将WGHS84地理坐标系转为CGSC2000地理坐标系，在投影栅格中选择  地理坐标系》Asia》GCS_China_Geodetic_Coordinate_System_2000
• 再将CGCS2000地理坐标系转为CGCS2000高斯克吕格投影坐标系，在投影栅格中选择  投影坐标系》Gauss Kruger》CGCS2000》CGCS2000_3_Degree_GK_CM_114E具体在哪个分度带根据经纬度坐标判断；

 3.剪裁并导出为ASCII文件在投影栅格中选择  地理坐标系》Asia》GCS_China_Geodetic_Coordinate_System_2000

（1）QGIS做法

• • • 加载图层：将.tif文件拖入左边即可看到图层，在右边图层样式里可以切换“单波段灰度”为“单波段伪彩色”，方便查看。
    • 在工具栏选择：栅格-提取-按范围剪裁栅格，在弹出的窗口中选择输入图层，剪裁范围点击右边▼，选择在地图画布中绘制，即可开始绘制，若想重新绘制在重复上述操作即可。同时设置[将指定的无数据值分配给输出波段]为-9999（若遗漏可通过下一步进行），在[保存为临时文件]处点击▼，选择保存到文件，输入文件名即可，然后点击运行，即可得到标准格式的asc文件。

 

• • • 在图层中选择栅格-转换-转换（格式转换）：设置[将指定的无数据值分配给输出波段]为-9999，再同上保存到文件。即可得到标准格式的asc文件，格式如下。

 

 注：有时候转出来nodata值不是-9999，而是其他值，这时候可以采用arcgis转，目前不会产生错误

（2）arcgis做法（推荐）

1. 1. 先剪裁，搜索“按掩膜提取”，在弹出的对话框中设置输入栅格为需要剪裁的文件，第二栏为想要剪裁的区域文件，第三栏保存路径。
   2. 在搜索中搜索“由栅格转出”→“栅格转ASCII”，在弹出的对话框中，将输入栅格设置为需要转换的文件，在第二栏设置保存位置，即可得到转换后的数据如下：

反之，asc转栅格，则搜索“转为栅格”→“ASCII转栅格”，同样设置输入文件、保存路径，第三栏不用变即可。但是转换完一般是没有坐标的需要“定义投影”。