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

GEODATA DAYS 2024

#1 Thu 29 November 2007 19:04

histoirede
Participant occasionnel
Date d'inscription: 10 May 2007
Messages: 14

[ArcGis 9.x] conversion polygones en grille définie

Bonjour !

Je m'excuse d'avance si je n'utilise pas le bon vocabulaire ... je ne suis pas spécialiste en SIG ...

Je possède une couche de polygones avec des émissions gazeuses (polygones=sols), et je souhaite "rasteriser" cette couche sur une grille régulière définie que j'ai construite à partir d'un outil dispo en OpenSource (Maillage). J'ai défini la taille de la grille selon les coordonnées longitude/latitude ainsi que le nombre de mailles. En fait, je ne sais pas quelle méthode utiliser ?
Pièce jointe : la couche grille et la couche polygones (émissions) de différente couleurs

Je pensais calculer une valeur moyenne pondérée des émissions par élément de grille à partir de la surface des polygones contenus dans chaque élément de grille et de la surface des éléments de grille (la même surface!!). Je dois d'abord calculer les surfaces d'éléments de grille, puis intersecter la couche grille et la couche de polygones (avec les émissions) ou faire une jointure spatiale via les centroides "grille-polygones" pour pouvoir calculer les surfaces de polygones redécoupés dans les éléments de grille. est ce correct ?

En fait, je pensais utiliser l'outil "polygons to raster" dans Spatial Analyst, mais il ne permet pas de "projeter" les polygones sur un maillage défini. On peut seulement choisir la taille de la cellule, ce qui ne me convient pas puisque en faisant ça, les cellules sont décallés par rapport aux cellules du maillage que j'ai défini.

je ne sais pas si ma méthode (moyenne pondérée, jointure spatiale ...) est  propre ou pas ?? Qu'en pensez vous ?

S'il manque des infos (ce qui est bien possible), n'hésitez aps à me demander ...

Tout commentaire est vraiment le bienvenu (ça presse un peu ...)

Merci d'avance !


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

Hors ligne

 

#2 Thu 29 November 2007 19:57

Sylvie
Membre
Date d'inscription: 5 Sep 2005
Messages: 3066

Re: [ArcGis 9.x] conversion polygones en grille définie

Bonsoir,

je pensais utiliser l'outil "polygons to raster" dans Spatial Analyst, mais il ne permet pas de "projeter" les polygones sur un maillage défini.On peut seulement choisir la taille de la cellule, ce qui ne me convient pas puisque en faisant ça, les cellules sont décallés par rapport aux cellules du maillage que j'ai défini.


Non, non : tu peux caler ton raster sur ton maillage. Dans Spatial Analyst/Options, tu définis un masque d'analyse = ta couche maillage

Bon courage
Sylvie

Hors ligne

 

#3 Fri 30 November 2007 09:34

histoirede
Participant occasionnel
Date d'inscription: 10 May 2007
Messages: 14

Re: [ArcGis 9.x] conversion polygones en grille définie

Merci beaucoup !
Effectivement ça marche ! la couche raster est calée sur la couche maillage.

Une autre chose, j'ai une couche vectorielle d'occupation des sols et toujours la même couche maillage. J'aimerais connaître le % d'une classe (d'occupation des sols) dans chaque élément de grille.

Les étapes de calcul pourraient être :
- calculer les surfaces des éléménts de grille
- jointure spatiale de la couche grille avec la couche occupation des sols
- calculer les surfaces d'intersection de la classe d'occupation des sols en question
- calculer les % 'doccupation = surfaces "intersectées"/surfaces des éléments de grille
- rasteriser la couche d'intersection sur la couche maillage (grace au masque d'analyse posé)

Qu'en pensez vous ? est ce correct ?

Hors ligne

 

#4 Mon 03 December 2007 14:44

histoirede
Participant occasionnel
Date d'inscription: 10 May 2007
Messages: 14

Re: [ArcGis 9.x] conversion polygones en grille définie

pas d'idée ? (humhumpfff)

