Re: [R-sig-Geo] SAS to R translation help wanted

2011-11-09 Thread Michael Sumner
These filters are easy to implement and are available in contributed R packages, but they are not very good. When you say satellitedo you mean Argos or GPS? Argos recently updated their proprietary algorithms to give better estimates but it it is still a rather black box.To me the issue is about

[R-sig-Geo] bbox

2011-11-09 Thread Raffaele Morelli
Hi, once a Postgis layer X is retrieved, how can I juggle with subsets of X? eg. I am working with bubble plots on a EU dataset but I want to do test it only using a 30x30Km in the nearby of Rome and I have no attribute data but the value I use for the bubbles. Should I use subset? Can I set a

Re: [R-sig-Geo] bbox

2011-11-09 Thread Tom Gottfried
assuming your geometry column is 'the_geom': while querying SELECT ... FROM ... WHERE the_geom @ ST_geomfromtext(...); with ST_geomfromtext() you can set up your custom bbox. But of course you can load everything into R and then use subset() if your data are in a data.frame or use over() with

Re: [R-sig-Geo] Calculating perimeter

2011-11-09 Thread Sam Rabin
Roger — Setting the repositories as suggested worked well, for these and one other package I'd been having trouble with. gLength worked fine once I got them working. I was only talking about single-ring polygon objects. Robert — The shapefiles were indeed in lat/long, so the fact that

Re: [R-sig-Geo] [raster] substraction between two rasterLayers: error Readbin

2011-11-09 Thread Carsten Neumann
Am 09.11.2011 15:39, schrieb Mathieu Rajerison: Hi, I have a RasterStack detm composed of a dem and a dtm. I wanted to apply the difference between dem and dtm so as to get the elevation. But I got the following message (in French): Erreur dans result[, i]- readBin(raster@file@con,

Re: [R-sig-Geo] Export raster data as *.grd file format in ArcMap

2011-11-09 Thread Robert J. Hijmans
Leni, If you have raster data in Arc and export as GRID it will create a database (folder, including a parallel folder called info) that you can treat like a file and open in R. Alternatively you can export to tif, img (erdas) and other formats (avoid ascii) that raster can read (most of them via

[R-sig-Geo] Stack, brick (or other data object) to query large amounts of raster data

2011-11-09 Thread pgalpern
Hello! Trying to decide the best way to store a large amount of data in a function I am writing using the raster package. The goal is speed. I need to flexibly and efficiently handles cases where 10 to 5000 copies of a raster are required. The raster that requires replication would

Re: [R-sig-Geo] [raster] substraction between two rasterLayers: error Readbin

2011-11-09 Thread Mathieu Rajerison
Ok, thanks for the answer. It must be a problem of disk space. I'll try what you say. Best, Mathieu 2011/11/9 Robert J. Hijmans r.hijm...@gmail.com Mathieu/Carsten, Carsten's suggestion is also good: elev-overlay(detm@layers[[1]],detm@layers[[2]], fun=function(a,b) {return(a-b)} )