Parámetros de superficie (Spatial Analyst)

Disponible con una licencia de Spatial Analyst.

Disponible con una licencia de 3D Analyst.

Resumen

Determina parámetros de una superficie ráster, como la orientación, la pendiente y las curvaturas.

Más información sobre cómo funciona Parámetros de superficie.

Uso

  • Los parámetros de salida se calculan celda por celda, ajustando una superficie local en torno a una celda de destino. Se pueden calcular cinco parámetros de superficie: Pendiente, Orientación, Curvatura media, Curvatura tangencial (curva de nivel normal) y Curvatura de perfil (línea de pendiente normal). Estas opciones están disponibles en el parámetro Tipo de parámetro (parameter_type en Python).

  • Todos los parámetros de salida se calculan utilizando coordenadas geodésicas y ecuaciones.

  • Cuando se especifica la opción Pendiente (SLOPE en Python) para el Tipo de parámetro, la salida representa la tasa de cambio de elevación para cada celda del modelo digital de elevación (DEM). Es la primera derivada de un DEM. El rango de valores de la salida de pendiente depende del tipo de unidades de medida.

  • Cuando se especifica la opción Orientación (ASPECT en Python) para el Tipo de parámetro, la salida identifica la dirección de la brújula a la que apunta la pendiente descendente en cada ubicación. Se expresa en grados positivos de 0 a 360, medidos en el sentido de las agujas del reloj partiendo del norte.

  • La curvatura se utiliza para describir la forma de una superficie. En ciencias de la Tierra, se utiliza para ayudar a comprender el impacto de la gravedad, la erosión y otras fuerzas en la superficie, y se utiliza junto con otros parámetros de superficie para identificar y clasificar formas de suelo.

    • Curvatura media (MEAN_CURVATURE en Python) mide la curvatura general de la superficie. Se calcula como la media de las curvaturas mínima y máxima. Si se especifica en Tipo de parámetro, la salida equivale al valor medio de las curvaturas de perfil (línea de pendiente normal) y tangencial (curva de nivel normal). Su signo, positivo o negativo, no es un indicador definitivo, excepto en valores extremos. Unos valores positivos altos indican áreas de denudación máxima, y unos valores negativos altos indican áreas de acumulación máxima (Minár et al., 2020).
    • Curvatura de perfil (línea de pendiente normal) (PROFILE_CURVATURE en Python) mide la curvatura normal geométrica a lo largo de la línea de pendiente. Unos valores positivos indican áreas de aceleración del flujo de la superficie y erosión. Una curvatura de perfil negativa indica áreas de deceleración del flujo de la superficie y deposición. Una curvatura de perfil (línea de pendiente normal) positiva indica que la superficie es convexa en esa celda en la dirección de la pendiente. Una curvatura negativa indica que la superficie es cóncava en esa celda y en esa misma dirección. Un valor 0 indica que la superficie es plana.
    • Curvatura tangencial (curva de nivel normal) (TANGENTIAL_CURVATURE en Python) mide la curvatura normal geométrica perpendicular a la línea de pendiente, tangente a la línea de curvas de nivel. Unos valores positivos indican áreas de flujo de superficie divergente. Unas curvaturas tangenciales negativas indican áreas de flujo de superficie convergente. Una curvatura tangencial (curva de nivel normal) positiva indica que la superficie es convexa en esa celda en perpendicular a la dirección de la pendiente. Una curvatura negativa indica que la superficie es cóncava en esa celda y en dirección perpendicular a la pendiente. Un valor 0 indica que la superficie es plana.
  • La opción Cuadrática (QUADRATIC en Python) de Tipo de superficie local (local_surface_type en Python) no ajusta las celdas de vecindad de forma exacta. Esta es la opción predeterminada, y la recomendada para la mayoría de aplicaciones y datos.

    • La superficie cuadrática minimiza el efecto de datos de superficie ruidosos, como una superficie LIDAR de alta resolución, lo cual es especialmente importante al calcular la curvatura.
    • La superficie cuadrática debería utilizarse al especificar un tamaño de vecindad mayor que el tamaño de celda, así como al utilizar la opción de vecindad adaptable.
  • La opción Bicuadrática (BIQUADRATIC en Python) de Tipo de superficie local ajusta los datos de las celdas de vecindad de forma exacta.

    • Esta opción es adecuada para una superficie de entrada altamente precisa.
    • Si la distancia de vecindad es mayor que el tamaño de celda ráster de entrada, se perderán las ventajas de precisión del tipo de superficie Bicuadrática. La distancia de vecindad se debería dejar en su valor predeterminado (igual que el tamaño de celda).
  • Distancia de vecindad (neighborhood_distance en Python) determina el tamaño de la vecindad y calcula el parámetro de superficie a lo largo de esta distancia desde el centro de la celda de destino.

    • No puede ser menor que el tamaño de celda ráster de entrada.
    • Una distancia de vecindad menor captura más variabilidad local en el paisaje, en las características de las entidades de paisaje más pequeñas. Si se dispone de datos de elevación de alta resolución, es más adecuado utilizar distancias mayores.
  • Si el parámetro Utilizar vecindad adaptable está activado (use_adaptive_neighborhood tiene el valor ADAPTIVE_NEIGHBORHOOD en Python), la distancia de vecindad cambiará con la variabilidad del terreno. La distancia de vecindad se reducirá si existe demasiada variabilidad en la ventana de cálculo.

  • Especificar la Unidad z (z-unit en Python) de la superficie es importante para garantizar el cálculo correcto de la salida de pendiente.

    Si hay una unidad z disponible en el sistema de coordenadas verticales del ráster de entrada, se aplicará automáticamente. Se recomienda que defina una unidad z para el ráster de entrada si falta esta unidad. Puede usar la herramienta Definir proyección para especificar una unidad z. Si no se ha definido, se utilizará el metro de manera predeterminada.

  • El rango de valores de la salida de pendiente depende de las unidades del tipo de Medida de pendiente de salida (output_slope_measurement en Python):

    • En Grados (DEGREE en Python), el rango de valores de pendiente va de 0 a 90.
    • En Elevación en porcentaje (PERCENT_RISE en Python), el rango va de 0 hasta básicamente el infinito. Una superficie plana es 0 por ciento, una superficie de 45 grados es 100 por ciento y, a medida que la superficie se vuelve más vertical, la elevación en porcentaje es cada vez mayor.
  • Si el parámetro Acimuts geodésicos del proyecto está activado (project_geodesic_azimuths tiene el valor PROJECT_GEODESIC_AZIMUTHS en Python), se cumple lo siguiente:

    • El norte siempre se representa mediante 360 grados.
    • Se proyectan acimuts para corregir la distorsión causada por un Sistema de coordenadas de salida no conforme. Estos ángulos se pueden utilizar para localizar con precisión puntos a lo largo de la pendiente descendente más pronunciada.

    Active el parámetro Acimuts geodésicos del proyecto si usa la salida de Parámetros de superficie como entrada de dirección hacia atrás del parámetro Ráster de dirección de flujo o dirección hacia atrás de entrada (in_back_direction_raster en Python) para una herramienta del conjunto de herramientas Distancia.

  • Si el parámetro Utilizar orientación ecuatorial está activado (project_geodesic_azimuths tiene el valor USE_EQUATORIAL_ASPECT en Python), la orientación se medirá desde un punto a lo largo del ecuador para corregir el sesgo de la dirección que sucede al aproximarse a los polos. Esta opción garantiza que los ejes norte-sur y este-oeste sean perpendiculares entre sí.

    Active el parámetro Utilizar orientación ecuatorial si su terreno está cerca de los polos norte o sur.

  • Consulte Entornos de análisis y Spatial Analyst para obtener detalles adicionales sobre los entornos de geoprocesamiento que se aplican a esta herramienta.

  • Referencias:

    • James D.E., M.D. Tomer, S.A. Porter. (2014). Trans-scalar landform segmentation from high-resolution digital elevation models. Poster presented at: ESRI Annual Users Conference; July 2014; San Diego, California.
    • Minár, J., Evans, I. S., & Jenčo, M. (2020). A comprehensive system of definitions of land surface (topographic) curvatures, with implications for their application in geoscience modelling and prediction. Earth-Science Reviews, 103414. https://doi.org/10.1016/j.earscirev.2020.103414

