Gostaria de obter a interseção de vários polígonos. Usando o shapelypacote do Python , posso encontrar a interseção de dois polígonos usando a intersectionfunção Existe uma função eficiente semelhante para obter a interseção de vários polígonos? Aqui está um trecho de código para entender o que...
Estou tentando descobrir onde várias linhas cruzam um polígono para dois geodataframes diferentes: 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 =...
Eu tenho um quadro de dados geopandas de um monte de cadeias de linhas que têm alguns dados associados a cada vértice / ponto: Point_x = (Lat, Lon, Time, ID, Data1, Data2, Data3) Os pontos são convertidos em cadeias de linhas com base no ID e ordenados por Hora. Eu quero quebrar as cadeias de...
Estou fazendo um fluxo de trabalho em pequena escala no notebook ipython com geopandas e bem torneadas, obtendo um monte de dados geoespaciais, às vezes de shapefiles, às vezes de postgis (onde o processamento mais caro é realizado). Agora, puxo as tabelas postgis para python usando sqlalchemy,...
Se tiver dois pontos, dos quais quero criar um LineStringobjeto direto : from shapely.geometry import Point, LineString A = Point(0,0) B = Point(1,1) O manual Shapely paraLineString estados: Uma sequência de Pointinstâncias não é um parâmetro de construtor válido. A LineStringé descrito por...
Estou fazendo um cálculo muito simples do comprimento de uma polilinha usando o formato: 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) Minhas coordenadas estão em WGS84. Não consigo encontrar...
Estou tentando criar um retângulo girado em coordenadas lat / lon. Basicamente, tenho as coordenadas do centro e o ângulo de rotação ao redor do centro. É possível usar usando, bem torneado ou gdal, ou algo
Para obter as coordenadas externas, preciso converter um MultiPolygon bem torneado em um polígono. Eu faço assim: 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 =...
Desejo cortar um conjunto de polilinhas (linhas pretas na imagem abaixo) no limite externo de um polígono. Quaisquer vazios dentro do polígono devem ser ignorados. Minha saída ideal são as linhas amarelas tracejadas. As linhas iniciais podem ou não ser retas. A imagem é um exemplo simplificado, na...
Estou tentando gerar polígonos para faixas orbitais de satélite. Até agora, tenho um método para gerar duas linhas que representam a borda de cada faixa em [lat, long]. Alguns dos trechos cruzam a linha de dados internacional e, por isso, envolvem-se: Consegui resolver isso com ogr2ogr...
Sei que esse problema foi resolvido antes com um software diferente, mas como faço com o Shapely? Eu tenho um polígono com uma lacuna muito pequena no interior da forma: Aqui está um close-up: O polígono é válido, portanto não posso usar o truque de buffer, e simplificar o polígono também...
Eu tenho algum código que estou usando para determinar quais polígonos esculturais / MultiPolygons se cruzam com um número de linhas de formas bem torneadas. Através das respostas a esta pergunta, o código passou disso: import fiona from shapely.geometry import LineString, Polygon, MultiPolygon,...
Eu tenho dois polígonos POLÍGONO ((- 5 -5, -5 5,5 5,5 -5, -5 -5)) POLÍGONO ((1 1, 1 7, 7 7, 7 1, 1 1)) e deseja dividi-los uniformemente. O resultado deve ser uma divisão ao longo da linha marcada vermelha. Estou usando o qgis / python shapely e desejo fazê-lo usando
Estou executando uma análise específica em que usamos shapely para criar buffers em torno de pontos (locais da loja) e, em seguida, verificamos se outros pontos (locais do usuário) estão presentes nesse valor do buffer. Quando eu verifico a distância usando formas, ele acaba sendo diferente da...
Estou trabalhando com Python, Shapely e Fiona ... Existe alguma maneira de fornecer uma entrada de um shapefile de linha e um shapefile de polígono e obter um resultado dos pontos de interseção e suas coordenadas? A ilustração fornece uma descrição precisa
Eu tenho um quadro de dados de coordenadas X, Y que representam pontos ao longo dos caminhos percorridos por várias entidades diferentes. Pseudo-dados aqui, mas é aproximadamente da forma: entity_id lat lon time 1001 34.5 14.2 4:55 pm 1001 34.7 14.5 4:58 pm 1001 35.0 14.6 5.03 pm 1002 27.1 19.2...
Estou tentando obter um arquivo KML ( wijken.kml ) no Shapely . O arquivo KML é validado contra o esquema XML adequado, então acho que a entrada está correta. Rotas que tentei: 1) Convertendo para o formato WKT ou WKB e lendo com funções incorporadas Conversão: ogr2ogr -f CSV wijken.csv...
Como faço para desenhar uma elipse