Your intuition was right. The problem is not present any more with:
> require(rgdal)
Loading required package: rgdal
Geospatial Data Abstraction Library extensions to R successfully loaded
Loaded GDAL runtime: GDAL 1.9.2, released 2012/10/08
Path to GDAL shared files: /usr/share/gdal/1.9
Loaded PR
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 different
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
repro
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 different results for the same file and objects.
Agus
O
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 d
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 wrote:
> Please provide a c
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:
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
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
>
