У меня есть файл со значениями 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 из командной строки?