grdhisteq
- 贡献者:
- 最近更新日期:
2025-10-17
- 官方文档:
- 简介:
对网格做直方图均衡
grdhisteq 允许用户找到能将给定网格文件划分为等面积区域的数据值。 其中一个常见用途是实现图像的直方图均衡化。 例如，假设用户有一个包含平坦地形并在中间有一座山的网格。 如果用 grdimage 或 grdview 按地形高度线性映射到灰度值， 则图像的大部分区域会非常暗，而山峰部分几乎全白。 使用 grdhisteq 可以将输入网格划分为 n_cells 个等面积的区间，并输出这些区间对应的数据值。 通过 awk 或 makecpt 可以用这些值生成一个 CPT 文件， 再用该 CPT 配合 grdimage 绘图，就能让图像中各个灰度级的分布更加均匀。 参照 grd2cpt。
第二个常见用途是根据某种累积分布函数生成统计意义上的等概率网格。 在这种情况下，输出网格在 (x, y) 位置上与输入网格相同，但其值被替换为在累积分布中的位置。 例如：
若要找出数据中最低的 10%，可以使用 grdhisteq 生成一个 n_cells = 10 的网格，
然后绘制结果中的 1 等值线，这条线将包围数据中最低 10% 的部分。
又如，为了对 grdgradient 的输出进行均衡，使其具有更平滑（如高斯型）的分布，
可以运行 grdhisteq 并使用 Gaussian 选项生成输出网格。 结果网格的值将服从均值为 0、方差为 1 的高斯分布，其空间位置与输入数据对应： 最小的输出值对应输入中最小的位置，依此类推。
语法
gmt grdhisteq ingrid [ -Goutgrid ] [ -Cn_cells ] [ -D[file] ] [ -N[norm] ] [ -Q ] [ -Rregion ] [ -V[level] ] [ -hheaders ] [ --PAR=value ]
注意：选项标志和相关参数之间不允许有任何空格。
必需参数
ingrid[=ID|?varname][+bband][+ddivisor][+ninvalid][+ooffset][+sscale]
输入网格名。通过追加 =ID 可指定 网格格式 [默认为 =nf]； 追加 ?varname 可指定 NetCDF 变量 [默认为 GMT 找到的第一个 2-D 网格]。 对网格追加子选项可进行额外设置：
+b 选取一个波段 band（仅用于图片）[默认为 0]
+d 将网格除以一个数 divisor [默认为 1]
+n 将网格中的 invalid 替换为 NaN
+o 将网格中的值进行一定的偏移，即加一个数 offset [默认为 0]
+s 将网格中的值缩放，即乘以 scale [默认为 1]
注 ：offset 操作位于 scale 之前。
可选选项
- -C n_cells
设置将数据范围划分为多少个区间（默认 16）。
- -D [file]
将层级信息输出到 file，若未指定文件则输出到标准输出。
-Goutgrid[=ID][+ddivisor][+ninvalid][+ooffset|a][+sscale|a][:driver[dataType][+coptions]]
输出网格文件名。通过追加 =ID 可指定 网格格式。 对网格追加子选项可进行额外设置：
+d 将网格除以一个数 divisor [默认为 1]
+n 将网格中的 invalid 替换为 NaN
+o 将网格中的值进行一定的偏移，即加一个数 offset，或使用 a 自动对值 进行调整以保证整数网格的精度 [默认为 0]
+s 将网格中的值缩放，即乘以 scale，或使用 a 自动对网格缩放以保证整数 网格的精度 [默认为 1]
注 ：offset 操作位于 scale 之前；+sa 将会同时设置 +oa。如果需 使用 GDAL 指定网格格式，则 ID 应设置为 gd，并指定 :driver 和可选的数据 类型 dataType，以及 +coptions 选项传递给 GDAL 的 -co 选项。
- -N [norm]
生成高斯输出。与 -G 一起使用可输出标准正态分布网格。 若附加 norm，则强制值落在 ±*norm* 范围内（默认为标准正态分布）。
- -Q
选择二次（非线性）直方图均衡化方式（默认是线性）。
- -Rxmin/xmax/ymin/ymax[+r][+uunit] (more …)
指定数据范围。
- -V[level] (more …)
设置 verbose 等级 [w]
- -h[i|o][n][+c][+d][+msegheader][+rremark][+ttitle] (more …)
跳过或生成指定数目的头段记录
- -^ 或 -
显示简短的帮助信息，包括模块简介和基本语法信息（Windows下只能使用 -）
- -+ 或 +
显示帮助信息，包括模块简介、基本语法以及模块特有选项的说明
- -? 或无参数
显示完整的帮助信息，包括模块简介、基本语法以及所有选项的说明
- --PAR=value
临时修改GMT参数的值，可重复多次使用。参数列表见 配置参数
示例
注: 下面是该模块的一些有效语法示例。使用远程文件（文件名以 @ 开头）的示例可以剪切并粘贴到您的终端中进行测试。其他需要输入文件的命令只是常见的使用类型的虚拟示例，但不能按编写的原样运行。
计算全球地形文件 earth_relief_05m 中能将选定区域划分为 16 个等面积区间的高度值:
gmt grdhisteq -C16 -D @earth_relief_05m -R0/20/0/20 > levels.txt
对文件 raw_intens.nc 中分布不均的亮度数据进行均衡化，使其适合 grdimage 或 grdview 使用:
gmt grdhisteq raw_intens.nc -Gsmooth_intens.nc -N -V
备注
对地理网格会执行加权的直方图均衡化，因为节点面积随纬度变化。
若使用 grdhisteq 生成高斯分布输出用于 grdimage 或 grdview 的光照着色，需注意： 输出值范围为 [-x, x]，其中 x 取决于输入网格点数 (nx × ny) 和高斯分布函数 F(x)。 令 N = nx × ny，则 x 会被调整，使 F(x) = (N - 1 + 0.5)/N。 标准正态分布中约 68% 的值位于 ±1 范围内，因此输出网格也将如此。 当 N 很大时，x 可能大于 4。 在 grdview 中，若光照梯度被裁剪到 [-1, 1]，则仅 68% 的数据会正确显示， 16% 会被裁剪为 -1，另 16% 为 +1。 若图像过亮或过暗，可使用 grdmath 乘以小于 1 的系数缩小值域，或使用 -N 指定归一化因子。