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,
本文介绍了如何在R中下载NASA的LandCover和Rainfall数据,处理HDF文件,将其转换为raster格式,然后进行地理空间分析和绘图,最后展示了在中国地图上应用这些数据的实例和两种不同的数据可视化方法。
订阅专栏 解锁全文
1273

被折叠的 条评论
为什么被折叠?



