将xyz文件转换为geotiff

Question

我有一个包含xyz值的文件，我想将其转换为geotiff，以便使用geoserver在web上提供服务。xyz文件本质上是一个水深测量，点（xy）位于5x5网格上（但由于没有测量，因此缺少一些点）

有几点：

53827.50 437437.50 -12.71
53832.50 437437.50 -12.71
53837.50 437437.50 -12.71
53842.50 437437.50 -12.71
53847.50 437437.50 -12.71
53852.50 437437.50 -12.71
53857.50 437437.50 -12.71
53862.50 437437.50 -12.71
53867.50 437437.50 -12.71
53872.50 437437.50 -12.71
53877.50 437437.50 -12.71
53882.50 437437.50 -12.71
53887.50 437437.50 -12.71
53892.50 437437.50 -12.71
53897.50 437437.50 -12.71
53902.50 437437.50 -12.71
53907.50 437437.50 -12.71
53912.50 437437.50 -12.71
53917.50 437437.50 -12.71
53922.50 437437.50 -12.71

我尝试用gdal转换它：

gdal_translate -a_srs EPSG:28992 input.txt output.tif

这确实给了我一个文件，但结果是一个只有几个黑点的白色图像。我要找的是一个图像，它显示的颜色越深（最好是红色）的值越高，而颜色越浅（例如从浅红色到白色）的值越低。我还希望5x5的缺失单元格是透明的

从命令行使用python或GDAL是否可以实现这一点

Answer 1

这在GDAL中是可能的。 使用以下命令：

  gdaldem color-relief output.tif colormap.txt output_colored.tif

颜色映射（colormap.txt）行具有以下格式：高度R G B。 您的水深图颜色图可能如下所示：

0 254 236 237
-100 198 54 38

当然，您必须将-100调整为最小值，必要时调整RGB值。如果需要，还可以添加其他值

有趣的例子：https://blog.mastermaps.com/2012/06/creating-color-relief-and-slope-shading.html