Estou usando arquivos de forma da história Limites estaduais dos EUA do NHGIS para fazer mapas em R. Esses arquivos de forma têm uma resolução muito maior do que eu preciso ao fazer mapas na escala do país: mais de 2 milhões de observações para um arquivo. Existe uma maneira de reduzir a resolução desses shapefiles no R? Se necessário, posso fazê-lo no QGIS.
qgis
shapefile
r
resolution
Lincoln Mullen
fonte
fonte
Respostas:
Você pode usar
gSimplify
orgeos
pacote e, se você adicionar otopologyPreserve=TRUE
sinalizador, ele preservará a topologia.Observe que você ainda pode acabar com linhas sobrepostas - precisamos da implementação desse robusto algoritmo DP em R:
http://www.sciencedirect.com/science/article/pii/S0098300413002380
[esse link possivelmente atrás de um paywall]
fonte
rgeos
documentos estão em cran.r-project.org/web/packages/rgeos/rgeos.pdf .Um método padrão para "afinar" formas lineares (incluindo limites de polígonos) é o algoritmo Douglas-Peucker . Pelo menos dois
R
pacotes implementam isso:dp
no pacote shapefiles ethinnedSpatialPoly
no pacote Guerry .Observe que o desbaste de polígonos adjacentes normalmente cria pequenas lacunas e sobreposições em seus limites: não respeita as relações topológicas entre eles. Não conheço nenhum
R
pacote que mantenha a topologia durante o desbaste. Para fins de mapeamento em pequena escala como esses, isso não deve ser um problema sério.Aliás, como ninguém sabe mais pronunciar alemão :-), Tom Poiker mudou a ortografia de seu nome há um tempo.
fonte