栅格单元是DHSVM运行的计算单元,也是河道数据、初始状态数据等驱动数据建立的基础,本文将简单介绍DHSVM栅格驱动数据的制作方法。
一、所需工具
- Arcmap 10.0版本以上
- 二进制转化工具 Convent.exe,DHSVM官网可下载
二、流域栅格数据制作
1.DEM数据
(1)将DEM数据导入到Arcmap,裁剪掉目标流域外多余的部分。利用重采样工具,采样至想要的分辨率。
2.坡向数据
(1)打开ArcToolBox,Spatial Analysis 工具——水文分析——填洼,将DEM数据的洼地填充好。
(2)Spatial Analysis 工具——水文分析——流向,输入填洼后的DEM,得到坡向数据。一共8个方向,取值从1-128。
3.MASK数据
(1)根据Arcmap水文分析中的流域提取方法,提取流域范围:填洼——流向——流量——捕捉倾泻点——分水岭。具体过程可以直接百度,这里就不详细说了。
(2)制作MASK数据,将得到的流域范围栅格重分类,流域内部区域为1,流域外部的Nodata区域重分类为0,得到MASK数据。
4.土地类型、植被类型数据
(1)导入下载好的土地类型、植被类型栅格数据,根据DEM数据裁剪、重采样,输出像元大小选择与DEM数据相同。注意类型数据重采样时需要选择最邻近法。
(2)重采样后的数据进行重分类,将每个类型的像元值改成1,2,3,,,这样的,方便输入到DHSVM中。
根据此步骤可获取DEM、MASK、坡向、植被类型、土壤类型 5个栅格数据,土壤深度数据后期在制作河道数据时获取。
三、栅格数据转二进制数据
DHSVM支持两种栅格数据输入格式:二进制和NetCDF。NetCDF制作比较麻烦,所以本文只介绍怎么制作二进制栅格数据。
1.导出栅格数据
(1)ArcToolBox——转换工具——由栅格转出——栅格转ASCII,将5个栅格数据转为文本数据。
(2)打开转换完成的ASCII文本,将最上方的描述信息删掉。描述信息中的行列数等最好记下来,后面会用到。
2.ASCII转二进制
(1)将二进制转换工具convent.exe放到和栅格数据同一个文件夹,打开命令控制符(cmd)或PowerShell,输入convent.exe,回车运行,可以看到convent.exe的用法。用法为:convent.exe 原始数据类型 目标类型 原始数据文件名 目标文件名 行号 列号。
(2)以DEM为例,原始类型为ascii,目标类型为float,输入convent.exe ascii float dem.txt dem.bin 324 468
可以获得一个名叫dem.bin的二进制文件,这就是DHSVM需要的栅格输入数据。
(3)根据这个方法将其他栅格数据也转为二进制格式,注意DEM的目标格式为float,其他数据的目标格式都为char
以上就是DHSVM驱动数据中栅格数据的制作方法,希望对大家有帮助。河道数据制作需要用到DEM和MASK数据,我们下回再说。