как установить диапазоны координат X и Y для геометрии окна наблюдения в пакете spatstat в R

Привет? У меня есть данные о распространении саженцев, которые содержат типы видов, координаты X и Y в UTM. Я хочу создать точечный узор по их координатам X и Y с помощью функции ppp () в пакете spatstat. Я попробовал это двумя способами:

 p.patt <- ppp(mydata$X, mydata$Y)
 p.patt <- ppp(mydata$X, mydata$Y, owin(c(100,131), c(100,130)))

Но есть «Предупреждение: 435 точек были отклонены как лежащие вне указанного окна» для них обоих.

Я предполагаю, что это связано с диапазонами координат X и Y, которые должны быть указаны в этом коде в c (…), c (…). Я проверил диапазон X и Y, и R дал мне следующие диапазоны:

for X: 368615 and 368746, 
for Y: 4587355 and 4587485

Когда я рисую данные, форма графика выглядит как «наклонный ромбо». Не знаю, поможет ли это.

Здесь я случайно выбрал несколько чисел: 100, 131 и 130. Я не мог найти никакой информации, как их установить в Интернете.

Итак, мой вопрос в том, как я могу использовать эти диапазоны координат для установки геометрии окна наблюдения для точки в пакете spatstat в R ?.

Заранее большое спасибо!


r coordinates plot spatstat ppp
person R starter    schedule 26.02.2019    source источник

Ответы (2)


arrow_upward
0
arrow_downward

Числа в вызове owin не являются шириной и высотой окна; это X и Y координаты углов окна.

Поскольку диапазон X значений координат точек данных составляет от 368615 до 368746, окно должно содержать как минимум этот диапазон. Точно так же в окне должен содержаться диапазон значений Y. Минимальное окно, в котором не будет предупреждений, -

p.patt <- ppp(mydata$X, mydata$Y, owin(c(368615,368746), c(4587355,4587485)))

или эквивалентно

p.patt <- ppp(mydata$X, mydata$Y, c(368615,368746), c(4587355,4587485))

Но это только минимальное допустимое окно; для правильного анализа необходима информация о регионе опроса. Если это не прямоугольник, то, как говорит Эге, вам нужно указать owin(poly=...), используя координаты вершин многоугольника.

person Adrian Baddeley    schedule 27.02.2019
comment
это работает хорошо. Теперь я понимаю. Рассмотрю Edge и ваше предложение по углам области съемки. Большое спасибо вам обоим! - person R starter; 27.02.2019

arrow_upward
0
arrow_downward

У вас нет информации о сюжете? Например. координаты углов многоугольной области, ограничивающей участок? Если у вас есть эти координаты, используйте их в качестве входных данных в аргументе poly из owin. Подробности см. В файле справки для owin. При отсутствии информации вы можете попробовать ripras оценить границы участка.

Что вы делаете прямо сейчас, так это то, что вы определяете точечный образец в прямоугольнике [0,131] × [0,130], а затем вы предоставляете набор точек с координатами за пределами этой области (гораздо большие значения координат), и все они отбрасываются.

person Ege Rubak    schedule 26.02.2019
comment
@ Эге Рубак. К сожалению, у меня нет информации о них. Можете ли вы сказать мне, что самое важное, что нужно учитывать при его настройке? - person R starter; 26.02.2019