Cómo funciona Pendiente

Disponible con una licencia de Spatial Analyst.

Disponible con una licencia de 3D Analyst.

La herramienta Pendiente identifica la inclinación en cada celda de una superficie de ráster. Mientras menor sea el valor de la pendiente, más plano será el terreno; mientras más alto sea el valor de la pendiente, más empinado será el terreno.

Nota:

La herramienta Parámetros de superficie proporciona una implementación de pendiente más reciente y se recomienda su uso en lugar de la herramienta Pendiente. La herramienta Pendiente ajusta un plano a las nueve celdas locales, pero un plano puede no ser un buen descriptor del paisaje y puede enmascarar o exagerar variaciones de interés naturales. La herramienta Parámetros de superficie ajusta una superficie a la vecindad de celdas, en lugar de un plano, lo que proporciona un ajuste más natural al terreno.

La herramienta Pendiente utiliza una ventana de celdas de 3 por 3 para calcular el valor, mientras que la herramienta Parámetros de superficie permite tamaños de ventana de 3 por 3 a 15 por 15 celdas. Los tamaños de ventana más grandes son útiles con datos de elevación de alta resolución para capturar procesos de superficie de suelo a una escala adecuada. Parámetros de superficie también proporciona una opción de ventana adaptable, que evalúa la variabilidad local del terreno e identifica el mayor tamaño de vecindad adecuado para cada celda. Puede resultar útil con terreno homogéneo gradual interrumpido por arroyos, carreteras o cortes agudos de pendiente.

Puede seguir utilizando el enfoque tradicional de la herramienta Pendiente si necesita que los resultados coincidan exactamente con las ejecuciones anteriores de la herramienta o si es más importante un tiempo de ejecución más corto que un algoritmo mejor.

El ráster de pendiente de salida se puede calcular en dos tipos de unidades: grados o porcentaje (elevación en porcentaje). La elevación en porcentaje puede comprenderse mejor si se considera como la elevación dividida entre el avance, multiplicada por 100. Considere el triángulo B más abajo. Cuando el ángulo es de 45 grados, la elevación es igual al avance, y la elevación en porcentaje es 100%. A medida que la pendiente alcanza la vertical (90 grados), como se puede ver en el triángulo C, la elevación en porcentaje comienza a acercarse al infinito.

Pendiente en grados frente a porcentaje
Comparación de los valores de la pendiente en grados frente al porcentaje.

La herramienta Pendiente se ejecuta con más frecuencia en un dataset de elevación, tal y como se muestra en las siguientes imágenes. Las pendientes más empinadas están sombreadas en marrón más oscuro en el ráster de pendiente de salida.

Ejemplo de salida de pendiente

La herramienta también se puede utilizar con otros tipos de datos continuos, como la población, para identificar cambios marcados en el valor.

Métodos de cálculo y el efecto del borde

Para calcular la pendiente hay dos métodos disponibles. Puede elegir entre realizar cálculos Planares o Geodésicos con el parámetro Método.

Para el método planar, la pendiente se mide como la tasa máxima de cambio del valor de una celda a sus vecinas inmediatas. El cálculo se realiza en un plano liso proyectado utilizando un sistema de coordenadas cartesianas 2D. El valor de la pendiente se calcula utilizando un estimador de diferencia finita de tercer orden.

Con el método geodésico, el cálculo se realizará en un sistema de coordenadas cartesianas 3D teniendo en cuenta la forma de la Tierra como un elipsoide. El valor de la pendiente se calcula midiendo el ángulo entre la superficie topográfica y el datum de referencia.

Tanto los cálculos planares como geodésicos se realizan utilizando una vecindad de celdas de 3 x 3 (ventana móvil). Para cada vecindad, si la celda de procesamiento (central) es NoData, la salida es no NoData. El cálculo también requiere que al menos siete celdas vecinas a la celda de procesamiento sean válidas. Si hay menos de siete celdas válidas, el cálculo no se realizará y la salida en esa celda de procesamiento será NoData.

Las celdas de las filas y columnas más exteriores del ráster de salida serán NoData. Esto se debe a que esas celdas no tienen suficientes vecinos válidos a lo largo del límite del dataset de entrada.

Método planar

La pendiente se calcula como la tasa de cambio (delta) de la superficie en las direcciones horizontal (dz/dx) y vertical (dz/dy) desde la celda central hasta cada celda adyacente. El algoritmo básico utilizado para calcular la pendiente es el siguiente:

slope_radians = ATAN ( √ ([dz/dx]2 + [dz/dy]2) )

