Jan,
If you have the boundaries of Germany and you want to make your own grid, you can use the overlay method to 'cut' the grid so that you only get the points inside Germany. You can use this object in gstat to make your predictions. Hope this helps, Virgilio El mar, 24-02-2009 a las 18:21 +0100, Jan Hackenberg escribió: > Hello > im doing a project with Interpolating ( Kriging ) Temperature Data from many > Stations in Germany. To do a Kriging i have to use a Grid. So I have now 2 > options. The better one is too download a prepared grid with the borders of > Germany. Then the interpolated map will look great, because i have only > Points inside German borders. So does anyone have got such a grid? > If no, the other is to make a grid inside a German Bounding Box. Then my map > will be a quare :). I have looked in the meuse Data set how such a grid > looks like : > > meuse.grid > Object of class SpatialPixelsDataFrame > Object of class SpatialPixels > Grid topology: > cellcentre.offset cellsize cells.dim > x 178460 40 78 > y 329620 40 104 > SpatialPoints: > x y > [1,] 181180 333740 > [2,] 181140 333700 > [3,] 181180 333700 > and so on. > So is there any Possibility to make such a grid, if i know : Coordinates of > my Bounding Box (upper left coodinates for example (200000,540000), upper > right(900000,650000)) and Cellsize for Example 10000? > I woult prefer too download an existing grid. > Sorry for my bad english and bad knowledge of R. :) > Greetings Jan Hackenberg > > [[alternative HTML version deleted]] > > _______________________________________________ > R-sig-Geo mailing list > R-sig-Geo@stat.math.ethz.ch > https://stat.ethz.ch/mailman/listinfo/r-sig-geo _______________________________________________ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo
