gdaldem

创建日期:2024-06-21
更新日期:2025-01-19

分析和可视化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用作输出节点值。

没有具体的选择。