Comment lire dans un raster et incorporer des données de référence dans .aux.xml, .tfw etc?

Question

J'ai un tif que j'essaie de lire avec la projection/datum/etc assignée. Ce sont des tifs exportés à partir d'ArcMap avec les fichiers .tif.xml et .tfw contenant des informations de projection. Est-il possible en R pour la coord.ref attribué à la .tif

Lire dans TIF

r<-'example.tif' 
r <- raster(r) 
r 

Sortie

class  : RasterLayer 
dimensions : 199, 695, 138305 (nrow, ncol, ncell) 
resolution : 50000, 50000 (x, y) 
extent  : -17367529, 17382471, -4692230, 5257770 (xmin, xmax,ymin, ymax) 
coord. ref. : NA 
data source : in memory 
names  : layer 
values  : 0.268264, 5.886104 (min, max) 

Je sais l'information de projection est contenue dans le la fichiers associés: .aux.xml, .tfw, .tif.xml.

Je cherche les meilleurs moyens pour l'assigner efficacement au tif?

Les noms des tifs et des fichiers de métadonnées associés sont les ensembles de conventions générés par l'exportation ArcMap. Le répertoire est également partagé.

Answer 1

La projection est normalement stockée directement dans le tif (comme un géotiff). Il semble qu'arcMAP le stocke également dans le fichier .tif.xml mais je ne sais pas vraiment pourquoi, car il n'est pas obligé de le faire. Quoi qu'il en soit, voici quelque chose que vous pouvez essayer:

1) trouvez votre projection. Soit le proj4string ou l'EPSG (http://spatialreference.org)

2) assigner à votre raster: proj4string(r) <- CRS("+init=EPSG:4326")

Ceci est ok si la projection de votre est la projection affectation réelle dans laquelle les données sont.