Re: [R-sig-Geo] boundaries for variogram

-- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 253 5773 http://intamap.geo.uu.nl/~paul http://nl.linkedin.com/pub/paul-hiemstra/20/30b/770

[R-sig-Geo] Adding a scale bar and north arrow to a ggplot

] ggplot2_0.8.7 digest_0.4.2 reshape_0.8.3 plyr_0.1.9proto_0.3-8 [6] sp_0.9-62 loaded via a namespace (and not attached): [1] lattice_0.19-13 hiems...@fg-113:~$ uname -a Linux fg-113 2.6.32-21-generic #32-Ubuntu SMP Fri Apr 16 08:10:02 UTC 2010 i686 GNU/Linux -- Drs. Paul Hiemstra Department

Re: [R-sig-Geo] sp.text background

from inkscape. At least in the version I was using. cheers, Paul -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 253 5773 http://intamap.geo.uu.nl/~paul http://nl.linkedin.com/pub/paul

Re: [R-sig-Geo] Plotting a tif file

and try to use in the aes() argument (col = col, where col is the name of the column with the color). cheers, Paul On 11/06/2010 03:03 PM, Paul Hiemstra wrote: Gpplot cannot plot SpatialPixelsDataFrame's in one go. You need to convert the grid to a data frame. An example using the meuse

Re: [R-sig-Geo] Plotting a tif file

don't understand which information should I provide to the new column. Best, Manuel On 07/11/2010 04:00 a.m., Paul Hiemstra wrote: Hi Manuel, Please do not respond to me alone, include the mailing list in the cc. You waste the option of someone else answering your question

Re: [R-sig-Geo] Plotting a tif file

So what is in the three bands? Paul On 11/07/2010 02:06 PM, Manuel Spínola wrote: Species richness (from 3 to 20 more or less) On 07/11/2010 06:46 a.m., Paul Hiemstra wrote: What is the data captured by the tif file? Paul On 11/07/2010 12:45 PM, Manuel Spínola wrote: Thank you very much

Re: [R-sig-Geo] Plotting a tif file

Honestly Paul, I don't know. The tif is a result from the kriging in Quantum GIS. Is very strange, right? Manuel On 07/11/2010 07:03 a.m., Paul Hiemstra wrote: So what is in the three bands? Paul On 11/07/2010 02:06 PM, Manuel Spínola wrote: Species richness (from 3 to 20 more or less

Re: [R-sig-Geo] Plotting a tif file

://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 253 5773 http://intamap.geo.uu.nl/~paul http://nl.linkedin.com/pub/paul-hiemstra/20/30b

Re: [R-sig-Geo] Plotting a tif file

much Paul. With plot I got a black polygon and with spplpot I got: spplot(riquezakrig) Error: cannot allocate vector of size 22.6 Mb with image show me what I want but no legend. Any specific geom to plot with ggplot? Best, Manuel On 06/11/2010 06:46 a.m., Paul Hiemstra wrote: On 11/06

Re: [R-sig-Geo] Plotting a tif file

specific geom to plot with ggplot? Best, Manuel On 06/11/2010 06:46 a.m., Paul Hiemstra wrote: On 11/06/2010 11:15 AM, Manuel Spínola wrote: Dear list members, I read a .tif file with readGDAL. What is the best way to display it?. spplot or ggplot I tried function image

Re: [R-sig-Geo] fix color scale

___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC

Re: [R-sig-Geo] Basic Info. requried

/mailman/listinfo/r-sig-geo -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 253 5773 http://intamap.geo.uu.nl/~paul http://nl.linkedin.com/pub/paul-hiemstra/20/30b/770

Re: [R-sig-Geo] number points interpolation

]] ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 253

Re: [R-sig-Geo] gstat: choice of inverse distance power in idw() function

. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 253 5773 http://intamap.geo.uu.nl/~paul http://nl.linkedin.com/pub/paul-hiemstra/20/30b/770 ___ R-sig

Re: [R-sig-Geo] Function for subsampling spatial data

