Как получить область для ограничивающей рамки для каждой точки в данных плоского точечного массива (ppp) с помощью Spatstat в R?

У меня есть объект точечного шаблона (ppp), то есть набор данных, содержащий географические координаты определенных позиций. Например:

> data(bei)
> df1<-as.data.frame(bei)
> View(df1)
> head(df1)
      x     y
1  11.7 151.1
2 998.9 430.5
3 980.1 433.5
4 986.5 425.8
5 944.1 415.1
6 940.5 410.4

Теперь я хотел бы установить небольшую ограничивающую рамку некоторого размера над каждой точкой объекта ppp и вычислить ее площадь. Как это сделать с помощью spatstat?

Или я должен пикселизировать объект ppp? В таком случае, как получить площадь пикселя?


r geospatial spatial spatstat
person Lesnar    schedule 16.12.2018    source источник
comment
Я не знаком с spatstat, но, возможно, все, что вам нужно, это найти наименьший прямоугольник, содержащий все точки, заданные, например, head(df1) и его площадью? Я спрашиваю, потому что кажется, что на ваш вопрос можно ответить, не зная о spatstat, и в этом случае у вас гораздо больше шансов получить быстрый ответ.   -  person Julius Vainora    schedule 16.12.2018
comment
Ну, я не хочу, чтобы прямоугольник содержал все точки, это можно легко сделать, создав окно. Я хочу создать прямоугольник или квадраты для каждой точки. Как это сделать в R (даже без использования spatstat)?   -  person Lesnar    schedule 16.12.2018
comment
neonscience.org/field-data-polygons-centroids   -  person G. Cocca    schedule 16.12.2018

Ответы (1)


arrow_upward
1
arrow_downward

В пакете spatstat функция boundingbox вычисляет наименьший прямоугольник, содержащий указанные пространственные объекты. В вашем примере вы можете ввести boundingbox(bei), чтобы получить наименьший прямоугольник, содержащий все точки данных в наборе данных точечного массива bei. Площадь рассчитывается с использованием area(boundingbox(bei)). Если вместо набора данных точечного шаблона у вас был только фрейм данных с координатами, например df <- as.data.frame(bei), затем аналогично boundingbox(df) и area(boundingbox(df)) предоставят вам ту же информацию.

Для поиска информации о spatstat вы можете ввести help(spatstat) для обзора доступных команд или использовать help.search для поиска по ключевым словам. Или посетите веб-сайт spatstat.org.

person Adrian Baddeley    schedule 17.12.2018