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

Question

Привет?У меня есть данные о распределении сеянцев, которые содержат виды видов, координаты 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 ?.

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

Answer 1

Числа в вызове 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=...), используя координаты вершин многоугольника.

Answer 2

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

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