Sintaxis

SurfaceParameters(in_raster, {parameter_type}, {local_surface_type}, {neighborhood_distance}, {use_adaptive_neighborhood}, {z_unit}, {output_slope_measurement}, {project_geodesic_azimuths}, {use_equatorial_aspect})
ParámetroExplicaciónTipo de datos
in_raster

Ráster de entrada de superficie.

Raster Layer
parameter_type
(Opcional)

Determina el tipo de parámetro de superficie de salida que se va a calcular.

  • SLOPELa tasa de cambio de la elevación. Esta es la opción predeterminada.
  • ASPECTLa dirección de pendiente descendente de la tasa máxima de cambio de cada celda.
  • MEAN_CURVATUREMide la curvatura general de la superficie. Se calcula como la media de las curvaturas mínima y máxima. Esta curvatura describe la convexidad o concavidad intrínseca de la superficie, independientemente de la dirección o la influencia de la gravedad.
  • TANGENTIAL_CURVATUREMide la curvatura normal geométrica perpendicular a la línea de pendiente, tangente a la línea de curvas de nivel. Esta curvatura se suele aplicar para caracterizar la convergencia o divergencia del flujo por la superficie.
  • PROFILE_CURVATUREMide la curvatura normal geométrica a lo largo de la línea de pendiente. Esta curvatura se suele aplicar para caracterizar la aceleración y deceleración del flujo por la superficie.
