Le portail francophone de la géomatique

MathieuB · Wed 17 November 2010 15:02

Bonjour à tous,

je souhaite sur certains sites d'études produire une couche vecteur, constituée de maille carrées de 100m de coté, renseignée avec les valeurs moyennes du MNT pour l'altitude, la pente et l'exposition.

J'ai créé ma grille avec GRASS :

Code:

v.mkgrid map=grille_lr_100m grid=2600,3000 --o position=coor coor=588600,6137100 box=100,100

Je souhaite maintenant intérroger le mnt pour l'altitude :

Code:

v.rast.stats -c -e vector=grille_lr_100m@PERMANENT raster=mnt11@PERMANENT colprefix=mnt_

L'erreur suivante est renvoyée, aprés un avertissement :

Code:

WARNING: This module has been replaced by r.category
ERROR: No categories found in raster map

Avez vous été confrontés à ce problème, avez vous une piste de contournement ?

Merci d'avance,

MathieuR · Wed 17 November 2010 15:57

Comment obtenir les valeurs moyennes avec ton MNT?
De plus, il faut qu'elles soient obtenues selon le même maillage

Sinon, avec QGIS, une fois que tu as ton raster MNT interpolé, tu peux obtenir les valeurs comme ceci
1-on génère la grille
Vecteur > Outils d'analyse > Grille Vecteur => grilleVecteur.shp

2-On rasterise
Raster > Rasteriser en prenant l'ID comme attribut => grilleVecteur.tiff

3-On croise les 2 rasters
Extension > Raster Calculator > on multiple les valeurs du MNT par celui du raster généré et on divise par cette même valeur [N_MNT_02B]@1 * [grilleVecteur]@1 / [grilleVecteur]@1=>grilleVecteur-MNT.tiff

4-On polygonise
Raster > Polygoniser=>grilleVecteur-MNT.shp

ppluvinet · Wed 17 November 2010 16:08

Salut Mathieu,

voilà ce que j'ai lu :

v.rast.stats still tries to use the old 'r.cats' module to gather raster categories. This should be updated to use the new 'r.category' module. This may also allow extended functionality for v.rast.stats

http://trac.osgeo.org/grass/ticket/161

Or ce problème a été souligné il y a plus de 3 ans.
Une solution serait peut-être de modifier la fonction v.rast.stats légèrement. Si je me souviens bien les fonctions sont accessibles dans des fichiers texte. Je me rappelle d'ailleurs avoir modifié cette fonction car je ne souhaitais n'avoir que la moyenne et pas toutes les autres statistiques qui ralentissait les calculs. Ceci dit, ca fait 2 ans que je n'ai quasiment pas touché à Grass donc je peux pas beaucoup t'aider.

MathieuB · Wed 17 November 2010 16:46

Merci à tous les deux,

baobazz, la fonction v.rast.stats semble faite pour cela :

Calculates univariate statistics from a GRASS raster map based on vector polygons and uploads statistics to new attribute columns

Je vais tester ta méthode avec QGIS.

pascal, merci pour le tuyau concernant les fonctions. Comme on en avait parlé, la piste wktraster de postgis permettra de répondre à ce genre de problème. Pour l'instant, je rencontre un problème avec le chargeur python ratser2pgsql, mais ça sera l'objet d'un autre post.

Si tu est de passage par Montpellier un de ces jours...

Je viendrai dire comment je m'en suis sorti.

MathieuR · Wed 17 November 2010 16:51

ok dans ce cas, cela me semble impossible avec QGIS...

MathieuR · Fri 19 November 2010 12:32

Il est marqué

WARNING: This module has been replaced by r.category

Il y a donc sûrement une méthode utilisant les catégories...

ce que je vois comme méthode est la suivante:
1-Donner une valeur à chaque carré de la maille vectorielle

2-Cette valeur sert à renseigner la catégorie de chaque pixel sur la base d'une jointure vecteur-raster
si pixel dans maille 1, on donne la catégorie 1

3-Ensuite, on affecte à tous les pixels d'une même catégorie la moyenne des valeurs des pixels de cette même catégorie.

Par contre, je ne vois pas trop comment faire sous GRASS, quels modules utiliser.

MathieuR · Fri 19 November 2010 14:30

En fait, je pense qu'il faudrait rasteriser la grille vecteur en choisissant des valeurs différentes pour chaque carré: v.to.rast. On l'appellerait grilleVecteur
http://grass.fbk.eu/gdp/html_grass63/r.average.html

Ensuite, on utilise r.average de cette façon:

Code:

r.average [-c] base=[grilleVecteur] cover=[MNT] output=[MNTGrille]

On peut vectoriser le raster par la suite avec r.to.vect feature=area

Dernière modification par baobazz (Fri 19 November 2010 14:58)

MathieuR · Mon 22 November 2010 11:12

Je viens d'effectuer un test et cela fonctionne très bien comme cela.
cf PJ
[img]http://data.imagup.com/1/1105093405.gif[/img]

MathieuB · Mon 22 November 2010 14:02

Bonjour baobazz,

merci beaucoup pour tes recherches. J'ai été pris par autre chose en fin de semaine et je me remettrai sur ce problème mercredi je pense.

Je viendrai faire un retour.

Bonne après-midi,

[EDIT]le 9/12/2010 : Je n'ai pas oublié de faire un retour, mais je n'ai pas encore pu me replonger sur le sujet.[/EDIT]

Dernière modification par Mathieu BOSSAERT (Thu 09 December 2010 10:21)

MathieuB · Fri 10 December 2010 09:11

Bonjour,

voici une page dédiée à la question sur le trés bon site http://datagistips.blogspot.com.

http://datagistips.blogspot.com/2010/11 … IS+tips%29

Pas encore eu le temps d'appliquer la méthode mais merci pour ces pistes et tutoriels.

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

#1 Wed 17 November 2010 15:02

QGIS/Grass : renseigner une grille vecteur depuis données d'un raster

Code:

Code:

Code:

#2 Wed 17 November 2010 15:57

Re: QGIS/Grass : renseigner une grille vecteur depuis données d'un raster

#3 Wed 17 November 2010 16:08

Re: QGIS/Grass : renseigner une grille vecteur depuis données d'un raster

#4 Wed 17 November 2010 16:46

Re: QGIS/Grass : renseigner une grille vecteur depuis données d'un raster

#5 Wed 17 November 2010 16:51

Re: QGIS/Grass : renseigner une grille vecteur depuis données d'un raster

#6 Fri 19 November 2010 12:32

Re: QGIS/Grass : renseigner une grille vecteur depuis données d'un raster

#7 Fri 19 November 2010 14:30

Re: QGIS/Grass : renseigner une grille vecteur depuis données d'un raster

Code:

#8 Mon 22 November 2010 11:12

Re: QGIS/Grass : renseigner une grille vecteur depuis données d'un raster

#9 Mon 22 November 2010 14:02

Re: QGIS/Grass : renseigner une grille vecteur depuis données d'un raster

#10 Fri 10 December 2010 09:11

Re: QGIS/Grass : renseigner une grille vecteur depuis données d'un raster

Pied de page des forums