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 _______________________________________________ R-sig-Geo mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-sig-geo
