Растр в Каппу

Question

Я пишу код для имитации лесных насаждений. Я смоделировал положение больших и средних деревьев с помощью Spatstats. Следующий шаг, который я хочу сделать, это моделировать меньшие деревья. Я пытался написать некоторый код, который позволяет меньшим деревьям располагаться дальше от больших и средних деревьев.

Я могу смоделировать правильное количество меньших деревьев, используя одно число каппа для модели сгущения Томаса в Spatstat, однако это дает сгущение близко к большим / средним деревьям. Я изо всех сил пытаюсь смоделировать правильное количество меньших деревьев, где родительские точки расположены в менее плотных частях древостоя.

Я попытался использовать функцию плотности в spatstat на больших и средних деревьях. Затем я попытался инвертировать растр, чтобы области с меньшим количеством деревьев имели более высокое значение пикселя, но я не могу найти способ изменить его на растр Каппа, который даст мне правильное количество меньших деревьев, которые будут сгенерированы.

Код ниже дает мне номер для правильного числа небольших деревьев, которые будут сгенерированы

Размер окна 100х100м

kappaS<-(Number_Small_Trees_to_be_Generated)/(mean(muS)*winarea)

Small.Trees1<-rThomas(kappa=kappaS,
                      scale=sigmaM,
                  mu=muS,
                  win=Window222)

Я надеюсь найти способ создания каппа-растра, который имитирует правильное количество точек, где родительские точки более удобны для создания на менее плотных участках стенда.

Большое спасибо за любую помощь!

Answer 1

Я рекомендую вам прочитать Раздел 12.3.4 "Модели неоднородных кластеров" в spatstat книга Пространственные точечные паттерны: методология и приложения с R. К сожалению, это не одна из бесплатных глав, поэтому вы должны надеяться, что в вашей местной библиотеке есть книга. (Отказ от ответственности: я соавтор.)

В этом разделе объясняется, как интенсивность предложенного вами процесса является сложным выражением (сверткой), и нелегко контролировать, сколько точек будет сгенерировано. С математической точки зрения гораздо легче справиться с неоднородным процессом потомства (т.е. mu быть пространственно изменяющимся). Это означает, что у вас есть родители повсюду, но они производят потомство в зависимости от местоположения, поэтому в областях со многими взрослыми / средними деревьями генерируется немного потомков, тогда как многие рождаются в пустых областях. Подробности в книге. Удачи!