光栅:下载worldclim 5分钟高度并计算地球表面面积
我想使用getData()函数下载worldclim 5分钟高度并计算表面积。对于这个品牌:光栅:下载worldclim 5分钟高度并计算地球表面面积,r,r-raster,R,R Raster,我想使用getData()函数下载worldclim 5分钟高度并计算表面积。对于这个品牌: 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 connecti
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
有什么想法吗 您可以这样做:
library(raster)
r <- getData('worldclim', var='alt', res=5)
a <- area(r)
b <- mask(a, r)
cellStats(b, sum)
#[1] 137881247
库(光栅)
R
[1]510072000
library(raster)
r <- getData('worldclim', var='alt', res=5)
a <- area(r)
b <- mask(a, r)
cellStats(b, sum)
#[1] 137881247