Hi Stuart-
Stuart Maclean wrote:
Hi Don, you are so right, my code was a comedy of errors. I should have
RTFM for GeoGrids. What gives me a nice contoured image from a nc
Variable is this:
String netcdfUrl = args[0];
GeoGridDataSource ggds = new GeoGridDataSource(null, netcdfUrl, null);
GridDataset gds = ggds.getDataset();
// get the geogrid by name, as Don says
GeoGrid gg = gds.findGridByName( "surface_temperature" );
GeoGridAdapter gga = new GeoGridAdapter( ggds, gg );
FieldImpl fi = gga.getData();
// sample 1 is 2nd time value, in my data set 1st time sample is
missing for this variable
FlatField ff = (FlatField)fi.getSample(1);
// now build the Displayable
Contour2DDisplayable cdd = new Contour2DDisplayable( "foo" );
// we have to set the realType, I was missing this
cdd.setC2DRealType(
(RealType)((FunctionType)ff.getType()).getRange() );
cdd.setData( fi );
cdd.setContourLevels( new RegularContourLevels( 5, 273, 273, 310 ) );
My biggest error was missing the line to setC2DRealType.
Hmmm, I didn't even know that method existed. ;-) Jeff must have
added that.
i was under
the impression that the IDV code could 'figure out' the RealType from
the FieldImpl passed in. Of ocurse this cannot be the case since the
FieldImpl could have many ranges, ie it could be realTupleType. Maybe
if the ((FunctionType)ff.getType()).getRange() call returns a RealType
that could be used automatically?
If you call loadData instead of setData, it will set the RealTypes
automagically. setData used to be a protected method (probably still
should be) so that each of the Displayables could do what is necessary
in a helper method like loadData() that would eventually call
setData. All setData does is set the data in the VisAD DataReference
(essentially).
I've updated the javadoc for Contour2DDisplayable to reflect this.
setGrid2D erroneously said to call setData in its deprecated tag.
Don
On Apr 4, 2008, at 7:22 AM, Don Murray wrote:
Hi Stuart-
Just a followup to this since I'm now looking at your first
message. In your code you have:
Stuart Maclean wrote:
I am now trying to contour a field from a netcdf grid file. The IDV
does this Ok, but my image is simply a white polygon of the grid
itself. My code is this:
String netcdfUrl = args[0];
GeoGridDataSource ggds = new GeoGridDataSource(null, netcdfUrl,
null);
GridDataset gds = ggds.getDataset();
GeoGrid gg = (GeoGrid)gds.getGrids().get(0);
This will just get the first grid in the GridDataset. What you want
to do is:
GeoGrid gg = gds.findGridByName("surface_temperature");
(or substitute whatever the variable name is for "surface_temperature").
GeoGridAdapter gga = new GeoGridAdapter( ggds, gg,
"surface_temperature" );
FieldImpl fi = gga.getData();
System.out.println( fi.getType().prettyString() );
which gives
(Time -> ((XC[unit:km], YC[unit:km]) -> surface_temperature[unit:Pa]))
Not sure why temp units are Pa??
The first grid in the GridDataset must be a pressure grid of
some sort.
FlatField ff = (FlatField)fi.getSample(1);
System.out.println( ff.getType().prettyString() );
which I presume extracts some time element, and gives
((XC[unit:km], YC[unit:km]) -> surface_temperature[unit:Pa])
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
*************************************************************
--
*************************************************************
Don Murray UCAR Unidata Program
dmurray@xxxxxxxxxxxxxxxx P.O. Box 3000
(303) 497-8628 Boulder, CO 80307
http://www.unidata.ucar.edu/staff/donm
*************************************************************
