ArcInfo ASC файл в трехколоночный XYZ ASCII-файл

[nickleb]

Какими средcтвами (про ArcGIS не говорим) конвертировать ArcInfo ASC grid-файл в трехколоночный (XYZ) ASCII-файл?

[lam]

Global Mapper

[Игорь Белов]

Открытый пакет GMT, команда grd2xyz -E

Проприетарный Golden Software Surfer ([Grid] → [Convert]).

[nickleb]

Открытый пакет GMT, команда grd2xyz -E

Проприетарный Golden Software Surfer ([Grid] → [Convert]).

Спасибо. Но в старенький Surfer - не загружается большой ASC.

[Игорь Белов]

в старенький Surfer - не загружается большой ASC.

Это смешно, честное слово. Чтобы последовательно читать значения Z и выводить X, Y, Z, не нужно загружать растр в память. Собственную программу можно написать за несколько минут, нужно только знать порядок записи узлов, чтобы выводить правильные X, Y. А этот порядок известен.

В конце-концов, поставьте GMT. Для целочисленных Z запускайте

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

grd2xyz myraster.asc -E > myraster.xyz

Для действительных чисел используйте

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

grd2xyz myraster.asc -Ef > myraster.xyz

[Донецков]

Что значит старенький? Помнится еще в 6 версии (это где-то 1996 - 1997 годы) в таблицу можно было загрузить до 2 млн. записей...

[nickleb]

записей более 134 млн. - 30''-ндный массив глубин Сев.Лед.Океана... В новом 12-ом Surfer'е можно, но у меня нет... спасибо всем... заPythоn'ю, наверное эту задачку... просто с ходу ни в Q, ни в SAGA, ни в DIVA этого не нашел... спасибо за отклики... извините за отнятое время... с праздником и победой на Олимпиаде!

[nickleb]

ErnieBoyd, lam, Донецков - спасибо Вам!

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

Как вариант:

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

gdal_translate -of "XYZ" infile.asc infile.xyz

[nickleb]

Как вариант:

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

gdal_translate -of "XYZ" infile.asc infile.xyz

Александр, спасибо. Запустил предложенный Вами вариант. Сработал. Из 1.2 Гб получил 7.5 Гб.

[gamm]

Из 1.2 Гб получил 7.5 Гб.

что-то мне подсказывает, что на этом создание трудностей и борьба с ними не закончится. Теперь обратно интерполировать будем, предварительно пересчитав проекцию ?

[nickleb]

да, gamm... здесь и прозорливости особой не надо, чтобы понять мыкания чайника к Surfer'у привыкшего, где имел дело с grid'ами небольших размеров и slice'овал оные... вот для начала в spatialite, что ли это положить или в PostGIS?... или, вообще, никуда не класть, а gdalwarp'ить в точки с lat/lon по мере надобности?

[nickleb]

...кто-нибудь подскажет ресурс, книгу с "множественными" примерами использования gdal/ogr?..

[gamm]

...кто-нибудь подскажет ресурс, книгу с "множественными" примерами использования gdal/ogr?..

у них же на сайте полный комплект документации + исходники, наше все. А чего делать-то собрались, можете описать?

[nickleb]

в старенький Surfer - не загружается большой ASC.

Это смешно, честное слово. Чтобы последовательно читать значения Z и выводить X, Y, Z, не нужно загружать растр в память. Собственную программу можно написать за несколько минут, нужно только знать порядок записи узлов, чтобы выводить правильные X, Y. А этот порядок известен.

В конце-концов, поставьте GMT. Для целочисленных Z запускайте

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

grd2xyz myraster.asc -E > myraster.xyz

Для действительных чисел используйте

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

grd2xyz myraster.asc -Ef > myraster.xyz

... в GMT 4.5 формат ASC ArcINFO-файла эти утилиты не распознают, а 5-рку не стал ставить... gdalytanslat'ом по совету Ал.-ра Мурого сделал... но у меня к Вам уже др. вопрос: а есть ли утилита в GMT по растру получить зн.-ия Z в наборе точек c XY нерегулярной сетки?... нечто аналогичное gdalwarp'у?.. или slicе'у Surfer'а?.. спасибо!