Hors ligne

 

#5 Mon 03 December 2007 15:15

alban
Participant assidu
Lieu: Chizé
Date d'inscription: 16 May 2007
Messages: 156
Site web

Re: [ArcGis 9.x] conversion polygones en grille définie

Bonjour,

Si j´ai bien compris ta requête, j´aurais une manipulation plus simple à te proposer :
- tu crées un identifiant unique pour chaque cellule de ton maillage (en vecteur)
- tu unies ta couche d´occupation du sol (union de 2 couches vectorielles)
- tu exportes le résultat dans une table pour excel
- dans Excel, tu crées un rapport de tableau croisé dynamique ("données / rapport de tableau...") avec les aires comme valeur, les occupations du sol comme entête de colonne et les Identifiants de parcelles comme titre de ligne.

En espérant avoir répondu à ta question,

Alban

Hors ligne

 

#6 Mon 03 December 2007 16:41

histoirede
Participant occasionnel
Date d'inscription: 10 May 2007
Messages: 14

Re: [ArcGis 9.x] conversion polygones en grille définie

merci !

ça ne me parle pas un rapport d tableau croisé dynamique ?? ça permet de constituer une matrice (lignes x colonnes) ? je vais regarder ce que tu veux dire ... ps : je travaille sous Open Office en tableur

Ok pourquoi pas ...
Je pourrai mettre les identifiants des cellules du maillage comme titre de ligne, puisque je cherche le % de plusieurs classes d'occupation des sols par élément de grille.
et après, je peux rebasculer sous ArcGis et projeter mes  % par classe d'occupation pour chacune des cellules ...

Hors ligne

 

#7 Mon 03 December 2007 16:56

histoirede
Participant occasionnel
Date d'inscription: 10 May 2007
Messages: 14

Re: [ArcGis 9.x] conversion polygones en grille définie

Une autre question, plutôt de méthodologie si quelqu'un a une idée ...

Lorsqu'on souhaite avoir des infos (surface de forêts via une base d'occupation des sols par exemple ...) sur un maillage régulier sachant que les données sont conetnus dans des polygones (ou vecteurs) :
- vaut il mieux rasteriser les polygones avec l'outil de Spatial Analyst : "polygones to raster"  ? si j'ai bien compris : l'outil prend la valeur du polygone qui remplit majoritairement la cellule ou prend le plus proche voisin entre les centroides (polygone - raster) mais une partie de l'information est perdue ??
- OU, n'est il pas mieux de garder les infos sur des vecteurs  et calculer les contributions des différents polygones contenus dans une grille (format vecteur) en utilisant différents outils (Union, calcul de surface ...)
Qu'en pensez vous ?? merci encore ...

Je ne sais pas si j'ai mis tous les éléments ... ça me semble clair mais je manque de recul ... (débutante SIG)
Demandez moi s'il manque quelquechose !

Hors ligne

 

#8 Tue 04 December 2007 08:48

histoirede
Participant occasionnel
Date d'inscription: 10 May 2007
Messages: 14

Re: [ArcGis 9.x] conversion polygones en grille définie

pas d'avis ?? svp svp svp, oh maitres du sig smile

Hors ligne

 

#9 Tue 04 December 2007 09:33

alban
Participant assidu
Lieu: Chizé
Date d'inscription: 16 May 2007
Messages: 156
Site web

Re: [ArcGis 9.x] conversion polygones en grille définie

Bonjour (je ne pense pas être 1 maitre du SIG mais je me lance quand même),

Selon moi, ta question amène un problème assez récurrent : celui de l´échelle.

J´explique : si tu travailles en raster, tu vas devoir choisir une taille de cellule. Cette taille aura en effet des répercussions sur toutes les analyses subséquentes. Avec le vecteur, c´est aussi vrai mais on le ressent moins facilement : il s´agit de l´échelle à laquelle a été digitalisée tes données (cf. BD CARTO vs. BD Topo par exemple).

