raster() and projection

Wed, Dec 2, 2009 1:38 PM

Hi!

Why is it that when we "read" a raster i.e.
dem <- raster("/media/Transcend/MONTSENY2008/MONTSENY_GEODATA/mde10mny.tif")

(where the tif file is a geotif file with correct projection informatio)
we actually do not get the projection information?

 > projection(dem)
[1] "NA"

By now I solve it thanks to a SPDF:
 > delme <- 
readOGR(dsn="/media/Transcend/MONTSENY2008/MONTSENY_GEODATA/Limitparcplaespecial", 
layer="AmbitPlaEspecial")
 > projection(dem) <- projection(delme)

Am I doing something wrong?

Thanks

Agus

raster() and projection

Thread (8 messages)