Je souhaite générer automatiquement des ID pendant que je crée de nouvelles fonctionnalités. Ainsi, chaque fois que je crée un nouveau point, son ID augmentera de 1 comme suit: Feature 1: id 1 Feature 2: id 2 Feature 3: id 3 Comment puis-je attaquer
Chaque fois que je crée un script pour manipuler un vecteur (point, ligne et polygone) en python, je rencontre ce doute. Quelle est la différence entre la boîte englobante (bbox), l'enveloppe, l'étendue et les
J'ai un projet QGIS qui contient un compositeur. Bien sûr, si je vais dans l'interface graphique, je peux l'exporter sous forme de fichier pdf, etc. De la lecture en ligne, voici ma tentative: #!/usr/bin/env python import sys from qgis.gui import QgsMapCanvas, QgsLayerTreeMapCanvasBridge from...
Je vois parfois srsName="EPSG:4326"dans des documents GML (par exemple https://postgis.net/docs/ST_GeomFromGML.html , bien que ma source particulière ne soit pas publique). L'intention (sauf peut-être pour l'ordre des axes qui n'est pas important ici) est claire et probablement interopérable sur...
Je travaille sur un projet qui nécessite un tampon de 30 m autour d'une entité linéaire. La caractéristique représente la côte et j'ai besoin du tampon pour montrer 30 m à l'intérieur des terres de la ligne représentant la côte. Existe-t-il un moyen de le faire en utilisant...
Comment afficher les étiquettes des points geojson dans une carte Leaflet? Il y a Leaflet.label qui est maintenant déconseillé en faveur de L.Tooltip , mais qui n'affiche que du texte au survol. Je souhaite afficher les étiquettes de texte directement sur la carte sans avoir besoin d'interaction...
J'essaie d'écrire un script de traitement en utilisant l' r.neighborsalgorithme Grass pour les rasters. Lorsque vous utilisez l'interface graphique pour cet algorithme, vous pouvez laisser le champ d'extension de couche GRASS blanc pour utiliser l'étendue minimale. Au contraire, lorsqu'il est...
J'essaie d'écrire un script qui enregistrera un rendu de plusieurs couches à l'aide du compositeur de carte. Le problème que je rencontre est que le script enregistre avant que qgis ait fini de rendre toutes les couches. Sur la base de plusieurs autres réponses ( 1 , 2 , 3 ), j'ai essayé d'utiliser...
Je dois effectuer un affûtage de Gram-Schmidt , je n'ai pas accès à ArcGIS ou ENVI, donc j'utilise QGIS . Existe-t-il un plug-in ou une capacité similaire disponible pour effectuer un affûtage de Gram-Schmidt dans QGIS? Sinon, est-ce possible dans tout autre logiciel open
J'essaie de trouver où plusieurs lignes coupent un polygone pour deux géodonnées différentes: from shapely.geometry import Polygon, LineString import geopandas as gpd polygon = Polygon([(0, 0), (1, 0), (1, 1), (0, 1), (0, 0)]) line1 = LineString([(0.5, 0.5), (0.7, 0.7)]) line2 = LineString([(0.9,...
J'ai un ensemble de données avec plus de 2000 fonctionnalités basées sur des adresses. J'ai besoin de zoomer sur chaque point individuellement et de les déplacer un peu (pas de décalage cohérent donc pas d'autre moyen de le faire). La sélection de l'entité suivante dans la table attributaire et le...
J'écris une fonction qui a une géométrie comme paramètre qui peut être dans différentes projections: CREATE OR REPLACE FUNCTION foo(in geometry) RETURNS boolean AS $$ DECLARE transformed_geom geometry := in; BEGIN IF ST_SRID(transformed_geom) != 32737 THEN transformed_geom :=...
J'ai une fonction qui crée des panneaux solaires photovoltaïques représentés comme des polygones. Essentiellement, il crée une grille rectangulaire où l'utilisateur peut spécifier les paramètres suivants: Longueur Largeur Distance horizontale Distance verticale Le code est basé sur le plugin...
J'ai besoin de créer une grille de secteur rayonnant à partir d'un point d'intérêt. Voici un exemple que j'ai laborieusement assemblé à la main: J'ai trouvé quelques solutions similaires dans Comment générer des fichiers de formes de fléchettes à l'aide d'ArcGIS for Desktop? et Draw rose des vents...
J'ai deux fichiers vectoriels dans le panneau des couches dans QGIS (V2.8.2, Wien). J'essaie de sélectionner des points dans les deux fichiers vectoriels simultanément pour les afficher dans les deux tables d'attributs - les deux tables sont définies pour afficher uniquement les entités...
Je cherche un moyen de fusionner les fichiers de bande Sentinel 2 jp2 ( B02, B03, B04 ) et de corriger les couleurs RVB. Tout doit être fait avec un script bash ou python. Pour mon exemple, je travaille sur ces images . Idéalement, la solution sera proche de ce tutoriel. Je peux fusionner les...
Je travaille avec des fichiers de formes de recensement (multipolygones) et j'ai un autre tableau avec des points. Nous enregistrons dans quel bloc de recensement chaque point se trouve. J'ai un point qui renvoie faux pour toutes les requêtes ST_Within et ST_Touches. Il retourne cependant vrai pour...
Je veux créer un calque temporaire à partir d'un calque existant et n'obtenir que les fonctionnalités du calque existant, qui sont sélectionnées. Je sais qu'ArcGIS possède cette fonctionnalité que vous pouvez créer une couche à partir d'entités sélectionnées. Cependant, je n'ai trouvé aucune...
J'essaie d'ouvrir un fichier de formes avec le bloc-notes Jupyter, mais je n'arrive pas à l'ouvrir. Le fichier de formes se trouve dans le même dossier que mon classeur et Python ne peut pas l'identifier même si j'utilise le chemin complet. import shapefile as shp import matplotlib.pyplot as plt sf...
J'ai un cadre de données géopandas dans EPSG: 4326 et je ferais un nouveau cadre de données composé de toutes les lignes qui se trouvent dans une certaine boîte englobante. D'abord, j'obtiens le cadre de délimitation qui m'importe (qui est en fait le cadre de délimitation d'une autre trame de...