Maintenant, comment définir ton échelle ? Il n´y a pas vraiment de réponse simple. Le mieux est de la définir en accord avec ton étude. Les limites de communes ne demandent pas la même échelle qu´une étude des bâtiments. Tu peux la définir en fonction de ton jeu de données (mes bâtiments sont en moyenne distants de x mètres donc je prends telle taille de cellule).

Un 2ème aspect dans la différence raster/vecteur. Avec le raster, tu te focalises sur la localisation, avec les vecteurs tu regardes surtout les entités. Dans le 1er cas, tu fais l´hypothèse que le milieu est continu (discret avec le vecteur).

Revenons désormais à ton étude. Si considérer ton milieu comme continu ne te pose pas de problème, je te conseillerai plutôt de tout faire en raster car beaucoup d´opérations s´en trouvent simplifiées. Second avantage : tous les traitements seront faits à la même échelle.

En espérant avoir été clair,

PS : existe-t´il un tableau croisé dans Open Office (j´ai pas encore fait le pas wink ) ?
PS2 : je crois avoir mieux compris ce que tu veux... peut-être tu peux gagner du temps en rasterisant à une échelle fine et en calculant ensuite la diversité ? (cf. "statistiques de voisinage" dans "Spatial analyst"

Hors ligne

 

#10 Sun 09 December 2007 19:32

histoirede
Participant occasionnel
Date d'inscription: 10 May 2007
Messages: 14

Re: [ArcGis 9.x] conversion polygones en grille définie

Merci beaucoup pour ton explication bien pédagogue !! C'est très clair ...

Mon étude s'est faite sur des vecteurs, car la plupart des bases utilisées (occupation des sols, base "type de sols" ...) était disponible à cette échelle et les infos des bases étaient vraiment importantes. L'échelle vectorielle a permis d'être bien précis. L'analyse a été réalisée à partir des  vecteurs également ...

Ensuite, le passage à une grille régulière ou raster est une 2ème étape de mon étude, et cette transformation est nécessaire pour caler les données sur le maillage "régulier" d'un modèle (pas d'analyse à priori, seulement une transformation). Les données géoréférencées seront des "données d'entrée" à un modèle... Du coup, effectivement il est important que la localisation des données soit précise et la taille des cellules correspond à celles du maillage du modèle.

J'ai une nouvelle question :
Je m'aperçois lorsqu'on veut rasteriser des entités (outil dans "Spatial Analyst" Entités to raster) un seul champ (ou colonne) peut être sélectionné.
J'explique comment sont définies mes couches de vecteurs (pour la rasterisation) :
à chaque entête de ligne, correspond un n° d'identifiant de polygone (surface de sol) et à chaque entête de colonne (ou champ), il s'agit d'une date. Les valeurs correspondent à des valeurs journalières d'émissions gazeuses. J'ai besoin de rasteriser tous les champs.

Donc ma question est : comment fait on une rasterisation "multi-colonne" ou "multi -champ" ?
Est ce que ce travail peut se faire sur une seule couche ?
j'ai besoin de constituer un même type de tableau d'émissions gazeuses (identifiants x détes) mais  les entêtes de ligne seront des identifiants de ma grille définie.

J'espère que mon explication est claire ... (j'en doute)

Merci d'avance

Hors ligne

 

#11 Mon 10 December 2007 11:49

Jeirhome
Membre
Lieu: Liverion
Date d'inscription: 22 Aug 2006
Messages: 4298
Site web

Re: [ArcGis 9.x] conversion polygones en grille définie

Bonjour,

