Suite

Exporter les résultats des comparaisons Select By Location de chaque ligne d'un fichier de formes avec des polygones d'un autre ?


J'utilise ArcMap 10 et j'essaie d'exporter plusieurs sélections (alias "Sélectionner par emplacement") d'un fichier à un autre.

Je commence par sélectionner la première ligne du shapefile A et je veux sélectionner tous les polygones du shapefile B qui se trouvent à l'intérieur du shapefile A (polygones B contenus dans A), pour ensuite les exporter en tant que shapefile_1.

Et ce processus doit être répété pour x enregistrements dans A pour produire plus de fichiers de formes (shapefile_1, shapefile_2,… shapefile_x).

Comment cela peut-il être fait?


À partir du commentaire, vous pouvez boucler sur chaque entité puis exécuter la sélection par emplacement sur une couche d'entités. Il existe différents types de boucles à cet effet 1) basées sur une sélection par attribut pour chaque FID dans une couche d'entités ou 2) directement basées sur la géométrie dans un curseur. Voici le code avec un curseur, qui est plus simple à mon humble avis. Vous obtenez la géométrie de chaque entité en tant que propriété de chaque ligne à l'aide de row.getValue(shapefieldname). Notez que pour les versions ultérieures d'ArcGIS, les curseurs arcpy.da sont meilleurs.

import arcpy fcA = "fcA.shp" fcB = "fcB.shp" arcpy.MakeFeatureLayer_management(fcB, "layerfcB") #la sélection est effectuée sur un calque desc = arcpy.Describe(fcA) shapefieldname = desc.ShapeFieldName # obtient le nom du champ contenant les informations géométriques i=0 avec arcpy.SearchCursor(fcA) comme curseur : pour la ligne dans le curseur : #boucle sur chaque entité de fcA i=i+1 arcpy.SelectByLocation_management("layerfcB", "within", row. getValue(shapefieldname)) #select basé sur la géométrie arcpy.CopyFeatures_management("layerfcB", "output_path" +str(i)+".shp") # exporter les entités sélectionnées dans une nouvelle classe d'entités

Faire une union de polygones dans GeoPandas, ou Shapely (en une seule géométrie)

J'essaie de trouver l'union de deux polygones dans GeoPandas et de produire une seule géométrie qui englobe les points des deux polygones en tant que sommets. La fonction geopandas.overlay me donne des polygones pour chaque union individuelle mais je voudrais un seul polygone.

Pour le contexte, j'utilise ceci pour combiner deux zones administratives en une seule zone (c'est-à-dire inclure un district municipal dans un pays).

L'exemple suivant provient du site geopandas et illustre ce que j'aimerais :

Aucune des géométries de sortie ne correspond à ce à quoi je m'attendais, c'est-à-dire la suivante :

Tout d'abord, comment puis-je sortir le polygone ci-dessus ( poly_union ) à partir des polygones d'entrée ( df1 , df2 ) à l'aide de GeoPandas ou de shapely ?

Deuxièmement, quelle est la nomenclature correcte associée à la géométrie ( poly_union ) que j'essaie de trouver ? Je l'appellerais un "union" mais chaque exemple que je trouve qui fait référence à des "unions" ne produit pas cette géométrie.


Si vous utilisez ArcGIS ArcCatalog :

Accédez à la table des matières, faites un clic droit sur .lyr, accédez à DONNÉES et sélectionnez l'option Exporter vers. Ici, vous pouvez l'enregistrer au format .shp.

Si vous utilisez ArcGIS ArcMap :

Si vous pouvez charger le fichier de couche dans ArcMap et que les données sont accessibles, vous devriez pouvoir cliquer avec le bouton droit sur la couche dans la table des matières et l'exporter en tant que fichier de formes. L'image ci-dessous vous aidera:

Cliquez sur « Exporter les données » et vous verrez cette boîte de dialogue. Cliquez sur le bouton mis en évidence dans l'image ci-dessous

Vous verrez maintenant une autre boîte de dialogue et vous pouvez choisir "Fichier de forme" dans le menu déroulant, comme indiqué dans l'image ci-dessous


Exporter les résultats des comparaisons Select By Location de chaque ligne dans un fichier de formes avec des polygones dans un autre ? - Systèmes d'information géographique

