主要内容如下: R-ggplot2 可视化绘制 R-rasterVis 可视化绘制 Arcgis 可视化结果展示 R-ggplot2 可视化绘制 由于对ggplot2的绘图体系还不是很了解,所以这一步花费很长时间 R-rasterVis 可视化绘制 在查阅相关资料时发现,绘制Raster数据时,也可采用R第三方拓展包 rasterVis 包进行快速绘制,绘制代码如下: library(raster) library (rasterVis) nlcd <- brick('E:\\Data_ning\\RStudio projects\\guanzhoulanduse\\gzlu.tif') names(nlcd) levels(nlcd2011)[[1]] test$landcover <- landclass levels(nlcd2011) <- test png(file = "F:\\DataCharm\\R-<em>rasterVis</em> 具体的其他图层属性设置,大家可以直接查看官网:<em>rasterVis</em>官网 Arcgis 可视化结果展示 前面介绍了两种代码可视化的绘制教程,其实最开始我是使用Arcgis进行可视化展示的,不得不说,Arcgis
1.读取数据 library(ncdf4) library(rgdal) library(gdalUtils) library(raster) library(rasterVis) library(sf image.png rasterVis::levelplot(hdf_raster, margin = NA, par.settings = RdBuTheme) + layer(sp.polygons # rotate the x and y axis s2 = t(flip(hdf_raster, direction='y' )) # plot rasterVis::levelplot(s2, margin ;rm(rasterNoProj) # crs(rasterDF)="+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 " # plot rasterVis
先加载所需R包及地图文件 library(ncdf4) library(rgdal) library(gdalUtils) library(raster) library(rasterVis) library names : X2011001 values : 0, 255 (min, max) 3.绘图 hdf_raster是我们提取到的25类 landcover,接下来就是绘图部分 rasterVis mask(hdf_raster, CHN) # take >5mins # plot mapTheme = rasterTheme(region=rev(brewer.pal(8,"Greens"))) rasterVis
0.2) image.png image.png library(ncdf4) library(rgdal) library(gdalUtils) library(raster) library(rasterVis
getwd() setwd("C:/Users/zxs/Desktop") getwd() 加载R包 package.list=c("geoviz","tidyverse","sf","terra","rasterVis
Spatial.html https://mgimond.github.io/Spatial/coordinate-systems-in-r.html https://oscarperpinan.github.io/rastervis