Comment créer des tracés de découpe 2D à partir de données 3D dans R?

Question

Voici le code R pour ma parcelle 3D ... Je veux générer des tranches 2d de ce tracé 3D. J'utilise RGL paquet dans R.

library(rgl) 

    > ncv1 
      [,1]   [,2]   [,3] 
[1,] 0.0087173253 0.0015347787 0.0007591226 
[2,] 0.0015347787 0.0066770936 -0.0001649119 
[3,] 0.0007591226 -0.0001649119 0.0011324989 

> ncv2 
      [,1]   [,2]   [,3] 
[1,] 0.0080515223 0.0.0006745652 
[2,] 0.0.0064364687 -0.0001695962 
[3,] 0.0006745652 -0.0001695962 0.0011687981 

> ntruemean 
[1] 4.775067 5.596896 1.142223 

> open3d() 
wgl 
10 

plot3d(ellipse3d(ncv2,ntruemean),col="red",type='shade',xlim=c(-2,2),ylim=c(-2,2),zlim=c(-2,2)) 

plot3d(ellipse3d(ncv1,ntruemean),col="blue",type='wire',add=TRUE,xlim=c(-2,2),ylim=c(-2,2),zlim=c(-2,2)) 

rgl.postscript("persp3dd1.eps","eps") 

Answer 1

Vous pouvez le couper en tranches assez facilement dans tout plan que vous aimez. Vous avez juste besoin de comprendre les coefficients de l'équation de l'avion au format a x + b y + c z + d = 0, qui est un problème d'algèbre linéaire de 1ère année.

Je ne peux pas faire cela pour vous, parce que je ne sais pas quelle ellipse vous avez l'intention de dessiner, ou ce que vous entendez par «droit au centre». Vous utilisez ntruemean comme nom de variable, mais passez le comme le paramètre scale, ce qui n'a pas vraiment de sens. Je soupçonne que vous vouliez le passer comme le paramètre centre.

En supposant que ce que vous vouliez dire, ce clip dans le centre:

library(rgl) 
ncv2 <- diag(1:3) # not your values 
ntruemean <- 4:6 # also not your values 
plot3d(ellipse3d(ncv2, centre = ntruemean), 
     col="red", type='shade') 
plot3d(ellipse3d(ncv2*1.01, centre = ntruemean), col = "blue", type="wire", add = TRUE) 
normal <- 7:9 # Choose your own normal! 
clipplanes3d(normal, d = -sum(normal * ntruemean))