为什么writeRaster会更改我的坐标参考系统?

时间:2019-03-28 09:37:43

标签: r raster r-raster

我正在使用r中的rasterrdgal软件包编辑栅格。

但是,在使用writeRaster函数进行编辑并保存了栅格之后,我注意到它已将坐标参考系从OSGB 1936更改为Airy 1830

为什么会发生这种情况,如何防止writeRaster这样做?

1 个答案:

答案 0 :(得分:1)

我看到了:

library(raster)
r <- raster(nrow=10, ncol=10)
values(r) <- 1:ncell(r)
crs(r) <- "+proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 +y_0=-100000 +ellps=airy +datum=OSGB36 +units=m +no_defs"

x <- writeRaster(r, "test.tif")
crs(x)
#CRS arguments:
#+proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 +y_0=-100000 +ellps=airy +towgs84=446.448,-125.157,542.06,0.1502,0.247,0.8421,-20.4894 +units=m +no_defs 

因此坐标参考系大部分保持不变(它仍然是横向墨卡托等),但是datum=OSGB36被删除了。

我不会为什么发生这种情况,但我认为这是同一回事(OSGB36基准面是由通风的1830椭球体定义的。)