[R-sig-Geo] How to produce a R2 for semivariance model fitting

2008-03-17 Thread Yong Li
Hi ALL, I try to use gstat to do spatial analyse. But I have difficulty to obtain the R2 for the fitted model of semivariance. Here I used R2 = 1 - attributes(model)$SSErr/sum((variogram$gamma-mean(variogram$gamma))^2), but the R2 is always close to 1.0. Any help would be greatly appreciated.

[R-sig-Geo] Finding closest match to GIS coordinates - R code

2008-03-17 Thread Mary.Barnes
Hi all, I recently joined, and have a rather basic question. I'd appreciate any help to find some R routines to find the closest match to GIS coordinates. I suspect simply using the Euclidean distance. I could write some code but I suspect there is some already. I describe three classes of difficu

Re: [R-sig-Geo] smoother contour maps

2008-03-17 Thread Michael Sumner
I once used Gabor Grothendieck's trick using apply across an array-redimmed matrix: http://tolstoy.newcastle.edu.au/R/help/06/07/31720.html You would have convert XYZ table data to image list form, and it performs a resize, but it might be worth pursuing. Cheers, Mike. ==Origina

[R-sig-Geo] Variogram Plot

2008-03-17 Thread PUJAN RAJ REGMI
Dear all, I am trying to plot semivarigram using the function "variog" but i am getting following error message "cannot allocate vector of size 757.6 Mb" I suppose my data set seems large?? it has 14090 samples. I will be thankful if anybody suggest me the alternative ways to overcome the probl

[R-sig-Geo] smoother contour maps

2008-03-17 Thread Fernando De Sales
Hello everyone, I was wondering if anyone knows how to smooth out a contour-line plot. I am currently using the following set of commands to create my plot filled.contour(xx,yy, PP, color=cm.colors, plot.axes = { contour(xx, yy, P2,add=T, col="grey",drawlabels=FALSE) axis(1,

Re: [R-sig-Geo] rgdal installation

2008-03-17 Thread Paul Hiemstra
Hi, I had a similar problem with reading HDF5 files. The standard gdal version compiled into rgdal probably does not support hdf4. One option is to compile rgdal yourself using a gdal and proj version that support hdf4. A good alternative is to translate the hdf4 file to tif (using shell() to

Re: [R-sig-Geo] spplot with two rasters

2008-03-17 Thread Dylan Beaudette
On Wednesday 12 March 2008, Roger Bivand wrote: > On Wed, 12 Mar 2008, Dylan Beaudette wrote: > > On Thursday 06 March 2008, Thomas Adams wrote: > >> Dylan, > >> > >> I think a solution using GRASS can be found on pages 110-111 of "Open > >> Source GIS: A GRASS GIS Approach", 3rd Ed. The same mater

Re: [R-sig-Geo] gdal error 1 (execution halted)

2008-03-17 Thread Roger Bivand
Dear Simon, On Mon, 17 Mar 2008, sj wrote: > Dear Roger > > Thanks for your detailed reply. > >> > I would like to use R from my bash console with the following command: >> > >> > /path/to/R --vanilla --slave --args < /data/myscript.R >> > >> > this script invokes a range of functions and tools

Re: [R-sig-Geo] rgdal installation

2008-03-17 Thread Michael Sumner
Hello, to find which drivers are available to you in rgdal use gdalDrivers(). These won't be as many as are available in FWTools, unless you build rgdal from source against FWTools itself. For HDF you must build rgdal against your own GDAL installation. I've done this with VC++ - see the notes

[R-sig-Geo] rgdal installation

2008-03-17 Thread David PINAUD
Dear all, I'm trying to import HDF (v4) files on R2.6.2 using the package rgdal and FWTools v2.1.0. I have the following message : > GDALinfo("T20060012006008.L3m_8D_SST_9.hdf") Erreur dans .local(.Object, ...) : GDAL Error 4: `T20060012006008.L3m_8D_SST_9.hdf' not recognised as a sup