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 !.
banniere

Le portail francophone de la géomatique


Toujours pas inscrit ? Mot de passe oublié ?
Nom d'utilisateur    Mot de passe              Toujours pas inscrit ?   Mot de passe oublié ?

Annonce

Suite à un problème technique intervenu entre le 22 et le 23 mars, nous avons du procéder dans la soirée du 25 mars, à la restauration de la base de données du 24 mars (matinée).

En clair, nous avons perdu vos contributions et inscriptions du dimanche 24 et du lundi 25 mars.
Nous vous prions de nous excuser.

#1 Wed 05 July 2017 12:22

CélineP
Juste Inscrit !
Date d'inscription: 15 Mar 2017
Messages: 8

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