tidyterra绘制具体空间数据地图

欢迎关注R语言数据分析指南

❝

昨天分享了如何使用R内的tidyterra包来进行空间数据分析绘制地图，本节来展示如何通过NASA下载各类的地理空间数据来进行绘图，非常的简单。对于颜色的选择tidyterra官方有很具体案例来介绍,更多详细内容请参考官方文档。

官方文档

❝

https://dieghernan.github.io/tidyterra/articles/tidyterra.html

NASA

❝

可以通过NASA下载各类地理空间数据
https://neo.gsfc.nasa.gov/view.php?datasetId=GISS_TA_Y&date=2023-12-01

数据格式转换

library(sf)
library(tidyverse)
library(terra)
library(tidyterra)


bathy_raster <- rast('GEBCO_BATHY_2002-01-01_gs_360x180.TIFF')
# 将无效值（99999）替换为 NA
bathy_raster[bathy_raster == 99999] <- NA
#以经度 -100° 为投影中心，纬度10为投影中心
crs_string <- "+proj=ortho +lon_0=-100 +lat_0=10 +x_0=0 +y_0=0 +a=6370000 +b=6370000 +units=m +no_defs"

案例1

ggplot() +
  geom_spatraster(data = bathy_raster) +
  scale_fill_princess_c(
    palette = "ella",
    labels = scales::label_number(),
    breaks = c(50,100,150,200,250),
    guide = guide_colorbar(reverse = TRUE)) +
  theme_void() +
  theme(plot.background = element_blank(),
        legend.position = 'none') 

案例2

# 将栅格数据投影到正射投影
bathy2 <- terra::project(bathy_raster, crs_string)

ggplot() +
  geom_spatraster(data = bathy2) +
  scale_fill_princess_c(
    palette = "cold",
    labels = scales::label_number(),
    breaks = c(50,100,150,200,250),
    guide = guide_colorbar(reverse = TRUE)) +
  theme_void() +
  theme(plot.background = element_blank(),
        legend.position = 'none') 

案例3

bathy_raster2 <- rast('RenderData.tiff')  

ggplot() +
  geom_spatraster(data = bathy_raster2) +
  scale_fill_hypso_c(palette = "utah_1") +
  theme_void() +
  theme(panel.background = element_blank(),
        plot.background = element_blank(),
        legend.position = 'none') 

案例4

bathy_raster3 <-  terra::project(bathy_raster2, crs_string)

ggplot() +
  geom_spatraster(data = bathy_raster3) +
  scale_fill_hypso_c(palette = "utah_1") +
  theme_void() +
  theme(panel.background = element_blank(),
        plot.background = element_blank(),
        legend.position = 'none') 

关注下方公众号下回更新不迷路

购买介绍

❝

本节介绍到此结束，有需要学习R数据可视化的朋友欢迎到淘宝店铺：R语言数据分析指南，购买小编的R语言可视化文档(2024版)，购买将赠送2023年的绘图文档内容。目前此文档(2023+2024)已经更新上传200案例文档，每个案例都附有相应的数据和代码，并配有对应的注释文档，方便大家学习和参考。

2024更新的绘图内容将同时包含数据+代码+注释文档+文档清单，2023无目录仅有数据文件夹，小编只分享案例文档，不额外回答问题，无答疑服务，零基础不推荐买。

案例特点

所选案例图均属于个性化分析图表完全适用于论文发表，内容异常丰富两年累计发布案例图200+，2024年6月起提供html版注释文档更加直观易学。文档累计上千人次购买拥有良好的社群交流体验。R代码结构清晰易懂，为防止中文乱码提供单独的注释文档

R代码结构清晰易懂，为防止中文乱码2024年6月起提供单独html注释文档

群友精彩评论

淘宝店铺

2024年已更新案例图展示

2023年案例图展示