Почему writeRaster меняет мою систему координат координат?

Question

Я редактирую растр, используя пакеты raster и rdgal в r.

Однако после внесения изменений и сохранения растра с помощью функции writeRaster я заметил, что система координат системы координат изменилась с OSGB 1936 на Airy 1830 .

Почему это происходит и как мне помешать writeRaster сделать это?

Answer 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 года).