You are right.
I was unjustifiably assuming that the projection info of the
raster was right just because I could display it in QGIS at
its right position, but...
$ gdalinfo -nomd mde10mny.tif
Driver: GTiff/GeoTIFF
Files: mde10mny.tif
mde10mny.aux
mde10mny.rrd
mde10mny.tif.aux.xml
Size is 2788, 2333
Coordinate System is:
LOCAL_CS[" Geocoding information not available Projection Name = Unknown Units =
other GeoTIFF Units = other",
UNIT["unknown",1]]
Origin = (435926.940000000002328,4636392.540000000037253)
Pixel Size = (10.000000000000000,-10.000000000000000)
Corner Coordinates:
Upper Left ( 435926.940, 4636392.540)
Lower Left ( 435926.940, 4613062.540)
Upper Right ( 463806.940, 4636392.540)
Lower Right ( 463806.940, 4613062.540)
Center ( 449866.940, 4624727.540)
Band 1 Block=64x64 Type=Float32, ColorInterp=Gray
Overviews: 697x584, 349x292, 175x146, 88x73, 44x37
Thanks
Agus
Robert J. Hijmans wrote:
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 <[email protected]> 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
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
_______________________________________________
R-sig-Geo mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
