On Fri, Apr 18, 2008 at 5:25 AM, Graham Carlyle
<[EMAIL PROTECTED]> wrote:
> Hi
>
> I'm wanting to be able to access the values at certain points in a
> raster file. Worldmill doesn't seem to support raster data? I did look
> at PCL-GDAL but it doesn't seem to provide access to the data. For the
> moment I guess I can use the slightly verbose python GDAL bindings.
>
> Also is PCL being developed further or is all new work in the other
> libraries (shapely, worldmill etc.) I'm slightly confused as to which
> libraries to use.
>
> thanks,
> Graham
>
> _______________________________________________
> Community mailing list
> [email protected]
> http://lists.gispython.org/mailman/listinfo/community
>
hi, have you tried gdal's ReadAsArray()? once you do something like:
>>> data = gdal.Open('raster')
>>> arr = data.ReadAsArray()
it's possible to index arr as a numpy array.
though you may want to use GetRasterBand() and send in some
offsets/sizes to ReadAsArray().
then you can do the normal numpy stuff--including plotting.
http://www.krugle.org/kse/files/svn/svn.sourceforge.net/matplotlib/toolkits/basemap/examples/testgdal.py
_______________________________________________
Community mailing list
[email protected]
http://lists.gispython.org/mailman/listinfo/community
