Область неправильного многоугольника как график на spatstat

Я впервые использую пакет spatstat, поэтому мне нужен совет. Я пытаюсь отобразить данные координат в области неправильного многоугольника (формат .shp), чтобы вычислить пространственный анализ, как К. Рипли. Как я могу добавить область неправильного многоугольника в качестве графика? Как я могу объединить данные .ppp из координат в область многоугольника? Я использовал следующие коды:

Преобразование координатных данных в формат .ppp

library(spatstat)
library(sp)
library(maptools)

tree.simu <- read.table("simulation.txt", h=T)
tree.simu.ppp <-ppp(x=tree.simu$X,y=tree.simu$Y,window=owin(c(min(tree.simu$X),max(tree.simu$X)),c(min(tree.simu$Y),max(tree.simu$Y))))
plot(tree.simu.ppp)

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

Построение площади неправильного многоугольника

area <- readShapePoly("Area/Fragment.shp")
plot(area)
plot(tree.simu.ppp, add=T)

or

points(tree.simu.ppp)

Пакет принимает последнюю функцию, но, когда я пытаюсь построить оба файла вместе, кажется, что файл .shp заполняет всю область. Я не могу визуализировать данные координат.

Спасибо, очень ценю вашу помощь!

ps: Если вы знаете какой-либо материал по этому вопросу, я буду рад взглянуть


r plot polygon spatial spatstat
person Renzo Ferreira    schedule 26.06.2014    source источник
comment
Должна быть проблема с вашим набором данных, так как это должно работать. Можете ли вы загрузить где-нибудь simulation.txt и набор файлов, составляющих Fragment.shp, и опубликовать ссылку в своем вопросе? В противном случае это просто догадки. Например, вы уверены, что ограничивающая рамка для tree.simu перекрывается с ограничительной рамкой для area ??   -  person jlhoward    schedule 26.06.2014
comment
Да, опубликуйте данные, чтобы мы могли увидеть, что пошло не так. Возможно, вы найдете вдохновение в виньетке spatstat с шейп-файлами: cran.r-project.org/web/packages/spatstat/vignettes/   -  person Ege Rubak    schedule 27.06.2014
comment
Привет, Jlhoward, я надеюсь, что это тренировка. Область фрагмента и данные моделирования координат находятся по следующей ссылке: drive.google.com/, если вам нужно что-то еще, дайте мне знать.   -  person Renzo Ferreira    schedule 27.06.2014

Ответы (1)


arrow_upward
0
arrow_downward

Это действительно связано с непоследовательными ограничивающими рамками, как предполагалось в комментарии @jlhoward. Ваши точки находятся в [273663.9, 275091.45] x [7718635, 7719267], а многоугольник содержится в [-41.17483, -41.15588] x [-20.619647, -20.610134].

Если предположить, что координаты действительно совпадают с окном, то правильный способ поместить их в объект ppp:

library(spatstat)
library(sp)
library(maptools)

area <- readShapePoly("Area/Fragment.shp")
area <- as(area, "owin")
tree.simu <- read.table("simulation.txt", h=T)
tree.simu.ppp <-ppp(x=tree.simu$X,y=tree.simu$Y,window=area)

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

person Ege Rubak    schedule 27.06.2014