R语言 sgeostat包 krige()函数中文帮助文档(中英文对照)

loveR

krige(sgeostat)
krige()所属R语言包：sgeostat

                                        Kriging
                                         克立格法

                                         译者：生物统计家园网 机器人LoveR

描述----------Description----------

Carry out spatial prediction (or kriging).
进行空间预测（克里格法）。


用法----------Usage----------


krige(s, point.obj, at, var.mod.obj, maxdist=NULL, extrap=FALSE, border)



参数----------Arguments----------

参数：s
a point object, generated by point(), at which prediction is carried out
点对象，所产生的point()，在预测进行


参数：point.obj
a point object, generated by point(), containing the sample points and data
所产生的point()，包含的样本点和数据点对象，


参数：at
the variable, contained in point.obj, for which prediction will be carried out
的变量，包含在point.obj，预测将进行


参数：var.mod.obj
variogram object
变异函数对象


参数：maxdist
an optional maximum distance. If entered, then only sample points (i.e, in point.obj) within maxdist of each prediction point will be used to do the prediction at that point. If not entered, then all n sample points will be used to make the prediction at each point.
一个可选的最大距离。如果输入，然后只样本点（即，在point.obj）内的预测点maxdist每个将被用来做什么的预测，在这一点上。如果没有输入，则所有的n个采样点将用于使在每个点的预测。


参数：extrap
logical, indicates if prediction outside the convex hull of data points should be done, default FALSE
逻辑，表示如果预测的数据点的凸包外应该做的，默认的FALSE


参数：border
optional polygon (list with two components x and y of same length) representing a (possibly non convex) region of interest to  be used instead of the convex hull. Needs extrap=TRUE.
可选的多边形（带有两个组件的列表x和y相同的长度），代表一个（可能的非凸）的感兴趣的区域的凸包，而不是使用。需要extrap=TRUE。


值----------Value----------

A point object which is a copy of the s object with two new variables, zhat and sigma2hat, which are,
A点对象，它是两个新的变量，s和zhat，这是sigma2hat对象的副本，


参考文献----------References----------



参见----------See Also----------

est.variogram,fit.variogram
est.variogram，fit.variogram


实例----------Examples----------



# a single point:[单点：]
prdpnt <- point(data.frame(list(x=180000,y=331000)))
prdpnt <- krige(prdpnt, maas.point, 'zinc', maas.vmod)
prdpnt

# kriging on a grid (slow!)[克里格法对网格（slow!的）]
grid <- list(x=seq(min(maas$x),max(maas$x),by=100),
             y=seq(min(maas$y),max(maas$y),by=100))
grid$xr <- range(grid$x)
grid$xs <- grid$xr[2] - grid$xr[1]
grid$yr <- range(grid$y)
grid$ys <- grid$yr[2] - grid$yr[1]
grid$max <- max(grid$xs, grid$ys)
grid$xy <- data.frame(cbind(c(matrix(grid$x, length(grid$x), length(grid$y))),
             c(matrix(grid$y, length(grid$x), length(grid$y), byrow=TRUE))))
colnames(grid$xy) <- c("x", "y")
grid$point <- point(grid$xy)
data(maas.bank)
grid$krige <- krige(grid$point,maas.point,'zinc',maas.vmod,
                    maxdist=1000,extrap=FALSE,border=maas.bank)
op <- par(no.readonly = TRUE)
par(pty="s")
plot(grid$xy, type="n", xlim=c(grid$xr[1], grid$xr[1]+grid$max),
                    ylim=c(grid$yr[1], grid$yr[1]+grid$max))
image(grid$x,grid$y,
      matrix(grid$krige$zhat,length(grid$x),length(grid$y)),
      add=TRUE)
contour(grid$x,grid$y,
        matrix(grid$krige$zhat,length(grid$x),length(grid$y)),
        add=TRUE)
data(maas.bank)
lines(maas.bank$x,maas.bank$y,col="blue")
par(op)

转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。


注：
注1：为了方便大家学习，本文档为生物统计家园网机器人LoveR翻译而成，仅供个人R语言学习参考使用，生物统计家园保留版权。
注2：由于是机器人自动翻译，难免有不准确之处，使用时仔细对照中、英文内容进行反复理解，可以帮助R语言的学习。
注3：如遇到不准确之处，请在本贴的后面进行回帖，我们会逐渐进行修订。