10

R の rgl パッケージを使用して、インタラクティブな 3D プロットに固定タイトルを追加しようとしていますが、これまでのところ、それを行うことができませんでした。また、メイン タイトルを 1 つと、メイン タイトルの下にサブタイトルを付けたいと考えています。

これがサンプルコードです。

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")

メイン タイトルを固定された場所に配置する (そして場合によってはサブ タイトルを下に追加する) 良い/簡単な方法があれば、それは素晴らしいことです。前もって感謝します、

4

4 に答える 4

14

昨年の更新により、タイトルを RGL プロットに入れることができるようになりました。

bgplot3dアイデアは、あたかも通常のプロットであるかのように、RGL ウィンドウの背景に出力する新しい関数を使用することです。

例えば:

# A minimal plot3d example
library(rgl)
a = matrix(runif(3000), ncol=3)
plot3d(a, col = rainbow(100)[cut(a[,1], breaks = seq(0,1,le=101), include.lowest = T)])

bgplot3d({
  plot.new()
  title(main = 'This is the main title', line = 3)
  mtext(side = 1, 'This is a subtitle', line = 4)
  # use here any other way you fancy to write your title
})

これにより、次のことが得られます。

                                                plot3d と bgplot3d を使用した RGL タイトル

ただし、 のヘルプ ページにメモとして記載されている、非常に重要な警告が 1 つありますbgplot3d

バックグラウンド プロットはビットマップとして描画されるため、適切にサイズ変更されません。最初にウィンドウのサイズを調整してから、そのサイズで背景を描画することをお勧めします。

于 2015-10-29T21:46:33.367 に答える
0

テストに関する提案と同様に、静止したままの凡例を追加して、それをタイトルとして使用できます。以下の凡例を rgl 3d プロットに追加するの例を参照してください。

legend3d("topright", legend = paste('Type', c('A', 'B', 'C')), pch = 16, col = rainbow(3), cex=1, inset=c(0.02))
于 2015-06-10T20:00:33.133 に答える