Mais

GToca erro sem erro de geometria?


Tenho um mapa que mostra os limites administrativos (disponível aqui). Eu carreguei no R e tentei executartoques(de rgeos) para determinar quais áreas são adjacentes a quais, mas receber um erro:

biblioteca (rgeos) library (rgdal)> Districts <- readOGR (…) fonte de dados OGR com driver: ESRI Shapefile com 9235 recursos e 16 campos Tipo de recurso: wkbPolygon com 2 dimensões> adj <- gTouches (Districts) Erro em RGEOSBinPredFunc (spgeom1 , spgeom2, byid, "rgeos_touches"): TopologyException: conflito de localização lateral em 34.353274826787519 -15.739538287828307

Depois de pesquisar esse erro no Google, encontrei muitos conselhos sobre "geometrias ruins". Em seguida, peguei o arquivo no QGIS (seguindo o conselho aqui do documento "Bad, bad polygon" do Faunalia) e no ArcMap (usando o código Python aqui) e executei as ferramentas de verificação / reparo de geometria. Nenhum programa encontrou um erro.

Aumentei meu mapa nas coordenadas ofensivas, mas não vejo nada que pareça errado - é um local onde duas linhas se juntam perto de outro vértice.

Como faço para reparar a geometria, ignorar essa região em particular ou alguma outra função R que pode funcionar melhor do quetoques?


Tentar:

adj <- gTouches (distritos, byid = TRUE)

Isso produziu uma matriz indicando quais regiões compartilham uma fronteira com outras regiões.


Assista o vídeo: Wilco en de leugendetector. Alles Kids. Afl. 3 (Outubro 2021).