Le portail francophone de la géomatique

Bonjour,

Je me tourne vers ce forum car je rencontre un problème avec la calculatrice raster de Gdal.

Pour résumé, je cherche à créer un modèle de traitement sous Qgis avec le modeleur graphique.
J'ai 4 raster avec lesquels je souhaite calculer des indices, dont les formules mathématiques sont plus ou moins complexes.
Parmi les algorithmes disponibles dans le modeleur, il y a la calculatrice raster de Gdal et celle de SAGA. J'ai donc choisi d'utiliser la première.

Après quelques recherches, j'ai remarqué que la syntaxe de cette calculatrice différait de celle que l'on trouve dans le menu raster de Qgis.
Gdal utilise le langage numpy (python). Je ne connais pas vraiment ce langage mais j'ai trouvé la plupart des "erreurs" que je pouvais faire dans mes formules : il faut utiliser un point et non pas une virgule et le signe ** pour une puissance par exemple.

De ce fait, lorsque je réalise un calcul simple : (A-B)/(A+B)
-> aucun problème, j’obtiens le raster calculé attendu

Mais lorsque j'essaye de réaliser un calcul plus compliqué (avec des racines et des puissances) -> j'obtiens un raster no data

Voici les deux formules qui me posent problème à l'heure actuelle :

(1.8*A-3.75*B+1.95*C) / sqrt ((2*A+1)**2-6*A+5* sqrt (B)-0.5)

et (1.2*A+1.5*B+0.3*C)/ sqrt (A-B**2+1.5)


Si vous avez des pistes ou que vous voyez où peut se situer mon erreur, je suis preneuse de toute suggestion !


Bonne journée à tous