"predict.krig.image" <-
function (obj, x, ...) 
{
    if (missing(x)) {
        temp1 <- c(obj$cov.function(obj$indexM, obj$indexM, Y = obj$delta, 
            cov.obj = obj$cov.obj))
        temp1 <- temp1[obj$rep.info]
        temp2 <- obj$beta[1] + obj$x %*% obj$beta[2:3]
        return(c(temp1 + temp2))
    }
    else {
        predict.interp.surface(obj$surface, x)
    }
}