Re: [R-sig-Geo] Non Square Pixels

2012-11-13 Thread Rinke
Last week we released a new version of the ARIS ASCII Grid Reader for ArcGIS 10.0: http://www.aris.nl/asciigridreader_arcgis http://www.aris.nl/asciigridreader_arcgis The ASCII Grid Reader version 3.0 will do a direct read of ESRI ASCII GRID's in ArcCatalog and ArcMap with small modifications to

[R-sig-Geo] raster CRS changes on tis own

2012-11-13 Thread Agustin Lobo
Hi! I modify the CRS of a raster object: vinedomosaicMCA03 - brick(paste(mosdir,mosaico_envi.dat,sep=/)) proj - CRS(+init=epsg:3042) proj CRS arguments: +init=epsg:3042 +proj=utm +zone=30 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs projection(vinedomosaicMCA03) - proj

Re: [R-sig-Geo] raster CRS changes on tis own

2012-11-13 Thread Agustin Lobo
Sorry I forgot my session info: sessionInfo() R version 2.15.2 (2012-10-26) Platform: x86_64-pc-linux-gnu (64-bit) locale: [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C LC_TIME=en_US.UTF-8LC_COLLATE=en_US.UTF-8 LC_MONETARY=en_US.UTF-8LC_MESSAGES=en_US.UTF-8 [7] LC_PAPER=C

Re: [R-sig-Geo] raster CRS changes on tis own

2012-11-13 Thread Roger Bivand
Please provide a completely reproducible example, otherwise nobody can say that this isn't simply a feature of your data. I can't see why helpers should be expected to create an example for you, which may not correctly represent your problem. Roger On Tue, 13 Nov 2012, Agustin Lobo wrote:

[R-sig-Geo] [plain text] Error: cannot allocate vector of size 274.7 Mb

2012-11-13 Thread Raffaele Morelli
Hi, sorry for double posting but I saw my messages were scrubbed from html garbage and don't know if nobody replies for that, but just in case... apologize if not. I am having trouble with this error when calculating a density on a point pattern dataset with 34480 points. My script exit when

Re: [R-sig-Geo] raster CRS changes on tis own

2012-11-13 Thread Agustin Lobo
http://dl.dropbox.com/u/3180464/mosaico_envi.dat http://dl.dropbox.com/u/3180464/mosaico_envi.hdr to be read in with vinedomosaicMCA03 - brick(mosaico_envi.dat) (both files must be present in the same directory) Agus On Tue, Nov 13, 2012 at 11:01 AM, Roger Bivand roger.biv...@nhh.no wrote:

Re: [R-sig-Geo] raster CRS changes on tis own

2012-11-13 Thread Roger Bivand
On Tue, 13 Nov 2012, Agustin Lobo wrote: http://dl.dropbox.com/u/3180464/mosaico_envi.dat http://dl.dropbox.com/u/3180464/mosaico_envi.hdr I do not think dumping 40MB on people is a helpful example. The idea is that you rather go through the trouble of re-creating the problem from a small

Re: [R-sig-Geo] [plain text] Error: cannot allocate vector of size 274.7 Mb

2012-11-13 Thread Roman Luštrik
Would it be at all possible if you made your problem reproducible? Cheers, Roman On Tue, Nov 13, 2012 at 11:20 AM, Raffaele Morelli raffaele.more...@gmail.com wrote: Hi, sorry for double posting but I saw my messages were scrubbed from html garbage and don't know if nobody replies for

Re: [R-sig-Geo] [plain text] Error: cannot allocate vector of size 274.7 Mb

2012-11-13 Thread Mauricio Zambrano-Bigiarini
On 13/11/12 11:20, Raffaele Morelli wrote: Hi, sorry for double posting but I saw my messages were scrubbed from html garbage and don't know if nobody replies for that, but just in case... apologize if not. I am having trouble with this error when calculating a density on a point pattern

Re: [R-sig-Geo] [plain text] Error: cannot allocate vector of size 274.7 Mb

2012-11-13 Thread Raffaele Morelli
2012/11/13 Roman Luštrik roman.lust...@gmail.com: Would it be at all possible if you made your problem reproducible? Cheers, Roman I can upload my .RData (~45Mb) on a server and the script if you're willing to try. In the meanwhile I am doing some test with my owin taken from the maps

Re: [R-sig-Geo] Azimuthal equidistant projection and rgdal

2012-11-13 Thread Ashton Shortridge
Dear Johan (and list), you are right in the first case: I think the result is expected. Using global data with this sort of projection will look kind of awful - the equatorial azimuthal equidistant illustrations in Snyder do not try to map data for the whole planet - they employ a hemisphere.

[R-sig-Geo] Strange error plotting SpatialPolygonsDataframe

2012-11-13 Thread Ralf Schäfer
Dear list, I encountered an odd behaviour. I have a R Script for a course I teach, which runs fine in my office. Now, when I attach the notebook to a beamer (i.e. resolution decreases), the plot function does not work (i.e. R does not respond anymore) for the following code: getData(GADM,

[R-sig-Geo] raster CRS changes on tis own

2012-11-13 Thread Michael Sumner
I think Roger showed remarkable patience with this, it really is asking too much if you cannot be bothered providing a single set of runnable code. That is after all what one must do to investigate the request, and it's simply impolite not to provide it. At the very least we cannot be sure we

Re: [R-sig-Geo] raster CRS changes on tis own

2012-11-13 Thread Roger Bivand
On Tue, 13 Nov 2012, Agustin Lobo wrote: Roger, I'm not dumping anything on anybody: those are links that you download if you want to, and that I posted only after you requested an example. Also, could you post your session info? According to what you say and what I reproduce we get

Re: [R-sig-Geo] Azimuthal equidistant projection and rgdal

2012-11-13 Thread Johan Mazel
Hi I don't think that Synder's book actually use a hemisphere for the azimuthal equidistant projection examples. For example, on figures 41.B and 41.C, you can see Australia (respectively left to top-left and bottom-left). I think that in both cases, you can see the whole world. This is quite