Eliminar los pequeños polígonos del shapefile

Question

He descargado GTOPO30 *.DEM y extrajo una línea de contorno a 2000m de altitud usando gdal_contour . En el shapefile resultante, tengo muchos polígonos pequeños en áreas donde la altitud está cerca de esa línea de contorno, y la imagen se ve bastante ruidosa.

¿Cómo puedo, preferiblemente usando gdal/ogr o QGis, de alguna manera "reducir" estos contornos, para que la imagen se vea más limpia? Básicamente quiero eliminar todas las pequeñas características.

PD: Una forma sería reducir la muestra de la cuadrícula antes de hacer la extracción del contorno; ¿es esta la forma de hacerlo, o debería reducir la muestra del archivo de formas?

Answer 1

En realidad, resolví esto usando una combinación de GDAL y QGis:

con GDAL:

Fusionar los diferentes azulejos de GTopo30 en un archivo

gdal_merge.py -o gtopo30_usa.tif *.DEM

Bajar la muestra de este archivo fusionado a 0.1°

gdalwarp gtopo30_usa.tif -tr 0.1 0.1 gtopo30_usa_downsampled.tif

Extraer la línea de contorno

gdal_contour -fl 2000 gtopo30_usa_downsampled.tif gtopo30_usa_contour2000.shp

en QGis:

• crear un nuevo proyecto
• añadir la capa de vectores, seleccionar la gtopo30_usa_contour2000.shp shapefile
• Vector->Geometría->Líneas a los polígonos
• Vector->Geometría->Añadir columnas de geometría
• Capa->Pregunta-> seleccione un área de umbral

Answer 2

Un método, utilizado para probar la Resultado de Paris-Harrington una declaración de la teoría de Ramsey que es independiente de la PA, funciona más o menos como sigue. La declaración a probar tiene la forma $ \forall n\ \exists k\ \varphi (n,k)$ . Por lo tanto, si se considera la función $f$ mapeando cada $n$ al menor testigo $k$ lo que la declaración equivale a la afirmación de que $f$ es una función total, que $f(n)$ se define para todos $n$ .

Ahora, la función que surge en el resultado de Paris-Harrington es una función de crecimiento enormemente rápido, de crecimiento alucinante. La forma en que el argumento funciona, muy a grandes rasgos, es trabajar en un modelo no estándar $M$ de PA, pero encontrar un segmento inicial $I \lt M$ un corte, tal que $f$ salta sobre el corte $I$ pero de tal manera que $I$ continúa satisfaciendo a la AP de todos modos. Así, en el modelo $I$ la función $f$ es no total, y así se tiene el modelo deseado donde la declaración es falsa.