如何读取栅格并将参考数据合并到.aux.xml，.tfw等中？

Question

我有一个tif我试图用指定的投影/基准/等读入。这些是从ArcMap导出的tif，其中包含包含投影信息的.tif.xml和.tfw文件。在R中是否有一种方法可以使用.tif

引入指定的coord.ref

阅读TIF

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

输出

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)

我知道投影信息包含在相关文件中：.aux.xml，.tfw，.tif.xml。

我正在寻找有效地将其分配给tif的最佳方法吗？

tif和相关元数据文件的名称是由ArcMap导出设置和生成的约定。该目录也是共享的。

Answer 1

投影通常直接存储在tif中（作为geotiff）。似乎arcMAP也将它存储在.tif.xml中，但我真的不知道为什么它没有。无论如何，这是你可以尝试的东西：

1）找到你的投影。 proj4string或EPSG（http://spatialreference.org）

2）将其分配给您的栅格：proj4string(r) <- CRS("+init=EPSG:4326")

如果您的分配投影是数据的实际投影，则可以。