Le portail francophone de la géomatique

jlavaure · Tue 25 July 2017 11:15

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)

Kadiri28 · Tue 25 July 2017 12:07

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

jlavaure · Tue 25 July 2017 12:17

@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 ?

Kadiri28 · Tue 25 July 2017 16:15

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.

jlavaure · Tue 25 July 2017 16:31

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 ^^

Thom25 · Tue 25 July 2017 17:16

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.

jlavaure · Tue 25 July 2017 17:59

@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.

Thom25 · Wed 26 July 2017 00:57

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. ^^

jlavaure · Wed 26 July 2017 09:10

@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.

Thom25 · Thu 27 July 2017 11:41

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

jlavaure · Thu 13 December 2018 17:21

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.

Communauté	Ressources	L'emploi	Le marché	L'association	Fils RSS
Forum Blogs Wiki Planet Sites amis	Agenda Biblio Fiches AMO-CNIG Paris PCGIS Téléchargements	Offres d'emploi Banque de CV Formations Démarche Métiers	Geo-Entreprises Geo-Communiqué Appels d'offres	L'association Qui sommes-nous ? Nous aider Bienfaiteurs Partenaires Historique Statuts	Tout le site Job Géomatique Services GeoBlogs ...	Faire un don Mentions légales

Le portail francophone de la géomatique

Annonce

#1 Tue 25 July 2017 11:15

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

#2 Tue 25 July 2017 12:07

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

#3 Tue 25 July 2017 12:17

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

#4 Tue 25 July 2017 16:15

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

#5 Tue 25 July 2017 16:31

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

#6 Tue 25 July 2017 17:16

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

#7 Tue 25 July 2017 17:59

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

#8 Wed 26 July 2017 00:57

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

#9 Wed 26 July 2017 09:10

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

#10 Thu 27 July 2017 11:41

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

#11 Thu 13 December 2018 17:21

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

Pied de page des forums