Parfois, ils y a plusieurs colonnes dans des fichiers RASTER. Allons savoir sous quelle condition exactement un tel fichier peut contenir plusieurs champs. Je vois deux solutions possibles :
- Créer un table multi-champs que l'on joint à la couche raster.
- Tricher et concaténer toutes les infos dans un seul champ (problème d'extraction des données ensuite).


Jérôme Cuinet
L'avantage de la Chine, c'est que le soleil se couche plus tard !

Hors ligne

 

#12 Mon 10 December 2007 14:45

histoirede
Participant occasionnel
Date d'inscription: 10 May 2007
Messages: 14

Re: [ArcGis 9.x] conversion polygones en grille définie

Merci pour ces infos.
Concernant la 1ère méthode, pourrais tu me dire comment crée t'on une table multi - champ ?

Sinon, je pense que je ne comprends pas bien ce qu'est un raster ? ou je fais mal le passage "entités to raster" pour un champ avec Spatial analyst. Je souhaitais obtenir en rasterisant ma couche vectorielle de données, un tableau similaire à la couche vectorielle (x lignes = x polygones) mais où chaque ligne correspond à une cellule (ou pixel ou raster ou maille). En réalité, la table du raster créé contient comme entête de ligne, la valeur du champ rasterisé et comme en tête de colonne, le nombre de cellules (count) ayant cette valeur.

J'imagine que la manipulation n'est pas bonne ... mais peut-être la jointure d'une table multi champ avec le raster permet de constituer un raster multi champs, et dont le nombre de lignes correpond au  nombre de cellules ???

Hors ligne

 

#13 Mon 10 December 2007 15:20

Jeirhome
Membre
Lieu: Liverion
Date d'inscription: 22 Aug 2006
Messages: 4298
Site web

Re: [ArcGis 9.x] conversion polygones en grille définie

La table multi-champs est une simple table .dbf. Lorsque tu affiches les informations attributaires de ta couche vectorielle, tu fais exporter Table... et tu as ce qu'il faut.

Un raster n'est une simple abréviation pour couche vectorielle maillée. Le raster est plutôt le format de dessin classique pour un ordinateur (jpeg, tiff, bmp...). Par définition, le nombre de maille est très vite très important, facilement des dizaines de milliers, des millions ou plus. Le traitement des rasters se fait le plus souvent par statistiques. Même si on peut consulter la valeur d'un pixel en particulier, il n'est pas envisageable de traiter les pixels un à un et donc d'afficher une table attributaire comme si le raster était une simple couche vectoriellle en forme de maille.

Si tu veux avoir un nombre de lignes correspondant au nombre de cellules, il faut rester en couche vectorielle, ne pas passer à du raster. Est-ce un peu plus clair ?


Jérôme Cuinet
L'avantage de la Chine, c'est que le soleil se couche plus tard !

Hors ligne

 

#14 Mon 10 December 2007 18:05

histoirede
Participant occasionnel
Date d'inscription: 10 May 2007
Messages: 14

Re: [ArcGis 9.x] conversion polygones en grille définie

Oui complément, merci beaucoup ! je pense que je m'étais bien mélangée les pinceaux avec la notion de raster.
Du coup, je vais rester en  couche vectorielle et passer à des vecteurs en forme de maille smile
merci bien .....

Hors ligne

 

#15 Mon 10 December 2007 18:16

histoirede
Participant occasionnel
Date d'inscription: 10 May 2007
Messages: 14

Re: [ArcGis 9.x] conversion polygones en grille définie

alban a écrit:

existe-t´il un tableau croisé dans Open Office (j´ai pas encore fait le pas wink ) ?


Finalement j'ai choisi ta méthode pour calculer les pourcentages de classes d'occuaption des sols.

Il existe sous Openoffice un outil simialire de tableau croisé dynamique. Il s'agit de "Pilote de données" dans le signet Données.

L'utilisation est assez facile.
1. Tu colles les 3 colonnes qui t'intéressent dans un tableau,
Dans mon cas, ce sont les identifiants de maille, les types d'occupation des sols (numéro de nomenclature) et les aires des types d'occupation.
2. Tu cliques sur le signet "pilote des données",  tu conserves ta sélection (les 3 colonnes),
3. tu précises "quoi et où" comme "identifiant de maille pour chaque ligne".
Le tour ets joué (OOo calcule pendant 10 bonnes minutes), on obtient un tableau 'lignes x colonnes) et perso, ça m'a gagné pas mal de temps !

Merci bien ...

Hors ligne

 

Pied de page des forums

Powered by FluxBB