[R-sig-Geo] block and output grid size in block kriging

[Edzer Pebesma]

On 02/09/2015 08:48 PM, subash wrote:
> Dear Group,
> 
> I was using bock kriging  for groundwater level data. I was wondering, if
> the output.grd size and the block size can be different.

Yes, they can be different. Blocks are centered over grid cell centers,
and may overlap, be identical to grid cells, or smaller than grid cells
and have space between them.

> 
> For example, In the example for the case given below, I need block kriged
> average of 5000x5000. The output grid size given is 1000m*1000m. Is it
> correct?
> 
> Or should the output.grd size also be 5000mx5000m as in block size?
> 
> 
> # IMP : There are the centre coordinates for block kriging and not corner
> coordiantes.
> Xmin = 581000 + 2500
> /Xmax = 711000 - 2500
> Ymin = 1266000 + 2500
> Ymax = 1371000 - 2500
> 
> output.grd = expand.grid(x=seq(from=Xmin, to=Xmax,by=1000),y=seq(from=Ymin,
> to=Ymax, by=1000))
> coordinates(output.grd) = c("x", "y")
> gridded(output.grd) = TRUE
> 
> krig =
> krige(borewell.piezo[,3]~1,location~Easting+Northing,data=borewell.piezo,newdata=output.grd,model
> = vario.fit,block=c(5000,5000))/
> 
> 
> Thanks & Regards,
> Subash
> 
> 
> 
> --
> View this message in context: http://r-sig-geo.2731867.n2.nabble.com/block-and-output-grid-size-in-block-kriging-tp7587770.html
> Sent from the R-sig-geo mailing list archive at Nabble.com.
> 
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
> 

-- 
Edzer Pebesma
Institute for Geoinformatics, University of Münster,
Heisenbergstraße 2, 48149 Münster, Germany; +49 251 83 33081
Journal of Statistical Software:   http://www.jstatsoft.org/
Computers & Geosciences: http://elsevier.com/locate/inca/398

-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 490 bytes
Desc: OpenPGP digital signature
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20150209/814ae8de/attachment.bin>