Форумы GIS-Lab.info

В 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) - выдаёт адекватную картинку.

Отвечу сам себе:

Код: Выделить всё

> S <- readShapePoints("/.../myPoints.shp")
> SP <- as(S, "SpatialPoints")
> P <- as(SP, "ppp")
> Z <- readShapeSpatial("/.../myPolyBoundary.shp")
> Z1 <- as(Z, "SpatialPolygons")
> W <- as.owin(Z1)
> P <- P[W]
> plot(P)

Форумы GIS-Lab.info

Как загрузить мультиполигон в качестве окна?

Как загрузить мультиполигон в качестве окна?

Re: Как загрузить мультиполигон в качестве окна?