La salida de GridMetrics no se alinea con el MDT base?

Question

Estoy tratando de derivar las métricas de la cuadrícula a partir de un archivo LiDAR (.las) utilizando la herramienta de FUSION GridMetrics herramienta.

Mi problema es que la trama de salida de GridMetrics no se alinea con el MDT base utilizado para ejecutar la herramienta.

Comando:

E:\FUSION\gridmetrics.exe E:\FUSION\testy\N-34-138-B-
b-1-1-4-3.dtm 3 0.5 E:\FUSION\testy\gridmetrics.csv E:\FUSION\testy\*.las

Ejemplo:

El punto se pone en el centro del píxel "desplazado".

Raster con "offset" en las coordenadas:

• Las coordenadas del centro del píxel son: (628933.5000, 496151.0000)
• Coordenadas de la esquina superior izquierda del píxel: (628933.2500, 496151.2500)

Raster con coordenadas "adecuadas" (y mi archivo MDT que va en GridMetrics ):

• Las coordenadas del centro del píxel son: (628933.7500, 496150.7500)
• Coordenadas de la esquina superior izquierda del píxel: (628933.5000, 496151.0000)

¿Cómo puedo entregar GridMetrics para dar salida a las coordenadas de píxeles adecuadas? También he intentado utilizar el interruptor /extent

E:\FUSION\gridmetrics.exe /extent:E:\FUSION\testy\N-34-138-B-b-1-1-4-3.dtm E:\FUSION\testy\N-34-138-B-
b-1-1-4-3.dtm 3 0.5 E:\FUSION\testy\gridmetrics.csv E:\FUSION\testy\*.las

y /align :

E:\FUSION\gridmetrics.exe /align:E:\FUSION\testy\N-34-138-B-b-1-1-4-3.dtm E:\FUSION\testy\N-34-138-B-
b-1-1-4-3.dtm 3 0.5 E:\FUSION\testy\gridmetrics.csv E:\FUSION\testy\*.las 

con mi archivo DTM en vano.

El tamaño de los píxeles de la trama es de 0,5 m. El MDT original fue creado con LAStools lasgrid .

Answer 1

Trabajar en torno a : utilizar interruptores /grid:X,Y,W,H o gridxy:X1,Y1,X2,Y2 para desplazar manualmente la malla de salida introduciendo sus coordenadas de extensión correctas.

---

Causa de la desalineación (hipótesis):

Tal vez la cabecera del DTM no se está leyendo correctamente. Eso explicaría tanto el desajuste con GridMetrics salida y también conmuta /align y /extent no funcionan, porque se basan en el MDT (también conocido como DEM).

Haga algunas pruebas para corroborar si el encabezado del DEM está relacionado con la causa raíz.

En primer lugar, asegúrese de que el MDE creado desde LAStools lasgrid está correctamente alineado con la nube de puntos. Puede superponer la nube de puntos sobre el MDE e inspeccionarlos visualmente, y puede comparar las extensiones de la nube de puntos con las del MDE (ejecutar lasinfo para comprobar la extensión de las nubes de puntos).

Si todo está bien con el DEM (llamémoslo DEM 1), emite a un archivo .asc (utilice LAStools las2las ) para comprobar su información de cabecera.

Si la cabecera tiene:

• XLLCENTER y YLLCENTER cambiarlos por XLLCORNER y YLLCORNER y menos coordenadas de desplazamiento con la mitad del valor del tamaño de la celda.

• En caso contrario, si tiene XLLCORNER y YLLCORNER cambiarlos por XLLCENTER y YLLCENTER y más coordenadas de desplazamiento con la mitad del valor del tamaño de la celda.

Convertirlo de nuevo a .dtm (LAStools las2las ) y esta es la DEM 2. Compruebe si el MDE 1 y el MDE 2 están correctamente alineados. A continuación, ejecute GridMetrics de nuevo utilizando el DEM 2 y verificar los resultados.