Hi Lionel,

On 9/20/07, Lionel Roubeyrie <[EMAIL PROTECTED]> wrote:
> Hi Jeff,
> thanks for the reply. Effectively, I saw the warpimage example, and based on
> that I just want to know if somebody has already used the PCL module to
> retrieve geographical informations (or another module).

With GDAL, it's trivial. I work with data from many sources, and I am
usually bypassing basemap altogether and going GDAL/OGR+MPL.

To show the first band of a raster map, I just do the following:
import pylab
import gdal

fname="/path/to/file"

g = gdal.Open ( fname)
data = g.GetRasterBand(1).ReadAsArray()
#data is a Numeric/numpy array. Can modify at will
pylab.imshow(data)
pylab.show()

Cheers,
Jose

-------------------------------------------------------------------------
This SF.net email is sponsored by: Microsoft
Defy all challenges. Microsoft(R) Visual Studio 2005.
http://clk.atdmt.com/MRT/go/vse0120000070mrt/direct/01/
_______________________________________________
Matplotlib-users mailing list
Matplotlib-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/matplotlib-users

Reply via email to