[gdal-dev] Null data and CRS. gdal_warp overkill???

[nicolas.cadieux at archeotec.ca]

Hi,
I am creating a Windows .bat files with a Python script to make DEMs 
with LiDAR points.  The idea is simple.
1) read input txt  file to get the cloud geometries

2)gdal_grid with nodata value marker set to -32768

3)gdal_warp with -dstnodata to set the value marker (of grid created in 
step 2)

4)gdal_fillnotdata.py to replace the no data values (of grid created in 
step 3)

5)gdal_warp again with  -dstnodata to set the value marker (again as 
some values may still be outside of spet 4 search radius)

6)gdaldem hillshade made with DEM of set 5

7)gdal_warp again with -t_srs srs_def: to set the CRS of the Hillshade 
created in step 6.

It looks like I am loosing time with gdal_warp (step 3,5 and7) as I am 
reading a file, adding novalue data or setting CRS, then writing the new 
file and then deleting the old file.

Question: Is there another way of doing this without steps 3, 5 and 7 or 
did i get it right?

Thanks,
Nicolas
P.S.
I am using win7_64 and gdal 1.11.2.
//
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20151125/43cd5f55/attachment.html>