Nous utilisons des cookies pour vous garantir la meilleure expérience sur notre site. Si vous continuez à utiliser ce dernier, nous considèrerons que vous acceptez l'utilisation des cookies. J'ai compris ! ou En savoir plus !.
Nom d'utilisateur    Mot de passe              Toujours pas inscrit ?   Mot de passe oublié ?

#1 mer. 05 juillet 2017 12:22

CélineP
Membre
Date d'inscription: 15 mars 2017
Messages: 4

Problème calculatrice raster Gdal (via modeleur Qgis)

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

Hors ligne

 

Pied de page des forums

Powered by FluxBB

Partagez  |