Eu tenho a seguinte camada usando SRID 27700 no postgis:
É todas as regiões administrativas do Reino Unido e (como você pode ver no agrupamento de cores) cada uma delas possui um campo de texto especificando o município em que se encontram.
O que eu gostaria de fazer é criar polígonos maiores do condado dos menores em um determinado condado, de modo que, na figura acima, todos os polígonos de cor verde-azulado formariam um polígono grande a partir do anel externo que contém todos os polis do cor, como todos os roxos, marrons, rosa, cinza etc, todos devem formar um polígono.
Eu já tentei o seguinte:
insert into parishesmerged (geometry)
select astext(multi(ST_Union(the_geom))) as the_geom from parishes
group by county_name
Mas ele continua gerando geometrias quebradas, que então eu tenho grandes problemas para processar ainda mais.
Estou tentando fazer um mapa mais simples no nível do condado com as principais áreas de produção de.
Quaisquer soluções não precisam estar no Postgis também, eu tenho a pilha OS4Geo completa instalada, a versão mais recente do QGis e mais utilitários do que posso usar.
As únicas coisas que não tenho são os meninos grandes como ArcGis (embora eu possa ter um Old Mapinfo por aí em algum lugar)
Para o registro, o conjunto de dados que estou tentando criar é acompanhar um livro GIS que estou escrevendo, destinado a programadores .NET que desejam escrever aplicativos GIS usando .NET
Depois de tentar as sugestões abaixo, a que melhor funcionou foi a solução 'Paul Ramseys'.
Agora tenho um bom arquivo simplificado de municípios e distritos, que é simples o suficiente para o meu livro, mas complexo o suficiente para permitir que eu demonstre algumas SQL geoespaciais interessantes.
Embora a solução de Paul tenha sido a que funcionou para mim, eu também recorri às outras respostas para simplificar o mapa de polígonos e reduzir ainda mais a complexidade.
Porém, quando observei isso, enquanto ST_Collect é realmente mais rápido que ST_Union, corrida por corrida também foi a principal responsável por geometrias quebradas. Meu palpite é que o aumento de velocidade é obtido à custa de menos precisão na função principal.
astext(multi())
peça? Estou saindo do que vejo em outros exemplos de dissolução do PostGIS.