r语言 gstat_R语言空间数据分析（五）：栅格数据处理

作者：黄天元，复旦大学博士在读，热爱数据科学与开源工具（R），致力于利用数据科学迅速积累行业经验优势和科学知识发现，涉猎内容包括但不限于信息计量、机器学习、数据可视化、应用统计建模、知识图谱等，著有《R语言高效数据处理指南》（《R语言数据高效处理指南》(黄天元)【摘要 书评 试读】- 京东图书）。知乎专栏：R语言数据挖掘。邮箱：huang.tian-yuan@qq.com.欢迎合作交流。

HopeR：R语言空间数据分析（零）：总目录​zhuanlan.zhihu.com

本帖子会简单介绍如何读入并处理栅格数据。首先，我们会用到一个矢量数据，数据来自：https://gadm.org/download_country_v3.html，用到的是澳洲的地图。读取方法如下：

# 获得数据的方法之一
# wget --no-check-certificate https://biogeo.ucdavis.edu/data/gadm3.6/gpkg/gadm36_AUS_gpkg.zip

library(pacman)
p_load(sf,raster,tidyverse)

# 查看有哪些图层
st_layers(
  "data/gadm36_AUS.gpkg"
)

# 读取特定图层
Ausoutline <- st_read("data/gadm36_AUS.gpkg", 
                      layer='gadm36_AUS_0')

可以对这个数据集进行观察：

# 观察
print(Ausoutline)

# 查看proj4坐标系
st_crs(Ausoutline)$proj4string

plot(Ausoutline)

然后，让我们读入栅格数据，这是一个全球平均气温数据，来源于：Historical climate data（tavg 5m）。

# 读入栅格数据
jan<-raster( "data/wc2.1_5m_tavg/wc2.1_5m_tavg_01.tif")
# have a look at the raster layer jan
jan
plot(jan)