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

[subash]

Dear Edzer,

Thanks for quick reply. A question still lingering in my mind is below

In a similar lines I was using block kriging of rainfall with block size of
25x25Km. The output grid  chosen is 1x1Km. This is resulting in kriging
estimate computed at 1x1Km . However, my interest is to get 25x25Km kriging
estimate.

The approach I tried was setting the output grid mesh as 25Km, but this is
wrong as the block average would be arrived from 1 point only.

output.grd = expand.grid(x=seq(from=tmp$min[1],
to=tmp$max[1],by=1),y=seq(from=tmp$min[2], to=tmp$max[2], by=1))
coordinates(output.grd) = c("x", "y")
gridded(output.grd) = TRUE

block_Krig = krige(stn.rain.data[,3]~1,location = ~easting+northing, data =
stn.rain.data,newdata=output.grd,model = clim.vrmod,block=c(25,25)) 

What is that I am missing.


Thanks in advance
Subash



--
View this message in context: http://r-sig-geo.2731867.n2.nabble.com/block-and-output-grid-size-in-block-kriging-tp7587770p7587773.html
Sent from the R-sig-geo mailing list archive at Nabble.com.