Re: [R-sig-Geo] generate simulation data for a theoretical spatial model

2010-01-31 Thread Tomislav Hengl
Dear rusers.sh, Here are few simple examples of how to simulate (not-normal) distributions and point processes using geoR and spatstat: http://spatial-analyst.net/book/node/388 See also: http://leg.ufpr.br/geoR/geoRdoc/vignette/geoRintro/geoRintrose8.html#x9-120008 I guess that covariates

Re: [R-sig-Geo] How to diplasy a compass rose on a map

2010-02-03 Thread Tomislav Hengl
Package StatDA has an option to generate a nort arrow (I am sure you could easily replace this with any symbol you prefer): Northarrow(Xbottom, Ybottom, Xtop, Ytop, Xtext, Ytext, Alength, Aangle, Alwd, Tcex) Here are some examples of publication ready maps:

Re: [R-sig-Geo] Making web-repository of gridded maps: NetCDF or WKT Raster?

2010-02-08 Thread Tomislav Hengl
-Original Message- From: Pierre Racine [mailto:pierre.rac...@sbf.ulaval.ca] Sent: Friday, February 05, 2010 4:49 PM To: Tomislav Hengl; r-sig-geo Subject: RE: [R-sig-Geo] Making web-repository of gridded maps: NetCDF or WKT Raster? Hi Tom, I understand that PostGIS WKT Raster

Re: [R-sig-Geo] Public reference dataset

2010-03-03 Thread Tomislav Hengl
Hi Christoph, I have prepared few datasets of interest for testing geostatistical algorithms: http://spatial-analyst.net/book/data-sets I personally like the BCI (soil samples) dataset very much (used in the spatstat package also). It has enough points and they are sampled at both short and

Re: [R-sig-Geo] Geocode for Bogota

2010-03-12 Thread Tomislav Hengl
-Original Message- From: Dan Putler [mailto:dan.put...@sauder.ubc.ca] Sent: Thursday, March 11, 2010 6:46 PM To: hadley wickham Cc: Tomislav Hengl; r-sig-geo@stat.math.ethz.ch Subject: Re: [R-sig-Geo] Geocode for Bogota I'm afraid Hadley is correct. Section 10.12 of Google Maps

[R-sig-Geo] Updated repository of worldmaps (300MB of data at 5km resolution)

2010-04-01 Thread Tomislav Hengl
Dear R-sig-geo, FYI: I have just updated the small repository of publicly available data sets of interest for global modeling/mapping that I have launched about a year ago. This now contains 62 layers at resolution of 0.05 arcdegrees and with a complete world coverage (it use to be 65S-65N

Re: [R-sig-Geo] saga grid_gridding problem

2010-04-05 Thread Tomislav Hengl
Quoting Francisco Alonso Sarria alons...@um.es: Hi: I'm begining to work with rsaga and saga compiled on ubuntu linux and I have problems when I want to generate a grid from a point file (by the way is the example of section 6.3.1 in Hengl's book A practical guide to geoestatistical

[R-sig-Geo] Derivation of intermediate points (great circle geometry) in the geosphere package

2010-04-15 Thread Tomislav Hengl
I have started using the geosphere package kindly provided by Robert Hijmans [http://cran.r-project.org/web/packages/geosphere/]. My intention is to derive a global map of the air routes density using the open access flight routes [http://openflights.org/data.html]. I have discovered

Re: [R-sig-Geo] Updated repository of worldmaps (300MB of data at 5km resolution)

2010-04-16 Thread Tomislav Hengl
-sig-geo- boun...@stat.math.ethz.ch] On Behalf Of Tomislav Hengl Sent: Thursday, April 01, 2010 3:54 PM To: r-sig-geo@stat.math.ethz.ch Cc: caspar.hallm...@sovon.nl; bob.macmil...@wur.nl; 'Batjes, Niels' Subject: [R-sig-Geo] Updated repository of worldmaps (300MB of data at 5km resolution

Re: [R-sig-Geo] polygon/raster overlay error

2010-04-19 Thread Tomislav Hengl
Evann Smith wrote: Hello, I am attempting to calculate the mean value in a raster (3432 x 8640) for each of 689 polygons (code at end of email). To that end, I have a few questions: First, the for loop runs fine for a while (2 hours, maybe) before terminating itself and returning this error:

Re: [R-sig-Geo] converting lat/long to North American lambert x/y

2010-04-21 Thread Tomislav Hengl
You should use the proj4 functionality implemented in the rgdal package. This is what I get: library(rgdal) Geospatial Data Abstraction Library extensions to R successfully loaded Loaded GDAL runtime: GDAL 1.7.0dev, released 2008/11/26 Path to GDAL shared files:

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

2010-04-24 Thread Tomislav Hengl
I think you want to use the kernel density funcionality of spatstat [http://spatstat.org]; this will allow you to atach population at polygon centers as weigths. Here are some examples: library(spatstat) ... wosCR.ppp - ppp(wos...@coords[,1], wos...@coords[,2], marks=wosmap$CR,

Re: [R-sig-Geo] Summary of multiple rasters

2010-05-03 Thread Tomislav Hengl
Dear Johannes, Looks like you are in the field of error propagation [http://spatial-accuracy.org/workshopSUP]. Both mean and median are valid statistical parameters to describe the central tendencies. What is often more interesting is the propagated uncertainty in the final estimates. I

Re: [R-sig-Geo] projection

2010-05-20 Thread Tomislav Hengl
It is a pity that people that have provided this data did not attach any metadata (see for example the simple metadata format I used to prepare worldmaps e.g. [http://spatial-analyst.net/worldmaps/globedem.rdc]). My favorite places to look for correct proj4 strings (i.e. projection

[R-sig-Geo] GEOSTAT 2010 summer school Live broadcast

2010-06-24 Thread Tomislav Hengl
Dear R-sig-geo, Next week, we (Roger Bivand, Edzer Pebesma, Gerard Heuvelink, Olaf Conrad, Markus Metz, Victor Olaya and I) will run the GEOSTAT 2010 summer school in Plasencia, Spain. We have decided to open a live broadcast channel (sound and screencasting + separate web-cams), and open

[R-sig-Geo] Revolution R Enterprise - solution to handling large (spatial) data?

2010-09-10 Thread Tomislav Hengl
(Maybe this is not really suited for the R-sig-geo, but since we often experience problems with loading and visualizing large data) Has anyone yet used the Revolution R Entrerprise Version of R? Apparently it fixes/reduces the problem of large data and has an excellent GUI.

Re: [R-sig-Geo] Step characteristics on raster maps

2010-10-05 Thread Tomislav Hengl
Hi Basile, If you work with large data, then you should definitively consider using SAGA GIS. Here are some examples: -- # download gridded data: download.file(http://geomorphometry.org/sites/default/files/volcano_maungawhau.zip;, destfile=paste(getwd(),

Re: [R-sig-Geo] KED trend coefficients in Gstat with R

2010-10-26 Thread Tomislav Hengl
Op 26-10-2010 9:17, Edzer Pebesma schreef: On 10/26/2010 09:02 AM, eric.ja...@ujf-grenoble.fr wrote: Hello, I would like to know if there is a way to get the trend coefficents when doing kriging with external drift. I would like to compare with the coefficients that I can have when doing

<    1   2