R语言：将栅格数据坐标系从27700转换为WGS84经纬度

使用以下代码将坐标系从27700转换为WGS84经纬度坐标系：

library(raster)
library(rgdal)

# 读取Raster数据
r <- raster('your_raster_file.tif')

# 定义源坐标系为27700
proj4string(r) <- CRS('+init=epsg:27700')

# 定义目标坐标系为WGS84经纬度坐标系
target_crs <- CRS('+init=epsg:4326')

# 投影转换
r_wgs84 <- projectRaster(r, crs = target_crs)

# 保存结果
writeRaster(r_wgs84, 'your_output_file.tif', format = 'GTiff')

在这个代码中，我们首先使用 raster 和 rgdal 库读取和处理数据。然后，我们定义源坐标系为27700，并将其赋值给Raster对象的 proj4string 属性。接下来，我们定义目标坐标系为WGS84经纬度坐标系，并使用 projectRaster 函数将数据从源坐标系转换为目标坐标系。最后，我们使用 writeRaster 函数保存结果。