Tengo un mosaico de un 1m LiDar derivados del DEM. Necesito a la salida de un subconjunto de los datos como de los 10m de la DEM. Actualmente estoy usando el agregado de la herramienta en ARCGIS 10 para producir un valor promedio para cada nuevo 10m de píxeles. Cualquier consejo sobre si esta es la mejor técnica para ese tipo de tarea? Es el valor de la media es el mejor enfoque con este tipo de datos?
Respuestas
¿Demasiados anuncios?
hernan43
Puntos
566
Un error común (que me he hecho demasiado) es a la baja muestra un mapa de bits mediante el remuestreo de la herramienta con la interpolación bilineal. Ver esta respuesta para una explicación de por qué esto no es bueno. Una trama puede ser muestreados en tres pasos.
El primer paso puede no ser necesario. Reproyectar la trama al destino de las extensiones. El uso de la interpolación bilineal, y mantener la celda de salida tamaño de la misma como la resolución de entrada (por ejemplo, 1 m). Utilice el punto de registro para "ajustar" la trama de las esquinas para la proyección. La salida de las extensiones puede ser especificado en los "Entornos", y sugiero que la especificación de las extensiones con un múltiplo de 10 m (o cualquiera que sea la resolución). Estos límites de control donde las estadísticas se determinan para el final de la trama.
Realizar el Bloque de Estadísticas (que se encuentra en Spatial Analyst Tools > Barrio). El uso de un rectángulo con 10 células por tanto la altura y la anchura, y elegir la opción "MEDIA" para un tipo de estadística. Pruebe diferentes formas y tipos, si quieres. El tamaño de la celda es la relación de muestreo.
Desde el bloque de estadísticas no cambia la resolución de la trama, el último paso es volver a muestrear (que se encuentra en Herramientas de Administración de Datos > Raster > Raster de Procesamiento). Elegir 10 m, y el uso de "más CERCANO" para elegir el bloque de estadística en el centro de la célula.
Los pasos 2 y 3 pueden ser reemplazados con la Curtvprice la sugerencia de utilizar el Agregado de la herramienta, lo que efectivamente va a producir los mismos resultados usando el rectángulo de medios.
Lazer
Puntos
3926
Usted puede cortar un paso del proceso mediante el Agregado de la herramienta, lo que evita tener que replicar todos aquellos valor de resumen de las células, y el volver a muestrear.
Dan
Puntos
16
Todos los puntos mencionados anteriormente son importantes para tomar en cuenta y estoy totalmente de acuerdo en que bilineal de remuestreo es bastante problemática. Aunque, tengo curiosidad por qué nadie está discutiendo la convolución cúbica? El problema con el uso de un bloque de función es que la media es bastante irrelevante cuando la distribución no es normal o multimodal, como era de esperar con un lidar derivados del DEM.
Si usted tiene acceso a los originales de los datos lidar, sólo interpolar los datos para la resolución deseada utilizando el "Topo a Ráster" de la herramienta de ArcGIS. Si sólo tiene acceso a 1m DEM trama, me parece la mejor, aunque de manera menos eficiente, el método sería convertir la trama de puntos y el uso de una placa delgada o Bi-spline cúbico. Esto permitiría la remuestrear barrio para que se ajuste de una curva no lineal a los datos.
Alternativamente, usted podría suavizar el 1m ráster utilizando un núcleo Gaussiano, aproximar el tamaño de su deseado volver a muestrear la resolución (10x10) y, a continuación, un bilineal remuestrear sería mucho más apropiado. Este enfoque permitiría tener un control directo sobre el parámetro de suavizado y resultaría en un "local" de la distribución normal, donde la media se convierte en relevante como indicador de tendencia central y de un ajuste lineal es compatible.