Por lo general, la pendiente se mide en grados, y se utiliza este algoritmo:

slope_degrees = ATAN ( √ ([dz/dx]2 + [dz/dy]2) ) * 57.29578

Nota:

El valor 57,29578 que se muestra aquí es una versión truncada del resultado de 180/pi.

El algoritmo de pendiente también se puede interpretar de la siguiente manera:

slope_degrees = ATAN (rise_run) * 57.29578
  • donde:

    rise_run = √ ([dz/dx]2 + [dz/dy]2]

Los valores de la celda central y sus ocho vecinas determinan las deltas horizontal y vertical. Las vecinas se identifican como letras, de la a a la i, con la letra e representando a la celda para la cual se calcula la pendiente.

Ventana de superficie
Ventana de escaneado de superficie

La tasa de cambio en la dirección x de la celda e se calcula con el siguiente algoritmo:

[dz/dx] = ((c + 2f + i)*4/wght1 - (a + 2d + g)*4/wght2) / (8 * x_cellsize)
  • donde:

    wght1 y wght2 son los recuentos ponderados horizontales de las celdas válidas.

    Por ejemplo, si:

    • c, f e i tienen valores válidos, wght1 = (1+2*1+1) = 4.
    • i es NoData, wght1 = (1+2*1+0) = 3.
    • f es NoData, wght1 = (1+2*0+1) = 2.

    Se aplica una lógica similar a wght2, salvo que las ubicaciones vecinas sean a, d y g.

La tasa de cambio en la dirección y de la celda e se calcula con el siguiente algoritmo:

[dz/dy] = ((g + 2h + i)*4/wght3 - (a + 2b + c)*4/wght4) / (8 * y_cellsize)
  • donde:

    wght3 y wght4 son el mismo concepto que en cálculo de [dz/dx].

Ejemplo de cálculo de pendiente planar

En este ejemplo, se calcula el valor de la pendiente de la celda central de la ventana móvil.

Entrada de ejemplo de pendiente
Entrada de ejemplo de pendiente

La tasa de cambio en la dirección x de la celda central e es:

[dz/dx] = ((c + 2f + i)*4/wght1 - (a + 2d + g)*4/wght2) / (8 * x_cellsize) = ((50 + 60 + 10)*4/(1+2+1) - (50 + 60 + 8)*4/(1+2+1)) / (8 * 5) = (120 - 118) / 40 = 0.05

La tasa de cambio en la dirección y de la celda e es:

[dz/dy] = ((g + 2h + i)*4/wght3 - (a + 2b + c)*4/wght4) / (8 * y_cellsize) = ((8 + 20 + 10)*4/(1+2+1) - (50 + 90 + 50)*4/(1+2+1)) / (8 * 5) = (38 - 190) / 40 = -3.8

Tomando la tasa de cambio en la dirección x e y, la pendiente de la celda central e se calcula utilizando:

rise_run = √ ([dz/dx]2 + [dz/dy]2) = √ ((0.05)2 + (-3.8)2) = √ (0.0025 + 14.44) = 3.80032
slope_degrees = ATAN (rise_run) * 57.29578 = ATAN (3.80032) * 57.29578 = 1.31349 * 57.29578 = 75.25762

El valor de la pendiente en número entero para la celda e es 75 grados.

Salida de ejemplo de pendiente
Salida de ejemplo de pendiente

Método geodésico

El método geodésico mide la pendiente en un sistema de coordenadas 3D geocéntricas, también llamado sistema de coordenadas centrado en la Tierra, Tierra fija (ECEF), teniendo en cuenta la forma de la Tierra como un elipsoide. La forma en que esté proyectado el dataset no afectará al resultado del cálculo. Se utilizarán las unidades z del ráster de entrada si están definidas en la referencia espacial. Si la referencia espacial de la entrada no define las unidades z, deberá hacerlo con el parámetro de unidad z. El método geodésico produce una pendiente más exacta que el método planar.

Transformación de coordenadas geodésicas

El sistema de coordenadas de ECEF es un sistema de coordenadas cartesianas dextrógiro 3D cuyo origen es el centro de la Tierra, donde cualquier ubicación se representa mediante coordenadas X, Y y Z. Consulte la figura siguiente para ver un ejemplo de una ubicación de destino T expresada con coordenadas geocéntricas.

El sistema de coordenadas de ECEF
El ráster de la superficie se transforma desde el sistema de coordenadas de entrada en un sistema de coordenadas geocéntricas 3D.