Tous les logiciels SIG doivent contenir des routines d'exportation de données. Pour certains logiciels SIG, tels qu'ArcInfo, qui a une position forte en tant qu'application de développement de données, les routines d'exportation de données ont été un élément précoce et important.

D'autres produits logiciels conçus à l'origine pour l'interrogation et l'affichage, tels qu'ArcView 3.x, ont généralement des capacités d'exportation limitées. Ces produits ont toujours été plus intéressés par l'obtention de données dans que d'obtenir des données en dehors.

ArcGIS propose de nombreuses options pour exporter des données.

Conversion de couches d'entités en fichiers de formes

Les couches d'entités peuvent être "exportées" vers des fichiers de formes. Cela signifie qu'une couverture ArcInfo, un dessin CAO, une couche d'événements ou même un autre fichier de formes peuvent être convertis en un fichier de formes. En fait, tout jeu de données vectorielles pris en charge (c'est-à-dire tout ce qui apparaît dans une vue sous forme de points, de lignes ou de polygones) peut être converti au format shapefile. Les jeux de données pris en charge sont décrits dans Intégration de données dans ArcGIS .

La conversion d'autres sources de données en fichiers de formes peut être pratique lorsque vous devez créer une source de données d'entités à déplacer sur un système de fichiers ou d'une machine à une autre. Si vous disposez d'un fichier de formes, vous pouvez copier les fichiers de composants sur des systèmes de fichiers ou des ordinateurs.

Comparez cela à la structure de données ArcInfo, dans laquelle une seule source de couche d'entités est constituée de fichiers de données répartis sur plusieurs répertoires, et il est impossible de gérer la couverture avec uniquement le système d'exploitation.

Lorsqu'un nouveau fichier de formes est créé à partir d'une couche existante, ArcGIS vous demande si vous souhaitez ajouter ce nouveau fichier de formes en tant que couche au bloc de données actuel. Si vous n'ajoutez pas la couche au bloc de données, le fichier de formes existera toujours sur le disque et pourra être utilisé ultérieurement.

Ici, l'original des stands est un jeu de données de couverture ArcInfo.

Ici, il sera converti en un fichier de forme.

Notez qu'une couche vectorielle peut être exportée vers l'un des trois types de données différents (classe d'entités File et Personal Geodatabase, Shapefile ou classe d'entités Spatial Database Engine) :

Vous pouvez voir qu'il y a maintenant deux couches représentant les peuplements, l'une provient d'une source de données de couverture et l'autre d'une source de données de fichier de formes.

Conversion des entités sélectionnées d'une couche en un fichier de formes

Parfois, vous souhaiterez peut-être afficher uniquement un sous-ensemble d'entités d'une couche donnée. Une façon d'y parvenir consiste à appliquer un filtre à la couche, de sorte que seules les entités sélectionnées s'affichent. Une autre méthode consiste à créer un nouveau jeu de données basé sur une sélection d'entités. Lorsqu'un ensemble d'entités (une ou plusieurs) est sélectionné et que la couche est convertie en un fichier de formes, seul l'ensemble sélectionné est écrit dans le nouveau fichier de formes.

Ici, seul le > 100 ans se dresse de la stands.shp shapefile sont convertis en leur propre shapefile :

Les peuplements sélectionnés sont surlignés en cyan.

Cet ensemble sélectionné est converti en un fichier de formes :

Noter: si vous souhaitez créer un fichier de formes composé de toutes les caractéristiques d'un calque, assurez-vous de n'avoir aucune sélection active sur le calque. Sinon, votre nouvelle couche ne contiendra que des entités sélectionnées !

Conversion de couches d'entités en géodatabase

Il est possible de convertir des couches d'entités et des sélections de couches d'entités en classes d'entités de géodatabase. Pour ce faire, vous devez disposer d'une géodatabase existante. La seule différence entre l'exportation vers un fichier de formes et une classe d'entités de géodatabase est que vous sélectionnez la géodatabase vers laquelle exporter et donnez un nom à la classe d'entités en sortie. Dans l'image suivante, la géodatabase personnelle est proj_mgmt.mdb et la classe d'entités est appelée des stands. La création de géodatabases sera abordée dans la leçon sur la création de couches d'entités.

Conversion de couches d'entités en rasters

