Описание тега kriging

Вопросы с тегом
Метод статистической интерполяции, также известный как регрессия гауссовского процесса, наиболее часто используемый в геостатистике. Цель состоит в том, чтобы нанести на карту поверхность с учетом ограниченных выборочных данных. Процесс оценивает изменчивость предоставленных данных, затем использует средневзвешенное значение соседних точек с учетом расстояния и направления для интерполяции желаемых точек карты.
1 ответ

Условное моделирование процесса Кригинга / Гаусса в Matlab

Я хотел бы выполнить условное моделирование для моделей гауссовского процесса (GP) в Matlab. Я нашел учебное пособие Мартина Коларжа ( http://mrmartin.net/?p=223). sigma_f = 1.1251; %parameter of the squared exponential kernel l = 0.90441; %paramete…
02 май '15 в 09:15
1 ответ

Зацикливание функции IDW в R

Я новичок в R и хотел бы протестировать различные параметры для выполнения интерполяции с обратным взвешиванием расстояния (IDW). data.idw.n <- idw(variable~1, data, data.grid, nmax=n) Я хотел бы повторить следующую функцию несколько раз, просто …
27 фев '15 в 14:55
0 ответов

Как провести 1-D кригинг с pykrige?

Я пытаюсь выполнить одномерный кригинг с помощью модуля python 'pykrige', но получаю несколько ошибок. Я попытался передать нулевой вектор в качестве аргумента для второй переменной. Мой код, как показано ниже hs = np.array(np.arange(0,1,bw)) hs0 = …
19 фев '18 в 09:01
1 ответ

Условное моделирование (с Кригингом) в R с распараллеливанием?

Я использую пакет gstat в R для генерации последовательных гауссовых симуляций. Мой компьютер имеет 4 ядра, и я попытался распараллелить функцию krige(), используя параллельный пакет, следуя сценарию, предоставленному Chris McKee, чтобы ответить на …
25 май '18 в 11:53
0 ответов

Интерполяция категориальных переменных в r

У меня есть набор данных, который состоит из шейп-файлов округов страны, и каждый округ имеет категориальное значение: ему назначено 1,2,3. Это приводит к карте, показывающей районы с одним из 3 цветов. Тем не менее, я считаю, что границы цветов сли…
29 ноя '18 в 00:38
1 ответ

R: плохая подгонка вариограммы, плохие результаты кригинга

Я пытаюсь сделать кригинг в Джакартском заливе. У меня есть набор точек измерения с соответствующими координатами и атрибутами (pH, соленость,...) Чтобы сделать кригинг, мне сначала нужно найти модель для моей вариограммы. Когда я использую функцию …
03 сен '15 в 04:19
1 ответ

R универсальный кригинг с автоКригом ()

Я пытаюсь использовать функцию autoKrige() в пакете automap для простого применения универсального кригинга. У меня нерегулярно расположенная сетка измерений, и я хочу интерполировать их в точном пространственном масштабе. Пример кода: library('auto…
07 сен '15 в 13:00
1 ответ

R gstat пространственно-временная вариограмма кригинга

Я пытаюсь использовать функцию variogramST из пакета R gstat для вычисления пространственно-временной вариограммы. Существует 12 лет данных с 20000 точек данных в нерегулярных точках в пространстве и времени (без полной или частичной сетки). Я долже…
02 авг '16 в 21:05
1 ответ

Как написать формулу в цикле для autoKrige {automap} в R

У меня большая матрица и имена столбцов следующие: colid=vector(length = 60) for(i in 1"60) { colid[i]=paste0("V",i) } Когда я использую autoKrige функция в automapформула должна быть записана в таком формате: V1~1. Когда я пытаюсь сделать это с пом…
10 апр '17 в 22:34
0 ответов

Обычный Кригинг в 2D и 3D: проблемы с дублированием?

Я хочу криге в 2-х мерной и 3-х мерной области (отдельно). У меня есть много моментов, которые были воспроизведены ниже. НОТА: Атрибут, который я хотел бы обозначить, - это Y. Они намеренно отрицательны, поскольку являются значениями экспонент.Домен…
01 май '15 в 16:36
1 ответ

Как перекрыть карту пространственного прогнозирования кригинга на определенной области карты страны в R?

У меня есть почасовой набор данных PM10 для 81 наблюдения под названием "seoul032823". Вы можете скачать здесь. Я выполнил обычный кригинг для этого набора данных, а также получил пространственную карту для прогнозирования кригинга. Я также могу пок…
29 мар '16 в 08:01
1 ответ

Создать вариограмму в пакете R's gstat

Предположим, у меня есть данные об осадках, полученные на четырех метеостанциях за период 2004-2016 гг. Я загрузил данные в базу данных для поиска в R. Моя цель состоит в том, чтобы брать данные за каждый день этого периода и повторно использовать э…
16 апр '17 в 05:31
2 ответа

Может ли кто-нибудь описать метод двумерной интерполяции, который лучше, чем билинейная интерполяция?

У меня есть сетка точек данных, на которой я в настоящее время использую билинейную интерполяцию, чтобы найти недостающие точки в сетке. Мне указали в направлениях Кригинг, или, как ты, лучшая линейная объективная оценка, но я не смог найти хороший …
21 ноя '14 в 14:58
1 ответ

Как извлечь конкретные значения с точечными координатами из интерполяций Кригинга, выполненных в R?

Используя R версии 3.4.2 и библиотеку "geoR", я сделал кригол интерполяции для разных переменных (ниже я привожу пример моего процесса). Я также сделал матрицу с координатами для 305 деревьев с различными отметками (виды, DBH, высота), которые наход…
29 май '18 в 23:53
1 ответ

ГАМ с "gp" более плавным: прогноз на новых локациях

Я использую следующую геоаддитивную модель library(gamair) library(mgcv) data(mack) mack$log.net.area <- log(mack$net.area) gm2 <- gam(egg.count ~ s(lon,lat,bs="gp",k=100,m=c(2,10,1)) + s(I(b.depth^.5)) + s(c.dist) + s(temp.20m) + offset(log.n…
01 авг '18 в 13:46
1 ответ

Получение теоретической вариограммы из параметризованной ковариационной функции (СТК)

Я использую набор инструментов STK в течение нескольких дней для кригинга полей параметров окружающей среды, то есть в геостатистическом контексте. Я нахожу набор инструментов очень хорошо реализованным и полезным (большое спасибо авторам!), И предс…
02 мар '18 в 20:37
0 ответов

Выполнить Kriging интерполяцию с использованием Arcpy

У меня есть список классов точечных объектов. Я пытаюсь написать скрипт Python для выполнения интерполяции Кригинга. Я получаю сообщение об ошибке массажа в этом коде "Point_Num" не определено, Ниже сценария я работаю import arcpy from arcpy import …
22 ноя '16 в 11:58
0 ответов

Пространственно-временное прогнозирование с помощью KrigeST в R

Я новичок в R и ST kriging, я создал объект STFDF с данными за первые 4 часа и хочу прогнозировать данные за 5 часов. Используя krigeST в R, я создал модель vgm. Но после запуска: pred = krigeST(formula=values~1, data=my_stfdf, newdata = my_pred, mo…
06 ноя '13 в 08:50
0 ответов

openmdao Kriging суррогатные претензии плохие данные, но это не так

Я пытаюсь обучить суррогат в openmdao 1.7.1 с помощью FloatKrigingSurrogate(). Я получаю сообщение об ошибке: /applications/Frameworks/Python_Frameworks/OpenMDAO/openmdao-1.7.1/lib/python2.7/site-packages/openmdao/surrogate_models/kriging.py:145: Ru…
30 авг '16 в 21:11
1 ответ

Существуют ли программы, в которых реализован процесс Гаусса с несколькими выходами?

Я пытаюсь реализовать байесовскую оптимизацию, используя регрессию Гаусса, и сначала хочу попробовать GP с несколькими выходами. Есть много программ, которые реализуют GP, например, fitrgp функция в MATLAB и панели инструментов ooDACE. Но я не нашел…
07 окт '16 в 07:05