ArcInfo ASC файл в трехколоночный XYZ ASCII-файл
Добавлено: 23 фев 2014, 12:45
Какими средcтвами (про ArcGIS не говорим) конвертировать ArcInfo ASC grid-файл в трехколоночный (XYZ) ASCII-файл?
Геоинформационные системы (ГИС) и Дистанционное зондирование Земли
https://gis-lab.info/forum/
Спасибо. Но в старенький Surfer - не загружается большой ASC.ErnieBoyd писал(а):Открытый пакет GMT, команда grd2xyz -E
Проприетарный Golden Software Surfer ([Grid] → [Convert]).
Это смешно, честное слово. Чтобы последовательно читать значения Z и выводить X, Y, Z, не нужно загружать растр в память. Собственную программу можно написать за несколько минут, нужно только знать порядок записи узлов, чтобы выводить правильные X, Y. А этот порядок известен.nickleb писал(а):в старенький Surfer - не загружается большой ASC.
Код: Выделить всё
grd2xyz myraster.asc -E > myraster.xyzКод: Выделить всё
grd2xyz myraster.asc -Ef > myraster.xyzКод: Выделить всё
gdal_translate -of "XYZ" infile.asc infile.xyzАлександр, спасибо. Запустил предложенный Вами вариант. Сработал. Из 1.2 Гб получил 7.5 Гб.Александр Мурый писал(а):Как вариант:Код: Выделить всё
gdal_translate -of "XYZ" infile.asc infile.xyz
что-то мне подсказывает, что на этом создание трудностей и борьба с ними не закончится. Теперь обратно интерполировать будем, предварительно пересчитав проекцию ?nickleb писал(а):Из 1.2 Гб получил 7.5 Гб.
у них же на сайте полный комплект документации + исходники, наше все. А чего делать-то собрались, можете описать?nickleb писал(а):...кто-нибудь подскажет ресурс, книгу с "множественными" примерами использования gdal/ogr?..
... в GMT 4.5 формат ASC ArcINFO-файла эти утилиты не распознают, а 5-рку не стал ставить... gdalytanslat'ом по совету Ал.-ра Мурого сделал... но у меня к Вам уже др. вопрос: а есть ли утилита в GMT по растру получить зн.-ия Z в наборе точек c XY нерегулярной сетки?... нечто аналогичное gdalwarp'у?.. или slicе'у Surfer'а?.. спасибо!ErnieBoyd писал(а):Это смешно, честное слово. Чтобы последовательно читать значения Z и выводить X, Y, Z, не нужно загружать растр в память. Собственную программу можно написать за несколько минут, нужно только знать порядок записи узлов, чтобы выводить правильные X, Y. А этот порядок известен.nickleb писал(а):в старенький Surfer - не загружается большой ASC.
В конце-концов, поставьте GMT. Для целочисленных Z запускайтеДля действительных чисел используйтеКод: Выделить всё
grd2xyz myraster.asc -E > myraster.xyzКод: Выделить всё
grd2xyz myraster.asc -Ef > myraster.xyz