Hola,
Estoy intentando hacer un kriging con deriva externa con la funci�n "krige()",
pero me da un error:
Error in eval(expr, envir, enclos) : object 'radPcp' not found
El c�digo principal es:
radar.mx <- rasterToPoints(radar.cu)
colnames(radar.mx) <- c("x","y","radPcp")
radar.spdf <- as.data.frame(radar.mx)
lluvia.pluv <- as.data.frame(t(rbind(pluvio.h[1,log.NA], coord.aux[,log.NA])))
colnames(lluvia.pluv) <- c("pluvPcp", "x", "y")
lluvia.pluv.spdf <- lluvia.pluv
coordinates(lluvia.pluv.spdf) <- ~x+y
#radarpx <- as(radar.spdf, "SpatialPixelsDataFrame")
plantilla.sp <- as.data.frame(radar.mx)[,1:2]
coordinates(plantilla.sp) <- ~x+y
class(plantilla.sp)
for (v in 1:length(pluvio.h[1,log.NA])) {
lluvia.pluv.df <- data.frame(lluvia.pluv.spdf) [-v,]
lluvia.pluv.spdf <- lluvia.pluv.df
coordinates(lluvia.pluv.spdf) <- ~x+y
# Hacemos el KED. Ver funci�n "krige()":
KED.rad <- krige(
formula=pluvPcp~radPcp, # covariable -> radar
locations=lluvia.pluv.spdf,
newdata=plantilla.sp, # podr�a ser cualquier objeto Spatial
model=v.fit, # modelo de semivariograma.
maxdist=5000
)
�Me podr�ais ayudar?
Un saludo,
[[alternative HTML version deleted]]
_______________________________________________
R-help-es mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-help-es
