Растр из текстовой матрицы

[bolotoved]

Можно ли в QGIS или напрямую в GDAL преобразовать матрицу с разделителями пробелами в tif? Может ее нужно перед импортом как-то специальным образом подготовить?

[Александр Мурый]

txt --> CSV --> VRT --> gdal_rasterize

[bolotoved]

Что-то застрял на стадии CSV --> VRT
В приложении собственно матрица, разделители -- запятые. 360*180.
Прописал команду:

Код: Выделить всё

gdal_rasterize -ts 360 180 matrix.csv matrix.vrt

Спрашивает про параметры -l, -a, как их правильно указать?

[Александр Мурый]

gdal_rasterize делает tif уже из готового vrt. Готовить VRT надо по этому рецепту.

Матрица довольно ужасная. Хотя нет, не матрица, а то, что это CSV.
Для перевода в VRT текст нужен в несколько полей, желательно XYZ, где Z -- значение, по к-й будем делать растр, а XY -- координаты точек.

[bolotoved]

Да уж
Полагаю, что должен быть какой-то иной способ преобразования матрицы m*n в растр, пусть и без привязки.

[Voltron]

Такой способ есть. Матрицу загоняем в массив NumPy, потом это безобразие сохраняется как канал растра. Все

[Максим Дубинин]

или превратить матрицу в ASCII растр убрав все переводы строки и добавив минимальный заголовок

[bolotoved]

или превратить матрицу в ASCII растр убрав все переводы строки и добавив минимальный заголовок

Подытожим.
Чтобы преобразовать матрицу m*n в растр нужно:
1. Представить ее в виде массива где в качестве разделителя используется пробел.
2. Прописать заголовок файла в формате Arc/Info ASCII Grid, который понимает GDAL через драйвер AAIGrid. Для матрицы 360*180, (представляющий массив с разрешением 1 градус на весь мир) он будет выглядеть так:

Код: Выделить всё

ncols        360
nrows        180
xllcorner    -180.000000000000
yllcorner    -90.000000000000
cellsize     1.000000000000
nodata_value 0

Где:
ncols -- количество колонок матрицы
nrows -- количество строк матрицы
xllcorner -- координаты западной границы участка
yllcorner -- координаты южной границы участка
cellsize -- разрешение растра (в градусах?)
nodata_value -- значение "нет данных".
Примечание: координаты могут указываться в десятичных градусах или в формате DD:MM:SS

3. Поменять расширение текстового файла на asc.
4. Полученный файл можно открыть в QGIS и экспортировать в любой другой растровый формат.

[ericsson]

То есть, говоря русским языком, добавить заголовок в формате Arc/Info ASCII Grid, который понимает GDAL через драйвер AAIGrid. Там еще всякие параметры имеются.

[bolotoved]

Спасибо, поправил.

[ericsson]

Есть "подводный камень".
Те параметры xllcorner, yllcorner, которые предложены, указывают на соответствующие края растра, если природа данных такова, что отсчеты находятся в центрах ячеек, и известны именно их координаты, то использовать надо параметры xllcenter, yllcenter