Nous aborderons les données raster et l'analyse plus tard dans le terme, mais il est possible, sous certaines contraintes, d'effectuer une conversion entre les couches d'entités et les rasters. Cela n'est généralement effectué que lorsque vous devez analyser des données de point, de ligne et de polygone avec d'autres données raster, ou si vous disposez de données raster catégorielles pouvant être converties en polygones pour une analyse vectorielle.

Conversion de calques d'images en grilles

Les jeux de données raster pris en charge, y compris les images, peuvent être convertis en jeux de données de grille. Les grilles de sortie auront des valeurs de cellule identiques aux valeurs de pixel dans l'image d'entrée. Les images multibandes peuvent être converties en grilles, une grille de sortie par bande d'image d'entrée. Alors que les sources de données d'image et de grille sont des rasters, les grilles sont un type spécial de raster qui peut avoir des tables attributaires avec plusieurs champs (alors que les images simples n'ont qu'une seule valeur par pixel)

Nous aborderons l'exportation de grilles plus tard dans le cours. Voir l'exercice sur la télédétection pour un exemple de conversion d'une image LandSat TM à 7 bandes en une série de 7 grilles. La conversion des images monobande s'effectue de la même manière, sauf qu'une seule grille de sortie est créée.

Exportation de couches de grille au format raster générique

Une couche de grille sélectionnée peut être convertie en un fichier raster générique en utilisant le Fichier > Exporter la source de données choix de menus. Cela permet à une couche de grille d'être convertie en un format raster générique.

Les deux formats raster en sortie pris en charge sont ASCII et binaire (format IEEE).

Le fichier raster ASCII générique est un fichier unique composé de 6 lignes de données d'en-tête (descriptives et de géoréférencement), suivies de valeurs de point de maillage de treillis (cellule de grille).

La routine générique d'exportation au format de fichier raster binaire crée deux fichiers, un fichier d'en-tête ASCII et un fichier binaire contenant les données de grille réelles. Lors de l'archivage ou de la distribution de fichiers raster génériques binaires, assurez-vous d'inclure à la fois le fichier d'en-tête et le fichier binaire.

Ces deux options d'exportation créent des fichiers autonomes qui peuvent être déplacés sur des systèmes de fichiers ou des plates-formes. Les fichiers peuvent être importés à tout moment, ce qui recrée une grille au format ArcInfo. D'autres logiciels SIG ou de traitement d'images peuvent lire ce type de formats de fichiers raster non propriétaires.

Exportation de couches d'entités au format de génération ArcInfo

Le format de génération ArcInfo est un fichier ASCII qui peut être utilisé pour créer une couverture ArcInfo (vecteur). Étant donné que la structure d'un fichier de génération est simple, les fichiers de génération sont également utilisés dans de nombreuses autres applications qui ne peuvent pas gérer les formats de données SIG propriétaires.

Un exemple de script existe dans ArcGIS qui vous permettra de créer des fichiers de génération à partir de fichiers de formes.

Les fichiers de génération sont formatés comme ces exemples simples :

format de fichier de génération de points

exemple de fichier

sens (de chaque ligne)

id (du point), x, y
identifiant, x, y
identifiant, x, y
identifiant, x, y
identifiant, x, y
d'autres points
fin de fichier)

format de fichier de génération de ligne

exemple de fichier

sens (de chaque ligne)

101
2,3
4,3
4,6
FINIR
101
2,3
4,3
4,6
FINIR
102
2,1
3,2
FINIR
.
FINIR
identifiant (de ligne)
x, y (sommet)
x, y
x, y
fin de ligne)
identifiant (de ligne)
x, y
x, y
x, y
fin de ligne)
identifiant (de ligne)
x, y
x, y
fin de ligne)
autres lignes
fin de fichier)

format de fichier de génération de polygone

goûter

sens (de chaque ligne)

701,7,7.5
3,9
4,5
8,4
9,7
6,9
3,9
FINIR
702,20,31
19,19
18,36
31,35
27,15
FINIR
.
FINIR
id (du poly), x, y (du point label)
x, y (sommet)
x, y (sommet)
x, y (sommet)
x, y (sommet)
x, y (sommet)
x, y (sommet)
fin (du polygone)
id (du poly), x, y (du point label)
x, y (sommet)
x, y (sommet)
x, y (sommet)
x, y (sommet)
fin (du polygone)
autres polygones
fin de fichier)

