xyzポイント形式で3Dサーフェス(つまり、地震断層面)を表すデータがあります。これらのサーフェスの3D表現を作成したいと思います。rglとakimaを使用してある程度の成功を収めましたが、それ自体で折り返されたり、同じx、yポイントに複数のz値を持つ可能性のあるジオメトリを実際に処理することはできません。または、ジオメトリ(qhullのconvhulln関数)を使用して、rglで見栄えのする凸包を作成できますが、これらは実際にはオブジェクトが開いている閉じたサーフェスです(ポイントセットを完全に囲まないでください)。これらのサーフェスを作成して、できればrglでレンダリングする方法はありますか?
編集
明確にするために、ポイントはサーフェスを定義するポイントクラウド内にあります。それらは、表面全体でさまざまなカバレッジ密度を持っています。ただし、主な問題は、サーフェスが片側で閉じられていないことです。より複雑なジオメトリで閉じられていないメッシュ/サーフェスを生成する方法がわかりません。
例として...
require(rgl)
require(akima)
faultdata<-cbind(c(1,1,1,2,2,2),c(1,1,1,2,2,2),c(10,20,-10,10,20,-10))
x <- faultdata[,1]
y <- faultdata[,2]
z <- faultdata[,3]
s <- interp(x,z,y,duplicate="strip")
surface3d(s$x,s$y,s$z,col=a,add=T)
これは一般的に私が欲しいものを作成します。ただし、より複雑な飛行機の場合、これは必ずしも機能するとは限りません。たとえば、データは次のとおりです。
faultdata<-cbind(c(2,2,2,2,2,2),c(1,1,1,2,2,2),c(10,20,-10,10,20,-10))
ポイントはすべて垂直方向に同一平面上にあるため、このアプローチを使用することはできません。私も同じ問題のためにconvhullnを使用できません、そして一般的に私は閉じた船体を望んでいません、私は表面が欲しいです。alphashape3dを調べたところ、有望に見えますが、この問題にどのように使用すればよいかわかりません。