Shapely é um pacote Python de código aberto para construção e análise de objetos geométricos planares 2D comuns aos sistemas de informações geográficas.
Estou procurando algumas sugestões sobre como tornar meu código python mais eficiente. Normalmente, a eficiência não importa para mim, mas agora estou trabalhando com um arquivo de texto de locais nos EUA com mais de 1,5 milhão de pontos. Com a configuração fornecida, leva cerca de 5 segundos para...
Existe uma maneira fácil de transformar objetos Shapely (polígonos e poligonais) de uma projeção para outra sem precisar cavar e extrair coordenadas manualmente? Na verdade, eu nem me importo se eles são objetos Shapely neste momento, só quero passar recursos e uma projeção e recuperar um...
Estou tentando executar o Shapely em uma nova máquina e tentei instalá-lo da fonte, usando o instalador em sua página da Web e usando pip install shapely. Toda referência que tento fazer a partir da biblioteca, exceto from shapely import *lança o erro WindowsError: [Error 126] The specified module...
Alguém pode demonstrar uma maneira simples de escrever estruturas de dados de geometria de bem torneados em arquivos de forma? Estou particularmente interessado em polígonos com furos e cordas de linha. Também seria benéfico ficar longe do arcpy (então osgeo, pyshp etc. seriam...
Estou procurando calcular estatísticas focais para cada célula de uma varredura, dentro de uma vizinhança de um critério especificado. Antecedentes - Eu tenho três rasters binários, cada um representando um único tipo de vegetação de interesse. Eu gostaria de calcular a porcentagem de cobertura...
fundo De um ponto conhecido, eu preciso estabelecer o "perímetro visível" ao redor mais próximo de uma tabela de MultiLineStrings, conforme mostrado no diagrama. Pesquisei este site com vários termos (por exemplo, borda mínima, perímetro mínimo, vizinho mais próximo, clipe, contendo polígono,...
Dado que geopandas GeoDataFramecontém uma série de polígonos, gostaria de obter a área em km² de cada elemento da minha lista. Este é um problema bastante comum, e a solução sugerida no passado foi usar shapelye pyprojdiretamente (por exemplo, aqui e aqui ). Existe uma maneira de fazer isso de...
Estou tentando capturar todos os polígonos não sobrepostos indicados abaixo usando Shapely (dados polígonos A, B e C). Além disso, espero fazê-lo sem iteração, teste de interseção etc. A resposta aceita a essa pergunta expressa o método PostGIS, mas parece que 'união' significa coisas diferentes...
Eu tenho tentado criar uma função que faz basicamente a mesma coisa que o QGIS "dissolve" a função. Eu pensei que seria super fácil, mas aparentemente não. Então, pelo que entendi, o uso de fiona com bem torneado deve ser a melhor opção aqui. Eu comecei a mexer com arquivos vetoriais, então esse...
Estou tentando encaixar duas linhas entre si usando Shapely / Geopandas, mas o resultado do encaixe é muito estranho. Eu tentei : import geopandas as gpd from shapely.geometry import * from shapely.ops import snap lines1 = gpd.GeoDataFrame.from_file('lines1.shp') lines1 = lines1.to_crs({'init':...
Eu tenho polígonos de coordenadas em (python shapely) que se parece com isso POLYGON ((24.8085317 46.8512821, 24.7986952 46.8574619, 24.8088238 46.8664741, 24.8155239 46.8576335, 24.8085317 46.8512821)) Gostaria de calcular a área deste polígono em km ^ 2. Qual seria a melhor maneira de fazer...
parece uma pergunta boba, mas não consigo encontrar uma resposta clara sobre isso: que unidades geopandas / bem torneadas usam o cálculo da distância / área entre os
Eu tentei vários exemplos de código usando bibliotecas como shapefile, fiona e ogr para tentar verificar se um ponto (x, y) está dentro dos limites de um multipolígono criado com o ArcMap (e, portanto, no formato shapefile). No entanto, nenhum dos exemplos funciona bem com multipolígonos, embora...
Estou usando Shapely em python e recebi MultiLineStringum monte de Linestringobjetos. Posso garantir que todos os LineStringobjetos são linhas simples com apenas 2 vértices e que todos fazem parte de uma única linha (sem ramificações). Eu quero "conectar os pontos" e criar um único LineString....
Estou procurando uma solução python de código aberto para converter raster em polígono (sem ArcPy). Eu conhecia a função GDAL para converter raster em polígono, e aqui está o manual: http://pcjericks.github.io/py-gdalogr-cookbook/raster_layers.html#polygonize-a-raster-band No entanto, espero que...
Estou tendo problemas para entender o uso de índices espaciais com o RTree. Exemplo: tenho 300 pontos de buffer e preciso conhecer a área de interseção de cada buffer com um shapefile de polígono. O arquivo de forma do polígono possui> 20.000 polígonos. Foi sugerido o uso de índices espaciais...
Eu tenho geopandas GeoDataFrame contendo centenas de formas Polygone MultiPolygonformas geométricas. Os polígonos se sobrepõem em muitos lugares. Gostaria de fazer uma nova geometria que contenha quantas delas se sobrepõem. Algo assim: Alguém tem alguma idéia de como abordar isso? Eu nem consigo...
Ao tentar instalar o Shapely no meu computador com Windows de 64 bits, não consigo fazer com que a biblioteca GEOS funcione. Até agora, eu executei o instalador OSGeo4W a partir do qual instalei o GDAL (acredito que a biblioteca geográfica esteja incluída nesse pacote). Depois disso, verifiquei e...
Estou lendo um shapefile usando Fiona: with fiona.open('data/boroughs/boroughs_n.shp') as source: mpolys = [p for p in source] candidate = polys[0]['geometry'] isso fornece um ditado contendo uma lista de listas de coordenadas e o tipo 'MultiPolygon'. No entanto, não sei como usar os dados de...
[EDIT: a solução para isso era simplesmente usar OGR para ler shapefiles. Veja o exemplo da geographika.] Em um shapefile da ESRI, não há distinção entre polígonos e multipolígonos. Além disso, não há distinção explícita entre orifícios internos e anéis externos (além da "mão" de um determinado...