Hi Jan,

Thanks.
You can also use
proj4string(r) <- "+proj=longlat +ellps=WGS84"

but the problem (the bug ?) is that it does not work anymore defining
directly the crs into the raster function.

Arnaud


2013/11/15 Verbesselt, Jan <jan.verbesselt@wur.nl>

> Hi Arnaud,
>
> Something like:
>
> r <- raster()
> projection(r)
> projection(r) <- "+proj=lcc +lat_1=48 +lat_2=33 +lon_0=-100 +ellps=WGS84"
> projection()
>
>
> Might help,
>
> Best,
> Jan
>
>
> On 15/11/13 17:52, "Arnaud Mosnier" <a.mosnier@gmail.com> wrote:
>
> >Hi all,
> >
> >I re-run a script written some times ago containing a line like the
> >following
> >
> >test <- raster("myraster.tif", crs = "+proj=longlat +ellps=WGS84")
> >
> >The file is loaded correctly, but the crs was not attributed.
> >
> >>proj4string(test)
> >
> >gives me
> >
> >[1] "NA"
> >
> >
> >When looking at the help file, nothing seems to have changed, so what the
> >problem ?
> >Any suggestion ?
> >
> >Thanks
> >
> >Arnaud
> >
> >       [[alternative HTML version deleted]]
> >
> >_______________________________________________
> >R-sig-Geo mailing list
> >R-sig-Geo@r-project.org
> >https://stat.ethz.ch/mailman/listinfo/r-sig-geo
> >
>
>
>

	[[alternative HTML version deleted]]

