Определить CRS в растровом файле

Question

Я хотел бы определить правильную систему координат для следующего ASCII растрового файла :

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)

Себастьян Рошет уже указал, что это, безусловно, lon / lat и что вы можете настроить 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