Растр в Каппу

Question

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

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

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

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

Размер окна составляет 100x100 м

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).Это означает, что у вас есть родители повсюду, но они производят потомство в зависимости от местоположения, поэтому в областях со многими взрослыми / средними деревьями генерируется немного потомков, тогда как многие рождаются в пустых областях.Подробности в книге.Удачи!