raster() and projection
Hi Agus,
How do you know the tif file has correct information. It likely is not
something that gdal understands
Can you show:
raster("/media/Transcend/MONTSENY2008/MONTSENY_GEODATA/mde10mny.tif")
GDALinfo("/media/Transcend/MONTSENY2008/MONTSENY_GEODATA/mde10mny.tif")
and a command line gdalinfo to the same file?
Thanks, Robert
On Wed, Dec 2, 2009 at 1:38 PM, Agustin Lobo <alobolistas at gmail.com> wrote:
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 R-sig-Geo at stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo