Kriging: Какой плагин QGIS или processing.runalg(....)

[nickleb]

Подскажите, пожалуйста, какой плагин QGIS или алгоритм из geoprocessing'а, доступного в инструментах QGIS 2.8, использовать для получения регулярного растра kriging'ом из ASCII-файла (Lon/Lot/Observation Value)? Надо в PyQGIS скрипт такую процедурку внедрить, т.к. исходных файлов много, а далее ещё и контуры надо строить и сохранять картинки...

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

Вы пробовали поискать?
Какой вид кригинга вам нужен, с какими параметрами? Или просто "в лоб", как обычно делают в Сёрфере?

[nickleb]

Вы пробовали поискать?
Какой вид кригинга вам нужен, с какими параметрами? Или просто "в лоб", как обычно делают в Сёрфере?

Александр, мне бы примерчик второй строчки фрагмента:

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

import processing
processing.runalg(....)

там бы по документации можно было бы и с возможными параметрами и опциями разобраться...
А, в Сёрфере, кстати, параметров гридирования для kriging'а хватает всяческих также... Сужу по старенькой 9-ой версии...
Спасибо, разберусь полегоньку...
В SAGA, по-моему, : и ординарный, и универсальный... с ними через QGIS и попробую, наверное...

[gamm]

для кригинга нужно сначала вариограмму вычислить. А если не вычислять, то лучше сплайн использовать, поскольку кригинг с линейной вариограммой без эффекта самородков, который стоит по умолчанию в SURFER - тоже сплайн

и делать все это лучше в R, там хоть вариограмму можно автоматом посчитать ...

P.S. да, и из опыта подшефных геологов - попытки угадать параметры ни к чему хорошему не приведут, там все нужно делать по науке.
P.P.S. с кригингом, если точек много, придется использовать локальное окно (его, кстати, тоже нужно по науке выбрать). А это приведет к негладкой модели, от чего подшефные геологи страшно гневались

[nickleb]

для кригинга нужно сначала вариограмму вычислить. А если не вычислять, то лучше сплайн использовать, поскольку кригинг с линейной вариограммой без эффекта самородков, который стоит по умолчанию в SURFER - тоже сплайн

и делать все это лучше в R, там хоть вариограмму можно автоматом посчитать ...

P.S. да, и из опыта подшефных геологов - попытки угадать параметры ни к чему хорошему не приведут, там все нужно делать по науке.

задался, м.б. дурацкой и блаженной, целью, но всё же...: в QGIS'е сделать одним py-скриптом: от гридирования (через построение вариограммы - почему бы и нет?) до готовой картинки... раз в нём можно через processing процедуры и из GRASS, и из SAGА, и GDAL/OGR... использовать... с R же QGIS тоже интегрируется...

[gamm]

я в R все делал, одним скриптом, включая построение картинок. Несколько суток считалось

но там хоть известно было, как модель примерно устроена, только параметры уточнялись. Так что вариограммный анализ лучше глазками контролировать и направлять (опять же личный опыт), тут дело такое ...

[gamm]

и для начала попробуйте МВА, который в SAGA есть (в R - тоже) - может, и кригинг не понадобится. Там хоть понятно, чем и как рулить ...

[nickleb]

но там хоть известно было, как модель примерно устроена, только параметры уточнялись. Так что вариограммный анализ лучше глазками контролировать и направлять (опять же личный опыт), тут дело такое ...

это да, конечно... это уже экспертный пилотаж...

[nickleb]

и для начала попробуйте МВА, который в SAGA есть (в R - тоже) - может, и кригинг не понадобится. Там хоть понятно, чем и как рулить ...

ежели по RBF Вы меня как-то просвещали, то по MBA пока нет... где можно почитать?

[gamm]

ежели по RBF Вы меня как-то просвещали, то по MBA пока нет... где можно почитать?

первоисточник
многа букв
тыц
тыц
тыц
тыц

[nickleb]

ежели по RBF Вы меня как-то просвещали, то по MBA пока нет... где можно почитать?

первоисточник
многа букв
тыц
тыц
тыц
тыц

спасибо за подборку!