Murakami's Memorandum

Home Page > My Memo               災害科学ホワイトボード   災害科学コース   理学部

(R) 3次元散布図・等高線マップ作成メモ

 3次元散布図(固定)

library(lattice)
cloud(Porosity~Depth*Inv_Vp,data=Dataset,groups=HCS)

 3次元散布図(マウスで回転)

library(rgl)
scatter3d(Dataset$Depth,Dataset$Inv_Vp,Dataset$Porosity,xlab="Depth",ylab="Inv_Vp",zlab="Porosity")

 観測点の分布を密度関数として等高線を作成

> require(MASS)
要求されたパッケージ MASS をロード中です
> x<-Dataset$Depth
> y<-Dataset$Inv_Vp
> d<-kde2d(x,y,n=50)
> image(d)
> contour(d,add=T)
> points(x,y,pch=20)

 不等間隔のデータを補完して等高線図を作成

> library(akima)
> dd<-interp(Dataset$Depth, Dataset$Inv_Vp, Dataset$Porosity, duplicate='mean')
> image(dd$x, dd$y, dd$z, xlab="Depth", ylab="Inv_Vp")
> contour(dd$x,dd$y, dd$z, add=TRUE)
> points(Dataset$Depth, Dataset$Inv_Vp, pch=20)

最終更新時間:2020年01月22日 18時34分45秒


HomePage > FrontPage