Discussion :

[viv_odc]

Bonjour,

Je cherche à retrouver l'altitude z(x,y) pour chaque point (x,y) d'une image.
J'ai pour données deux matrices, appelées cartes de pente, representant la derivée partielle de z par rapport à x pour l'une, par rapport à y pour la seconde.
Je note z_x la derivée partielle premiere de z par rapport à x.

Je n'ai pas de conditions au bords spécifiées. Dans l'idéal, j'aimerais ne fixer arbitrairement qu'un seul point. On peut rêver ...

Merci d'avance de votre aide (très) rapide

[pseudocode]

Tu as essayé de simplement intégrer z_x et z_y (Runge-Kutta) et en faire la somme ?

[viv_odc]

bonjour et merci de vous interesser à mon cas

Effectivement, ma première idée a été de simplement integrer en X d'une part, et en Y d'autre part et de faire la somme.
Je n'ai pas été cherché loin pour l'intégration : cumtrapz de Matlab.
Malheureusement je finsi par diverger (le bruit s'accumule certainement).
D'ailleurs je me demande si c'est vraiment réglo de faire la somme ...

Pour prendre une formulation EDP, je veux résoudre grad(z) = F, avec F = (f1,f2) connue.

[FR119492]

Salut!

J'ai pour données deux matrices

Tout d'abord, ce ne sont pas des matrices mais des tableaux.

Ensuite, une question: je suppose que tu as une grille à mailles carrées. Les valeurs de z_x et z_y sont-elles données sur les noeuds de la grille ou sur les milieux des côtés?
Jean-Marc Blanc

[viv_odc]

Des tableaux et non des matrices, je plaide coupable d'abus de langage idiot.
J'ai effectivement des mailles carrées. Mes tableaux font 3600*3248 (pixels).
Je ne saurais pas répondre avec assurance à votre question mais j'aurais une toute petite tendance à dire que les valeurs sont sur les noeuds. A choisir, quel serait le mieux ?
Je reste en ligne et réactive car c'est urgent

merci !

[FR119492]

Salut!

A choisir, quel serait le mieux ?

Ce n'est pas toi qui choisis: ça dépend de la nature de tes données.
Jean-Marc Blanc

[viv_odc]

Je me doute bien que ce n'est pas à moi de choisir malheureusement ...
A quelles pistes pensez vous dans l'un ou l'autre des cas ?

[FR119492]

Salut!
Comme ta question est formulée, z_x et z_y sont les composantes selon x et y du gradient de z(x,y). Malheureusement, le fait qu'un champ vectoriel soit un gradient n'est qu'un cas particulier. En conséquence, si le rotationnel du champ vectoriel de composantes z_x et z_y n'est pas nul, ton problème n'a pas de solution. d'une manière plus terre à terre, si, partant d'un point, tu fais un pas selon x puis un pas selon y, ou un pas selon y puis un pas selon x, tu n'arrives pas à la même altitude.
Jean-Marc Blanc