Как построить L Рипли, используя функцию Lest () из spatstat - PullRequest
Новая
       21

Как построить L Рипли, используя функцию Lest () из spatstat

0 голосов
/ 12 июня 2019

У меня есть файл csv с количеством точек.Данные включают в себя долготу, широту и дату.Я классифицировал свои данные в зависимости от сезона (весна, лето, осень и зима).Мне нужно рассчитать L Рипли для каждого сезона.Я использовал Lest() из пакета spatstat в R. Я мог построить L Рипли, используя эту функцию.Но диапазон оси x (расстояние) для каждого участка (сезона) различен.Как я могу сделать все четыре участка (сезона) в одном диапазоне (расстоянии)?

sa  <- readOGR(".", "boundary")  # Don't add the .shp extension

W   <- as(sa, "owin")

Summer.ppp <- as(Sum_UTM, "ppp")

Autumn.ppp <- as(Autumn_UTM, "ppp")

Winter.ppp <- as(Winter_UTM, "ppp")

Spring.ppp <- as(Spring_UTM, "ppp")


Summer <- Lest(Summer.ppp, correction="Ripley")
plot(Summer , xlab="distance (m)", ylab="K(d)", cex.lab= 2, cex.axis=2)
plot(Summer, . -r ~ r, ylab=expression(hat("L")), xlab="distance (m)",cex.lab= 2, cex.axis=2)


Autumn <- Lest(Autumn.ppp, correction="Ripley")
plot(Autumn, xlab="distance (m)", ylab="K(d)", cex.lab= 2, cex.axis=2)
plot(Autumn, . -r ~ r, ylab=expression(hat("L")), xlab="distance (m)",cex.lab= 2, cex.axis=2)

Winter <- Lest(Winter.ppp, correction="Ripley")
plot(Winter, xlab="distance (m)", ylab="K(d)", cex.lab= 2, cex.axis=2)
plot(Winter, . -r ~ r, ylab=expression(hat("L")), xlab="distance (m)",cex.lab= 2, cex.axis=2)

Spring <- Lest(Spring.ppp, correction="Ripley")
plot(Spring , xlab="distance (m)", ylab="K(d)", cex.lab= 2, cex.axis=2)
plot(Spring, . -r ~ r, ylab=expression(hat("L")), xlab="distance (m)",cex.lab= 2, cex.axis=2)

1 Ответ

0 голосов
/ 12 июня 2019

Вы пишете, что у вас есть широта и долгота, но в то же время ваши наборы данных называются как-то с "UTM", что будет означать, что у вас есть плоские координаты x, y. Что он? В spatstat вы можете работать только с плоскими координатами. Вы не объявляете, какие пакеты вы используете. Приведенная выше функция as для преобразования в ppp, вероятно, из maptools? В этом случае я думаю, что окно ppp просто установлено в ограничивающий прямоугольник координат. У вас есть объект с именем W, который вы, вероятно, захотите использовать в качестве окна ppp. Вы можете сделать что-то вроде: 

Window(Spring.ppp) <- W

Как уже упоминалось в комментариях, вы можете просто указать xlim = c(lower,upper) в вызове plot, где lower и upper - числовые значения желаемых пределов по оси x.

...