![935131028a7cf3dcbe2b227931a77de9.png](https://i-blog.csdnimg.cn/blog_migrate/43fb0cfbfcd0c6672454ddc4067e7985.jpeg)
文章参考:https://www.tylermw.com/a-step-by-step-guide-to-making-3d-maps-with-satellite-imagery-in-r/
数据准备:建议下载30米以下dem数据,以及30米以下的遥感影像,两者具有相重叠的区域(推荐利用GEE平台下载更方便)
1.相关包准备
# install.packages(c("rayshader", "raster", "sp", "magick"))
library(rayshader)
library(sp)
library(raster)
library(scales)
2.高程数据输入(替换数据路径即可)
zion_elevation = raster::raster("C:/Users/Administrator/Google 云端硬盘/caohai_dem/caohai_dem.tif")
print(zion_elevation)
height_shade(raster_to_matrix(zion_elevation)) %>%
plot_map()
3.卫星数据输入(选择合