3D 표면 렌더링

Question

xyz 포인트 형식으로 3D 표면 (즉, 지진 단층면)을 나타내는 데이터가 있습니다. 이 표면들에 대한 3D 표현을 만들고 싶습니다. 나는 rgl과 akima를 사용하여 약간의 성공을 거두었지만 동일한 x, y 지점에서 다중 z 값을 가질 수있는 지오메트리는 실제로 처리 할 수 ​​없습니다. 또는 기하학 (qhull의 convhulln 함수)을 사용하여 rgl에 멋지게 표시되는 볼록 선체를 만들 수 있지만 실제로는 객체가 열려있는 경우 (점 집합을 완전히 둘러싸 지 않음) 닫힌 서피스입니다. 이러한 서페이스를 만들고이를 rgl로 렌더링하는 방법이 있습니까?

명확히하기 위해 EDIT는 점은 표면을 정의하는 포인트 클라우드에 있습니다. 그들은 표면 전체에 걸쳐서 커버리지의 밀도를 변화시킨다. 그러나 가장 중요한 문제는 표면이 닫히지 않은 일방적 인 것이고 더 복잡한 형상에 대해 닫히지 않은 메쉬/표면을 생성하는 방법을 모르겠습니다. 예를 들어

...

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를 보았는데 유망 해 보였지만이 문제를 해결하는 방법을 모르겠습니다.

Answer 1

포인트가 표면으로 연결되는 방법을 어떻게 결정합니까? 거리별로? 이는 한 가지 방법 일 수 있으며 alphashape3d 패키지가 유용 할 수 있습니다. 그렇지 않으면 연결 방법을 정확하게 알고 있다면 rgl 구조로 직접 시각화 할 수 있습니다.