Calcul de volume sur une ellipsoïde

[mathieu_t]

Bonjour,

Bon j'ai un peu honte de poser la question, mais je patauge un peu :
Comment calculer le volume d'une ellipsoïde, ou plutôt une partie de ce volume ?

Je m'explique : J'ai un volume 3D discrétisé en voxels (éléments de volume en petits cubes), et je voudrais connaître pour chaque élément de volume à la surface de l'ellipsoïde la proportion d'ellipsoïde qu'il y a dedans (cette proportion allant de 0 s'il l'ellipsoïde ne passe pas dedans, à 1 si le voxel est totalement compris dedans).

Il m'est donc naturellement venu l'idée de calculer le volume de la portion d'ellipsoïde traversant un voxel.

Comme je travaille en coordonnées cartésiennes (les voxels étant des petits cubes de 1*1*1 mm^3), mon raisonnement est dans ces coordonnées :

L'équation de l'ellispoïde est avec a, b, c les 3 demis axes, réels strictement positifs.
J'ai alors pris deux cas:
+ si z >=0,
Dans notre cas, on peut définir un voxel dans le demi-plan positif (z>=0) par sa projection sur le plan Oxy (à un voxel, il y a une unique projection) : ses coordonnées sont le bord inférieur gauche (u,v) et le bord supérieur droit (u+1,v+1), et le volume que l'on souhaite calculer correspond alors à :
où

... Et c'est là que je coince, car je n'arrive pas à intégrer cette fonction (si j'intègre par rapport à x par exemple, je vais me retrouver avec une formule compliquée pleine de racines carrées et de ln, que je n'arrive pas à intégrer)...

+ si z < 0 : même raisonnement, mais avec -c au lieu de c...

++ J'ai pensé à passer en coordonnées sphériques, mais je patauge...

Suis-je clair ? Mon raisonnement est-il bon ? Une idée ?

[jose_latino]

La méthode est, avant d'utiliser coordonées sphériques, transformer l'ellipsoïde à une sphère, avec une transformation linéaire:



Si tu connais le théorème du changement de variables pour les intégrals triples, tu n'auras pas de problème. Bonne chance! :zen:

[jose_latino]

Je crois que tu sais utiliser seulement des intégrals doubles. C'est la même idée. Pour intégrer une fonction, il faut savoir décrire l'ensemble où tu vas intégrer, dans ce cas:

Donc, tu vas intégrer , sur l'ensemble (ellipse et son intérieur)
On peut faire des petit changement de variables:


Ça veut dire qu'on va utiliser la transformation
Alors: si (cercle), alors
Par le théorème du changement de variables:

Ici, tu peux appliquer les coordonées polaires, car, est un cercle. Bon courage! :id:
P.S: Tu peux remarquer que est la moitié du volume de la sphère unitaire.

[mathieu_t]

Salut !

Tout d'abord merci pour ta réponse.

Ensuite, quelques éclaircissements, et surtout quelques questions :
D'abord je ne souhaite pas intégrer sur toute l'ellipse (ou sphère par changement de variables), mais uniquement sur une petite partie.

Ensuite, ça n'est pas que je ne sais utiliser que des intégrales doubles, c'est surtout que je ne sais pas comment m'y prendre autrement (et oui, les maths commencent à cruellement dater pour moi, et pourtant j'en ai de plus en plus besoin en ce moment) :
D'où la question :
Euh ? C'est possible d'intégrer directement avec l'équation ?? (Le souci c'est que je ne vois pas comment... ça n'est pas une fonction, mais une équation...).
Parce que si ça existe, je suis preneur...

En attendant, j'ai trouvé ma méthode "bidouille" qui consiste à exprimer z en fonction des 2 autres variables puis d'intégrer...

Bon, en partant du principe que je n'arrive à utiliser que cette bidouille pour le moment, je suis arrivé au résultat suivant (je mets ça ici histoire de voir si j'ai pas fait une c... dans mes calculs) :
Avec le changement de variable, ça donne :


A présent, passons en coordonnées polaires :
Je pose et .
Le changement de variable me donne :
.

Si je ne m'abuse,


Finalement,

Avec





J'ai bon ?

Merci au fait !

[jose_latino]

C'est pas aconseillé intégrer directement, mais je te donne la formule:
(*)

[jose_latino]

C'est en (*) où tu peux utiliser les coordonées polaires

À partir d'ici, c'est pas difficil arriver à la réponse. Bon courage!

[jose_latino]

D'où la question :
Euh ? C'est possible d'intégrer directement avec l'équation ?? (Le souci c'est que je ne vois pas comment... ça n'est pas une fonction, mais une équation...)

Avec l'équation donnée, tu as seulement la frontière de cet ensemble. Tu as que la région est décrite de la façon suivante:

Pour intégrer il faut décrire la région à intégrer de la façon indiqué. Pour faire ça, j'ai utilisé l'équation que tu as donné. Une autre façon de décrire la région est:

[mathieu_t]

Salut !

Pour ce qui est de la solution, c'est ok.

Par contre,

Avec l'équation donnée, tu as seulement la frontière de cet ensemble. Tu as que la région est décrite de la façon suivante:

Pour intégrer il faut décrire la région à intégrer de la façon indiqué. Pour faire ça, j'ai utilisé l'équation que tu as donné. Une autre façon de décrire la région est:

Ok, mais j'intègre quoi ???

[jose_latino]

Si une région est décrite de la manière:

alors, on a que son volume est:

[pizzouille]

Bonsoir,

Soit S la surface de niveau 0 de F : U R -> R².

On demande théorème de fonction implicite au point A=(a,b,c) S en le détaillant le plus précisement possible où g(a,b) = c est la fonction implicite utilisée.
Puis de montrez que la surface paramétrée : V x I -> R, (x,y,z) -> (x,y,g(y,z)) est partout régulière et de donnez l'équation de son plan tangent en (a,b).

Merci à l'avance de votre aide.