Question

我想为以下ASCII raster file确定正确的坐标参考系统：

class       : RasterLayer 
dimensions  : 2160, 4320, 9331200  (nrow, ncol, ncell)
resolution  : 0.0833333, 0.0833333  (x, y)
extent      : -180, 179.9999, -90, 89.99993  (xmin, xmax, ymin, ymax)
coord. ref. : NA 
data source : C:/popc_0AD.asc 
names       : popc_0AD

我尝试通过将CRS设置为某些常见格式并对其进行绘制来猜测正确的投影，如相关文章中所建议的那样。但是我仍然不确定正确的设置。就我而言，raster和相关软件包并不包含任何能够估算丢失的CRS信息的功能。您是否知道此栅格文件的CRS可能是什么或如何找到？

Answer 1

范围表明未投影坐标。这似乎是地球的度数。然后，您可能要使用EPSG 4326，它也是crs="+proj=longlat +datum=WGS84 +no_defs"：

library(raster)

r <- raster("0AD_lu/cropland0AD.asc")
projection(r) <- "+proj=longlat +datum=WGS84 +no_defs"

但是，最好使用使用坐标参考系统正确构建的数据集。永远不建议猜测它...但是我知道拥有干净的元数据并不总是可能的...

Answer 2

您有

r <- raster(nrow=2160, ncol=4320, xmn=-180, xmx=179.9999, ymn=-90, ymx=89.99993, crs=NA)

塞巴斯蒂安·罗切特（SébastienRochette）已经指出，这肯定是长期存在的问题，您可以设置CRS来解决此问题

crs(r) <- "+proj=longlat +datum=WGS84"

在我看来，这种程度有点令人怀疑。看起来它应该是全局栅格，但是精度有所损失。如果是这样，您可以这样纠正：

extent(r) <- c(-180, 180, -90, 90)

获取

r
#class      : RasterLayer 
#dimensions : 2160, 4320, 9331200  (nrow, ncol, ncell)
#resolution : 0.08333333, 0.08333333  (x, y)
#extent     : -180, 180, -90, 90  (xmin, xmax, ymin, ymax)
#crs        : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0

识别栅格文件中的CRS

2 个答案: