Consertando milhares de pequenos polígonos sobrepostos?

Estou usando o ArcMap 10.3.

Eu tenho uma camada de diferentes tipos de uso da terra. Descobri que muitos polígonos se sobrepõem e, como resultado, o mesmo polígono tem dois tipos diferentes de uso da terra.

Usando a resposta de @ radouxju , construí uma topologia para a camada de uso da terra (com a regra de não se sobrepor). Em seguida, usei a mesclagem para corrigir os polígonos grandes sobrepostos.

No entanto, tenho milhares de pequenos polígonos sobrepostos. Eu posso consertar todos eles manualmente, o que levará tempo.

Existe uma maneira mais rápida de corrigir todos esses erros?

Atualizei o mapa para mostrar todos os tipos de uso da terra. Polígonos estão se sobrepondo; não lascas.

A preocupação levantada por @Adam Cara é muito válida. No entanto, quando olhei para a auto-interseção, parecia que a maioria deles é um pedaço fino na fronteira entre dois polígonos maciços ou pedaços / ilhas desconectados de um grande polígono sentado dentro de outro. Nesse caso, apagar sobreposições e atribuir sobreposições dissolvidas podem funcionar como atributo do vizinho mais próximo (o que ocorrer primeiro).

Processo:

arcpy.Intersect_analysis("Landuse2b #","D:/Scratch/mpart.shp","ALL","#","INPUT")
arcpy.MultipartToSinglepart_management("mpart","D:/Scratch/spart.shp")
arcpy.Dissolve_management("spart","D:/Scratch/dissolved_1.shp","#","#","SINGLE_PART","DISSOLVE_LINES")
arcpy.Erase_analysis("Landuse2b","dissolved_1","D:/Scratch/erased.shp","#")
arcpy.DeleteField_management("spart","FID_Landus;OBJECTID;Lnduse2013;Shape_Leng;Shape_Area")
arcpy.SpatialJoin_analysis("spart","erased","D:/Scratch/sj_1.shp","JOIN_ONE_TO_ONE","KEEP_ALL"…..)
arcpy.Merge_management("sj;erased","D:/Scratch/merge_1.shp"….)
arcpy.Dissolve_management("merge_1","D:/Scratch/dissolved.shp","Lnduse2013","#","SINGLE_PART","DISSOLVE_LINES")

Teste:

arcpy.Intersect_analysis("dissolved #","D:/Scratch/test.shp","ALL","#","INPUT")

Mensagem:

• Executando: Intersecione "dissolveu" D: \ Scratch \ test.shp ALL # INPUT
• Hora de início: Sex. 22 de abr. 10:16:07 2016
• Recursos de leitura ...
• Recursos de craqueamento ...
• Recursos de montagem ...
• AVISO 000117: Aviso de saída vazia gerada.
• Êxito em 22 de abril - 10:16:09 2016 (Tempo decorrido: 2,34 segundos)

Observe que esta mensagem significa que não há interseções automáticas na saída do Dissolve final

Resultado:

https://drive.google.com/file/d/0B6B33FrGNitPT0FvbVdHQm0xMDg/view?usp=sharing

Se seus polígonos são lascas, o comando eliminar funciona bem para mesclá-las no polígono da área maior ou na borda mais longa.

Se os polígonos forem sobrepostos, pode haver uma maneira mais fácil, mas eu selecionaria as sobreposições em uma camada separada e as uniria novamente, criando as lascas e usando o comando delete.

1. Selecione uma área pequena e converta em dois arquivos de forma para ver se esse processo está funcionando.
2. Crie um novo campo chamado teste (número) em s1 e s2.
3. Interseção de dois arquivos de forma S1 e S2.
4. As áreas sobrepostas serão identificadas em um arquivo separado, por exemplo, denominado Common1 (Área hachurada na figura).
5. Calcule o teste de campo como 1 em Comum.
6. União S1 e Comum
7. Selecione 1 no campo de teste e exclua a área hachurada.
8. Agora Mesclar s2 e comum. O comando Apagar pode ser usado para apagar a área comum.

