Calculer l'angle formé par deux gradients

[invitec4d2e6d6]

Bonjour,

Je travaille actuellement avec un modèle du champ magnétique terrestre. Le modèle me produit en sortie deux rasters (= matrices) contenant les valeurs d'inclinaison et d'intensité du champ magnétique terrestre. Ce sont donc des données discrètes.

Comme vous pouvez le voir sur l'image suivante inclinaison et intensité forment chacun un gradient. Je souhaiterais savoir comment, à partir des deux matrices, calculer l'angle que forment les gradients d'inclinaison et d'intensité. En autre terme pour chaque point de la carte, je cherche l'angle formé par l'intersection des isolignes d'intensité et d'inclinaison. Est t'il possible d'utiliser les valeurs des matrices directement au lieu de chercher a extraire les isolignes séparément pour calculer les angles?



Je tiens à préciser que je ne peux pas modifier le code source du modèle magnétique, je dois donc travailler avec les sorties qu’il produit.

Merci Beaucoup
-----

[invite9dc7b526]

Si tu utilises pour calculer le gradient en un pixel les 8 pixels voisins ça va être peu précis. Ce que je ferais c'est ajuster un plan en utilisant mettons les 20 pixels les plus proches. A partir de l'équation du plan tu peux calculer le gradient.

[invitec4d2e6d6]

Oui je vois ce que tu veux dire, j'ai essayé et en effet il y a pas mal de bruit dans les résultats. Je ne sais pas encore si cela sera vraiment problématique ou non.

Je me demandais aussi quelles pondérations donnerais tu aux 20 pixels pour construire ton filtre?