Eliminación de polígonos (definidos) en astilla mediante PostGIS?

Question

¿Cómo escribo SQL (para PostGIS en PostgreSQL) para eliminar polígonos de astillas (incluyendo astillas largas y estrechas y algunos pequeños polígonos irregulares que eran "huecos" erróneos en el conjunto de datos original) de mi conjunto de datos fusionándolos con un polígono adyacente con el borde más largo?

Tengo un conjunto de datos (tablename: "Vegetation") que contiene varios miles de polígonos identificados como "sliver" o "real" en una columna de descripción (nombre de columna: "desc") como "sliver" o "main". Los polígonos no se superponen a otros polígonos y son mayores que la tolerancia que puedo aceptar en las coordenadas.

[Editar: Antes de importar mi conjunto de datos a PostGIS, lo limpié importándolo a GRASS (así es como convertí todos los huecos en pequeños polígonos, que posteriormente fusioné con el conjunto de datos principal). Descubrí que si aumentaba la tolerancia de ajuste y el área mínima más allá de '1' en la importación a GRASS, empezaba a perder las características deseadas].

Entiendo que puedo necesitar los siguientes pasos:

1. Identifique el polígono (Vegetation.desc="main") adyacente a cada astilla (Vegetation.desc="sliver") que comparta el borde más largo.
2. Identifique la arista entre cada astilla y su polígono principal adyacente que comparta la arista más larga identificada en el Paso 1. (Posiblemente utilizando ST_Relate)
3. Elimine la arista entre cada astilla y su polígono principal adyacente que comparta la arista más larga utilizando el identificador de arista devuelto en el paso 2. (Posiblemente utilizando ST_RemEdgeModFace)

Sólo tengo un conocimiento muy básico de SQL y el proceso anterior está más allá de mi capacidad.

Rem : Sé que esto es posible en ArcGIS o GRASS con una sola función ("eliminar" o "v.clean rmarea"), pero quiero una solución para PostGIS.

Answer 1

Lo utilizo cuando agrego las formas del conjunto de datos GADM2 en SQL Server:

SET @g2 = @g1.STBuffer(1).Reduce(1).STBuffer(-1).Reduce(1)

En reduce() elimina los artefactos extendidos y acelera las cosas por un factor de 100x - un poco duro, pero bueno para aproximaciones.