Hi Wesley<br><br>I don&#39;t know about numpy <br>but here&#39;s an example in which a  gdal command for interpolating point .csv values using inverse distance is called from python<br>the raster layer is created using vrt <br>
------------------------------------------------------------------------------------------------------------------------------------<br><meta http-equiv="CONTENT-TYPE" content="text/html; charset=utf-8"><title></title><meta name="GENERATOR" content="OpenOffice.org 3.0  (Win32)"><style type="text/css">
        &lt;!--
                @page { margin: 2cm }
                P { margin-bottom: 0.21cm }
        --&gt;
        </style>
<p style="margin-bottom: 0cm; text-align: left;"><i><font face="Times New Roman, serif"><font style="font-size: 11pt;" size="2">gdal_grid
-a invdist:power=1.0:smoothing=50.0 -txe&quot;+str(xmin)+str(xmax)+&quot;-tye&quot;</font></font></i></p><div>
</div><p style="margin-bottom: 0cm; text-align: left;"><i><font face="Times New Roman, serif"><font style="font-size: 11pt;" size="2">+str(ymin)+str(ymax)+&quot;
-of GTiff -ot Float64 -l input  input.vrt output.tif&quot;).readlines()</font></font></i></p><p style="margin-bottom: 0cm; text-align: left;"><i><font face="Times New Roman, serif"><font style="font-size: 11pt;" size="2">--------------------------------------------------------------------------------------------------------------------<br>
</font></font></i></p><p style="margin-bottom: 0cm; text-align: left;">you can reach other gdal_grid interpolation techniques here</p><p style="margin-bottom: 0cm; text-align: left;"><a href="http://www.gdal.org/gdal_grid.html#gdal_grid_algorithms">http://www.gdal.org/gdal_grid.html#gdal_grid_algorithms</a></p>
<p style="margin-bottom: 0cm; text-align: left;">i don&#39;t know if it these algorithms are still maintained <br></p><p style="margin-bottom: 0cm; text-align: left;">Regards</p><p style="margin-bottom: 0cm; text-align: left;">
Nicolas BOZON<br></p><p style="margin-bottom: 0cm; text-align: left;"><br></p>
<br>