Re: Elevation surfaces -- Cutting planes

Hi Mohamad and Ugo,

> The idea is to be able to depict a movable plane that sweeps the
> DEM and sends the data of the intersection curve to the 2D display.

If the DEM is "FlatField elevation" with MathType ((x, y) -> z) then
to make a vertical cutting plane, construct "Gridded2DSet cutting_set"
with MathType Set(x, y) and manifold dimension = 1, whose samples
lie along the line in the (x, y) plane that is the intersection of
the (x, y) plane with the vertical cutting plane.  Then:

  FlatField cut_elevation = (FlatField)
    elevation.resample(cutting_set, Data.WEIGHTED_AVERAGE, Data.NO_ERRORS);

and cut_elevation will be displayed in your 3-D DEM display as the
curve that is the intersection of your DEM with the vertical cutting
plane.  For your 2-D display, construct another RealType w that
represents location along the line in the (x, y) plane and construct
"Gridded1DSet dem2d_samples" with MathType Set(w) and whose samples are
the w locations of the samples of cutting_set.  Then:

  FunctionType dem2d_type = new FunctionType(w, z);
  FlatField dem_2d = newFlatField(dem2d_type, dem2d_samples);
  double[][] values = cut_elevation.getValues();

Now you can display dem_2d in 2-D with the ScalarMaps:

  ScalarMap(w, Display.XAxis);
  ScalarMap(z, Display.YAxis);

If your DEM display includes 3-D FlatFields [i.e., MathType
((x, y, z) -> attribute)] representing sub-surface data, then you can do
something similar by defining a "Gridded3DSet cutting_plane" with MathType
Set(x, y, z) and manifold dimension = 2, whose samples are a 2-D grid lying
on the cutting plane.  Then resample the 3-D FlatFields to cutting_plane,

Bill Hibbard, SSEC, 1225 W. Dayton St., Madison, WI  53706
hibbard@xxxxxxxxxxxxxxxxx  608-263-4427  fax: 608-263-6738