分析和可视化DEM的工具。
# 概要
```
- 从任何GDAL支持的高程栅格生成阴影浮雕图:
gdaldem hillshade input_dem output_hillshade
[-z ZFactor (default=1)] [-s scale* (default=1)]"
[-az Azimuth (default=315)] [-alt Altitude (default=45)]
[-alg ZevenbergenThorne] [-combined | -multidirectional]
[-compute_edges] [-b Band (default=1)] [-of format] [-co "NAME=VALUE"]* [-q]
- 从任何GDAL支持的高程栅格生成斜率图:
gdaldem slope input_dem output_slope_map"
[-p use percent slope (default=degrees)] [-s scale* (default=1)]
[-alg ZevenbergenThorne]
[-compute_edges] [-b Band (default=1)] [-of format] [-co "NAME=VALUE"]* [-q]
- 从任何GDAL支持的高程栅格生成方面图
输出一个32位浮点栅格,像素值为0-360,表示方位角:
gdaldem aspect input_dem output_aspect_map"
[-trigonometric] [-zero_for_flat]
[-alg ZevenbergenThorne]
[-compute_edges] [-b Band (default=1)] [-of format] [-co "NAME=VALUE"]* [-q]
- 从任何GDAL支持的高程栅格生成颜色浮雕图
gdaldem color-relief input_dem color_text_file output_color_relief_map
[-alpha] [-exact_color_entry | -nearest_color_entry]
[-b Band (default=1)] [-of format] [-co "NAME=VALUE"]* [-q]
where color_text_file contains lines of the format "elevation_value red green blue"
- 从任何GDAL支持的高程栅格生成地形坚固性指数(TRI)地图:
gdaldem TRI input_dem output_TRI_map
[-compute_edges] [-b Band (default=1)] [-of format] [-q]
- 从任何GDAL支持的高程栅格生成地形位置索引(TPI)图:
gdaldem TPI input_dem output_TPI_map
[-compute_edges] [-b Band (default=1)] [-of format] [-q]
- 从任何GDAL支持的高程栅格生成粗糙度图:
gdaldem roughness input_dem output_roughness_map
[-compute_edges] [-b Band (default=1)] [-of format] [-q]
注意:
gdaldem通常假设x,y和z单位是相同的。 如果x(东西)
和y(南北)单位是相同的,但z(高)单位是不同的
scale(-s)选项可用于设置垂直单位与水平单位的比率。 对于
在赤道附近的Latlong投影,纬度和单位的单位
经度相似,海拔(z)单位可以转换为兼容性
通过使用scale
370400(如果海拔以英尺为单位)或scale
111120(如果海拔高度为米)。 对于不靠近赤道的地方,最好重新装备
网格使用gdalwarp之前使用gdaldem。
```
该实用程序有7种不同的模式:
~*~*hillshade~*~*
从任何GDAL支持的高程栅格生成阴影浮雕图
~*~*slope~*~*
从任何GDAL支持的高程栅格生成斜率图
~*~*aspect~*~*
从任何GDAL支持的高程栅格生成方面图
~*~*color-relief~*~*
从任何GDAL支持的高程栅格生成颜色浮雕图
~*~*TRI~*~*
从任何GDAL支持的高程栅格生成地形粗糙度索引的地图
~*~*TPI~*~*
从任何GDAL支持的高程栅格生成地形位置索引的地图
~*~*roughness~*~*
从任何GDAL支持的高程栅格生成粗糙度图
以下一般选项可用:
~*~*input_dem:~*~*
要处理的输入DEM栅格
~*~*output_xxx_map:~*~*
产生输出栅格
~*~*-of format:~*~*
选择输出格式。 默认值为GeoTIFF(GTiff)。 使用短格式名称。
~*~*-compute_edges:~*~*
(GDAL> = 1.8.0)在光栅边缘和近点数值处进行计算
~*~*-alg ZevenbergenThorne:~*~*
(GDAL> = 1.8.0)使用Zevenbergen&Thorne公式,而不是Horn的公式来计算斜率和方差。 文学研究表明,Zevenbergen&Thorne更适合平滑景观,而霍恩(Horn)在更粗糙的地形上表现更好。
~*~*-b band:~*~*
选择要处理的输入频带。 频带的编号从1。
~*~*-co "NAME=VALUE":~*~*
将创建选项传递给输出格式驱动程序。 可能会列出多个-co选项。 有关每种格式的合法创建选项,请参阅格式化的文档
~*~*-q:~*~*
抑制进度监视器和其他无错误输出。
对于所有算法,除了色彩浮雕外,如果在以每个源像素为中心的3x3窗口中发现至少一个设置为节点值的像素,则将发射目标数据集中的节点值。 结果是每个图像集周围将有一个1像素的边界值, 从GDAL 1.8.0,if -compute_edges被指定,gdaldem将计算图像边缘的值,或者如果在3x3窗口中找到一个nodata值,则通过内插缺失值。
~#~# 模式
~#~## hillshade
该命令输出8位光栅,具有很好的阴影浮雕效果。 这对于可视化地形非常有用。 您可以选择指定光源的方位角和高度,垂直夸张因子和缩放因子,以考虑垂直和水平单位之间的差异。
值0用作输出节点值。
以下具体选项可用:
~*~*-z zFactor:~*~*
垂直夸张用于预先乘以高程
~*~*-s scale:~*~*
垂直单位与水平的比例。 如果源DEM的水平单位是度数(例如Lat / Long WGS84投影),则如果垂直单位是米,则可以使用scale
111120(如果尺寸为英尺,则为scale
370400)~*~*-az azimuth:~*~*
光的方位角,以度数表示。 0,如果它来自栅格的顶部,90从东部,...默认值315很少会被更改,因为它是通常用于生成阴影图的值。
~*~*-alt altitude:~*~*
光线的高度,以度为单位。 90如果光从DEM上方开始,如果是耙光,则为0。
~*~*-combined:~*~*
(从GDAL 1.10开始)组合阴影,斜率和斜率阴影的组合。
~*~*-multidirectional:~*~*
(从GDAL 2.2开始)多向着色,从225度,270度,315度和360度方位角照明的山体阴影的组合。
多向山体阴影应用http:~/~/pubs.usgs.gov/of/1992/of92-422/of92-422.pdf的公式。
~#~## slope
该命令将采用DEM光栅并输出具有斜率值的32位浮点栅格。 您可以选择所需的斜率值的类型:度或百分比斜率。 在水平单位不同于垂直单位的情况下,您还可以提供比例因子。
值-9999用作输出节点值。
以下具体选项可用:
~*~*-p :~*~*
如果指定,斜率将以百分比斜率表示。 否则,它表示为度数
~*~*-s scale:~*~*
垂直单位与水平的比例。 如果源DEM的水平单位是度数(例如Lat / Long WGS84投影),则如果垂直单位是米,则可以使用scale
111120(如果尺寸为英尺,则为scale
370400)~#~## aspect
该命令输出一个32位浮点栅格,值为0°至360°,表示斜率面向的方位角。 方位角的定义是:0°意味着斜率面向北,90°面向东,180°面向南,270°面向西(前提是输入栅格的顶部为 北面)。 方面值-9999用作nodata值,以指示斜率= 0的平坦区域中的未定义方面。
以下具体选项可用:
~*~*-trigonometric:~*~*
返回三角角而不是方位角。 因此0°指东,北90°,西180°,南270°
~*~*-zero_for_flat:~*~*
对于斜率为0的平坦区域,而不是-9999,返回0
通过使用这两个选项,gdaldem方面返回的方面应与GRASS r.slope.aspect中的相同。 否则,它与Matthew Perry的aspect.cpp实用程序相同。
~*~*color-relief~*~*
该命令输出一个3波段(RGB)或4波段(RGBA)光栅,其值由高程和基于文本的颜色配置文件计算,其中包含各种高程值和相应希望颜色之间的关联。 默认情况下,给定高程值之间的颜色是平滑混合的,结果是一个很好的着色DEM。 可以使用-exact_color_entry或-nearest_color_entry选项来避免与颜色配置文件的索引不匹配的值的线性插值。
以下具体选项可用:
~*~*color_text_file:~*~*
基于文本的颜色配置文件
~*~*-alpha :~*~*
向输出栅格添加Alpha通道
~*~*-exact_color_entry :~*~*
在颜色配置文件中搜索时使用严格匹配。 如果没有找到匹配的颜色条目,将使用“0,0,0,0”RGBA四元组
~*~*-nearest_color_entry :~*~*
使用与颜色配置文件中最接近的条目相对应的RGBA quadruplet。
色彩浮雕模式是支持VRT作为输出格式的唯一模式。 在这种情况下,它会将颜色配置文件转换为适当的LUT元素。 请注意,指定为百分比的高程将被转换为绝对值,当源栅格的统计信息与构建VRT时使用的统计信息不同时,必须将其考虑在内。
基于文本的颜色配置文件通常包含每行4列:高程值和相应的红色,绿色,蓝色组件(0到255之间)。 高程值可以是任何浮点值,也可以是nodata值的nv关键字。高程也可以表示为百分比:0%是光栅中找到的最小值,最大值为100%。
可以为alpha组件添加额外的列。 如果未指定,则假定完全不透明度(255)。
接受各种字段分隔符:逗号,制表,空格,':'。
GRASS使用的常用颜色也可以通过使用其名称来指定,而不是RGB三元组。 支持的列表是:白色,黑色,红色,绿色,蓝色,黄色,品红色,青色,水色,灰色/灰色,橙色,棕色,紫色/紫色和靛蓝色。
由于GDAL 1.8.0,也支持GMT .cpt调色板文件(仅COLOR_MODEL = RGB)。
注意:颜色配置文件的语法派生自GRASS r.colors实用程序支持的语法。 ESRI HDR颜色表格文件(.clr)也符合该语法。 alpha组件和tab和逗号作为分隔符的支持是GDAL特定的扩展。
例如 :
```
3500 white
2500 235:220:175
50% 190 185 135
700 240 250 150
0 50 180 50
nv 0 0 0 0
```
~#~## TRI
此命令将输出具有从高程计算出的值的单波段栅格。 TRI代表地形坚固性指数,其定义为中心像素及其周围细胞之间的平均差异(参见Wilson等人2007,Marine Geodesy 30:3-35)。
值-9999用作输出节点值。
没有具体的选择。
~#~## TPI
此命令将输出具有从高程计算出的值的单波段栅格。 TPI代表地形位置指数,其定义为中心像素与其周围细胞平均值之间的差异(参见Wilson等人2007,Marine Geodesy 30:3-35)。
值-9999用作输出节点值。
没有具体的选择。
~#~## roughness
此命令将输出具有从高程计算出的值的单波段栅格。 粗糙度是Wilson等人(2007年,海洋大地测量学30:3-35)中定义的中心像素及其周围细胞的最大的细胞间差异。
值-9999用作输出节点值。
没有具体的选择。