[gdal-dev] The slowness of gdal_grid

[Ari Jolma]

I'm doing some laser scanning data -> geotiff interpolations and I'm 
finding that they take some time - the current run has been going on for 
~20 hours now.

The current run is for a 3 by 2 kilometer area, which I'm interpolating 
into a raster with 1 m x 1 m cells. The data is in PostGIS but it 
shouldn't matter:

gdal_grid -a invdist:radius1=10:radius2=10 -txe 3350505 3353786 -tye 
6692068 6694105 -outsize 3281 2037 -of GTiff -ot Float32 -l nummela 
PG:dbname=nummela nummela.tiff
Grid data type is "Float32"
Grid size = (3281 2037).
Corner coordinates = (3350504.500000 6694105.500000)-(3353786.500000 
6692067.500000).
Grid cell size = (1.000000 1.000000).
Source point count = 2588809.
Algorithm name: "invdist".
Options are 
"power=2.000000:smoothing=0.000000:radius1=10.000000:radius2=10.000000:angle=0.000000:max_points=0:min_points=0:nodata=0.000000"

Any ideas how to make things faster? I did a interpolation to 2 m x 2 m 
raster earlier, which was faster but not very fast either.

Best regards,

Ari