R语言,rayshader包，render_points函数，打点，画散点。

三维散点图
添加阴影
虚线，紧贴地面,高质量渲染
# Thu Oct 15 16:41:18 2020 -

# 字符编码：UTF-8
# R 版本：R x64 4.0.2 for window 10
# cgh163email@163.com
# 个人笔记不负责任
# —— 拎了个梨🍐
.rs.restartR()
rm(list=ls());gc()
render_points函数，打点，画散点。
require(rayshader)
# 使用范围对象定义的参考系统中的纬度/经度或坐标将3D数据点添加到当前场景。 如果未提供高度，则将这些点升高到高度图上方一定的偏移量。 如果这些点偏离边缘，则将使用高度图上最近的高度。
#Starting at Moss Landing in Monterey Bay, we are going to simulate a flight of a bird going
#out to sea and diving for food.

#First, create simulated lat/long data
set.seed(2009)
moss_landing_coord = c(36.806807, -121.793332)
x_vel_out = -0.001 + rnorm(1000)[1:300]/1000
y_vel_out = rnorm(1000)[1:300]/200
z_out = c(seq(0,2000,length.out = 180), seq(2000,0,length.out=10),
          seq(0,2000,length.out = 100), seq(2000,0,length.out=10))

bird_track_lat = list()
bird_track_long = list()
bird_track_lat[[1]] = moss_landing_coord[1]
bird_track_long[[1]] = moss_landing_coord[2]
for(i in 2:300) {
  bird_track_lat[[i]] = bird_track_lat[[i-1]] + y_vel_out[i]
  bird_track_long[[i]] = bird_track_long[[i-1]] + x_vel_out[i]
}


#Render the 3D map
montereybay %>%
  sphere_shade() %>%
  plot_3d(montereybay,zscale=50,water=TRUE,
          shadowcolor="#40310a", watercolor="#233aa1", background = "tan",
          theta=210,  phi=22, zoom=0.20, fov=55)

#传递基础栅格的范围（存储在montereybay #dataset的属性中）以及轨道的纬度，经度和海拔高度。
render_points(extent = attr(montereybay,"extent"),
              lat = unlist(bird_track_lat), long = unlist(bird_track_long),
              altitude = z_out, zscale=50,color="white")
render_snapshot()
dev.copy(png,'三维散点图.jpg');dev.off()


#我们将海拔高度设置为零，以使轨道在水面上具有“阴影”。
render_points(extent = attr(montereybay,"extent"),
              lat = unlist(bird_track_lat), long = unlist(bird_track_long),
              altitude = 0, zscale=50, color="black")
render_camera(theta=30,phi=35,zoom=0.45,fov=70)
render_snapshot()
dev.copy(png,'添加阴影.jpg');dev.off()

#移除作业数据
render_points(clear_previous=TRUE)

# 最后，我们还可以通过离开海拔高度NULL来绘制从地面偏移的GPS坐标。＃在这里，我们绘制了一个围绕Moss Landing的值的圆。 这需要原始的高度图。
t = seq(0,2*pi,length.out=100)
circle_coords_lat = moss_landing_coord[1] + 0.3 * sin(t)
circle_coords_long = moss_landing_coord[2] + 0.3 * cos(t)
render_points(extent = attr(montereybay,"extent"), heightmap = montereybay,
              lat = unlist(circle_coords_lat), long = unlist(circle_coords_long),
              zscale=50, color="red", offset=100, size=5)
render_camera(theta = 160, phi=33, zoom=0.4, fov=55)
render_snapshot()
dev.copy(png,'虚线，紧贴地面.jpg');dev.off()

# 高质量渲染：
render_highquality(point_radius = 3, clamp_value=10)
rgl::rgl.close() #  关闭3D设备