# \donttest{
data("URI1_tooth")
Enamel<-URI1_tooth$mesh1
Dentin<-URI1_tooth$mesh2
Pulp<-URI1_tooth$mesh3
outline<-URI1_tooth$outline
set<-URI1_tooth$set
#Example on the root
AlignMeshesR<-ToothAlignment(mesh1=Enamel,mesh2=Dentin,mesh3=Pulp,set,outline,analyse = "r")
require(rgl)
open3d()
shade3d(AlignMeshesR$almesh1$mesh,col="white",alpha=0.5)
shade3d(AlignMeshesR$almesh2$mesh,col="pink",alpha=0.5,add=TRUE)
shade3d(AlignMeshesR$almesh3$mesh,col="orange",alpha=0.5,add=TRUE)
spheres3d(AlignMeshesR$alset,radius=0.25)
spheres3d(AlignMeshesR$outline,radius=0.1,col="blue")
lines3d(AlignMeshesR$outline)
text3d(rbind(AlignMeshesR$outline[AlignMeshesR$margins_sel,],AlignMeshesR$alset[4,]),
texts=1:5,cex=4)
spheres3d(AlignMeshesR$alset[2,],radius=0.3,col="red")
arrow3d(colMeans(AlignMeshesR$aloutline),AlignMeshesR$alset[4,],lwd=3,col="green",
type="lines",s=1/10)
lines3d(rbind(AlignMeshesR$alset[2,],AlignMeshesR$alset[1,]),lwd=3,col="green")
axes3d()
# }
#Example on the crown
# \donttest{
AlignMeshesC<-ToothAlignment(mesh1=Enamel,mesh2=Dentin,mesh3=Pulp,set,outline,
analyse = "c")
require(rgl)
open3d()
shade3d(AlignMeshesC$almesh1$mesh,col="white",alpha=0.5)
shade3d(AlignMeshesC$almesh2$mesh,col="pink",alpha=0.5,add=TRUE)
shade3d(AlignMeshesC$almesh3$mesh,col="orange",alpha=0.5,add=TRUE)
spheres3d(AlignMeshesC$alset,radius=0.25)
spheres3d(AlignMeshesC$outline,radius=0.1,col="blue")
lines3d(AlignMeshesC$outline)
text3d(rbind(AlignMeshesC$outline[AlignMeshesC$margins_sel,],AlignMeshesC$alset[4,]),
texts=1:5,cex=4)
spheres3d(AlignMeshesC$alset[2,],radius=0.3,col="red")
arrow3d(colMeans(AlignMeshesC$aloutline),AlignMeshesC$alset[4,],lwd=3,col="green",
type="lines",s=1/10)
lines3d(rbind(AlignMeshesC$alset[2,],AlignMeshesC$alset[1,]),lwd=3,col="green")
lines3d(rbind(AlignMeshesC$alset[2,],AlignMeshesC$alset[1,]),lwd=3,col="green")
axes3d()
# }
Run the code above in your browser using DataLab