Due to the current gap in continued funding from the U.S. National Science Foundation (NSF), the NSF Unidata Program Center has temporarily paused most operations. See NSF Unidata Pause in Most Operations for details.
Due to the current gap in continued funding from the U.S. National Science Foundation (NSF), the NSF Unidata Program Center has temporarily paused most operations. See NSF Unidata Pause in Most Operations for details.
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
visad archives: