3 votos

¿Añadir una capa de shapefile de puntos a una ventana de shapefile de polígonos en R?

Soy nuevo en Spatstat para R y estoy intentando añadir una capa de puntos a una ventana que he creado a partir de un polígono. Veo que ya se han hecho preguntas como esta pero sigo sin conseguir lo que necesito con las respuestas dadas. Hasta ahora tengo:

##Run libraries: 
library(“spatstat”)
library(“maptools”)
library(“sp”)
library(“rgdal”)

##set the working directory
setwd(“Z:/CP/spatstat/data/”)

##read in the shapefiles
Points <- readOGR(getwd(), "MyPointShapefile")
StudyArea <- readOGR(getwd(), "MyPolygonShapefile")

##create the window from polygon shp
w <- as.owin(StudyArea)

¿Cómo se suman los dos para que los puntos se proyecten en la ventana?

3voto

Dan Puntos 16

Sólo hay que coaccionar "Puntos" a un objeto ppp usando "w" como ventana.

w <- as.owin(StudyArea)
  pts.ppp <- as.ppp(coordinates(Points), w)
    plot(pts.ppp)

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X