#1 Tue 02 June 2009 14:19
- Eddy Rivers
- Participant actif
- Lieu: France
- Date d'inscription: 6 Aug 2007
- Messages: 53
- Site web
[ArcGis 9.1] Extraire des données par maille
Bonjour à tous.
J'espère que le weekend fût bon !
J'ai une question concernant l'extraction de données (je mets le terme en italique car je ne suis pas sûr qu'il faille employer celui-ci...).
Je m'explique :
- j'ai une table_1 qui contient des données (flux de gaz) sous forme de polygones (en rouge sur le belle figure ci-joint) ;
- j'ai un raster_1 régulier (en noir sur la fig.) ;
- j'aimerais extraire (répartir, ou encore intersecter : quel est le terme exact ?) les données de la table_1 dans les cases du raster_1.
Donc prendre en compte la part de surface des polygones présents dans chaque maille de mon raster, puis multiplier cette part de surface par le flux du polygone correspondant, pour connaitre mon flux par case de raster...
Clair ou pas clair ?
Et donc ma question : comment procéder (sachant que l'interpolation IdW ou le krigeage ne me satisfont pas) ?
Dans tous les cas, d'avance merci !
Franck.
NB : j'ai Spatial Analyst !
Dernière modification par Eddy Rivers (Tue 02 June 2009 15:04)
Franck (aménagement du territoire, écologie, écotourisme)
Hors ligne
#2 Tue 02 June 2009 19:11
- carteq
- Participant assidu
- Date d'inscription: 8 Mar 2008
- Messages: 829
Re: [ArcGis 9.1] Extraire des données par maille
Je ne te donnerai pas une démarche complète, mais plutôt mon avis: je procéderais plutôt en faisant une grille vectorielle (fishnet) plutôt qu'avec un raster. Ça sera plus facile de compiler une fois l'intersection de la grille et des flux réalisée.
Dernière modification par carteq (Tue 02 June 2009 19:11)
Hors ligne
#3 Wed 03 June 2009 15:40
- Eddy Rivers
- Participant actif
- Lieu: France
- Date d'inscription: 6 Aug 2007
- Messages: 53
- Site web
Re: [ArcGis 9.1] Extraire des données par maille
Je ne te donnerai pas une démarche complète, ...
Snif
... mais plutôt mon avis :
Ah !
je procéderais plutôt en faisant une grille vectorielle (fishnet) plutôt qu'avec un raster. Ça sera plus facile de compiler une fois l'intersection de la grille et des flux réalisées.
D'accord.
Donc par contre vous procéderiez avec la fonction "Intersect" ?
Parce que je n'obtiens pas de résultats concluants...
A l'instant j'étais en train de tester avec la fonction "Union", qui me permet de "découper" chaque partie de polygone à l'intérieur de chaque maille ; c'est un début dans ce que j'attends !
Mais par contre dans ma table_union finale, j'ai mes surfaces de polygones par nouvel identifiant, et en fait ce n'est pas la valeur de la surface découpée qui apparait, mais la valeur de la surface de mon polygone original.
En plus simple (j'espère) :
- j'ai un polygone A, qui fait 100m2 ;
- il est présent dans 3 carrés de ma maille, à resp. 10, 20 et 70%.
Je n'obtiens pas un tableau du genre :
Fid_union Pol Area
1 A 10
2 A 20
3 A 70
Mais plutôt :
Fid_union Pol Area
1 A 100
2 A 100
3 A 100
Avez vous une idée quant à comment récupérer une nouvelle surface après un une union ?
Merci d'avance !
Franck (aménagement du territoire, écologie, écotourisme)
Hors ligne
#4 Wed 03 June 2009 16:20
- Eddy Rivers
- Participant actif
- Lieu: France
- Date d'inscription: 6 Aug 2007
- Messages: 53
- Site web
Re: [ArcGis 9.1] Extraire des données par maille
Après recherches, il semblerait que je puisse utiliser l'outil Analyse par mailles, disponible ici !
Franck (aménagement du territoire, écologie, écotourisme)
Hors ligne
#5 Wed 03 June 2009 16:56
- pbader
- Participant actif
- Date d'inscription: 21 Sep 2007
- Messages: 79
Re: [ArcGis 9.1] Extraire des données par maille
Bonjour.
Pourquoi ne pas recalculer la surface des nouveaux polygones dans un nouveau champ 'Area2' et faire le ratio '(Area2 * 100) / Area' ?
P
Hors ligne
#6 Wed 03 June 2009 17:36
- Eddy Rivers
- Participant actif
- Lieu: France
- Date d'inscription: 6 Aug 2007
- Messages: 53
- Site web
Re: [ArcGis 9.1] Extraire des données par maille
Pourquoi ne pas recalculer la surface des nouveaux polygones dans un nouveau champ 'Area2' et faire le ratio '(Area2 * 100) / Area' ?
Tout simplement parce que je ne sais pas faire !
Mais au vu de votre commentaire, c'est qu'il doit y avoir un outil pour le faire ; je vais le chercher de ce pas !
Merci.
Edit : je pense qu'il faut passer par la calculatrice de valeurs de champs, mais pour calculer area_2 je ne vois pas comment procéder... C'est la surface de mon polygone de départ dans ma maille d'arrivée ? Si oui, et bien c'est tout à fait ce que je cherche à obtenir....
Dernière modification par Eddy Rivers (Wed 03 June 2009 17:46)
Franck (aménagement du territoire, écologie, écotourisme)
Hors ligne
#7 Wed 03 June 2009 17:43
- carteq
- Participant assidu
- Date d'inscription: 8 Mar 2008
- Messages: 829
Re: [ArcGis 9.1] Extraire des données par maille
Est-ce possible de poster un fichier avec un échantillon ?
Hors ligne
#8 Thu 04 June 2009 10:13
- Eddy Rivers
- Participant actif
- Lieu: France
- Date d'inscription: 6 Aug 2007
- Messages: 53
- Site web
Re: [ArcGis 9.1] Extraire des données par maille
Est-ce possible de poster un fichier avec un échantillon ?
Oui !
En PJ (entre 5 et 6 Mo, j'espère que ça passe**), mon fichier *.shp avec l'union de ma grille régulière et de mes polygones pour lesquels je voudrais connaitre leur proportion dans chaque pixel de ma grille.
Merci.
Franck.
** D'ailleurs, j'utilise Free pour uploader et échanger des fichiers, mais souvent chez mes contacts, cela rame.
Qu'utilisent les gens du forum pour partager de gros fichiers (à partir de 100 Mo), en hébergeur gratuit ?
Franck (aménagement du territoire, écologie, écotourisme)
Hors ligne
#9 Thu 04 June 2009 12:51
- Eddy Rivers
- Participant actif
- Lieu: France
- Date d'inscription: 6 Aug 2007
- Messages: 53
- Site web
Re: [ArcGis 9.1] Extraire des données par maille
Pour le suivi : il semblerait que j'ai trouvé la bonne manière de procéder (enfin j'espère !).
Pour le calcul des surfaces de polygones :
- ouvrir la table attributaire contenant mes données unies ;
- créer un nouveau champ double ;
- cliquer droit sur le nouveau champ créé, "Calculer les valeurs" ;
- cocher la case "Avancé" ;
- taper ce bout de code :
Dim dblArea as double
Dim pArea as IArea
Set pArea = [shape]
dblArea = pArea.area
- enfin, mettre "dblArea" dans la dernière zone de saisie.
Cette manip. ayant été vue ici.
Normalement je devrais pouvoir récupérer les pourcentages d'occupation de mes polygones dans mes zones.
Dernière modification par Eddy Rivers (Thu 04 June 2009 12:52)
Franck (aménagement du territoire, écologie, écotourisme)
Hors ligne
#10 Wed 26 May 2010 17:55
- loloti
- Juste Inscrit !
- Date d'inscription: 26 May 2010
- Messages: 7
Re: [ArcGis 9.1] Extraire des données par maille
Bonjour,
Je voudrais relancer ce sujet où je rencontre exactement le même problème. J'ai une couche de polygones avec une valeur de variable et j'ai effectué l'analyse par maille avec des pixels de 250*250mètres découpant la couche. Les trois colonnes (EFFECTIF, S_BRUTE et S_RATIO) ont été remplies mais je ne sais pas comment obtenir la nouvelle valeur de ma variable à l'échelle des pixels??? Je suppose comme le définit le tutorial de l'analyse par maille, le S_RATIO intervient dans ce calcul mais je ne sais pas comment? :-(
Si quelqu'un aurait une piste, merci d'avance.
Bye
Laurent
Hors ligne
#11 Thu 27 May 2010 07:41
Re: [ArcGis 9.1] Extraire des données par maille
Bonjour,
il faut effectivement utiliser le champs S_RATIO.
Si la maille a été créée avec des points S_RATIO = S_BRUTE
Si la maille a été créée avec des polygones c'est bien le champs S_RATIO qui
doit être utilisé et dont la somme doit correspondre au chiffre de départ
Par exemple pour moi dans mes calculs de densités le champs S_RATIO
correspond au nombre de logements par maille.
Bon travail
AlineC
Hors ligne
#12 Mon 31 May 2010 21:19
- loloti
- Juste Inscrit !
- Date d'inscription: 26 May 2010
- Messages: 7
Re: [ArcGis 9.1] Extraire des données par maille
Bonjour,
Merci pour la précision pour le champ S_RATIO. Celui-ci contient des valeurs comprises entre 0,1 et 0,9 pour mes 72000 pixels. Mon problème est comment calculer la nouvelle valeur lorsqu'un pixel considéré coupe deux ou plus de polygones (colonne EFFECTIF), la colonne S_BRUTE me donne le total des valeurs des polygones compris dans le pixel, jusqu'à là aucun souci. Mais comment à partir du S_RATIO qu'on m'indique définir la proportion de la valeur de chaque polygone qui concerne le pixel en question pour avoir une valeur estimée proche de la réalité.
J'espère j'ai été clair dans mes précisions. Merci d'avance.
Laurent
Hors ligne