Les fichiers de génération peuvent être utilisés pour recréer des couvertures ArcInfo et, s'ils sont utilisés conjointement avec des tables exportées, peuvent également recréer les données attributaires des couvertures. Vous n'aurez probablement jamais besoin d'utiliser des fichiers de génération, mais il s'agit d'un format de fichier très simple si vous avez besoin de forcer une solution pour importer des coordonnées dans une autre application.

Les tables sont facilement exportées vers n'importe quel format source (table de géodatabase personnelle, dBASE, texte délimité, INFO et table SDE). Comme pour la création de tout autre fichier, le type de fichier, le nom et l'emplacement sur le système de fichiers doivent être fournis pendant le processus d'exportation. L'exportation des tableaux s'effectue à l'aide de la Options > Exporter contrôle pendant qu'une table est ouverte et active.

S'il y a des sélections actives sur les tables, seuls les enregistrements sélectionnés seront exportés. Si vous souhaitez exporter tous les enregistrements de la table, veillez à effacer toutes les sélections actives.

Un mot d'avertissement : les fichiers dBASE qu'ArcGIS utilise pour stocker les données ont une limite de 10 caractères pour les noms de champ. Si vous avez une table source contenant des champs avec des noms de plus de 10 caractères, les tables exportées ArcGIS ne sortiront que les 10 premiers caractères des noms de champ.

Si vous exportez vers du texte délimité, les noms de champ complets seront exportés, même si les noms de champ comportent plus de 10 caractères.

"Exportation" des blocs de données et des mises en page

Les mises en page et les blocs de données peuvent être enregistrés en tant que fichiers graphiques en sélectionnant Fichier > Exporter la carte du menu. Les détails de ceci sont couverts dans l'exercice de laboratoire de cartographie.

Il existe un certain nombre de formats d'image de sortie différents. Les fichiers créés à partir de l'exportation de blocs de données et de mises en page peuvent être utilisés dans n'importe quelle application prenant en charge ces types de fichiers.

Exportation de scènes 3D au format VRML

Les vues ArcScene ou ArcGlobe, qui seront couvertes par l'analyse 3D et la modélisation de surface, peuvent être exportées sous forme de fichiers VRML (Virtual Reality Markup Language), qui peuvent être chargés dans Netscape ou d'autres applications prenant en charge VRML. Le fichier VRML ressemblera exactement à la scène 3D, y compris toutes les caractéristiques et symboles, les décalages et les extrusions. Un plug-in VRML gratuit pour Netscape et Microsoft Internet Explorer est disponible chez Cosmo Software. Plusieurs autres visualiseurs VRML sont disponibles, tels que GLView.

La symbologie de la scène VRML ne peut pas être modifiée dans le navigateur Web, mais les utilisateurs n'ont pas besoin d'ArcGIS ou d'un autre logiciel spécialisé pour afficher et naviguer dans le modèle de surface.

Les fichiers VRML peuvent prendre beaucoup de temps à générer et les fichiers peuvent être très volumineux. Les TIN sont généralement beaucoup plus compacts que les grilles lors de l'exportation vers VRML. Voici une image d'une scène 3D dans ArcGIS.

Et la même scène exportée sous forme de fichier VRML affiché dans un navigateur Web :

Vous pouvez créer des documents ArcScene ou ArcGlobe, puis les publier sur des pages Web. Il existe un certain nombre de visionneuses VRML gratuites.


SIG à mi-parcours

(3) Petite échelle : 1 : 100 000 000 (le monde entier)
pourrait montrer toute l'Asie avec moins de détails spécifiques.
les images paraissent plus petites avec moins de détails.
une carte d'un continent serait une carte à petite échelle
Les pieds ou les mètres sont l'unité cartographique pour la plupart des applications SIG du gouvernement local

