On Oct 8, 8:20 am, Timmie <[email protected]> wrote:
> Week 5: Map algebra and writing ra
> -> use numpy.ma
>
> Week 7: Misc stuff
> us subprocess instad of os.system
>
> Week 1: Reading and writing vector
> * use PythonXY or EPD gives you a full distribution for science.
>
> I question would like to see tackled in a future course:
> If I import vector data as points in X-Y-Z format,
> create an array with numpy or any more advanced
> geo statistics tools (R-Project), how do I set the
> properties of the ds.GeoTransform?
> Or independantly from gdal:
> How do I best save an array as a georeferenced image.
i agree, these examples are a great resource.
here's an example of how, given some data in a 2d numpy array/grid to
set the geotransform, save to a.tif, then reload and plot with MPL.
http://gist.github.com/205115
there's a few interpolation resources in scipy/mpl, griddata in mpl is
pretty easy to use:
http://www.scipy.org/Cookbook/Matplotlib/Gridding_irregularly_spaced_data
-brentp
