Auxiliary variables that are used to explain the trend-part of variation need
to be available also
at all new prediction locations.
see ?krige:
newdata - data frame or Spatial object with prediction/simulation locations;
should contain
attribute columns with the independent variables (if
sorry, a type error:
cellsize - 10 # see also http://spatial-analyst.net/pixel.php
xLL - round(slot(aSHAPE,bbox)[1])+cellsize/2
xUR - round(slot(aSHAPE,bbox)[3])-cellsize/2
yLL - round(slot(aSHAPE,bbox)[2])+cellsize/2
yUR - round(slot(aSHAPE,bbox)[4])-cellsize/2
xN -