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 Tue 25 July 2017 11:15

jlavaure
Participant assidu
Date d'inscription: 20 Apr 2013
Messages: 201

QGIS 2.14 + GRASS: HELP - MNE à partir BDALTI 75m et BDTOPO

Bonjour à tous,

Alors voilà, je maîtrise plutôt bien les SIG en mode vecteur, mais dès qu'il y a une histoire de raster de suis loin d'être à l'aise.
De toute façon il faut bien s'y confronter, je dois réaliser un MNE (Relief + Bâtiment + Végétation).

J'ai en ma possession la BDALTI 75m, la BDTOPO BATI et VEGETATION.

J'utilise QGIS avec les modules GRASS.

Première étape : Je ré-échantillonne la BD ALTI 75m en 5m afin d'avoir une résolution plus "logique et précise par rapport a la BDTOPO" --->  r.resamp.interp bicubique (je ne connais pas la différence avec les autres méthodes de ré-échantillonnage ?

Deuxième étape : Je rastérise mon bâti avec h comme paramètre d'altitude et 5m de résolution

Troisième étape : J'estime une hauteur par type de végétation et rastérise comme le bâti

Quatrième étape, le blocage : Comment additionner avec la calculatrice raster mes deux rasters issus de la BDTOPO ? Si j'ai une habitation de 5m au milieu d'une forêt de 10m je ne veux pas obtenir une parcelle de 15m mais bien de 10m. A contrario si j'ai une batisse de 10m au milieu des vignes (1m) il faudrait un résultat de 10m et non 11m.

Je précise que la projection est la même pour l'ensemble des mes couches.

Du coup une idée pour m'aiguiller, un pro des raster et QGIS/GRASS ?

Merci d'avance. J.

Dernière modification par jlavaure (Tue 25 July 2017 11:25)

Hors ligne

 

#2 Tue 25 July 2017 12:07

Kadiri28
Participant actif
Lieu: Toulouse
Date d'inscription: 10 Mar 2010
Messages: 121

Re: QGIS 2.14 + GRASS: HELP - MNE à partir BDALTI 75m et BDTOPO

Bonjour,
Il faut combiner tous les shp avant de rasteriser en mettant la hauteur dans la table attributaire, par ex découpez la couche forêt avec la couche bâti puis assemblez les et puis rasterisez en mettant la colonne hauteur comme valeur.

Bon courage

Hors ligne

 

#3 Tue 25 July 2017 12:17

jlavaure
Participant assidu
Date d'inscription: 20 Apr 2013
Messages: 201

Re: QGIS 2.14 + GRASS: HELP - MNE à partir BDALTI 75m et BDTOPO

@Kadiri28 :

Merci pour ce conseil, mais en soit je ne vois pas comment cela peut fonctionner si par exemple si un bout de fôret de 10m est découpé et lorsque j'assemble avec le bâti il y a la maison de 5m, alors mon MNT a un "trou" ... non ?

Hors ligne

 

#4 Tue 25 July 2017 16:15

Kadiri28
Participant actif
Lieu: Toulouse
Date d'inscription: 10 Mar 2010
Messages: 121

Re: QGIS 2.14 + GRASS: HELP - MNE à partir BDALTI 75m et BDTOPO

Je n'ai pas bien compris la question!
sinon pour avoir le MNE, il ya deux raster : un "MNS" et un"MNH", le MNH s'ajoute au MNH, si le bâtiment fait 5m il est additionné au 75m d'altitude du MNS donc 80m, la forêt autour aura 85m.

Hors ligne

 

#5 Tue 25 July 2017 16:31

jlavaure
Participant assidu
Date d'inscription: 20 Apr 2013
Messages: 201

Re: QGIS 2.14 + GRASS: HELP - MNE à partir BDALTI 75m et BDTOPO

En soit, si j'ai une maison dans de la forêt je ne souhaite pas avoir comme valeur la maison + la forêt mais la plus haute des deux. Donc mon MNE sera le MNT + les plus hautes valeurs à chaque fois entre mon raster bâti et végétation.

Je ne sais pas si c'est plus compréhensible ^^

Hors ligne

 

#6 Tue 25 July 2017 17:16

Thom25
Participant occasionnel
Date d'inscription: 16 Sep 2013
Messages: 31

Re: QGIS 2.14 + GRASS: HELP - MNE à partir BDALTI 75m et BDTOPO

Salut,

Je n'ai pas testé, mais essaye peut être avec la fonction de la calculatrice raster de grass : r.mapcalc . Tu devrais pouvoir additionner les valeurs de chaque cellule des différents raster.

Hors ligne

 

#7 Tue 25 July 2017 17:59

jlavaure
Participant assidu
Date d'inscription: 20 Apr 2013
Messages: 201

