[gdal-dev] Gdal + numpy + interpolation

[Nicolas B.]

Hi Wesley

I don't know about numpy
but here's an example in which a  gdal command for interpolating point .csv
values using inverse distance is called from python
the raster layer is created using vrt
------------------------------------------------------------------------------------------------------------------------------------

*gdal_grid -a invdist:power=1.0:smoothing=50.0
-txe"+str(xmin)+str(xmax)+"-tye"*

*+str(ymin)+str(ymax)+" -of GTiff -ot Float64 -l input  input.vrt
output.tif").readlines()*

*
--------------------------------------------------------------------------------------------------------------------
*

you can reach other gdal_grid interpolation techniques here

http://www.gdal.org/gdal_grid.html#gdal_grid_algorithms

i don't know if it these algorithms are still maintained

Regards

Nicolas BOZON
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20090501/80d35a86/attachment.html