Irregular XYZ de datos de punto a ráster (TIF) - LASTools alternativa?

Question

Estoy tratando de convertir mi irregulares espaciado de punto de datos (xyz) a DEM y luego, en un ráster (geotiff). Debido a que el espaciado irregular necesito una triangulación en primer lugar, a la derecha (así, sin espacios están en el final de la trama)?

Así que la forma es: XYZ -> Triangular -> Interpolar A la Trama. Tengo un montón de puntos (mínimo 150.000), mi raster tamaño debe ser de 0,25 x 0,25 m o menor. Necesito una manera de hacer esto con el open source / free herramientas!

Lo que yo he probado hasta ahora:

• LASTools (las2dem) - funciona muy bien, es muy rápido <- eso sería exactamente lo que necesito, pero no es gratis.

• FUSIÓN (TINSurfaceCreate entonces DTM2ASCII) - funciona bien, pero ralentiza extremas, como la trama de tamaño se reduce (1 x 1 m y 0,5 x 0,5 m, está bien, lo que está debajo es muy lento).

• La SAGA de Herramientas / gdal_grid (a través de QGIS) -> muy lento así.

Cualquier otra sugerencia que debería probar?

Answer 1

Creo que el quid de lidar es si usted necesita el filtro de la no-tierra devuelve - supongo que tiene este hecho ya.

Usted no necesita triangulación necesariamente, también hay kriging, IDW, thin-plate spline y otros. Ver el gstat y de trama de paquetes en R para más.

Aún así, el RTriangle paquete en R ofrece una rápida triangulación (con restricciones) y así me gustaría empezar por ahí. Otros ofrecen paquetes de triangulaciones de Delaunay y que sería igual de bueno en este caso.

He hecho una muy aproximada de las rebajas documento aquí: https://github.com/mdsumner/xyztrisurf/blob/master/xyztrisurf.md

A escala de este a sus datos, sería posible, pero podría necesitar algo de bastante cuidado en su manipulación - no he hecho nada para ver cómo se hace aquí.

FWIW, este "interpolar a la red de triángulos" es, esencialmente, lo que de Matlab griddata por defecto (afaict), y lo Manifold GIS hace con la "triangulación" método para crear una superficie a partir de puntos.

Answer 2

PDAL puede hacer esto con points2grid. Ver http://www.pdal.io/stages/writers.p2g.html para más detalles sobre cómo funciona, y ver https://github.com/Applied-GeoSolutions/lidar2dems para una canalización de procesamiento de solución si usted tiene que hacer un montón de esto.

Nota no hay un ASCII lector para PDAL sin embargo, y se supone que está utilizando algún tipo de LiDAR formato.

Answer 3

Así que mi solución fue la siguiente: estoy usando FUSION de triangulación y de trama de cálculo. se pone un *.dtm archivo (binario) -> la conversión con el FUSION herramientas DTM2ASCII es el cuello de botella.

Así que me puse mis manos en el binario lector de C#, lea la totalidad de la *.dtm archivo y escribí mi propio tif-salida con libtiff.net (http://bitmiracle.com/libtiff/).

No es tan rápido como LAStools, pero para mis necesidades es bueno y funciona bastante bien :).