HTML version deleted]] ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O

Re: [R-sig-Geo] gstat fit.variogram, attr(fit.1,SSErr)

in SSerr is caused by a different fit by fit.variogram, which in turn is caused by a different initial variogram model guess. regards, Paul -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone

Re: [R-sig-Geo] Mapping multiple attributes at once

version deleted]] ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O

Re: [R-sig-Geo] spRbind for spplot

and not even more cryptical... Best regards, Klaus On Tue, 17 Aug 2010 20:21:43 +0200 Paul Hiemstra p.hiems...@geo.uu.nl wrote: Dear Klaus, You need to be more specific as to what you want exactly. From what I understand both your runoff and predictions are attributes that you want to use to fill

Re: [R-sig-Geo] spRbind for spplot

advice, Klaus -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 253 5773 http://intamap.geo.uu.nl/~paul http://nl.linkedin.com/pub/paul-hiemstra/20/30b/770

Re: [R-sig-Geo] Error autoKrige

can't run the script myself to see what goes wrong. cheers, Paul -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 253 5773 http://intamap.geo.uu.nl/~paul http://nl.linkedin.com/pub/paul

Re: [R-sig-Geo] Error autoKrige

On 08/11/2010 05:42 PM, kad wrote: I use automap 1.0-7 with version R 2.11.0. My data http://r-sig-geo.2731867.n2.nabble.com/file/n5412713/k.txt k.txt That's your data, not an example of R code using the data, resulting in your error. regards, Paul -- Drs. Paul Hiemstra Department

Re: [R-sig-Geo] how to do randomly sampling in raster layer

___ 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 -- Drs. Paul Hiemstra Department

Re: [R-sig-Geo] scale semivariance befor fit.variogram

]] ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130

Re: [R-sig-Geo] plot.autoKrige and standard breaks for coloring

list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo Hi Mark, Use the at parameter to set the breaks. cheers, Paul -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone

Re: [R-sig-Geo] plot.autoKrige and standard breaks for coloring

AM, Paul Hiemstra wrote: On 07/07/2010 11:55 PM, Mark Connolly wrote: plot.autoKrige does a nice job of laying out the predicted surface, the standard error and the fitted variogram (over the sample variogram). Is there a way to set the breaks used for coloring the prediction surface so I can

Re: [R-sig-Geo] spatial data

, USA ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508

Re: [R-sig-Geo] How to locate local minima points on a Krige-fitted surface?

with security by Microsoft. https://signup.live.com/signup.aspx?id=60969 -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue Phone: +3130 253 5773 Wed-Fri http

Re: [R-sig-Geo] How to extract coordinates values from a shapefile?

On 06/11/2010 11:43 AM, Paul Hiemstra wrote: On 06/09/2010 11:01 PM, Dylan Beaudette wrote: On Wednesday 09 June 2010, Hadley Wickham wrote: On Wed, Jun 9, 2010 at 3:24 PM, Edzer Pebesma edzer.pebe...@uni-muenster.de wrote: The example provided by Matt assumes that each polygon consists

Re: [R-sig-Geo] How to extract coordinates values from a shapefile?

other attached packages: [1] maptools_0.7-34 lattice_0.18-5 sp_0.9-62 foreign_0.8-40 loaded via a namespace (and not attached): [1] grid_2.11.0 Cheers, Dylan -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O

Re: [R-sig-Geo] How to locate local minima points on a Krige-fitted surface?

-- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue Phone: +3130 253 5773 Wed-Fri http://intamap.geo.uu.nl/~paul http://nl.linkedin.com/pub/paul-hiemstra/20/30b/770

Re: [R-sig-Geo] test for CSR

-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue Phone: +3130 253 5773 Wed-Fri http

Re: [R-sig-Geo] How to extract coordinates values from a shapefile?

also works well in terms of performance. cheers, Paul Cheers, Dylan -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue Phone: +3130 253 5773 Wed-Fri http

Re: [R-sig-Geo] from R to ArcMap

