Hello,
I am having problems with kriging an ascii file DTM from ArcGIS which has many
-9999 (no data) values. Do these normally have to be removed first? I want to
make the elevation DTM the predictor base for kriging with external drift.
My code is the following:
elev <- read.asciigrid("elev.asc", as.image=FALSE, plot.image=TRUE)
str(elev)
elev_ked <- krige(meanRain~Zloc, locations=meanRain, newdata=elev,
model=fitted_vario)
Error in model.frame.default(terms.f, newdata, na.action = na.action) :
object is not a matrix
In addition: Warning message:
'newdata' had 550448 rows but variable(s) found have 1 rows
My data file meanRain has X Y Z coordinates and Precipitation values (4
columns).
Thank you so much for your time and help!
Cara
[[alternative HTML version deleted]]
_______________________________________________
R-sig-Geo mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
