1.下载NASA数据
这里不在赘述,参考如何获取NASA
数据,下面的例子根据下载的LandCover
与Rainfall
数据进行展示,如何利用R语音进行读取,然后绘图。先加载所需R包及地图文件
library(ncdf4)
library(rgdal)
library(gdalUtils)
library(raster)
library(rasterVis)
library(sf)
library(exactextractr)
library(maptools)
library(cleangeo)
library(rworldmap)
rm(list=ls())
# continental shapefile
cont = getMap()
cont = clgeo_Clean(cont)
cont = sapply(levels(cont$continent),
FUN = function(i){
poly = gUnionCascaded(subset(cont, continent == i))
poly = spChFIDs(poly, i)
SpatialPolygonsDataFrame(poly,