String
local_surface_type
(Opcional)

Determina el tipo de función de superficie ajustada alrededor de la celda de destino.

  • QUADRATICAjusta una función de superficie cuadrática a las celdas de la vecindad. Esta es la opción predeterminada.
  • BIQUADRATICAjusta una función de superficie bicuadrática a las celdas de la vecindad.
String
neighborhood_distance
(Opcional)

La salida se calcula a lo largo de esta distancia desde el centro de la celda de destino. Determina el tamaño de la vecindad.

El valor predeterminado es el tamaño de celda ráster de entrada, que da lugar a una vecindad de 3 x 3.

Linear Unit
use_adaptive_neighborhood
(Opcional)

Permite que la distancia de vecindad varíe con los cambios del paisaje. La distancia máxima viene determinada por la distancia de vecindad. La distancia mínima es el tamaño de celda ráster de entrada.

  • FIXED_NEIGHBORHOODUtilizar una única distancia de vecindad (fija) en todas las ubicaciones. Esta es la opción predeterminada.
  • ADAPTIVE_NEIGHBORHOODUtilizar una distancia de vecindad adaptable en todas las ubicaciones.
Boolean
z_unit
(Opcional)

La unidad lineal de los valores z verticales.

Se define mediante un sistema de coordenadas verticales si existe. Si no existe un sistema de coordenadas verticales, la unidad z debe definirse en una lista de unidades para garantizar el cómputo geodésico correcto. El valor predeterminado es metros.

  • INCHLa unidad lineal será pulgadas.
  • FOOTLa unidad lineal será pies.
  • YARDLa unidad lineal será yardas.
  • MILE_USLa unidad lineal será millas.
  • NAUTICAL_MILELa unidad lineal será millas náuticas.
  • MILLIMETERLa unidad lineal será milímetros.
  • CENTIMETERLa unidad lineal será centímetros.
  • METERLa unidad lineal será metros.
  • KILOMETERLa unidad lineal será kilómetros.
  • DECIMETERLa unidad lineal será decímetros.
String
output_slope_measurement
(Opcional)

Cuando parameter_type es SLOPE, determina las unidades de medida (grados o porcentajes) del ráster de pendiente de salida.

  • DEGREELa inclinación de la pendiente se calculará en grados.
  • PERCENT_RISELa inclinación de la pendiente se calculará como elevación en porcentaje, que también se conoce como pendiente en porcentaje.
String
project_geodesic_azimuths
(Opcional)