/marinescotland [[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 -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences

Re: [R-sig-Geo] Parameters in fitting variogram automatically in gstat

On 05/10/2010 06:52 AM, Kai Zhang wrote: Hi folks, I am working on a problem dealing with fitting variogram automatically due to a large dateset. Gstat R package is used to conduct this work. Part of this work is to compare performance of different kriging methods. The comparison criteria

Re: [R-sig-Geo] negative range from fit.variogram through autofitVariogram

=rescale(sparse_df$theta_percent)) ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht

Re: [R-sig-Geo] negative range from fit.variogram through autofitVariogram

I just uploaded a new version of automap (1.07) that fixes this problem. It deletes fitted variogram models with a negative sill/range/nugget. cheers, Paul Paul Hiemstra wrote: Hi Mark, Thanks for the reproducible example. The problem is that when I look at the sample variogram

Re: [R-sig-Geo] Calculate area of raster cells

) = TRUE bla = gridparameters(meuse.grid) area = bla$cellsize[1] * bla$cellsize[2] area cheers, Paul -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue Phone: +3130 253 5773

[R-sig-Geo] Plotting SpatialPolygonsDataFrame using ggplot2

to the problem? cheers, Paul -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue Phone: +3130 253 5773 Wed-Fri http://intamap.geo.uu.nl/~paul http://nl.linkedin.com/pub/paul

Re: [R-sig-Geo] Plotting SpatialPolygonsDataFrame using ggplot2

Hi Hadley, That worked like a charm. Thanks. cheers, Paul hadley wickham wrote: Hi Paul, You're missing the group specification: ggplot(dum, aes(x = long, y = lat)) + geom_path(aes(group = group)) Hadley On Tue, Apr 27, 2010 at 4:17 AM, Paul Hiemstra p.hiems...@geo.uu.nl wrote: Hi

Re: [R-sig-Geo] areal data to pixel image

version deleted]] ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115

Re: [R-sig-Geo] automap model defaults

https://stat.ethz.ch/mailman/listinfo/r-sig-geo Hi Mark, These models are the models that are most often used, so it seemed logical to use these ones. But why do you want to include the linear model specifically? cheers, Paul -- Drs. Paul Hiemstra Department of Physical Geography Faculty

Re: [R-sig-Geo] automap model defaults

handle this. cheers, Paul -- Edzer Mark Connolly wrote: On 04/15/2010 07:42 AM, Paul Hiemstra wrote: Mark Connolly wrote: I was wondering what is behind the selection of the default models tried by autofitVariogram (model=c(Sph,Exp,Gau,Ste)). For example, why was linear model

Re: [R-sig-Geo] a function to calculate distance

at the spDists and spDistsN1 functions from the sp-package. cheers, Paul -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue Phone: +3130 253 5773 Wed-Fri http

Re: [R-sig-Geo] Reversed raster after readGDAL()

___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue Phone

Re: [R-sig-Geo] help with areal interpolation I think

://www.scopus.com/inward/record.url?eid=2-s2.0-2497904partnerID=40md5=88b8634524719043b9fad1bb71f2c9ca} } You could take a look at the papers and see what you can use. cheers, Paul -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan

Re: [R-sig-Geo] ordinary kriging and elevation with GSTAT

. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue Phone: +3130 253 5773 Wed-Fri http://intamap.geo.uu.nl/~paul ___ R-sig-Geo mailing

Re: [R-sig-Geo] Tweak legend in image.plot()

any suggestions? Muhammad -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue Phone: +3130 253 5773 Wed-Fri http://intamap.geo.uu.nl/~paul

Re: [R-sig-Geo] Overlay-ing two grids

All the best, Renaud -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue Phone: +3130 253 5773 Wed-Fri http://intamap.geo.uu.nl/~paul

Re: [R-sig-Geo] Intersect polygon with grid

]] ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone

Re: [R-sig-Geo] How to automatically fit variogram in gstat

-- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue Phone: +3130 253 5773 Wed-Fri http://intamap.geo.uu.nl/~paul ��i

Re: [R-sig-Geo] How to automatically fit variogram in gstat

and best regards, Kai __ 赶快注册雅虎超大容量免费邮箱? ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Drs. Paul Hiemstra Department

Re: [R-sig-Geo] plotting 100 variograms on the same plot

with the xyplot function? Thank you very much! Laura 2009/12/4 Paul Hiemstra p.hiems...@geo.uu.nl mailto:p.hiems...@geo.uu.nl Paul Hiemstra wrote: Laura Poggio wrote: Dear all, I would like to plot 100 variograms on the same plot, possibly using

Re: [R-sig-Geo] problem saving spplot in a loop

Hi, This is a byproduct of the fact that under the hood spplot uses lattice commands (levelplot and xyplot) to make the plots. This behavior is a FAQ for lattice commands, and thus spplot. Even shorter than Henk's suggestion is to use: png(filename=...) print(spplot(...)) dev.off() cheers,

Re: [R-sig-Geo] plotting 100 variograms on the same plot

(...) { panel.xyplot(...) # points of the sample var ret = variogramLine(model_list[[panel.number()]], maxdist = max(sv_list$dist)) llines(ret$dist, ret$gamma, col = red) }) cheers, Paul -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht

Re: [R-sig-Geo] plotting 100 variograms on the same plot

Paul Hiemstra wrote: Laura Poggio wrote: Dear all, I would like to plot 100 variograms on the same plot, possibly using a loop. At the moment I am using autofitVariogram in the package Automap to fit the variograms automatically. I found some code calling the lattice package ( http

Re: [R-sig-Geo] Still error exporting ascii grid generated by DTM grid, please any suggestions?

:22337 NA's :22337.000 ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O

Re: [R-sig-Geo] 'Customizing' spsample

: namespace:sp -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue Phone: +3130 253 5773 Wed-Fri http://intamap.geo.uu.nl/~paul

Re: [R-sig-Geo] Newbie question on shape files and mapping

no liability for any damage caused by any such code transmitted by or accompanying this e-mail or any attachment. ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Drs. Paul Hiemstra Department

Re: [R-sig-Geo] Writing an ASCII grid with nodata values

to output floating point data to an ESRI ASCII grid file and specify a nodata value. Thanks, Ned ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Drs. Paul Hiemstra Department of Physical Geography

Re: [R-sig-Geo] 'LDLfactor' error in 'krige' function

and if it could be overcome in some way. Thanks in advance for any help. Mauricio -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue Phone: +3130 253 5773 Wed-Fri

Re: [R-sig-Geo] newbie questions

saldanha.plan...@gmail.com [[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 -- Drs. Paul Hiemstra Department of Physical Geography Faculty

Re: [R-sig-Geo] newbie questions - maps

-free GIS software, by the way. My cartographer colleagues who use ArcGIS generally haul mapping output from there into other packages for font and other effect rendering. Yours, Ashton -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht

Re: [R-sig-Geo] sub-selecting spatial data by distance from a different point(set)

___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130

Re: [R-sig-Geo] sub-selecting spatial data by distance from a different point(set)

Sorry forgot: spDistsN1 is from the sp package. cheers, Paul Paul Hiemstra wrote: Hi Sean, Take a look at the spDistsN1 function to calculate the distances from a single point to all other points. You can then make the subset based on the distance. What could also be an options is to build

Re: [R-sig-Geo] Variance unchanging with krige in gstat with iterations

[[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 -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan

Re: [R-sig-Geo] Various Variograms of gstat package together in a single page

=distance, ylab=semivariance) } In advance, thank you very much. Eduardo. -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue Phone: +3130 253 5773 Wed

Re: [R-sig-Geo] Spatial Panel data models

https://stat.ethz.ch/mailman/listinfo/r-sig-geo Hi, I'm not sure what spatial panel data models are, but that may be my problem :). A start to look for packages to do anything spatial in R is the Spatial Task View. http://cran.r-project.org/web/views/Spatial.html cheers, Paul -- Drs. Paul

Re: [R-sig-Geo] Spatial Panel data models

of applications of spatial panel data models in R. Thanks Juan ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Drs. Paul

Re: [R-sig-Geo] Export of an irregular grid

50103 Cherbourg Cedex Phone (Direct) : (00 33)(0)233887342 Fax : (00 33)(0)233887339 ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Drs. Paul

Re: [R-sig-Geo] variogram estimation with adaptive sampling weights

by it. __ The Shaw Group Inc. http://www.shawgrp.com [[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 -- Drs. Paul Hiemstra Department of Physical Geography

Re: [R-sig-Geo] plot points in grid-plot with sp.layout

://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue Phone: +3130 253 5773 Wed-Fri http://intamap.geo.uu.nl/~paul

Re: [R-sig-Geo] plot points in grid-plot with sp.layout

Paul Hiemstra wrote: Tom Gottfried wrote: Hi list, the last plot in the following code ends up with error messages in the plot: data(meuse.grid) data(meuse) coordinates(meuse.grid) - ~x+y coordinates(meuse) - ~x+y meuse.grid - as(as(meuse.grid, SpatialPixelsDataFrame

Re: [R-sig-Geo] maps with R

. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue Phone: +3130 253 5773 Wed-Fri http://intamap.geo.uu.nl/~paul ___ R-sig-Geo mailing

Re: [R-sig-Geo] Very basic figure question

-- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue Phone: +3130 253 5773 Wed-Fri http://intamap.geo.uu.nl/~paul ___ R-sig-Geo

[R-sig-Geo] Strange behavior of spsample when sampling from Sgrid

, dist_from_meuse2, dist_from_meuse2_pt2)) But this is ofcourse not ideal. Is there a way to force spsample to always sample from top left to bottom right? cheers, Paul -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht

Re: [R-sig-Geo] Strange behavior of spsample when sampling from Sgrid

;LC_IDENTIFICATION=C attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] sp_0.9-44 loaded via a namespace (and not attached): [1] grid_2.9.2 lattice_0.17-20 tools_2.9.2 Paul Hiemstra wrote: Dear list, I spotted

Re: [R-sig-Geo] Average kriging surface over region

per polygon. Or alternatively, look at the overlay() function. cheers, Paul -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue Phone: +3130 253 5773 Wed-Fri http

Re: [R-sig-Geo] Problem with gstat variogram estimation

-- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue Phone: +3130 253 5773 Wed-Fri http://intamap.geo.uu.nl/~paul

Re: [R-sig-Geo] Finding to which Voronoi/Dirichlet a point belongs to

[[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 -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht

Re: [R-sig-Geo] Determine the location of the actual R-file

@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue Phone: +3130 253 5773 Wed-Fri http://intamap.geo.uu.nl/~paul

Re: [R-sig-Geo] Importing ESRI Shapefile into RSAGA

]] ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue

Re: [R-sig-Geo] How to make meuse data generates meuse.grid data?

SpatialPixels( Spatialpoints(makegrid(meuse))).I found using this method, the results do not the same as meuse.grid. Could you define not the same. cheers, Paul could you tell me the details? Sincerely Jessica -- Drs. Paul Hiemstra Department of Physical Geography Faculty

Re: [R-sig-Geo] [R-sig-geo] Interpolation of measures with measurement errors

://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue Phone: +3130 253 5773 Wed-Fri http://intamap.geo.uu.nl/~paul

Re: [R-sig-Geo] Zooming into a specific area

it but it did not zoom into the area but just blocked out the areas around xlim and ylim. Are there other ways? --- On Wed, 9/23/09, Paul Hiemstra p.hiems...@geo.uu.nl wrote: From: Paul Hiemstra p.hiems...@geo.uu.nl Subject: Re: [R-sig-Geo] Zooming into a specific area To: Julius Tesoro

Re: [R-sig-Geo] Zooming into a specific area

, You can use the xlim and ylim parameters in spplot to zoom in. cheers, Paul -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue Phone: +3130 253 5773 Wed-Fri http

Re: [R-sig-Geo] Convert SpatialPointsDataFrame to a SpatialGridDataFrame

this. Cheers, Brendan Malone PhD candidate University of Sydney [[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 -- Drs. Paul Hiemstra

Re: [R-sig-Geo] spplot without legend

Tom Gottfried schreef: Hi list, is there any way to draw a map with spplot() WITHOUT a legend? By the way: is there any documentation available for the function mapLegendGrob()? I mean besides the usage and arguments section in the relevant help page. I'm sorry I can't find it. Thanks in

Re: [R-sig-Geo] Fw: Help

Hi, Loading the meuse.grid dataset is done using: data(meuse.grid) to generate meuse.grid from meuse use something like: data(meuse) coordinates(meuse) = ~x+y meuse_grid = spsample(meuse, type = regular, cellsize = c(40,40)) gridded(meuse_grid) = TRUE cheers and good luck, Paul Danlin Yu

Re: [R-sig-Geo] spplot question

Roger Bivand schreef: On Fri, 11 Sep 2009, Eric Morway wrote: Edzer, I apologize for the confusion. Is there an option in spplot that fills in blank/white cells (I'm assuming they are blank because they have a value less than the minimum cut) with the minimum color on the legend bar (in

Re: [R-sig-Geo] Multiple variograms in one graph

HTML version deleted]] ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Drs. Paul Hiemstra Department

Re: [R-sig-Geo] Raster to points

Cutberto Uriel Paredes Hernández schreef: Dear list, Is there a way to convert a raster to points in R (SpatialGridDataFrame to SpatialPointsDataFrame)? I know I can do this with ArcGIS but I would like not to use propietary software. Sorry for posting this basic question but it couldn't find

Re: [R-sig-Geo] majority filter for grids

-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue Phone: +3130 253 5773 Wed-Fri http

Re: [R-sig-Geo] Creating grid from xyz file

Mehdi Khan schreef: Hello everyone, I have a table with 2 columns containing geographical coordinates (lat, long) and one column containing attribute data. This is a data frame object which I would like to convert into a grid object. It is the 7.3 mb california file on this page:

Re: [R-sig-Geo] Creating a Grid out of a matrix of coordinates

, tolerance, round, fuzz.tol) : dimension 1 : coordinate intervals are not constant What do I do? Thanks a lot! On Fri, Aug 21, 2009 at 12:34 PM, Paul Hiemstra p.hiems...@geo.uu.nl mailto:p.hiems...@geo.uu.nl wrote: Hi, If your object's points are regularly spaced and in a SpatialPoints

Re: [R-sig-Geo] Creating a Grid out of a matrix of coordinates

Hi, If your object's points are regularly spaced and in a SpatialPoints object, use: gridded(myobj) = TRUE Or you can from a data.frame with the coordinates and attributes to a grid by: coordinates(myobj) = ~x+y gridded(myobj) = TRUE or shorthand: gridded(myobj) = ~x+y cheers, Paul

Re: [R-sig-Geo] daily precipitation idw kriging and beginner questions

mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue Phone: +3130 253 5773 Wed

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

software. Department of Agriculture and Food WA [[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 -- Drs. Paul Hiemstra Department of Physical Geography

Re: [R-sig-Geo] [SoC v.autokrige2.py] autoKrige() and projected data

] http://casoilresource.lawr.ucdavis.edu/drupal/node/438 -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue Phone: +3130 253 5773 Wed-Fri http

Re: [R-sig-Geo] convert long/lati to UTM

Hi, More in general you need to take a look at the sp-pacakge for Spatial classes in R. Once your data is in these classes you can use spTransform from the rgdal pacakge to perform the projection. spTransform uses the proj4 library. You need to find the so called proj4 strings describing

Re: [R-sig-Geo] can spplot add features ?

ChinaBorder to statSpat with the 'panel' argument of spplot. Any idea about how to proceed ? Patrick ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Drs. Paul Hiemstra Department of Physical

  1   2   >