Añadir un punto en un spplot en R

Question

Tengo que añadir un punto sobre un existente spplot en R. Lo necesito sólo para un propósito visual. Tengo las coordenadas del punto, pero no sé cómo cambiar el código para añadir el punto.

Este es el código del spplot :

spplot(pioggiaTP[1])

Y este es el mapa que se obtiene al ejecutar el comando:

Las coordenadas del punto son X=720731,8704 e Y=4861300,795 en el sistema de referencia WGS 84 / UTM zona 32N EPSG:32632, el mismo sistema de la parcela.

¿Alguien puede ayudarme?

Answer 1

spplot devuelve un celosía parcela de la clase "trellis". Puede utilizar el layer de la función latticeExtra para actualizar este objeto.

Ejemplo:

library("sp")
library("latticeExtra")

data(meuse)
data(meuse.grid)

coordinates(meuse.grid) <- ~x+y
proj4string(meuse.grid) <- CRS("+init=epsg:28992")

p <- spplot(meuse.grid, "dist")
print(p)

p + layer(panel.points(x, y, col="green", pch=19), data=meuse)

Answer 2

Otra opción sería utilizar el argumento 'sp.layout' que viene con spplot . Tenga en cuenta que, a diferencia del enfoque sugerido por @rcs, es necesario suministrar un objeto "sp" (que puede tener ya disponible, ya que está trabajando con spplot ). Utilizando el conjunto de datos "meuse.grid", su código podría ser aproximadamente el siguiente.

## sample point
test <- data.frame(x = 180000, y = 331500)
coordinates(test) <- ~ x + y
proj4string(test) <- "+init=epsg:28992"

## display gridded data including sample point
spplot(meuse.grid, "dist", scales = list(draw = TRUE),
       sp.layout = list("sp.points", test, pch = 16, cex = 2, col = "black"))