我正在尝试使用 R 中的 rgl 包为交互式 3d 绘图添加固定标题,但到目前为止我还无法做到。我还想在主标题下有一个主标题和一个副标题。
这是一个示例代码,
library(rgl)
data<-read.table(text=" X Y Z
1 147.0883 -18.69122 -13.90000
2 147.0894 -18.69455 -10.97250
3 147.0883 -18.69122 -17.45000
4 147.0883 -18.69122 -15.44000
5 147.0883 -18.69122 -13.45000
6 147.0909 -18.69922 -12.25000
7 147.0883 -18.69122 -17.30000
8 147.0883 -18.69122 -16.40000
9 147.0883 -18.69122 -14.30000
10 147.0883 -18.69122 -18.50000
11 147.0883 -18.69122 -15.67606
12 147.0883 -18.69122 -17.25780
13 147.0883 -18.69122 -3.64000
14 147.1164 -18.68133 -22.13000
15 147.0883 -18.69122 -18.54778
16 147.0883 -18.69122 -15.50000
17 147.1185 -18.68691 -14.55500
18 147.0883 -18.69122 -18.12500
19 147.0901 -18.69670 -14.39767",header=T)
data
a<-as.matrix(data)
# call the plug-in bandwidth estimator
H.pi<-Hpi(a,binned=TRUE)*3 ## a is a matrix of x,y,z points
# calculate the kernel densities
fhat<-kde(a,H=H.pi)
# produce a 3D plot
open3d()
title<-"Individual 1 (male)" # main title
title2<-"Site A" # subtitle
plot(fhat,cont=c(50,95),colors=c("green","red"),drawpoints=TRUE,
xlab="",ylab="", zlab="",size=1.5,ptcol="black",
axes=F,box=T,aspect=c(1,1,0.5),zlim=c(-40,0))
axes3d(c("x--","y--","z-+"),cex=0.8)
mtext3d("Longitude",edge="x--",line=3,las=2,cex=1.1)
mtext3d("Latitude",edge="y--",line=3,las=2,cex=1.1)
mtext3d("Depth (m)",edge="z-+",line=3,las=2,cex=1.1)
title3d(main=title,sub=title2,line=4)
我还没有找到固定标题(和副标题)的好方法,因此当我将情节导出为 3D 电影时,主标题不会随情节旋转......
movie3d(spin3d(axis=c(0,0,1), rpm=4),dir="~/Dropbox/R/Sample plots", duration=15, fps=10, movie="Plot 1")
如果您有一种很好/简单的方法将主要标题放在固定位置(并可能在下面添加次要标题),那就太好了。提前致谢,