r.fillnullsgrass: en línea en la nube

PROGRAMA:

NOMBRE

r.fillnulls - Rellena áreas sin datos en mapas ráster mediante interpolación de splines.

PALABRAS CLAVE

ráster, elevación, interpolación

SINOPSIS

r.fillnulls
r.fillnulls --ayuda
r.fillnulls Las opciones de entrada=nombre salida=nombre Método=cadena [tensión=flotar] [sencillo.=flotar]
[Edge=entero] [npmin=entero] [segmax=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

Método=cadena [requerido]
Método de interpolación a utilizar
Opciones: bilineal bicúbico rst
Por defecto: rst

tensión=flotar
Parámetro de tensión de spline
Por defecto: 40.

sencillo.=flotar
Parámetro de suavizado de spline
Por defecto: 0.1

Edge=entero
Ancho del borde del orificio utilizado para la interpolación (en celdas)
Opciones: 2 - 100
Por defecto: 3

npmin=entero
Número mínimo de puntos para la aproximación en un segmento (> segmax)
Opciones: 2 - 10000
Por defecto: 600

segmax=entero
Número máximo de puntos en un segmento
Opciones: 2 - 10000
Por defecto: 300

DESCRIPCIÓN

r.fillnulls llena píxeles NULOS (sin áreas de datos) en el mapa ráster de entrada y almacena datos completos
a un nuevo mapa ráster de salida. Las áreas de relleno se interpolan desde el área sin datos
búfer de límites usando v.surf.primero interpolación spline regularizada con tensión
(método = primero) o r.resamp.bspline interpolación spline cúbica o lineal con Tykhonov
regularización.

NOTAS

Cada zona de influencia de límite de área se establece en tres veces la resolución del mapa para obtener nominalmente tres
puntos alrededor del borde. De esta forma, el algoritmo se interpola en el agujero con un
pendiente y curvatura en los bordes, con el fin de evitar que dicho plano plano se genere en
un agujero. El ancho del área del borde se puede ajustar cambiando el parámetro del borde.

Durante la interpolación, es posible que se produzca la siguiente advertencia al utilizar el método RST:

Advertencia: la tira existe con datos insuficientes
Advertencia: tarda demasiado en encontrar puntos para la interpolación; cambie la región a
área donde están tus puntos

Esta advertencia se genera si existen grandes agujeros de datos dentro de la superficie. Como la idea de
r.fillnulls es llenar esos agujeros, el usuario puede ignorar la advertencia. La interpolación
continuará. Sin embargo, el usuario puede prestar atención a las notas siguientes.

Si falla la interpolación, los mapas ráster y vectoriales temporales se dejan en su lugar para permitir
Identificación del orificio del mapa sin rellenar (área NULL) y reparación manual.

Cuando se utiliza el método RST predeterminado, el algoritmo se basa en v.surf.primero regularizado
splines con módulo de interpolación de tensión que interpola los valores de la celda ráster para
Áreas de datos NULL de los valores límite del área de datos NULL. Una MÁSCARA raster eventual es
respetado durante el llenado de área (s) de datos NULL. Los valores interpolados se parchean en
las áreas de datos NULL del mapa de entrada y se guardan en un nuevo mapa ráster. De lo contrario, ya sea
se puede seleccionar la interpolación spline lineal o cúbica con regularización Tykhonov
(Residencia en r.resamp.bspline).

ADVERTENCIA

Dependiendo de la forma de las áreas de datos NULL, pueden ocurrir problemas debido a una
número de valores de celda de entrada para el proceso de interpolación. La mayoría de los problemas ocurrirán si un
El área de datos NULL alcanza una gran parte del límite del mapa. El usuario deberá tener cuidado
comprobar el resultado usando r.mapcalc (generando un mapa de diferencias para el mapa de entrada y
aplicando la tabla de colores "diferencias" con r.colores) y / o para consultar una celda individual
valores.

El método RST almacena mapas temporales en el disco duro. Requerirá al menos la misma cantidad de espacio libre
como toma un mapa ráster de entrada adicional.

EJEMPLO

En este ejemplo, el mapa de elevación SRTM en la ubicación del conjunto de datos de muestra de Carolina del Norte es
filtrado para valores de elevación atípicos; Los píxeles que faltan se vuelven a interpolar para obtener una
mapa de elevación completo:
g.region raster = elev_srtm_30m -p
d. mon wx0
d.histograma elev_srtm_30m
# eliminar elevaciones demasiado bajas (especialmente lagos)
r.mapcalc "elev_srtm_30m_filt = if (elev_srtm_30m <50.0, null (), elev_srtm_30m)"
d.histograma elev_srtm_30m_filt
d. rast elev_srtm_30m_filt
r.fillnulls input = elev_srtm_30m_filt output = elev_srtm_30m_complete tension = 20
d.histograma elev_srtm_30m_completo
d. rast elev_srtm_30m_complete

Referencias

· Mitas, L., Mitasova, H., 1999, Interpolación espacial. En: P. Longley, MF
Goodchild, DJ Maguire, DWRhind (Eds.), Sistemas de información geográfica:
Principios, técnicas, gestión y aplicaciones, Wiley, págs. 481-492

· Mitasova H., Mitas L., Brown WM, DP Gerdes, I. Kosinovsky, Baker, T.1995,
Modelado de fenómenos distribuidos espacial y temporalmente: nuevos métodos y herramientas para
SIG DE HIERBA. Internacional Actualidad of GIS, 9 (4), número especial sobre Integración de SIG
y Modelado ambiental, 433-446.

· Mitasova H. y Mitas L. 1993: Interpolación por spline regularizado con tensión:
I. Teoría e implementación, Matemático Geología 25, 641-655.

· Mitasova H. y Hofierka L. 1993: Interpolación por spline regularizado con
Tensión: II. Aplicación al modelado del terreno y al análisis de geometría de superficies,
Matemático Geología 25, 657-667.

Utilice r.fillnullsgrass en línea utilizando los servicios de onworks.net