(1) Points : les informations d'attribut seraient la hauteur du poteau, la dernière date à laquelle le poteau a été entretenu n'ont pas de dimension
utilisé pour définir les emplacements de petits objets tels que des puits, des bâtiments ou des étangs. (pages 33-34)
utilise une seule paire de coordonnées pour représenter l'emplacement d'une entité considérée comme n'ayant aucune dimension. Seul l'emplacement est important.
(2)Lignes/entités linéaires/ou arcs : a 2 paires de coordonnées, rivières EX, routes, limites (p. 33-34)
Une longue ligne droite peut être représentée par deux paires de coordonnées alors qu'une ligne courbe peut être représentée par trois coordonnées ou plus.
Point de départ et point d'arrivée (appelés nœuds)
les points intermédiaires d'une ligne sont appelés sommets.
Les attributs peuvent être attachés aux lignes entières, aux segments, aux nœuds ou aux sommets.
(3) Les polygones (entités surfaciques) sont définis par un ensemble fermé de lignes formées par un ensemble de lignes connectées, soit une ligne avec un point de fin qui se connecte au point de départ, soit un ensemble de lignes connectées de début à finir. (p. 36) Parcs EX, plaques d'égout, données d'attributs telles que la zone, le comté et peuvent inclure d'autres polygones.
peuvent être adjacents à d'autres polygones et ainsi partager des lignes de bord avec d'autres polygones.

(1) Intervalle égal : EX 100, 200, 300, etc., n'est pas adapté aux distributions de données à pics élevés ou asymétriques.

(2) Ruptures naturelles (Jenks) : les points de rupture utilisant une méthode de clustering, les plages de classes sont spécifiques à l'ensemble de données individuel, il est donc difficile de comparer les cartes (au lieu de 0,14 - 15,05, arrondissez et faites 0,14 % - 20 %). Les classes sont basées sur des regroupements naturels inhérents aux données.

(3) Manuel : vous permet de saisir manuellement ce que vous voulez pour vos intervalles

(5) L'intervalle défini vous permet de spécifier une taille d'intervalle utilisée pour définir une série de classes avec la même plage de valeurs. Par exemple, chaque intervalle s'étendra sur 75 unités

(6) Quantile : Chaque classe contient un nombre égal de caractéristiques. Différent de l'intervalle égal.

(7) Géométrique : crée des interruptions de classe basées sur des intervalles de classe qui ont une série géométrique

Les sorties:
Raster :
(1) .bmp, .tif, .jpg
(2) taille/résolution
(3) options de couleur

Vecteur:
(1) .emf (métafichier amélioré)
(2) .eps (postscript encapsulé)
(3) . ai (illustrateur Adobe)
(4) .lpk (paquets de couches)

(Récit condensé : dans le panneau de droite d'ArcCatalog, cliquez avec le bouton droit sur la géodatabase dans laquelle vous souhaitez importer le fichier de formes. Sélectionnez Importer. Classe d'entités (unique). Dans la boîte de dialogue Classe d'entités vers géodatabase (unique), cliquez sur le bouton Parcourir juste pour à droite du champ Entités en entrée. Accédez au fichier de formes que vous souhaitez importer. Cliquez sur Ajouter. Cliquez sur OK.)

Problèmes/problèmes courants avec les jointures :
(1) les types de champs sont différents (par exemple, l'un est numérique et l'autre est du texte)
EX. les valeurs de texte "left align" tandis que les valeurs numériques "right align" dans les colonnes EX STFID et STFID_NUM)
Solution : créez un nouveau champ de type numérique et utilisez la calculatrice de champ pour voler les données de l'ancienne valeur de texte, nommez votre nouveau champ puis choisissez le type que vous souhaitez, cliquez sur ok. Sous les champs, double-cliquez sur l'ancien champ qui avait une valeur de texte sous les champs, double-cliquez sur l'ancien champ qui avait une valeur de texte sous Types, choisissez Nombre. cliquez sur OK
(2) Le format des données varie (il faut supprimer les " tirets " entre les chiffres) EX MAPNO 48-S-138 et LOTBLOCK 28G249

DE KYLE :
Principes fondamentaux du livre : dans la table des matières, cliquez avec le bouton droit sur une couche. Cliquez sur Jointures et relations. Choisissez le champ de cette couche sur lequel la jointure sera basée. Choisissez le champ de la table sur lequel baser la jointure. Cliquez sur OK.

EX de Numérique : en calculant le nombre d'Hispaniques 0-5 et 5-17 ans des colonnes de la table attributaire, le calculateur de champ vous donnerait un total : [AGE_UNDER5] + [AGE_5_17]

(2) Projection sécante : Une projection dont la surface coupe la surface d'un globe. Une projection sécante conique ou cylindrique, par exemple, est encastrée dans un globe