Você já pensou em dividir as classes individuais de uso da terra em arquivos de forma separados para trabalhar com eles dessa maneira. Primeiro pensei em fazer um cruzamento, mas vi que havia falhado acima. Se você tiver sete shapefiles separados, poderá manipulá-los mais facilmente.

A principal pergunta que tenho é se houver sobreposição, você tem uma hierarquia para qual classe vence e deve aparecer no final? Por exemplo: se o Recorte e o Laticínio se sobrepuserem, quais você mantém? Depois de estabelecer essa hierarquia, cruze duas classes por vez, começando pela mais baixa e terminando com a mais alta. Se sua saída estiver vazia em um determinado cruzamento, isso significa que eles não se cruzam.

Depois de criar uma topologia, você identificou os erros.

Na lista que mostra os erros, você pode selecionar todos os erros do mesmo tipo, clicar com o botão direito do mouse nessa lista e escolher uma "correção". No seu caso, a correção seria "criar recurso", que criará um novo polígono onde dois polígonos se sobrepõem. Esta é a primeira etapa (observe que você pode obter um resultado semelhante usando a ferramenta de união com a camada de uso da terra como entrada única e, em seguida, com várias peças em uma única peça (mas as sobreposições ainda se sobrepõem, é necessário excluir uma das o valor do atributo) ou converta seu polígono em linhas e depois em polígono (mas, se você perder o atributo, precisará de junção espacial para recuperar os atributos))

O segundo passo consistiria na remoção dos pequenos polígonos. O problema é que você realmente não tem uma regra para decidir qual valor de atributo você gostaria de manter. Minha ferramenta favorita é a ferramenta de eliminação , onde você pode escolher em qual polígono vizinho os pequenos polígonos serão mesclados (com base no maior limite compartilhado ou na maior área)

Última observação: se você possui uma regra para a classe de uso prioritário do solo, selecione iterativamente os pequenos polígonos (que não possuem rótulo) por localidade para atribuir cada classe de uso do solo por ordem decrescente de prioridade e use a calculadora de campo para atribuir seus rótulo.

Enquanto eu concordo com os outros sobre se são polígonos sobrepostos ou lascas, pelo que estou vendo aqui nesta foto, você está mostrando a todos. Eu recomendo que você altere suas cores de simbologia para uso no solo. Desta forma, você pode ver qual deles deve ser ou não o outro. A única cor que você nos mostrou em rosa é realmente difícil de ler e pode levar muito tempo para você.

Na segunda parte, eu selecionaria todos os polígonos de sobreposições ou lascas e os exportaria para um novo arquivo de forma. Dessa forma, você pode identificar o que eles são no uso da terra.

Uma solução possível, se você tiver o Spatial Analyst ou 3D Analyst, converter os polígonos de uso da terra em uma imagem rasterizada usando a ferramenta (Conversão de polígono para varredura) para garantir que nenhum polígono sobreposto saia, pois não é possível obter pixels sobrepostos. Em seguida, converta a varredura novamente em polígonos usando a ferramenta (Raster to Polygon), se você ainda precisar dela em formato vetorial. Nesse caso, escolha o tamanho da célula de alta resolução (tamanho detalhado da célula) para obter polígonos de boa qualidade, mas é necessário considerar a memória. Pode ser necessário tentar várias vezes alterando o tamanho da célula até ficar satisfeito com os resultados. Talvez seja necessário suavizar os polígonos de saída, se não for perfeito, mas isso depende se você tiver pelo menos Licença Padrão para ArcGIS Desktop. Sei que não é uma tarefa fácil, mas pode resolver seu problema.

Verifique as poucas folgas ou sobreposições entre polígonos para encontrar tolerância e use a ferramenta Integrar do Arctoolbox com os valores de tolerância adequados.