Thanks for your quick response. I'm sorry, though- but I'm still
confused. I think perhaps I did not make it clear what I want to do. I
just want to take example P2_03 on the visAD tutorial page:
http://www.ssec.wisc.edu/~billh/tutorial/s2/Section2.html
and draw the sample points as bars instead of dots. The top of each bar
would be at the y (height) value of the point and the bottom of each bar
would be at y = 0. All I meant about the "midpoint" stuff was that if
barwidth = 1 then 0.5 of the bar would be to the left of the x (time)
position and 0.5 of the bar would be to the right of the x position. In
other words, the upper left corner of each bar would be at coordinate (x
- 0.5, y) and the lower right corner of the bar shape would be at (x +
0.5, 0). In this case, I don't think we need the third RealType ("dif")
in your example. Right?
I also didn't understand the section of code where you set up the
shapes. Why do you have to generate the Linear1DSet of samples when you
already have all the sample data you need? Which leads to a more general
question- is VisAD not a good program for plotting/charting discrete
data? Is it much better suited for continuous functions?
Anyway, here is my modified version of example P2_03 as I attempt to
draw the points as bars. The only section that's different from the one
on the web page is the "set up shapes" section (and the fact that I made
the heightMap range start at 0 instead of -10, which shouldn't matter
much). I guess my biggest problem is that I don't know how to get the
bar shapes to draw. I know I'm missing code to do that, but I'm not
quite sure what it should be. I don't even understand how the example in
the tutorial knew to draw dots as opposed to a line connecting the
points in the previous tutorial example. Any help is greatly
appreciated! Here is my code, followed by your prev. email for
reference. You may want to just scroll down to my "set up shapes"
section. THANK YOU!!
/*
VisAD Tutorial
Copyright (C) 2000 Ugo Taddei
*/
// Import needed classes
import visad.*;
import visad.java2d.DisplayImplJ2D;
import java.rmi.RemoteException;
import java.awt.*;
import javax.swing.*;
/**
VisAD Tutorial example 2_03
Data are organized as MathType ( index -> ( time, height ) ) and
represent some points from the parabola of the previous example
Data are indexed (time, height) points and get depicted as such.
Run program with java P2_03
*/
public class P2_03Bar{
// Declare variables
// The quantities to be displayed in x- and y-axes: time and height,
respectively
// Our index is alos a RealType
private RealType time, height, index;
// A Tuple, to pack time and height together
private RealTupleType t_h_tuple;
// The function ( time(i), height(i) ), where i = index,
// represented by ( index -> ( time, height) )
// ( time, height) are a Tuple, so we have a FunctionType
// from index to a tuple
private FunctionType func_i_tuple;
// Our Data values, the points, are now indexed by the Set
private Set index_set;
// The Data class FlatField, which will hold time and height data.
// time data are implicitely given by the Set time_set
private FlatField vals_ff;
// The DataReference from the data to display
private DataReferenceImpl data_ref;
// The 2D display, and its the maps
private DisplayImpl display;
private ScalarMap timeMap, heightMap;
public P2_03Bar (String []args)
throws RemoteException, VisADException {
// Create the quantities
// x and y are measured in SI meters
// Use RealType(String name, Unit u, Set set), set is null
time = new RealType("time", SI.second, null);
height = new RealType("height", SI.meter, null);
// Organize time and height in a Tuple
t_h_tuple = new RealTupleType( time, height);
// Index has no unit, just a name
index = new RealType("index");
// Create a FunctionType ( index -> ( time, height) )
// Use FunctionType(MathType domain, MathType range)
func_i_tuple = new FunctionType( index, t_h_tuple);
// Create the x_set, with 5 values, but this time using a
// Integer1DSet(MathType type, int length)
index_set = new Integer1DSet(index, 5);
// These are our actual data values for time and height
// Note that these values correspond to the parabola of the
// previous examples. The y (height) values are the same, but the x
(time)
// are now given given.
float[][] point_vals = new float[][]{{-3.0f, -1.5f, 0.0f, 1.5f,
3.0f,},
{0.0f, 33.75f, 45.0f, 33.75f, 0.0f,} };
// Create a FlatField, that is the Data class for the samples
// Use FlatField(FunctionType type, Set domain_set)
vals_ff = new FlatField( func_i_tuple, index_set);
// and put the height values above in it
vals_ff.setSamples( point_vals );
// Create Display and its maps
// A 2D display
display = new DisplayImplJ2D("display1");
// Get display's graphics mode control and draw scales
GraphicsModeControl dispGMC = (GraphicsModeControl)
display.getGraphicsModeControl();
dispGMC.setScaleEnable(true);
// Create the ScalarMaps: quantity time is to be displayed along
XAxis
// and height along YAxis
// Use ScalarMap(ScalarType scalar, DisplayRealType display_scalar)
timeMap = new ScalarMap( time, Display.XAxis );
heightMap = new ScalarMap( height, Display.YAxis );
// Add maps to display
display.addMap( timeMap );
display.addMap( heightMap );
// Scale heightMap. This will scale the y-axis, because heightMap
has DisplayRealType YAXIS
// We simply choose the range from -4 to 4 for the x-axis
// and -10.0 to 50.0 for
timeMap.setRange( -4.0, 4.0);
heightMap.setRange( 0.0, 50.0);
// Create a data reference and set the FlatField as our data
data_ref = new DataReferenceImpl("data_ref");
data_ref.setData( vals_ff );
///////////////////// set up shapes /////////////////////
// can't do this- get class cast exception because it's not a
// ShapeControl, it's a visad.java2d.ProjectionControlJ2D.
// So how do I get the control?
//ShapeControl control = (ShapeControl) heightMap.getControl();
float barWidth = 1f;
int nPts = point_vals[1].length;
VisADQuadArray[] shapes = new VisADQuadArray[nPts];
for (int i=0; i<nPts; i++) {
shapes[i] = new VisADQuadArray();
shapes[i].vertexCount = 4;
float hw = 0.5f * barWidth;
float x = point_vals[0][i];
float y = point_vals[1][i];
shapes[i].coordinates = new float[]
{x - hw, y, 0.0f, x + hw, y, 0.0f,
x + hw, y, 0.0f, x - hw, y, 0.0f};
}
// how do I draw these shapes? Do I just call this once I get
// a valid ShapeControl object? :
// control.setShapes(shape)
///////////////////// end set up shapes /////////////////////
// Add reference to display
display.addReference( data_ref );
// Create application window, put display into it
JFrame jframe = new JFrame("VisAD Tutorial example 2_03Bar");
jframe.getContentPane().add(display.getComponent());
// Set window size and make it visible
jframe.setSize(300, 300);
jframe.setVisible(true);
}
public static void main(String[] args)
throws RemoteException, VisADException
{
new P2_03Bar(args);
}
} //end of Visad Tutorial example 2_03
"From: Bill Hibbard " wrote:
>
> Hi Jen,
>
> > Hi. I am new to VisAD but it seems really great. Problem is, it's so
> > powerful that I'm having trouble figuring out the simple stuff. I've
> > been looking at the (very helpful) tutorial at:
> >
> > http://www.ssec.wisc.edu/~billh/tutorial/s2/Section2.html
> >
> > I want to do a plot similar to the one in example P2_03 or P2_04, but
> > insead of points to represent the data I would like bars. I want to make
> > a bar graph, where the bottom of each bar is at a certain point on the x
> > axis (in this example, maybe the midpoint of each bar would be at its
> > corresponding x value) and the top of each bar is at the height of the
> > corresponding y value. Do I have to make a VisADQuadArray for every bar
> > I want or is there a simpler solution? Does anybody have/know of a
> > simple Bar Chart example in 2d or 3d using VisAD?
>
> You can probably do this using ScalarMaps to Display.Shape.
> Say you have a FlatField 'field' with MathType
> (index -> (mid, top)) and you want the midpoint of bars at
> mid values and tops of bars at top values.
>
> You can probably do this using ScalarMaps to Display.Shape.
> Say you have a FlatField 'field' with MathType
> (index -> (mid, top)) and you want the midpoint of bars at
> mid values and tops of bars at top values. Then do
> something like:
>
> // create new_field with (top - mid) values
> RealType dif = new RealType("dif");
> RealTupleType range = new RealTupleType(mid, top, dif);
> FunctionType new_func = new FunctionType(index, range);
> Set set = field.getDomainSet();
> int length = set.getLength();
> float[][] old_values = field.getFloats();
> float[][] new_values = new float[3][length];
> new_values[0] = old_values[0];
> new_values[1] = old_values[1];
> float min_mid = Float.MAXIMUM_VALUE;
> float max_mid = Float.MINIMUM_VALUE;
> float max_dif = Float.MINIMUM_VALUE;
> for (int i=0; i<length; i++) {
> if (new_values[0][i] < min_mid) min_mid = new_values[0][i];
> if (new_values[0][i] > max_mid) max_mid = new_values[0][i];
> new_values[2][i] = new_values[1][i] - new_values[0][i];
> if (new_values[2][i] > max_mid) max_mid = new_values[2][i];
> }
> FlatField new_field = new FlatField(new_func, set);
> new_field.setSamples(new_values);
>
> // some tuning values for display
> float dif_resolution = ...; // desired resolution of dif values
> int ndifs = 1 + (int) (max_dif / dif_resolution); // # of dif samples
> float width = ...; // width of bar (e.g., 0.02f)
>
> // set up Scalarmaps for display
> display.addMap(new ScalarMap(index, Display.XAxis));
> ScalarMap mid_map = new ScalarMap(mid, Display.YAxis);
> display.addMap(mid_map);
> mid_map.setRange(min_mid, max_mid);
> display.addMap(new ScalarMap(mid, Display.YAxis));
> ScalarMap dif_map = new ScalarMap(dif, Display.Shape);
> display.addMap(dif_map);
>
> // set up shapes
> ShapeControl control = (ShapeControl) dif_map.getControl();
> Linear1DSet dif_set = new Linear1DSet(dif, 0.0, max_dif, ndifs);
> float[][] difs = dif_set.getSamples();
> VisADQuadArray[] dif_shapes = new VisADQuadArray[ndifs];
> for (int i=0; i<ndifs; i++) {
> dif_shapes[i] = new VisADQuadArray();
> dif_shapes[i].vertexCount = 4;
> float hw = 0.5f * width;
> float hh = 0.5f * difs[0][i] * (2.0f / (max_mid - min_mid));
> dif_shapes[i].coordinates = new float[]
> {hw, hh, 0.0f, -hw, hh, 0.0f, -hw, -hh, 0.0f, hw, -hh, 0.0f};
> }
> control.setShapeSet(dif_set);
> control.setShapes(dif_shapes);
>
> I haven't tested this, so it may have bugs, but its the general
> idea.
>
> Cheers,
> Bill
> ----------------------------------------------------------
> Bill Hibbard, SSEC, 1225 W. Dayton St., Madison, WI 53706
> hibbard@xxxxxxxxxxxxxxxxx 608-263-4427 fax: 608-263-6738
> http://www.ssec.wisc.edu/~billh/vis.html
