Le portail francophone de la géomatique

jack tippit · Mon 18 April 2011 14:24

Bonjour à tous,
Je travaille sur des zones viticoles de terrasse à partir d'un MNT à 1m. Je voulais pouvoir isoler les zones de terrasse, et je ne savais pas trop comment m'y prendre. J'ai donc utilisé r.mapcalc pour tester chaque pixel de mon MNT en fonction du voisin en reprenant une technique trouver ici http://montagnasicura.it/j/index.php?op … Itemid=69.
Ce qui me donne :

Code:

r.mapcalc

rupturePente50 = if(abs(SlopeDEM@etienne - SlopeDEM@etienne[1,1]>50)||(abs(SlopeDEM@etienne - SlopeDEM@etienne[0,1])>50) || (abs(SlopeDEM@etienne - SlopeDEM@etienne[-1,1])>50)||(abs(SlopeDEM@etienne - SlopeDEM@etienne[-1,0])>50) || (abs(SlopeDEM@etienne - SlopeDEM@etienne[0,-1])>50)||(abs(SlopeDEM@etienne - SlopeDEM@etienne[1,-1])>50) || (abs(SlopeDEM@etienne - SlopeDEM@etienne[-1,-1])>50)||(abs(SlopeDEM@etienne - SlopeDEM@etienne[1,0])>50))

end

et qui donne un résultat intéressant, mais où il reste encore du bruit... Donc je voulais savoir si l'un d'entre vous avait une idée à me donner pour que je puisse isoler mes petits murs de pierre du bruit ambiant...

Mille mercis

Etienne

PS je vous met une image du rendu pour vous faire une idée

Fichier(s) joint(s) :
Pour accéder aux fichiers vous devez vous inscrire.

jack tippit · Tue 26 April 2011 13:06

Bonjour a tous,
Je reviens sur ce poste parce que j'ai trouvé une solution qui me semble acceptable.
J'ai utilisé successivement r.clump pour sur mon MNT pour rassembler les pixels en classe, puis j'ai ulisé r.li.shannon pour classer ces pixels et enfin, histoire de discrétiser encore et d'éliminer le "bruit" j'ai utilisé r.neighbors en utilisant la méthode mode à 10 pixels. Le résultat me semble très intéressant.
Voilà, pour moi le problème est résolu...

MathieuR · Fri 06 May 2011 18:37

merci bcp pour l'auto-réponse!
en quoi consiste l'utilisation de r.li.shannon?

MathieuR · Mon 16 May 2011 14:52

J'avoue que la méthode m'intéresse bcp!
Est-il possible d'avoir une vision du rendu final? après r.li.shannon voire r.neighbors?

jack tippit · Wed 18 May 2011 18:05

Bonjour Baobazz
Je suis désolé pour cette réponse tardive, entre les déplacements et les fluctuations d'internet je n'avais pas regardé depuis un moment ce poste. Donc pour reprendre dans l'ordre, r.li.shannon permet d'appliquer un index de diversité sur des variables (pour moi les altitudes). Ce qui permet de discrétiser les valeurs en sortie de r.clump (pour moi, mais on peut l'appliquer dans bien d'autres problèmes).
Je peux bien sûr t'envoyer mon document de travail, je vais juste essayer de le relire un peu pour l'orthographie... :-S

MathieuR · Thu 19 May 2011 07:16

Il n'y pas à t'excuser jack tippit!

Dans le premier rendu, tu obtiens un raster avec deux classes dominantes possédant chacune une rangée de valeurs (d'où le r.clump d'après): les zones de forte pente: les limites de terrasse, et le reste.

si je comprends bien, r.clump, t'a permis d'affecter une valeur unique de cellule à tes deux catégories.

Dans ton rendu final, je crois que tu souhaitais délimiter les différentes terrasses par une couleur unique, un peu à la manière des bassins versants générés par r.watershed, non? C'est pourquoi tu as utilisé r.li.shannon?

Mais je n'arrive pas trop à comprendre comment un indice de richesse paysagère / ou de biodiversité peut arriver à capturer les différentes terrasses. C'est pourquoi un rendu post-r.li.shannon m'intéresse!

Aussi, comme r.li.shannon se base sur une fenêtre mobile pour le calcul des valeurs, sur quoi t'es-tu basé pour régler la taille de la fenêtre?

Que de questions! ;-)

Prends ton temps

--
Et pour ceux qui souhaitent se familiariser avec les indicateurs de richesse paysagère:
http://ec.europa.eu/agriculture/publi/landscape/ch1.htm
étude basée sur corine land cover

Dernière modification par baobazz (Thu 19 May 2011 07:17)

jack tippit · Thu 19 May 2011 10:00

Bonjour Baobazz
pour reprendre les choses dans l'ordre, j'ai essayé de combiner deux approches pour déterminer les zones de terrasse. La première chose qui m'est venu a l'esprit est de rechercher les zones de rupture de pente c'est ce qui est illustrer par l'image que j'ai mise dans le premier poste. Donc en partant d'un MNT a mètre (oui oui bcp de point :-) ), j'ai utiliser la commande r.mapcalc que j'ai évoquée (toujours dans le premier poste) et j'ai essayer aussi avec r.param.scale et l'option profc qui donne un résultat relativement proche.
La seconde voie que j'ai utilisée en repartant de zéro (de mon MNT) j'ai fait un un r.clump pour voire s'il pouvait classer les altitudes entre elles... et c'est sur r.clump que j'ai effectué un r.li.shannon histoire de discrétiser la multitude de classes issues de shannon. Enfin avec r.neighbors j'ai viré les valeurs toutes seul dans mes plages de cellules pour aboutir à des grosses "patates". Et enfin en croisant les deux je me suis rendu compte que les terrasses trouver grace a r.param et le résultat de r.neighbors concordait... et une petite journée de terrain ont confirmé que les zones de terrasses étaient bien déterminées..

MathieuR · Thu 19 May 2011 10:30

ok, très intéressant tout ça!

donc, si je récapitule bien, il y a deux méthodes utilisées:
1-r.mapcalc sur les pentes ou r.param.scale
2-r.li.shannon sur fichier altitude puis r.clump puis r.neighbors

1 et 2 donnant des résultats similaires

la méthode 2 n'utilise pas les pentes mais plutôt la diversité des altitudes autour de chaque cellule. Ainsi, si on est en bordure de terrasse, l'indice de shannon sera élevé.

Ai-je juste?

Merci en tout cas! Cela montre toute la puissance de GRASS dans l'analyse paysagère!

jack tippit · Thu 19 May 2011 12:08

Oui c'est un tres bon résumé

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 Mon 18 April 2011 14:24

[GRASS] déterminaton des zones de terrasses viticole

Code:

#2 Tue 26 April 2011 13:06

Re: [GRASS] déterminaton des zones de terrasses viticole

#3 Fri 06 May 2011 18:37

Re: [GRASS] déterminaton des zones de terrasses viticole

#4 Mon 16 May 2011 14:52

Re: [GRASS] déterminaton des zones de terrasses viticole

#5 Wed 18 May 2011 18:05

Re: [GRASS] déterminaton des zones de terrasses viticole

#6 Thu 19 May 2011 07:16

Re: [GRASS] déterminaton des zones de terrasses viticole

#7 Thu 19 May 2011 10:00

Re: [GRASS] déterminaton des zones de terrasses viticole

#8 Thu 19 May 2011 10:30

Re: [GRASS] déterminaton des zones de terrasses viticole

#9 Thu 19 May 2011 12:08

Re: [GRASS] déterminaton des zones de terrasses viticole

Pied de page des forums