Shapely est un package Python open source pour la construction et l'analyse d'objets géométriques plans 2D communs aux systèmes d'information géographique.
cela ressemble à une question idiote, mais je ne trouve pas de réponse claire à cela: quelles unités géopandas / galbées utilisent le calcul de la distance / surface entre les objets?
J'ai essayé plusieurs exemples de code utilisant des bibliothèques telles que shapefile, fiona et ogr pour tenter de vérifier si un point (x, y) se situe dans les limites d'un multipolygone créé avec ArcMap (et donc au format shapefile). Cependant, aucun des exemples ne fonctionne bien avec les multipolygones, bien …
Lorsque j'essaie d'installer Shapely sur mon ordinateur Windows 64 bits, je ne parviens pas à faire fonctionner la bibliothèque GEOS. Jusqu'à présent, j'ai exécuté le programme d'installation OSGeo4W à partir duquel j'ai installé GDAL (je pense que la bibliothèque geos est incluse dans ce package). Après cela, j'ai vérifié et …
Je voudrais obtenir l'intersection de plusieurs polygones. En utilisant le shapelypackage de Python , je peux trouver l'intersection de deux polygones en utilisant la intersectionfonction. Existe-t-il une fonction efficace similaire pour obtenir l'intersection de plusieurs polygones? Voici un extrait de code pour comprendre ce que je veux dire: from shapely.geometry …
Je lis un fichier de formes avec Fiona: with fiona.open('data/boroughs/boroughs_n.shp') as source: mpolys = [p for p in source] candidate = polys[0]['geometry'] cela me donne un dict contenant une liste de listes de coordonnées, et tapez "MultiPolygon" Cependant, je ne sais pas comment utiliser les données de coordonnées pour créer …
[EDIT: la solution à cela était simplement d'utiliser OGR pour lire les fichiers de formes. Voir l'exemple de geographika.] Dans un fichier de formes ESRI, il n'y a pas de distinction entre Polygones et MultiPolygons. De plus, il n'y a pas de distinction explicite entre les trous intérieurs et les …
J'ai un cadre de données géopandas d'un tas de chaînes de lignes qui ont des données associées à chaque sommet / point: Point_x = (Lat, Lon, Time, ID, Data1, Data2, Data3) Les points sont convertis en chaînes de lignes en fonction de l'ID et classés par heure. Je veux briser …
Je fais un workflow à petite échelle dans un cahier ipython avec des géopandas et bien fait, en tirant un tas de données géospatiales, parfois à partir de fichiers de formes, parfois à partir de postgis (où un traitement plus coûteux est effectué). Maintenant, je tire ensuite les tables postgis …
Si j'ai deux points, à partir desquels je veux créer un LineStringobjet droit : from shapely.geometry import Point, LineString A = Point(0,0) B = Point(1,1) Le manuel Shapely pour lesLineString États: Une séquence d' Pointinstances n'est pas un paramètre constructeur valide. A LineStringest décrit par des points, mais n'est pas …
Je fais un calcul très simple de la longueur d'une polyligne en utilisant galbée: from shapely.geometry import LineString ... xy_list = [map(float,e) for e in xy_intm] line = LineString(xy_list) s = '%s,%s,%s' % (fr,to,line.length) Mes coordonnées sont en WGS84. Je n'arrive pas à trouver d'informations sur l'attribut de longueur galbée. …
J'essaie de créer un rectangle tourné en coordonnées lat / lon. Fondamentalement, j'ai les coordonnées du centre et l'angle de rotation autour du centre. Est-il possible de faire en utilisant, galbé ou gdal, ou autre chose?
J'essaie de générer des polygones pour des andains orbitaux satellites. Jusqu'à présent, j'ai une méthode pour générer deux lignes qui représentent le bord de chaque bande en [lat, long]. Certains des andains traversent la ligne de données internationale et s'enroulent donc: J'ai pu résoudre ce problème avec ogr2ogr -wrapdateline: ogr2ogr …
Afin d'obtenir les coordonnées extérieures, j'ai besoin de convertir un MultiPolygon bien fait en Polygone. Je le fais comme ça: if poly.geometry.type == 'Polygon': x, y = poly.geometry.exterior.xy elif poly.geometry.type == 'MultiPolygon': allparts = [p.buffer(0) for p in poly.geometry] poly.geometry = shapely.ops.cascaded_union(allparts) x, y = poly.geometry.exterior.xy # here happens the …
Je souhaite couper un ensemble de polylignes (lignes noires dans l'image ci-dessous) à la limite extérieure d'un polygone. Tout vide dans le polygone doit être ignoré. Ma sortie idéale sont les lignes jaunes en pointillés. Les lignes initiales peuvent être droites ou non. L'image est un exemple simplifié, en réalité …
We use cookies and other tracking technologies to improve your browsing experience on our website,
to show you personalized content and targeted ads, to analyze our website traffic,
and to understand where our visitors are coming from.
By continuing, you consent to our use of cookies and other tracking technologies and
affirm you're at least 16 years old or have consent from a parent or guardian.