Automap - настройка выходного разрешения и использование нескольких ядер?

Question

У меня есть два вопроса:

Я использую функцию autoKrige в automap для генерации некоторых растровых данных и пытаюсь выяснить, как увеличить разрешение выходных данных. Вывод слишком гранулированный, когда я открываю его в QGIS.

Пока что все, что я понял, это то, что это, вероятно, связано с аргументом "newdata". Я создал сетку размером 3000x3000 метров и использовал ее как "новые данные", но этот процесс занимает вечность, и я не уверен, что он закончится до сбоя.

Q1) Как вы устанавливаете разрешение вывода autoKrige?

Q2) Многоядерный пакет R будет работать с Autop? Я понял, что он использует только одно из двух моих ядер.

Answer 1

Вы уже используете правильные способы установки выходного разрешения: увеличение разрешения сетки, которую вы передаете newdata, Что касается производительности, узким местом часто является не выходное разрешение, а размер входного набора данных. Посмотрите на maxdist а также nmax аргументы krige, который может быть передан autoKrige и будет передано. Они ограничивают количество точек, используемых при прогнозировании, то есть ограничивают размер ковариационной матрицы, которую необходимо инвертировать. настройка nmax скажем, 100 должны помочь с производительностью.

Что касается вашего второго вопроса. Вы можете использовать несколько ядер только в том смысле, чтобы иметь возможность запускать несколько autoKrige анализирует рядом друг с другом. automap не поддерживает параллельную обработку в течение одного вызова autoKrigeи это было бы совсем немного работы, чтобы заставить это работать. Посмотрите на parallel пакет, который входит в стандартную комплектацию R, и mclapply функция в этом пакете.