Este es el comando r.thingrass que se puede ejecutar en el proveedor de alojamiento gratuito de OnWorks utilizando una de nuestras múltiples estaciones de trabajo en línea gratuitas, como Ubuntu Online, Fedora Online, emulador en línea de Windows o emulador en línea de MAC OS.
PROGRAMA:
NOMBRE
r.delgado - Adelgaza las celdas no nulas que denotan entidades lineales en una capa de mapa ráster.
PALABRAS CLAVE
trama, geometría
SINOPSIS
r.delgado
r.delgado --ayuda
r.delgado Las opciones de entrada=nombre salida=nombre [iteraciones=entero] [-exagerar] [-ayuda] [-verboso]
[-tranquilo] [-ui]
Banderas
--Sobrescribir
Permitir que los archivos de salida sobrescriban los archivos existentes
--ayuda
Resumen de uso de impresión
--verboso
Salida del módulo detallado
--tranquilo
Salida de módulo silencioso
--ui
Forzar el inicio del cuadro de diálogo GUI
parámetros:
Las opciones de entrada=nombre [requerido]
Nombre del mapa ráster de entrada
salida=nombre [requerido]
Nombre del mapa ráster de salida
iteraciones=entero
Número máximo de iteraciones
Por defecto: 200
DESCRIPCIÓN
r.delgado escanea el nombrado Las opciones de entrada capa de mapa ráster y adelgaza las celdas no NULL que denotan lineal
entidades en entidades lineales que tienen un solo ancho de celda. Las líneas ráster a menudo necesitan ser
adelgazado (esqueletizando entidades ráster) a un solo ancho de píxel antes de que puedan ser
transformado a datos vectoriales.
r.delgado adelgazará solo las celdas ráster no NULL (sin datos) del nombre Las opciones de entrada mapa raster
capa dentro de la configuración actual de la región geográfica. El ancho de la celda del adelgazado salida
La capa de mapa ráster será igual a la resolución de celda de la ubicación geográfica establecida actualmente.
región. Todas las entidades lineales adelgazadas tendrán el ancho de una sola celda.
r.delgado creará un nuevo salida archivo de datos ráster que contiene las entidades lineales adelgazadas.
r.delgado asume que las características lineales están codificadas con valores positivos sobre un fondo de
NULL está en el Las opciones de entrada archivo de datos ráster, por lo que crea un mapa de salida NULL / 1.
NOTAS
r.delgado solo crea capas de mapa ráster. Para crear un mapa vectorial, el usuario necesitará
para ejecutar r.a.vect en el mapa ráster resultante.
r.delgado puede crear pequeñas espuelas o "líneas colgantes" durante el proceso de adelgazamiento. Estas
las espuelas se pueden eliminar (después de crear una capa de mapa vectorial) v.limpio (rmdangle herramienta).
Este código implementa el algoritmo de adelgazamiento descrito en "Análisis de algoritmos de adelgazamiento
Usando morfología matemática "por Ben-Kwei Jang y Ronlad T. Chin en Transacciones on
Patrón de Costura Analisis y Máquina Intelligence, vol. 12, núm. 6, junio de 1990. La definición Jang
y Chin da del proceso de adelgazamiento es "eliminación sucesiva de capas externas de píxeles
de un objeto conservando los píxeles cuya eliminación alteraría la conectividad o
acortar las patas del sceleton ".
El sceleton finalmente se adelgaza cuando el proceso de adelgazamiento converge; es decir, "no más
los píxeles se pueden eliminar sin alterar la conectividad ni acortar las patas del sceleton "
(pág.541). Los autores prueban que el proceso de aclareo descrito siempre converge y
produce sceletons de un píxel de grosor. El número de iteraciones depende del original
espesor del objeto. Cada iteración despega los píxeles externos del objeto.
Por lo tanto, si el objeto tiene un grosor de <= n píxeles, el algoritmo debería converger en <=
iteraciones.
EJEMPLO
Para vectorizar el mapa ráster flujos_derivados en el conjunto de datos de muestra de Carolina del Norte que
representa la red de transmisión derivada del DEM de resolución de 10 m por cuenca r., correr:
g.region ráster = elevación -p
# crear mapa de acumulación de flujo
r. elevación de la cuenca = elevación acumulación = acum_50K umbral = 50000
# extraer corrientes del mapa de acumulación de flujo
r.mapcalc "streams_from_flow = if (abs (acum_50K)> 1000, 1, null ())"
# esqueletizar mapa
r.thin streams_from_flow out = streams_thin
d. mon wx0
d.rast streams_from_flow
d. borrar
d. rast streams_thin
Adelgazamiento de características ráster (esqueletización)
El mapa resultante cabe opcionalmente vectorizado:
r.to.vect streams_thin salida = streams_thin tipo = línea
# visualizar
d.rast acum_50K
d.vect streams_thin color = rojo ancho = 2
Red de arroyos vectorizada después del adelgazamiento extraída del mapa de acumulación de flujo
Utilice r.thingrass en línea utilizando los servicios de onworks.net