Re: QGIS 2.14 + GRASS: HELP - MNE à partir BDALTI 75m et BDTOPO

@Thom25 :

Je dois décidément mal m'exprimer ^^

L'idée n'étant pas justement de les additionner, mais de sortir un MNE le plus proches de la réalité. La finalité de se projet est de faire des études d’influences visuelles.

Concrètement :

-  j'ai ma base avec le MNT de la BDTOPO 75m que j'ai ré-échantillonner avec un pas de 5m.

- Dessus je souhaiterais y ajouter la végétation et le bâti qui sont tous les deux dans 2 rasters différents avec également un pas de 5m.

- Si j'additionne "classiquement" je vais avoir le MNT + Bâti + Végétation ... Donc dans une forêt de 15m avec un château de 20m je vais obtenir 35m
... alors qu'en réalité la hauteur qui sera perçu du sol sera 20m (la plus haute). Cas contraire, j'ai une vigne (1m) avec un mas en plein milieu (10m), la hauteur réelle qui doit être ajouté au MNT est 10m et non (11+1).

Et c'est exactement ça que je cherche à faire.

Mais comment faire en restant dans QGIS + Modules Grass ?

J.

Hors ligne

 

#8 Wed 26 July 2017 00:57

Thom25
Participant occasionnel
Date d'inscription: 16 Sep 2013
Messages: 31

Re: QGIS 2.14 + GRASS: HELP - MNE à partir BDALTI 75m et BDTOPO

D'accord, mais lorsque que tu rasterises tes données bâti et végétation, il devrait y avoir des cellules nulles dans tes rasters bâti et végétation (correspondant soit à une absence de bâtiments soit à une absence de végétation).

Je ne sais pas avec quelles données sources tu as créé les modèles numériques de bâti et de végétation, mais si les données de base sont bien faites, il ne devrait pas y avoir d'entités bâti qui intersectent des entités végétation (car je pense qu'aucun château n'est construit sur la cime des arbres d'une forêt et inversement, on ne fait généralement pas pousser d'arbre sur un château)

Si les données sont mal fichues alors il faut que tu passes par une (ou des) étape(s) de découpe avant de rasteriser. (outil "Différenciation" sur QGIS ou "v.overlay" sur grass) Ceci afin supprimer les hauteurs qui peuvent se superposer entre le bâti et la végétation

J'espère que j'ai mieux compris cette fois-ci. ^^

Hors ligne

 

#9 Wed 26 July 2017 09:10

jlavaure
Participant assidu
Date d'inscription: 20 Apr 2013
Messages: 201

Re: QGIS 2.14 + GRASS: HELP - MNE à partir BDALTI 75m et BDTOPO

@Thom25 :

... pas con ! C'est moi qui suis à côté de la plaque depuis le début.

Les données sources sont issues de la BDTOPO de l'IGN.
En regardant de plus près dans ma zone d'étude il y a bien des bâtiments au milieu de "zone arborée" par exemple. Du coup je vais différencier en amont la végétation avec le bâti. Puis rastériser mes deux couches vecteurs.

Dans la calculatrice raster il suffit de faire "a + b+ c" ?

J.

Hors ligne

 

#10 Thu 27 July 2017 11:41

Thom25
Participant occasionnel
Date d'inscription: 16 Sep 2013
Messages: 31

Re: QGIS 2.14 + GRASS: HELP - MNE à partir BDALTI 75m et BDTOPO

Pour la calculatrice raster, il me semble que la formule que tu as écrite doit fonctionner dans notre cas.

Hors ligne

 

#11 Thu 13 December 2018 17:21

jlavaure
Participant assidu
Date d'inscription: 20 Apr 2013
Messages: 201

Re: QGIS 2.14 + GRASS: HELP - MNE à partir BDALTI 75m et BDTOPO

Bonjour à tous,

Je remonte le sujet car je rencontre une nouvelle difficulté sur ce sujet :

J'ai décidé de ré-échantillonner mes raster à 1mx1m (au lieu de 5mx5mx).
Résultat, j'ai un raster de 4Go (MNT) + 2Go (bati) + 120Mo (Végétation).

Lorsque je lance la formule d'addition classique dans la calculatrice raster j'ai droit à un minidump ou au message d'erreur suivant : Raster calculator: Insufficient memory available for operation

Après avoir un peu fouillé sur le net, QGIS semble avoir du mal avec les traitements des "gros" raster.
J'ai donc tenté de passer par GRASS avec r.mapcalc. Disons que ça tourne mais ne produit rien, même pas un raster vide.

Je précise que j'utilise QGIS 2.15 et GRASS 7 dans QGIS.

Une solution ?

Merci d'avance

J.

Hors ligne

 

Pied de page des forums

Powered by FluxBB