NOTICE: This version of the NSF Unidata web site (archive.unidata.ucar.edu) is no longer being updated.
Current content can be found at unidata.ucar.edu.
To learn about what's going on, see About the Archive Site.
NOTICE: This version of the NSF Unidata web site (archive.unidata.ucar.edu) is no longer being updated.
Current content can be found at unidata.ucar.edu.
To learn about what's going on, see About the Archive Site.
Silvere-
I have been at http://www.ssec.wisc.edu/~billh/globe.gif and found exactly what I am looking for. From the refering page from Visad "See the visad/examples/SatDisplay.java program for an example of how to use this capability." I haven't had the time yet, I admit. What I would like to do is a "best resolution" real time global cloud cover with Java3D. I can't use the Global Montage http://www.ssec.wisc.edu/data/comp/latest_moll.gif since the projection system seems difficult to use in Java3D and resolution as well as continental contour don't fit my needs. But it appears you are using the kind of data I am looking for.
The globe.gif image is from an application I wrote a long time ago so I guess I'm best qualified to answer. The latest_moll.gif image is generated from a satellite composite in McIDAS AREA file format with a Mollweide projection that has been displayed in McIDAS. The globe.gif image is a display in VisAD of an AREA file similar to the one in latest_moll.gif. VisAD supports most of the McIDAS navigation modules, so if your data were stored in a McIDAS AREA file, you could use SatDisplay to map it to a globe. If all you have is a GIF or JPEG, mapping to a sphere is going to be a little difficult. If your image is rectillinear and you know the upper left corner of the image, you could use something like the attached program to do this. This image will add a CoordinateSystem to the domain of the image that will transform image coordinates to lat/lon. You can test this program by compiling it and running: java -Xmx128m WorldImageAdapter ftp://ftp.unidata.ucar.edu/pub/dmurray/bigmap_541x271.jpg (you'll need visad.jar and . in your classpath to compile and run) The image is a global topography image that covers -180 to 180 lon and 90 to -90 lat. You could adjust the code accordingly if your domain is smaller. Any non-rectillinear image would require a different CoordinateSystem. Good luck. Don ************************************************************* Don Murray UCAR Unidata Program dmurray@xxxxxxxxxxxxxxxx P.O. Box 3000 (303) 497-8628 Boulder, CO 80307 http://www.unidata.ucar.edu/staff/donm *************************************************************
import java.awt.event.*;
import javax.swing.*;
import visad.*;
import visad.bom.*;
import visad.data.gif.*;
import visad.data.jai.*;
import visad.data.tiff.*;
import visad.data.DefaultFamily;
import visad.data.mcidas.BaseMapAdapter;
import visad.data.netcdf.Plain;
import visad.georef.*;
import visad.java2d.*;
import visad.java3d.*;
import visad.jmet.GRIBCoordinateSystem;
public class WorldImageAdapter {
FlatField newField = null;
public WorldImageAdapter(String source)
throws Exception {
DefaultFamily df = new DefaultFamily("images", true);
FlatField image = (FlatField) df.open(source);
FunctionType fType = (FunctionType) image.getType();
Linear2DSet domainSet = (Linear2DSet) image.getDomainSet();
int neles = domainSet.getLengths()[0];
int nlines = domainSet.getLengths()[1];
double nlat = 180./(nlines-1);
double nlon = 360./(neles-1);
MapProjection mp =
new GRIBCoordinateSystem(0,
neles,
nlines,
-90.0,
-180.0,
90.0,
180.0,
nlon,
nlat);
RealTupleType domainType = fType.getDomain();
RealTupleType rangeType = fType.getFlatRange();
RealTupleType newDomainType =
new RealTupleType((RealType)domainType.getComponent(0),
(RealType)domainType.getComponent(1), mp, null);
Linear2DSet newDomain =
new Linear2DSet(newDomainType, 0, neles-1, neles,
nlines - 1, 0.0, nlines);
FunctionType newFieldType =
new FunctionType(newDomainType, rangeType);
newField = new FlatField(newFieldType, newDomain);
newField.setSamples(image.getFloats(false), false);
}
public FlatField getData() {
return newField;
}
public static void main(String[] args)
throws Exception {
if (args.length == 0) {
System.out.println("Must supply a filename");
System.exit(1);
}
WorldImageAdapter wia = new WorldImageAdapter(args[0]);
final FlatField data = wia.getData();
RealTupleType rangeType =
((FunctionType)data.getType()).getFlatRange();
DisplayImpl display =
/*
new DisplayImplJ2D("display");
new DisplayImplJ3D("display", new TwoDDisplayRendererJ3D());
*/
new DisplayImplJ3D("display");
/*
GraphicsModeControl gmc = display.getGraphicsModeControl();
gmc.setTextureEnable(false);
*/
//display.addMap(new ScalarMap(RealType.Latitude, Display.YAxis));
//display.addMap(new ScalarMap(RealType.Longitude, Display.XAxis));
display.addMap(new ScalarMap(RealType.Latitude, Display.Latitude));
display.addMap(new ScalarMap(RealType.Longitude, Display.Longitude));
display.addMap(
new ScalarMap((RealType)rangeType.getComponent(0), Display.Red));
display.addMap(
new ScalarMap((RealType)rangeType.getComponent(1), Display.Green));
display.addMap(
new ScalarMap((RealType)rangeType.getComponent(2), Display.Blue));
display.getGraphicsModeControl().setScaleEnable(true);
DataReference ref = new DataReferenceImpl("ref");
ref.setData(data);
display.addReference(ref);
/*
try {
BaseMapAdapter bma = new BaseMapAdapter("OUTLSUPW");
DataReference mapRef = new DataReferenceImpl("map");
mapRef.setData(bma.getData());
display.addReference(mapRef);
} catch (Exception e) {
System.out.println("Map OUTLSUPW not available");
}
*/
JFrame frame = new JFrame("adapter test");
frame.addWindowListener(new WindowAdapter() {
public void windowClosing(WindowEvent e) {
System.exit(0);
}
});
frame.getContentPane().add(display.getComponent());
JButton save = new JButton("Save");
save.addActionListener(new ActionListener() {
public void actionPerformed(ActionEvent e) {
try {
new Plain().save("world.nc", data, true);
} catch (Exception exp) {
System.out.println("couldn't save image");
}
}
});
frame.getContentPane().add("South", save);
frame.pack();
frame.show();
}
}
visad archives: