Another option is to use VisAD's Jython connection to do these. Here's a
sample script that will take a GIF image, assume that it is a rectangular
projection between longitudes -110:-70 and latitudes 33:53, recast it into a
new domain, and display it along with a basemap. It does, however, extract the
number of rows and columns from the GIF file.
You will, of course, have to substitute your own files and coordinates, but I
only offer this as an example:
from visad.python.JPythonMethods import *
from visad import *
clearplot()
a=load("../data/mcidas/goes.gif")
xsize = a.getDomainSet().getX().getLength()
ysize = a.getDomainSet().getY().getLength()
dom=Linear2DSet(RealTupleType.SpatialEarth2DTuple,-110,-70,xsize, 53,33,ysize)
typ = MathType.stringToType("((Longitude,Latitude)->(red,green,blue))")
ff = FlatField(typ,dom)
ff.setSamples(a.getFloats())
plot(ff)
map=load("../data/mcidas/OUTLUSAM")
plot(map)
You can download the latest Jython release from:
http://jython.sourceforge.net/
tom
--
Tom Whittaker (tomw@xxxxxxxxxxxxx)
University of Wisconsin-Madison
Space Science and Engineering Center
Phone/VoiceMail: 608/262-2759
Fax: 608/262-5974
