В R разбираюсь пока плохо; пытаюсь сделать так:
Код: Выделить всё
S <- readShapePoints("/../myPoints.shp")
SP <- as(S, "SpatialPoints")
P <- as(SP, "ppp")
S1 <- readShapePoly("/../myPolyBoundary.shp")
owin(poly = S1)
вылезают ошибки: "Отсутствует метод для конверсии данного класса S4 в вектор" либо (если owin(poly = S1) выполнять отдельно) "'x' and 'y' lengths differ"
некоторые точки выходят за границы мультиполигона, однако, если использовать мультиполигон, покрывающий все точки, ситуация не изменяется.
при этом > plot(S1) - выдаёт адекватную картинку.
Последний раз редактировалось
SS_Rebelious 27 июн 2011, 22:25, всего редактировалось 1 раз.
Look for something long enough, and you will find it. Look for something without understanding, and it will find you...
"All paid jobs absorb and degrade the mind." Aristotle
If you take 1 step towards freedom it'll take 2 steps towards you!