Форумы GIS-Lab.info

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

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

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

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

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

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

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

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

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

Подытожим.
Чтобы преобразовать матрицу m*n в растр нужно:
1. Представить ее в виде массива где в качестве разделителя используется пробел.
2. Прописать заголовок файла в формате Arc/Info ASCII Grid, который понимает GDAL через драйвер AAIGrid. Для матрицы 360*180, (представляющий массив с разрешением 1 градус на весь мир) он будет выглядеть так: Где:
ncols -- количество колонок матрицы
nrows -- количество строк матрицы
xllcorner -- координаты западной границы участка
yllcorner -- координаты южной границы участка
cellsize -- разрешение растра (в градусах?)
nodata_value -- значение "нет данных".
Примечание: координаты могут указываться в десятичных градусах или в формате DD:MM:SS

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

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

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

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