растр: скачать worldclim 5 минут высоты и рассчитать площадь поверхности земли

Question

Мне нравится загружать worldclim 5 минут над уровнем моря с помощью функции getData () и вычислять площадь поверхности. Для этого сделайте:

library(raster)

#Download worldclim 5 minutes altitude for whole world 

elevation_world <- getData('worldclim', var='alt', res=5, lon=180, lat=85)

Error in file(con, "r") : cannot open the connection
In addition: Warning messages:
1: In utils::unzip(zipfile, exdir = dirname(zipfile)) :
  error 1 in extracting from zip file
2: In file(con, "r") :
  cannot open file 'C:/Users/fores/Documents/wc5/alt.hdr': No such file or directory

elevation_world <- getData('worldclim', var='alt', res=5, lon=-180, lat=-85)

Error in file(con, "r") : cannot open the connection
In addition: Warning messages:
1: In utils::unzip(zipfile, exdir = dirname(zipfile)) :
  error 1 in extracting from zip file
2: In file(con, "r") :
  cannot open file 'C:/Users/fores/Documents/wc5/alt.hdr': No such file or directory

Не работает с моим лат / длинным выбором, почему?

Во-вторых, мне тоже нравится, когда объект elevation_world рассчитывает площадь поверхности земли:

grids <- as(elevation_world, "SpatialGridDataFrame")
proj4string(grids) <- CRS(wgs.84)
areaSpatialGrid(grids)

и я ожидал, что на выходе в км2 будет что-то вроде:

[1]510072000

Есть идеи?

Answer 1

Вы можете сделать это:

library(raster)
r <- getData('worldclim', var='alt', res=5)
a <- area(r)
b <- mask(a, r)
cellStats(b, sum)
#[1] 137881247

Но учтите, что это исключает Антарктику.Вместо этого вы можете сделать

Что касается вашего первого вопроса: эти аргументы не применяются.