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]]

