Конвертировать xyz файл в geotiff - PullRequest
0 голосов
/ 15 апреля 2020

У меня есть файл со значениями xyz, который я хочу преобразовать в геотиф, чтобы я мог использовать его в Интернете с помощью геосервера. Файл 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

I Я попытался преобразовать его с помощью gdal:

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

Это дает мне файл, но в результате получается белое изображение с несколькими черными точками. То, что я ищу, - это изображение, которое показывает более темный цвет (предпочтительно красный) для более высоких значений и более светлый цвет (например, от светло-красного до белого) для более низких значений. Я также хотел бы, чтобы отсутствующие ячейки 5x5 были прозрачными.

Возможно ли это с помощью python или GDAL из командной строки?

1 Ответ

0 голосов
/ 18 апреля 2020

Это возможно с GDAL. Используйте следующую команду:

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

Строка цветовой карты (colormap.txt) имеет следующий формат: высота RG B. Ваша цветовая карта c карты может выглядеть следующим образом:

0 254 236 237
-100 198 54 38

Конечно, вы должны настроить минимальное значение -100, при необходимости отрегулируйте значения RGB. Также вы можете добавить другие значения при необходимости.

Интересный пример: https://blog.mastermaps.com/2012/06/creating-color-relief-and-slope-shading.html

...