[R-sig-Geo] Defining CRS when loading a raster doesn't seems to work anymore

[Verbesselt, Jan]

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 at 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 at r-project.org
>https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>