[R-sig-Geo] creating irregularly-shaped grids for kriging

Pierre Roudier pierre.roudier at gmail.com
Mon Jul 20 10:53:08 CEST 2009


Hi Craig,

I have a workaround for this using the getpoly() inout()

2009/7/20, Scanlan, Craig <craig.scanlan at agric.wa.gov.au>:
> Hi all,
>
> I need to krige data that is irregularly spaced within irregularly
> shaped areas (agricultural fields).
>
> At present, I am creating a regular grid for newdata in
> krige(data,newdata....) by
>
> x.range <- as.integer(range(a2 at coords[,1]))
> y.range <- as.integer(range(a2 at coords[,2]))
> grd <- expand.grid(x=seq(from=x.range[1], to=x.range[2], by=10),
> y=seq(from=y.range[1], to=y.range[2], by=10) )
>
> I have not been able to find a method to trim the newdata grid to the
> same shape as the original data - is this possible?
>
> The task I am trying to achieve is the same as creating meuse.grid from
> meuse data.
>
> Can anyone help?
>
> Regards,
>
> Craig.
>
>
>
> This e-mail and files transmitted with it are privileged and confidential
> information
> intended for the use of the addressee. The confidentiality and/or privilege
> in this e-mail is
> not waived, lost or destroyed if it has been transmitted to you in error. If
> you received this
> e-mail in error you must
> (a) not disseminate, copy or take any action in reliance on it;
> (b) please notify the Department of Agriculture and Food, WA immediately by
> return e-mail to the sender;
> (c) please delete the original e-mail.
>
> This email has been successfully scanned by
> McAfee Anti-Virus software.
> Department of Agriculture and Food WA
>
> 	[[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>



More information about the R-sig-Geo mailing list