作者:黄天元,复旦大学博士在读,热爱数据科学与开源工具(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')
![bdee60b3624f86ea6d07e9ab5a7fb7ec.png](https://i-blog.csdnimg.cn/blog_migrate/1b50d932cdaf878aa51f30df3ad35a71.jpeg)
可以对这个数据集进行观察:
# 观察
print(Ausoutline)
# 查看proj4坐标系
st_crs(Ausoutline)$proj4string
![b2ad8495f698c2b3960e70230a656805.png](https://i-blog.csdnimg.cn/blog_migrate/6343bb4b8eb6fcddea6c77083f5d940d.png)
plot(Ausoutline)
![6f94253ee13f409f6e8e270e2d1a3783.png](https://i-blog.csdnimg.cn/blog_migrate/54f8894ddc6a60d48a63a5783de4dca8.png)
然后,让我们读入栅格数据,这是一个全球平均气温数据,来源于: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)