Etiqueta | Explicación | Tipo de datos |
Ráster de superficie de entrada | Ráster de entrada de superficie. | Raster Layer |
Ráster de salida | El ráster de salida. Será de tipo punto flotante. | Raster Dataset |
Tipo de parámetro (Opcional) | Especifica el tipo de parámetro de superficie de salida que se calculará.
| String |
Tipo de superficie local (Opcional) | Especifica el tipo de función de superficie que se ajustará alrededor de la celda de destino.
| String |
Distancia de vecindad (Opcional) | La salida se calculará 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 por 3. | Linear Unit |
Utilizar vecindad adaptable (Opcional) | Especifica si la distancia de vecindad variará con los cambios del paisaje (adaptable). 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.
| Boolean |
Unidad Z (Opcional) | Especifica la unidad lineal que se utilizará para los valores z verticales. Se define mediante un sistema de coordenadas verticales si existe. Si no existe ningún sistema de coordenadas verticales, defina la unidad z mediante la lista de unidades para garantizar el cómputo geodésico correcto. El valor predeterminado es metros.
| String |
Medida de pendiente (Opcional) | Unidades de medida (grados o porcentajes) que se utilizarán para el ráster de pendiente de salida. Este parámetro solo está disponible si el parámetro Tipo de parámetro está establecido en Pendiente.
| String |
Acimuts geodésicos del proyecto (Opcional) |
Especifica si se proyectan acimuts geodésicos para corregir la distorsión de ángulo causada por la referencia espacial de salida.
Este parámetro solo está disponible si el parámetro Tipo de parámetro está establecido en Orientación. | Boolean |
Utilizar orientación ecuatorial (Opcional) |
Especifica si la orientación se medirá desde un punto del ecuador o desde el polo norte.
Este parámetro solo está disponible si el parámetro Tipo de parámetro está establecido en Orientación. | Boolean |
Máscara de análisis de entrada (Opcional) | Los datos de entrada definen las ubicaciones en las que se producirá el análisis. Puede ser un ráster o un dataset de entidades. Si la entrada es un ráster, puede ser de tipo entero o de punto flotante. Si la entrada son datos de entidad, puede ser de tipo de punto, línea o polígono. Cuando los datos de máscara de entrada son un ráster, el análisis se realizará en ubicaciones que tengan un valor válido, incluido cero. Las celdas que son NoData en la entrada de la máscara serán NoData en la salida. | Composite Geodataset |
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. Las opciones de parámetros de superficie disponibles para el parámetro Tipo de parámetro (parameter_type en Python) son Pendiente, Orientación, Curvatura media, Curvatura tangencial (curva de nivel normal), Curvatura de perfil (línea de pendiente normal), Curvatura del plano (curva de nivel proyectada),Torsión geodésica de curvas de nivel, Curvatura gaussiana y Curvatura de Casorati.
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 otros factores 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): 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): 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): 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.
- Curvatura del plano (curva de nivel proyectada) (CONTOUR_CURVATURE en Python): curvatura a lo largo de las líneas de curvas de nivel.
- Torsión geodésica de curvas de nivel ( CONTOUR_GEODESIC_TORSION en Python): tasa de cambio en el ángulo de pendiente a lo largo de las líneas de curvas de nivel.
- Curvatura gaussiana (GAUSSIAN_CURVATURE en Python): curvatura general de una superficie. Se calcula como el producto de la curvatura mínima y máxima y puede tomar valores negativos y positivos. Los valores positivos indican que la superficie es convexa en esa celda, mientras que los valores negativos indican que es cóncava. Un valor 0 indica que la superficie es plana.
- Curvatura de Casorati ( CASORATI_CURVATURE en Python): curvatura general de la superficie. Puede ser cero o positivo. Los valores positivos altos indican áreas de curva pronunciada en varias direcciones.
Las unidades de todas las salidas de tipo de curvatura serán el recíproco (el cuadrado del recíproco en el caso de Curvatura gaussiana) de las unidades x,y del entorno Sistema de coordenadas de salida.
La opción Cuadrática del parámetro Tipo de superficie local (local_surface_type = "QUADRATIC" 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.
- Utilice la superficie cuadrática 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 del parámetro Tipo de superficie local (local_surface_type = "BIQUADRATIC" en Python) 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. Deje la distancia de vecindad como predeterminada (igual al tamaño de celda).
El parámetro 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, por ejemplo, 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 = "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.
Al especificar el valor del parámetro Unidad Z (z_unit en Python) de la superficie, se garantiza un cómputo 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 definir 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 de medida del parámetro Medida de pendiente (output_slope_measurement en Python):
- Grados (DEGREE en Python): el rango de valores de pendiente abarca de 0 a 90.
- 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 = "PROJECT_GEODESIC_AZIMUTHS" en Python), se cumple lo siguiente:
- El norte se representa mediante 360 grados.
- Se proyectan acimuts para corregir la distorsión causada por un valor no conforme del entorno Sistema de coordenadas de salida. Estos ángulos se pueden utilizar para localizar con precisión puntos a lo largo de la pendiente descendente más pronunciada.
Si el parámetro Utilizar orientación ecuatorial está activado (project_geodesic_azimuths = "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 ocurre al aproximarse a los polos. Este parámetro garantiza que los ejes norte-sur y este-oeste sean perpendiculares entre sí.
Active el parámetro Utilizar orientación ecuatorial si el terreno está cerca de los polos norte o sur.
Utilice el parámetro Máscara de análisis de entrada (in_analysis_mask en Python) para limitar el análisis a ubicaciones específicas de interés dentro del ráster de superficie de entrada. Las ubicaciones se pueden definir mediante datos ráster o de entidad. El parámetro Máscara de análisis de entrada tendrá prioridad sobre la configuración del entorno Máscara.
Si el valor del Ráster de superficie de entrada (in_raster en Python) y los datos ráster de Máscara de análisis de entrada (in_analysis_mask en Python) tienen el mismo tamaño de celda y las celdas están alineadas, se utilizarán directamente en la herramienta. No se remuestrearán internamente mientras la herramienta esté en funcionamiento.
Si el tamaño de celda es diferente, el tamaño de celda de salida será el máximo de las entradas y el valor de Ráster de superficie de entrada se utilizará internamente como el ráster de alineación. Si el tamaño de celda es el mismo, pero las celdas no están alineadas, el valor del Ráster de superficie de entrada se utilizará internamente como ráster de alineación. En ambos casos, se desencadenará un remuestreo interno antes de realizar la operación de extracción.
Para obtener más información, consulte los temas de los entornos Tamaño de celda y Ráster de alineación.
Referencias:
- James, D. E., Tomer, M. D. y Porter, S. A. 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. y Jenčo, M. A comprehensive system of definitions of land surface (topographic) curvatures, with implications for their application in geoscience modelling and prediction. Earth-Science Reviews, 103414, 2020. https://doi.org/10.1016/j.earscirev.2020.103414
Parámetros
arcpy.ddd.SurfaceParameters(in_raster, out_raster, {parameter_type}, {local_surface_type}, {neighborhood_distance}, {use_adaptive_neighborhood}, {z_unit}, {output_slope_measurement}, {project_geodesic_azimuths}, {use_equatorial_aspect}, {in_analysis_mask})
Nombre | Explicación | Tipo de datos |
in_raster | Ráster de entrada de superficie. | Raster Layer |
out_raster | El ráster de salida. Será de tipo punto flotante. | Raster Dataset |
parameter_type (Opcional) | Especifica el tipo de parámetro de superficie de salida que se calculará.
| String |
local_surface_type (Opcional) | Especifica el tipo de función de superficie que se ajustará alrededor de la celda de destino.
| String |
neighborhood_distance (Opcional) | La salida se calculará 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 por 3. | Linear Unit |
use_adaptive_neighborhood (Opcional) | Especifica si la distancia de vecindad variará con los cambios del paisaje (adaptable). 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.
| Boolean |
z_unit (Opcional) | Especifica la unidad lineal que se utilizará para los valores z verticales. Se define mediante un sistema de coordenadas verticales si existe. Si no existe ningún sistema de coordenadas verticales, defina la unidad z mediante la lista de unidades para garantizar el cómputo geodésico correcto. El valor predeterminado es metros.
| String |
output_slope_measurement (Opcional) | Unidades de medida (grados o porcentajes) que se utilizarán para el ráster de pendiente de salida.
Este parámetro solo está disponible si el parámetro parameter_type está establecido en SLOPE. | 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.
Este parámetro solo está disponible si el parámetro parameter_type está establecido en ASPECT. | Boolean |
use_equatorial_aspect (Opcional) |
Especifica si la orientación se medirá desde un punto del ecuador o desde el polo norte.
Este parámetro solo está disponible si el parámetro parameter_type está establecido en ASPECT. | Boolean |
in_analysis_mask (Opcional) | Los datos de entrada definen las ubicaciones en las que se producirá el análisis. Puede ser un ráster o un dataset de entidades. Si la entrada es un ráster, puede ser de tipo entero o de punto flotante. Si la entrada son datos de entidad, puede ser de tipo de punto, línea o polígono. Cuando los datos de máscara de entrada son un ráster, el análisis se realizará en ubicaciones que tengan un valor válido, incluido cero. Las celdas que son NoData en la entrada de la máscara serán NoData en la salida. | Composite Geodataset |
Muestra de código
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.ddd import *
SurfaceParameters("elevation_1m.tif", "C:/data/output/outsurfaceparameters01.tif",
"", "", "5 METERS", "ADAPTIVE_NEIGHBORHOOD", "", "PERCENT_RISE")
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: 3D Analyst Extension
# Import system modules
import arcpy
from arcpy.ddd import *
# Set environment settings
arcpy.env.workspace = "C:/data"
# Check out the ArcGIS 3D Analyst extension license
arcpy.CheckOutExtension("3D")
# Set local variables
inRaster = "elevation_1m.tif"
outRaster = "C:/data/output/outsurfaceparameters02.tif"
inParameterType = "PROFILE_CURVATURE"
inNeighborhoodDistance = "10 METERS"
inUseAdaptiveNeighborhood = "ADAPTIVE_NEIGHBORHOOD"
# Execute the tool
SurfaceParameters(inRaster, outRaster, inParameterType, "",
inNeighborhoodDistance, inUseAdaptiveNeighborhood)
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: 3D Analyst Extension
# Import system modules
import arcpy
from arcpy.ddd import *
# Set environment settings
arcpy.env.workspace = "C:/data"
# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("3D")
# Set local variables
inRaster = "elevation_1m.tif"
outRaster = "C:/data/output/outsurfaceparameters03.tif"
inParameterType = "ASPECT"
inNeighborhoodDistance = "5 METERS"
inProjectGeodesicAzimuths = "PROJECT_GEODESIC_AZIMUTHS"
# Execute the tool
SurfaceParameters(inRaster, outRaster, inParameterType, "", inNeighborhoodDistance,
"", "", "", inProjectGeodesicAzimuths)
Entornos
Información de licenciamiento
- Basic: Requiere 3D Analyst o Spatial Analyst
- Standard: Requiere 3D Analyst o Spatial Analyst
- Advanced: Requiere 3D Analyst o Spatial Analyst