(2) Mercator transversal (cylindre, puis "déroulé" sur une surface plane) sont les types de projection les plus courants utilisés pour les données spatiales en Amérique du Nord et dans une grande partie du monde p 07

(3) L'azimut est l'angle formé entre une direction de référence (Nord) et une ligne allant de l'observateur à un point d'intérêt projeté sur le même plan que la direction de référence (un point de référence A, et deux autres points B et C)

(2) Système de coordonnées universel transversal de Mercator : le système de coordonnées (UTM) est un système de coordonnées global, basé sur la projection transversale de Mercator (utilisée aux États-Unis, en Amérique du Nord et dans d'autres pays113 Le système UTM divise la terre en zones de 6 degrés de large en longitude (numérotées de 1 à 60), et s'étendent
Mercator transverse universel : Le système de coordonnées Mercator transverse universel est une application spécialisée de la projection Mercator transverse. Le globe est divisé en 60 zones, chacune s'étendant sur six degrés de longitude. Ex. Lorsque vous entendez "UTM Zone 11", c'est là que l'on peut trouver l'Idaho.

(3) Mercator transversal : Ceci est similaire au Mercator sauf que le cylindre est tangent le long d'un méridien au lieu de l'équateur. Le résultat est une projection conforme qui ne maintient pas les vraies directions.

(4) Projections continentales et mondiales : inclure des variantes des projections de Mercator, Goode, Mollweide et Miller, entre autres p115
Aux États-Unis, les systèmes de coordonnées UTM et State Plane définissent un ensemble standard de projections cartographiques largement utilisées p123

(5) Système de coordonnées géographiques : le système de coordonnées géographiques n'est pas une projection cartographique. La terre est modélisée comme une sphère ou un sphéroïde.


Modification des détails dans Illustrator

Si vous souhaitez utiliser illustrator pour sélectionner et modifier la géométrie provenant d'ArcMap, vous constaterez que bon nombre des couches individuelles provenant d'ArcMap sont en fait des groupes de couches qui utilisent des masques d'écrêtage. Ceux-ci peuvent être libérés en sélectionnant la géométrie en question et en choisissant Objet-> Libérer le masque d'écrêtage. De plus, une partie de la géométrie avec une couche peut être regroupée. Et bien sûr, ceux-ci peuvent être dégroupés en cliquant avec le bouton droit sur le calque et en choisissant Dégrouper.


Option 2 : Exportation d'objets vers KML/KMZ via l'assistant

La deuxième option que vous pouvez choisir dans ‘Spatial Manager’ est l'utilisation du Assistant d'exportation général (commande SPMEXPORT dans les versions CAO), qui fonctionne de manière similaire pour exporter des objets de votre dessin/carte vers n'importe quel fichier spatial ou conteneur de données (Shapefile, SDF, GML, PostGIS, SQL Server, etc., en plus de KML/KMZ )

Cette fonction n'est pas aussi directe et intuitive que la précédente mais elle vous permettra de affiner beaucoup plus les données et propriétés à exporter et l'apparence des objets exportés dans Google Earth. Vous pouvez choisir les données à exporter à partir des tableaux de données, ajouter des propriétés de dessin de données alphanumériques (dans les versions CAO, les détails des hachures, les longueurs ou les zones, etc.) et définir les couleurs de ligne et de remplissage ainsi que les épaisseurs de ligne ou les transparences dans le KML /Fichier KMZ

Vous pouvez aussi définir une structure de dossier pour les objets sortants (même en les ajoutant à un fichier KML/KMZ existant), ou sélectionnez comment considérer le Position Z et altitude des objets résultants dans Google Earth

Les remarques sur les potentiels des fichiers Google Earth et KML/KMZ qui sont mentionnées dans la première option, sont également valables ici

Dans la deuxième section des vidéos ci-dessus, vous verrez comment cet assistant fonctionne lors de l'exportation vers des fichiers KML/KMZ. Merci de les jouer

‘Spatial Manager’ Exportation (vers KML/KMZ) fonction et paramètres – Fichiers KML/KMZ structurés

Ouvrir le KML/KMZ exporté résultant dans Google Earth

En savoir plus sur la publication vers Google Earth et/ou l'exportation vers KML/KMZ dans ‘Spatial Manager’ :


Possibilités d'accès

Acheter un seul article

Accès instantané à l'intégralité de l'article PDF.

Le calcul des taxes sera finalisé lors du paiement.

Abonnez-vous à la revue

Accès en ligne immédiat à tous les numéros à partir de 2019. L'abonnement se renouvellera automatiquement chaque année.

Le calcul des taxes sera finalisé lors du paiement.

La source: U.S. Energy Information Administration, 2019

La source: Laboratoire national des énergies renouvelables, 2019

La source: Laboratoire national des énergies renouvelables, 2019


Recherche de polygones voisins - requête postgis

Le problème : j'ai une table pleine de polygones (districts d'un pays). Certains de ces polygones ont un certain attribut qui peut ou non être vérifié. Dans ce cas, l'attribut est appelé "spread" et la valeur "cochée" est 1.