Especifica si se proyectan acimuts geodésicos para corregir la distorsión de ángulo causada por la referencia espacial de salida.

  • GEODESIC_AZIMUTHSNo se proyectan acimuts geodésicos. Esta es la opción predeterminada.
  • PROJECT_GEODESIC_AZIMUTHSSe proyectan acimuts geodésicos.
Boolean
use_equatorial_aspect
(Opcional)

Medir la orientación desde un punto del ecuador.

  • NORTH_POLE_ASPECTMedir la orientación desde el polo norte. Esta es la opción predeterminada.
  • EQUATORIAL_ASPECTMedir la orientación desde un punto del ecuador.
Boolean

Valor de retorno

NombreExplicaciónTipo de datos
out_raster

El ráster de salida.

Raster

Muestra de código

Ejemplo 1 de SurfaceParameters (ventana de Python)

En el siguiente ejemplo se muestra cómo usar esta herramienta en la ventana de Python.

Este ejemplo genera un ráster de pendiente con valores de salida en porcentajes, utilizando un método de vecindad adaptable. La distancia de vecindad máxima es 5 metros.

from arcpy.sa import *
outSurfaceParameters = SurfaceParameters("elevation_1m.tif", "", "", "5 METERS",
                                         "ADAPTIVE_NEIGHBORHOOD", "", "PERCENT_RISE")
outSurfaceParameters.save("C:/sapyexamples/output/outsurfaceparameters01.tif")
Ejemplo 2 de SurfaceParameters (script independiente)

En el siguiente ejemplo se muestra cómo usar esta herramienta en un script independiente de Python.

Este ejemplo genera un ráster de curvatura de perfil (línea de pendiente normal), utilizando un método de vecindad adaptable. La distancia de vecindad máxima es 10 metros.

# Name: SurfaceParameters_Ex_02.py
# Description: Derive profile (normal slope line) curvature for a 1m resolution
# elevation raster over an adaptive neighborhood distance of maximum 10m. 
# Requirements: Spatial Analyst Extension

# Import system modules
import arcpy
from arcpy.sa import *

# Set environment settings
arcpy.env.workspace = "C:/sapyexamples/data"

# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")

# Set local variables
inRaster = "elevation_1m.tif"
inParameterType = "PROFILE_CURVATURE"
inNeighborhoodDistance = "10 METERS"
inUseAdaptiveNeighborhood = "ADAPTIVE_NEIGHBORHOOD"

# Execute the tool
outSurfaceParameters = SurfaceParameters(inRaster, inParameterType, "",
                                         inNeighborhoodDistance, inUseAdaptiveNeighborhood)

# Save the output 
outSurfaceParameters.save("C:/sapyexamples/output/outsurfaceparameters02.tif")
Ejemplo 3 de SurfaceParameters (script independiente)

En el siguiente ejemplo se muestra cómo usar esta herramienta en un script independiente de Python.

Este ejemplo genera un ráster de orientación utilizando una distancia de vecindad de 5 metros. Corregir las distorsiones de dirección derivadas del uso de una proyección no conforme.

# Name: SurfaceParameters_Ex_03.py
# Description: Derive aspect for an elevation surface over a distance of 5m, correct
# for direction distortion from non-conformal projection system. 
# Requirements: Spatial Analyst Extension

# Import system modules
import arcpy
from arcpy.sa import *

# Set environment settings
arcpy.env.workspace = "C:/sapyexamples/data"

# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")

# Set local variables
inRaster = "elevation_1m.tif"
inParameterType = "ASPECT"
inNeighborhoodDistance = "5 METERS"
inProjectGeodesicAzimuths = "PROJECT_GEODESIC_AZIMUTHS"

# Execute the tool
outSurfaceParameters = SurfaceParameters(inRaster, inParameterType, "",
                                         inNeighborhoodDistance, "", "", "",
                                         inProjectGeodesicAzimuths)

# Save the output 
outSurfaceParameters.save("C:/sapyexamples/output/outsurfaceparameters03.tif")

Información de licenciamiento

  • Basic: Requiere Spatial Analyst o 3D Analyst
  • Standard: Requiere Spatial Analyst o 3D Analyst
  • Advanced: Requiere Spatial Analyst o 3D Analyst

Temas relacionados