Re: [R-sig-Geo] problem with NY8 <- readOGR(".", "NY8_utm18")

2014-11-24 Thread Pascal Oettli
Masha, Are you running the command in the correct directory? You should replace "." by the path to your file. Regards, Pascal On Tue, Nov 25, 2014 at 2:09 PM, Masha Pautrel wrote: > Dear all, > > I try to use NY8 data set, but I can't open the file. > > Here the errors message > >> library(sp)

[R-sig-Geo] problem with NY8 <- readOGR(".", "NY8_utm18")

2014-11-24 Thread Masha Pautrel
Dear all, I try to use NY8 data set, but I can't open the file. Here the errors message > library(sp) > library(rgdal) > NY8 <- readOGR(".", "NY8_utm18") Erreur dans ogrInfo(dsn = dsn, layer = layer, encoding = encoding, use_iconv = use_iconv, : Cannot open file > cities <- readOGR(".", "N

Re: [R-sig-Geo] Variable search radius for interpolation

2014-11-24 Thread Tomislav Hengl
I've implemented recently a function that helps speeding up kriging with spatially clustered data (for it fits kernel density function first and then sub-samples regular points proportionally to the density of observations). It comes together with the GSIF package: http://gsif.r-forge.r-proj

Re: [R-sig-Geo] questions on spatial interpolation exercise(NL)

2014-11-24 Thread Tomislav Hengl
You need to provide access to "monv_201310.txt" so that we can see what is the structure of that file. The original file with rainfall data from that URL is simply a text file from the KNMI website (without any tabular structure) so that e.g.: ... write.table(t(matrix(unlist(strsplit(tmp, ",

Re: [R-sig-Geo] subplots of maps in R using hydroTSM package

2014-11-24 Thread Zilefac Elvis
Many thanks, Jon. Just what I wanted. Zilefac. On Monday, November 24, 2014 2:48 AM, Jon Skoien wrote: Hi Zilefac, Your problem is not really related to hydroTSM, but to the fact that hydrokrige uses spplot for plotting, which is again based on the lattice package. This does not use the par(

[R-sig-Geo] Variable search radius for interpolation

2014-11-24 Thread Moshood Agba Bakare
Hi All, Could anyone guide me on how to define variable search radius in GSTAT while interpolating into a grid cell? The GSTAT documentation used MAXDIST to define fixed search radius. (2) Is there a function in GSTAT to determine the point density while interpolating into different grids size? Th

Re: [R-sig-Geo] shifting points in plot window

2014-11-24 Thread Edzer Pebesma
I believe plotting routines in raster were built from (or inspired by) those in package fields. The issue arises when the plotting area is split in two parts, one for the map and one for the legend, and the need to know the aspect ratio of the first when adding to the plot -- the legend width is fi

Re: [R-sig-Geo] shifting points in plot window

2014-11-24 Thread Michael Krabbe Borregaard
Thanks, This was very helpful. I have a package that depends on raster, and have included a note on this in the documentation. From: Michael Sumner: This is a known problem with no known fix. I've been trying to isolate it but without success. If anyone can see how to fix it we can help get i

Re: [R-sig-Geo] subplots of maps in R using hydroTSM package

2014-11-24 Thread Jon Skoien
Hi Zilefac, Your problem is not really related to hydroTSM, but to the fact that hydrokrige uses spplot for plotting, which is again based on the lattice package. This does not use the par()-arguments. Instead there are two other ways for producing multiple plots in one figure: - Save each of

Re: [R-sig-Geo] questions on spatial interpolation exercise(NL)

2014-11-24 Thread Edzer Pebesma
Mengxi, please provide a minimal script that reproduces your error. On 11/23/2014 11:51 PM, Mengxi Yang wrote: > Dear list, > I have a question on write.table. I found an nice exercise on spatial > interpolation from > http://spatial-analyst.net/wiki/index.php?title=Spatial_interpolation_exercises

Re: [R-sig-Geo] subplots of maps in R using hydroTSM package

2014-11-24 Thread Edzer Pebesma
It looks like hydroTSM uses sp::spplot to plot maps; this function uses lattice which builds on grid graphics, and hence mfrow does not work. As hydrokrige does not return the trellis plot object but merely plots it, it's not easy to arrange plots using lattice::print.trellis It seems relatively

Re: [R-sig-Geo] R crashes performing gstat IDW with omax argument

2014-11-24 Thread Edzer Pebesma
Thanks for the reproducible and small example! I found the bug and committed a fix to r-forge. For the the version that causes this problem, you can use omax when in addition you specify nmax to the same value (but pls check that the output makes sense). On 11/24/2014 09:20 AM, Tim Appelhans wro