Maintenant, j'aimerais exécuter une requête qui trouve tous les polygones "cochés". Quelque chose comme:

puis j'aimerais définir l'attribut "spread" sur "1" dans chaque polygone voisin qui n'est pas déjà coché. (Je veux aussi définir un deuxième attribut, mais ce n'est qu'un ajout mineur)

Tout d'abord, commençons par une requête qui sélectionne tous les polygones voisins d'un polygone avec une valeur de propagation = 1

Cette requête renvoie tous les polygones voisins d'un polygone avec spread = 1

Maintenant, je veux mettre à jour la table en fonction des résultats de cette sous-requête. Ceci est réalisé avec ce morceau de code fourni par John Powell aka Barca (voir la réponse et aussi les commentaires ci-dessous):

Exécutez cette requête et elle définira les attributs répartis sur 1 et l'heure sur 2 des polygones adjacents sans toucher les polygones d'origine avec propagation = 1. Par conséquent, cela pose la réponse parfaite à ma question.


Exporter les résultats des comparaisons Select By Location de chaque ligne d'un fichier de formes avec des polygones d'un autre ? - Systèmes d'information géographique

Plus de traitement des données spatiales

n Dissoudre des entités en fonction d'un attribut.

n Découpez un calque en fonction d'un autre.

Exemple d'analyse spatiale - Localisation d'un centre commercial

n Critères d'un « bon » emplacement :

n À moins de trois kilomètres d'une autoroute principale, mais pas à moins d'un quart de mille d'une autoroute principale.

n Au moins un mile d'un plan d'eau (par exemple, des rivières, des lacs).

n A proximité d'une zone résidentielle importante (au moins 2500 ménages dans un rayon de 5 miles).

n Nécessite un grand terrain vacant (au moins 10 acres)

n Analyse classique de « conformité du site »

n Décomposez la tâche en plusieurs étapes d'analyse ArcGIS

n Couche autoroute (entité linéaire),

n Couche de terrain vacant (entité polygonale),

n Couche Rivers (entité linéaire ou surfacique),

n Couche des lacs (entité polygonale),

n Couche d'occupation des sols (entité polygonale),

n Couche de zone inondable (entité polygonale),

n Couche de recensement (niveau de groupe de blocs, entité surfacique).

1. Sélectionnez des parcelles vacantes de plus de 10 acres, résultat : vacantlg10.

2. Couche d'autoroute tampon avec des tampons de 2 miles et 0,25 mile, résultat : hw2 et hwquartier .

3. Couches tampons Rivers and Lakes avec des tampons de 1 mile, résultat : rivière1 et lac1.

4. Intersecter vacantlg10 avec hw2 = vacanthw1, terrain utile.

5. Union hwquarter avec river1, lake1 et flood (il faut faire deux couches à la fois), résultat : non utilisable , (terrain ne répondant pas aux critères)

6. supprimer notusable de vacanthw1, résultat : utilisable.

7. intersection utilisable avec la couche de groupes de blocs de recensement (touchez la limite de), résultat : bloc utilisable .

8. calculez et resélectionnez le bloc utilisable dont le nombre total de ménages est supérieur à 2500, résultat : grosbloc utilisable.

n Considérez un critère supplémentaire : « Le revenu médian des ménages de la zone résidentielle proche devrait être supérieur à 30 000 $ » ?

n Affiner le critère 'zone résidentielle majeure'

Exemple ArcGIS du devoir n° 2 : recherchez quels groupes de blocs se trouvent dans quelle ville