Re: [R] ascii-grid export

2013-11-06 Thread Barry Rowlingson
On Mon, Nov 4, 2013 at 7:27 AM, Enzo Cocca wrote: > yes barry I really need this. > > I tried to use raster or rgdal but with poor results. > > I have this function: > > VGM_PARAM_A3 <- gstat(id="bos_bison", > formula=combusto~1,locations=~coord_x+coord_y, data=archezoology_table, > nmax = 10) > >

Re: [R] ascii-grid export

2013-11-03 Thread Enzo Cocca
yes barry I really need this. I tried to use raster or rgdal but with poor results. I have this function: VGM_PARAM_A3 <- gstat(id="bos_bison", formula=combusto~1,locations=~coord_x+coord_y, data=archezoology_table, nmax = 10) VGM_PARAM_A3 <- gstat(VGM_PARAM_A3, model=vgm(1, "Sph", 5, 0), fill.

Re: [R] ascii-grid export

2013-11-02 Thread Barry Rowlingson
Or do you mean you want to write the gridded output of an interpolation you did (perhaps using kriging) in gstat as an ESRI ASCII Grid file for reading into a GIS? If so, you can probably do it with writeGDAL from the rgdal package, or writeRaster from the raster package. I don't really know what

Re: [R] ascii-grid export

2013-11-02 Thread Gergely Daróczi
If I understand what you are up to, my "pander" package might help. You mean something like this? I used the example form the package manual. > library(pander) > data(meuse) > coordinates(meuse) = ~x+y > pander(variogram(log(zinc)~1, meuse)) np dist