El cálculo geodésico utiliza una coordenada X, Y, Z que se calcula según sus coordenadas geodésicas (latitud φ, longitud λ, altura h). Si el sistema de coordenadas del ráster de superficie de entrada es un sistema de coordenadas proyectadas (PCS), el ráster se reproyecta primero a un sistema de coordenadas geográficas (GCS) donde cada ubicación tiene una coordenada geodésica y, a continuación, se transforma en el sistema de coordenadas de ECEF. La altura h (valor z) es la altura de elipsoide asociada a la superficie de elipsoide. Consulte el gráfico de la ilustración a continuación.

Altura de elipsoide
Altura de elipsoide

Para transformar a coordenadas de ECEF desde una coordenada geodésica (latitud φ, longitud λ, altura h), utilice las fórmulas siguientes:

X = (N(φ)+h)cosφcosλ
Y = (N(φ)+h)cosφsinλ
Z = (b2/a2*N(φ)+h)sinφ
  • donde:
    • N( φ ) = a2/ √(a2cosφ2+b2sinφ2)
    • φ = latitud
    • λ = longitud
    • h = altura de elipsoide
    • a = eje mayor del elipsoide
    • b = eje menor del elipsoide

La altura h de elipsoide se expresa en metros en las fórmulas anteriores. Si la unidad z de su ráster de entrada se especifica en otra unidad, se transformará internamente a metros.

Cálculo de pendiente

La pendiente geodésica es el ángulo formado entre la superficie topográfica y la superficie del elipsoide. Toda superficie paralela a la superficie del elipsoide tiene una pendiente de 0. Para calcular la pendiente en cada ubicación, se ajusta un plano vecino con celdas de 3 x 3 alrededor de cada celda de procesamiento utilizando el método por mínimos cuadrados. El mejor ajuste en el método por mínimos cuadrados minimiza la suma de la diferencia cuadrada (dzi) entre el valor z real y el valor z ajustado. Consulte la ilustración a continuación para ver un ejemplo.

Ejemplo de ajuste por mínimos cuadrados
Ejemplo de ajuste por mínimos cuadrados

Aquí, el plano se representa como z = Ax + By + C. Para cada centro de celda, dzi es la diferencia entre el valor z real y el valor z ajustado.

El plano se ajusta mejor cuando ∑9i=1dzi2 se minimiza.

Una vez ajustado el plano, se calcula una superficie normal en la ubicación de la celda. En la misma ubicación, también se calcula una normal a la elipsoide en perpendicular al plano tangente de la superficie de elipsoide.

Cómputo de pendiente geodésica
Cómputo de pendiente geodésica

La pendiente, en grados, se calcula desde el ángulo entre la elipsoide normal y la superficie normal topográfica, representada aquí como β. Según la ilustración de arriba, el ángulo α es la pendiente geodésica, que es igual que el ángulo β, según la ley de geometría congruente.

Para calcular la pendiente como elevación en porcentaje, se utiliza la fórmula siguiente:

Slope_PercentRise = ATAN(β) * 100%

¿Debería utilizar la herramienta Parámetros de superficie?

Si el valor del parámetro Ráster de entrada (in_raster en Python) es de alta resolución con un tamaño de celda inferior a unos pocos metros o especialmente ruidoso, plantéese usar la herramienta Parámetros de superficie y su opción de distancia de vecindad definida por el usuario en lugar de la vecindad 3 por 3 inmediata de esta herramienta. Utilizar una vecindad más grande puede minimizar el efecto de las superficies ruidosas. Utilizar una vecindad más grande también puede representar mejor formas de suelo y características de superficie cuando se utilizan superficies de alta resolución.

Utilización de una GPU

Para el método geodésico, esta herramienta es capaz de aumentar el rendimiento si tiene determinado hardware de GPU instalado en su sistema. Consulte la sección Procesamiento de GPU con Spatial Analyst para obtener más información sobre la compatibilidad, cómo configurarla y cómo activarla.

Referencias

Marcin Ligas, and Piotr Banasik, 2011. Conversion between Cartesian and geodetic coordinates on a rotational ellipsoid by solving a system of nonlinear equations (GEODESY AND CARTOGRAPHY), Vol. 60, No 2, 2011, pp. 145-159

B. Hofmann-Wellenhof, H. Lichtenegger and J. Collins, 2001. GPS - theory and practice. Sección 10.2.1. p. 282.

David Eberly 1999. Least Squares Fitting of Data (Geometric Tools, LLC), pp. 3.

Temas relacionados