If I'm understanding you right, you want the surface of the hexahedra that
contain the "isovalue". If the data is connections dependent, then Include
will do it for you - max = min = "isovalue" to get the hexes, then
ShowBoundary to get their surfaces. Note that you want to use FaceNormals
to assign surface normals so it doesn't try to interpret the hexes as poor
approximations for spheres. If the data is positions dependent, then its a
little harder - "show me all the hex that have at least one vertex on
either side of the "isovalue". You can use Compute and Post to do that, as
in the following example. The three modules on the top left fake some
positions dependent data. The next Compute assigns 1 to each vertex >
"isovalue" and 0 to each <= "isovalue". The Post averages the values at
each vertex of a hex to get a value for the hex itself, which will be > 0
and < 1 for those hex with at least one vertex on either side of the
"isovalue". The next Compute assigns 1 to them, and 0 otherwise. The
Include selects the 1 hexes, and you have your "surface".
-------------------------------- CUT HERE
--------------------------------------
//
// time: Tue Feb 26 10:35:28 2002
//
// version: 3.1.2 (format), 4.1.3 (DX)
//
//
// MODULE main
// workspace: width = 426, height = 559
// layout: snap = 0, width = 50, height = 50, align = NN
//
macro main(
) -> (
) {
//
// node Construct[1]: x = 0, y = 6, inputs = 4, label = Construct
// input[1]: defaulting = 0, visible = 1, type = 16777224, value = {[-1 -1
-1]}
// input[2]: defaulting = 0, visible = 1, type = 16777224, value = {[.1 .1
.1]}
// input[3]: defaulting = 0, visible = 1, type = 8, value = [21 21 21]
//
main_Construct_1_out_1 =
Construct(
main_Construct_1_in_1,
main_Construct_1_in_2,
main_Construct_1_in_3,
main_Construct_1_in_4
) [instance: 1, cache: 1];
//
// node Replace[1]: x = 1, y = 76, inputs = 4, label = Replace
// input[3]: defaulting = 0, visible = 1, type = 32, value = "positions"
//
main_Replace_1_out_1 =
Replace(
main_Construct_1_out_1,
main_Construct_1_out_1,
main_Replace_1_in_3,
main_Replace_1_in_4
) [instance: 1, cache: 1];
//
// node Compute[1]: x = 10, y = 152, inputs = 2, label = Compute
// input[1]: defaulting = 0, visible = 0, type = 32, value =
"mag(abs($0)-0.5)"
// expression: value = mag(abs(a)-0.5)
// name[2]: value = a
//
main_Compute_1_out_1 =
Compute(
main_Compute_1_in_1,
main_Replace_1_out_1
) [instance: 1, cache: 1];
//
// node ShowBox[1]: x = 0, y = 276, inputs = 1, label = ShowBox
//
main_ShowBox_1_out_1,
main_ShowBox_1_out_2 =
ShowBox(
main_Compute_1_out_1
) [instance: 1, cache: 1];
//
// node Scalar[1]: x = 307, y = 16, inputs = 11, label = Scalar
// input[1]: defaulting = 0, visible = 0, type = 32, value = "Scalar_1"
// input[3]: defaulting = 0, visible = 0, type = 5, value = 0.7
// input[5]: defaulting = 1, visible = 0, type = 5, value = -1000000.0
// input[6]: defaulting = 1, visible = 0, type = 5, value = 1000000.0
// input[7]: defaulting = 1, visible = 0, type = 5, value = 1.0
// input[9]: defaulting = 1, visible = 0, type = 1, value = 5
// output[1]: visible = 1, type = 5, value = 0.7
//
//
// node Compute[3]: x = 175, y = 104, inputs = 3, label = Compute
// input[1]: defaulting = 0, visible = 0, type = 32, value = "$0 > $1 ? 1.0
: 0.0"
// expression: value = a > b ? 1.0 : 0.0
// name[2]: value = a
// name[3]: value = b
//
main_Compute_3_out_1 =
Compute(
main_Compute_3_in_1,
main_Compute_1_out_1,
main_Scalar_1_out_1
) [instance: 3, cache: 1];
//
// node Post[1]: x = 275, y = 118, inputs = 2, label = Post
// input[2]: defaulting = 0, visible = 1, type = 32, value = "connections"
//
main_Post_1_out_1 =
Post(
main_Compute_3_out_1,
main_Post_1_in_2
) [instance: 1, cache: 1];
//
// node Compute[4]: x = 337, y = 165, inputs = 3, label = Compute
// input[1]: defaulting = 0, visible = 0, type = 32, value = "$0 >= 1 || $0
<= 0 ? 0 : 1"
// expression: value = a >= 1 || a <= 0 ? 0 : 1
// name[2]: value = a
// name[3]: value = b
//
main_Compute_4_out_1 =
Compute(
main_Compute_4_in_1,
main_Post_1_out_1,
main_Compute_4_in_3
) [instance: 4, cache: 1];
//
// node Include[1]: x = 282, y = 252, inputs = 6, label = Include
// input[2]: defaulting = 0, visible = 1, type = 5, value = 1.0
// input[3]: defaulting = 0, visible = 1, type = 5, value = 1.0
// input[4]: defaulting = 1, visible = 1, type = 3, value = NULL
//
main_Include_1_out_1 =
Include(
main_Compute_4_out_1,
main_Include_1_in_2,
main_Include_1_in_3,
main_Include_1_in_4,
main_Include_1_in_5,
main_Include_1_in_6
) [instance: 1, cache: 1];
//
// node ShowBoundary[1]: x = 298, y = 335, inputs = 2, label = ShowBoundary
//
main_ShowBoundary_1_out_1 =
ShowBoundary(
main_Include_1_out_1,
main_ShowBoundary_1_in_2
) [instance: 1, cache: 1];
//
// node FaceNormals[1]: x = 264, y = 415, inputs = 1, label = FaceNormals
//
main_FaceNormals_1_out_1 =
FaceNormals(
main_ShowBoundary_1_out_1
) [instance: 1, cache: 1];
//
// node Collect[1]: x = 27, y = 392, inputs = 2, label = Collect
//
main_Collect_1_out_1 =
Collect(
main_ShowBox_1_out_1,
main_FaceNormals_1_out_1
) [instance: 1, cache: 1];
//
// node Image[1]: x = 38, y = 497, inputs = 49, label = Image
// input[1]: defaulting = 0, visible = 0, type = 67108863, value = "Image_1"
// input[4]: defaulting = 0, visible = 0, type = 1, value = 1
// input[5]: defaulting = 0, visible = 0, type = 8, value = [0 0 0]
// input[6]: defaulting = 0, visible = 0, type = 8, value = [6.50456
4.48337 3.65602]
// input[7]: defaulting = 0, visible = 0, type = 5, value = 4.66499
// input[8]: defaulting = 0, visible = 0, type = 1, value = 640
// input[9]: defaulting = 0, visible = 0, type = 5, value = 0.751
// input[10]: defaulting = 0, visible = 0, type = 8, value = [-0.159773
0.752661 -0.638728]
// input[11]: defaulting = 1, visible = 0, type = 5, value = 30.0001
// input[12]: defaulting = 0, visible = 0, type = 1, value = 0
// input[14]: defaulting = 0, visible = 0, type = 1, value = 1
// input[15]: defaulting = 1, visible = 0, type = 32, value = "none"
// input[16]: defaulting = 1, visible = 0, type = 32, value = "none"
// input[17]: defaulting = 1, visible = 0, type = 1, value = 1
// input[18]: defaulting = 1, visible = 0, type = 1, value = 1
// input[19]: defaulting = 0, visible = 0, type = 1, value = 0
// input[29]: defaulting = 1, visible = 0, type = 3, value = 0
// input[41]: defaulting = 0, visible = 0, type = 32, value = "rotate"
// depth: value = 24
// window: position = (0.4609,0.4219), size = 0.5109x0.5117
// internal caching: 1
//
main_Image_1_out_1,
main_Image_1_out_2,
main_Image_1_out_3 =
Image(
main_Image_1_in_1,
main_Collect_1_out_1,
main_Image_1_in_3,
main_Image_1_in_4,
main_Image_1_in_5,
main_Image_1_in_6,
main_Image_1_in_7,
main_Image_1_in_8,
main_Image_1_in_9,
main_Image_1_in_10,
main_Image_1_in_11,
main_Image_1_in_12,
main_Image_1_in_13,
main_Image_1_in_14,
main_Image_1_in_15,
main_Image_1_in_16,
main_Image_1_in_17,
main_Image_1_in_18,
main_Image_1_in_19,
main_Image_1_in_20,
main_Image_1_in_21,
main_Image_1_in_22,
main_Image_1_in_23,
main_Image_1_in_24,
main_Image_1_in_25,
main_Image_1_in_26,
main_Image_1_in_27,
main_Image_1_in_28,
main_Image_1_in_29,
main_Image_1_in_30,
main_Image_1_in_31,
main_Image_1_in_32,
main_Image_1_in_33,
main_Image_1_in_34,
main_Image_1_in_35,
main_Image_1_in_36,
main_Image_1_in_37,
main_Image_1_in_38,
main_Image_1_in_39,
main_Image_1_in_40,
main_Image_1_in_41,
main_Image_1_in_42,
main_Image_1_in_43,
main_Image_1_in_44,
main_Image_1_in_45,
main_Image_1_in_46,
main_Image_1_in_47,
main_Image_1_in_48,
main_Image_1_in_49
) [instance: 1, cache: 1];
// network: end of macro body
CacheScene(main_Image_1_in_1, main_Image_1_out_1, main_Image_1_out_2);
}
main_Construct_1_in_1 = {[-1 -1 -1]};
main_Construct_1_in_2 = {[.1 .1 .1]};
main_Construct_1_in_3 = [21 21 21];
main_Construct_1_in_4 = NULL;
main_Construct_1_out_1 = NULL;
main_Replace_1_in_3 = "positions";
main_Replace_1_in_4 = NULL;
main_Replace_1_out_1 = NULL;
main_Compute_1_in_1 = "mag(abs($0)-0.5)";
main_Compute_1_out_1 = NULL;
main_ShowBox_1_out_1 = NULL;
main_Scalar_1_in_1 = "Scalar_1";
main_Scalar_1_in_2 = NULL;
main_Scalar_1_in_3 = 0.7 ;
main_Scalar_1_in_4 = NULL;
main_Scalar_1_in_5 = NULL;
main_Scalar_1_in_6 = NULL;
main_Scalar_1_in_7 = NULL;
main_Scalar_1_in_8 = NULL;
main_Scalar_1_in_9 = NULL;
main_Scalar_1_in_10 = NULL;
main_Scalar_1_in_11 = NULL;
main_Scalar_1_out_1 = 0.7 ;
main_Compute_3_in_1 = "$0 > $1 ? 1.0 : 0.0";
main_Compute_3_out_1 = NULL;
main_Post_1_in_2 = "connections";
main_Post_1_out_1 = NULL;
main_Compute_4_in_1 = "$0 >= 1 || $0 <= 0 ? 0 : 1";
main_Compute_4_in_3 = NULL;
main_Compute_4_out_1 = NULL;
main_Include_1_in_2 = 1.0;
main_Include_1_in_3 = 1.0;
main_Include_1_in_4 = NULL;
main_Include_1_in_5 = NULL;
main_Include_1_in_6 = NULL;
main_Include_1_out_1 = NULL;
main_ShowBoundary_1_in_2 = NULL;
main_ShowBoundary_1_out_1 = NULL;
main_FaceNormals_1_out_1 = NULL;
main_Collect_1_out_1 = NULL;
macro Image(
id,
object,
where,
useVector,
to,
from,
width,
resolution,
aspect,
up,
viewAngle,
perspective,
options,
buttonState = 1,
buttonUpApprox = "none",
buttonDownApprox = "none",
buttonUpDensity = 1,
buttonDownDensity = 1,
renderMode = 0,
defaultCamera,
reset,
backgroundColor,
throttle,
RECenable = 0,
RECfile,
RECformat,
RECresolution,
RECaspect,
AAenable = 0,
AAlabels,
AAticks,
AAcorners,
AAframe,
AAadjust,
AAcursor,
AAgrid,
AAcolors,
AAannotation,
AAlabelscale,
AAfont,
interactionMode,
title,
AAxTickLocs,
AAyTickLocs,
AAzTickLocs,
AAxTickLabels,
AAyTickLabels,
AAzTickLabels,
webOptions) -> (
object,
camera,
where)
{
ImageMessage(
id,
backgroundColor,
throttle,
RECenable,
RECfile,
RECformat,
RECresolution,
RECaspect,
AAenable,
AAlabels,
AAticks,
AAcorners,
AAframe,
AAadjust,
AAcursor,
AAgrid,
AAcolors,
AAannotation,
AAlabelscale,
AAfont,
AAxTickLocs,
AAyTickLocs,
AAzTickLocs,
AAxTickLabels,
AAyTickLabels,
AAzTickLabels,
interactionMode,
title,
renderMode,
buttonUpApprox,
buttonDownApprox,
buttonUpDensity,
buttonDownDensity) [instance: 1, cache: 1];
autoCamera =
AutoCamera(
object,
"front",
object,
resolution,
aspect,
[0,1,0],
perspective,
viewAngle,
backgroundColor) [instance: 1, cache: 1];
realCamera =
Camera(
to,
from,
width,
resolution,
aspect,
up,
perspective,
viewAngle,
backgroundColor) [instance: 1, cache: 1];
coloredDefaultCamera =
UpdateCamera(defaultCamera,
background=backgroundColor) [instance: 1, cache: 1];
nullDefaultCamera =
Inquire(defaultCamera,
"is null + 1") [instance: 1, cache: 1];
resetCamera =
Switch(
nullDefaultCamera,
coloredDefaultCamera,
autoCamera) [instance: 1, cache: 1];
resetNull =
Inquire(
reset,
"is null + 1") [instance: 2, cache: 1];
reset =
Switch(
resetNull,
reset,
0) [instance: 2, cache: 1];
whichCamera =
Compute(
"($0 != 0 || $1 == 0) ? 1 : 2",
reset,
useVector) [instance: 1, cache: 1];
camera = Switch(
whichCamera,
resetCamera,
realCamera) [instance: 3, cache: 1];
AAobject =
AutoAxes(
object,
camera,
AAlabels,
AAticks,
AAcorners,
AAframe,
AAadjust,
AAcursor,
AAgrid,
AAcolors,
AAannotation,
AAlabelscale,
AAfont,
AAxTickLocs,
AAyTickLocs,
AAzTickLocs,
AAxTickLabels,
AAyTickLabels,
AAzTickLabels) [instance: 1, cache: 1];
switchAAenable = Compute("$0+1",
AAenable) [instance: 2, cache: 1];
object = Switch(
switchAAenable,
object,
AAobject) [instance:4, cache: 1];
SWapproximation_options =
Switch(
buttonState,
buttonUpApprox,
buttonDownApprox) [instance: 5, cache: 1];
SWdensity_options =
Switch(
buttonState,
buttonUpDensity,
buttonDownDensity) [instance: 6, cache: 1];
HWapproximation_options =
Format(
"%s,%s",
buttonDownApprox,
buttonUpApprox) [instance: 1, cache: 1];
HWdensity_options =
Format(
"%d,%d",
buttonDownDensity,
buttonUpDensity) [instance: 2, cache: 1];
switchRenderMode = Compute(
"$0+1",
renderMode) [instance: 3, cache: 1];
approximation_options = Switch(
switchRenderMode,
SWapproximation_options,
HWapproximation_options) [instance: 7, cache: 1];
density_options = Switch(
switchRenderMode,
SWdensity_options,
HWdensity_options) [instance: 8, cache: 1];
renderModeString = Switch(
switchRenderMode,
"software",
"hardware")[instance: 9, cache: 1];
object_tag = Inquire(
object,
"object tag")[instance: 3, cache: 1];
annoted_object =
Options(
object,
"send boxes",
0,
"cache",
1,
"object tag",
object_tag,
"ddcamera",
whichCamera,
"rendering approximation",
approximation_options,
"render every",
density_options,
"button state",
buttonState,
"rendering mode",
renderModeString) [instance: 1, cache: 1];
RECresNull =
Inquire(
RECresolution,
"is null + 1") [instance: 4, cache: 1];
ImageResolution =
Inquire(
camera,
"camera resolution") [instance: 5, cache: 1];
RECresolution =
Switch(
RECresNull,
RECresolution,
ImageResolution) [instance: 10, cache: 1];
RECaspectNull =
Inquire(
RECaspect,
"is null + 1") [instance: 6, cache: 1];
ImageAspect =
Inquire(
camera,
"camera aspect") [instance: 7, cache: 1];
RECaspect =
Switch(
RECaspectNull,
RECaspect,
ImageAspect) [instance: 11, cache: 1];
switchRECenable = Compute(
"$0 == 0 ? 1 : (($2 == $3) && ($4 == $5)) ? ($1 == 1 ? 2 : 3) : 4",
RECenable,
switchRenderMode,
RECresolution,
ImageResolution,
RECaspect,
ImageAspect) [instance: 4, cache: 1];
NoRECobject, RECNoRerenderObject, RECNoRerHW, RECRerenderObject =
Route(switchRECenable, annoted_object);
Display(
NoRECobject,
camera,
where,
throttle) [instance: 1, cache: 1];
image =
Render(
RECNoRerenderObject,
camera) [instance: 1, cache: 1];
Display(
image,
NULL,
where,
throttle) [instance: 2, cache: 1];
WriteImage(
image,
RECfile,
RECformat) [instance: 1, cache: 1];
rec_where = Display(
RECNoRerHW,
camera,
where,
throttle) [instance: 1, cache: 0];
rec_image = ReadImageWindow(
rec_where) [instance: 1, cache: 1];
WriteImage(
rec_image,
RECfile,
RECformat) [instance: 1, cache: 1];
RECupdateCamera =
UpdateCamera(
camera,
resolution=RECresolution,
aspect=RECaspect) [instance: 2, cache: 1];
Display(
RECRerenderObject,
camera,
where,
throttle) [instance: 1, cache: 1];
RECRerenderObject =
ScaleScreen(
RECRerenderObject,
NULL,
RECresolution,
camera) [instance: 1, cache: 1];
image =
Render(
RECRerenderObject,
RECupdateCamera) [instance: 2, cache: 1];
WriteImage(
image,
RECfile,
RECformat) [instance: 2, cache: 1];
}
main_Image_1_in_1 = "Image_1";
main_Image_1_in_3 = "X24,,";
main_Image_1_in_4 = 1;
main_Image_1_in_5 = [0 0 0];
main_Image_1_in_6 = [6.50456 4.48337 3.65602];
main_Image_1_in_7 = 4.66499;
main_Image_1_in_8 = 640;
main_Image_1_in_9 = 0.751;
main_Image_1_in_10 = [-0.159773 0.752661 -0.638728];
main_Image_1_in_11 = NULL;
main_Image_1_in_12 = 0;
main_Image_1_in_13 = NULL;
main_Image_1_in_14 = 1;
main_Image_1_in_15 = NULL;
main_Image_1_in_16 = NULL;
main_Image_1_in_17 = NULL;
main_Image_1_in_18 = NULL;
main_Image_1_in_19 = 0;
main_Image_1_in_20 = NULL;
main_Image_1_in_21 = NULL;
main_Image_1_in_22 = NULL;
main_Image_1_in_23 = NULL;
main_Image_1_in_25 = NULL;
main_Image_1_in_26 = NULL;
main_Image_1_in_27 = NULL;
main_Image_1_in_28 = NULL;
main_Image_1_in_29 = NULL;
main_Image_1_in_30 = NULL;
main_Image_1_in_31 = NULL;
main_Image_1_in_32 = NULL;
main_Image_1_in_33 = NULL;
main_Image_1_in_34 = NULL;
main_Image_1_in_35 = NULL;
main_Image_1_in_36 = NULL;
main_Image_1_in_37 = NULL;
main_Image_1_in_38 = NULL;
main_Image_1_in_39 = NULL;
main_Image_1_in_40 = NULL;
main_Image_1_in_41 = "rotate";
main_Image_1_in_42 = NULL;
main_Image_1_in_43 = NULL;
main_Image_1_in_44 = NULL;
main_Image_1_in_45 = NULL;
main_Image_1_in_46 = NULL;
main_Image_1_in_47 = NULL;
main_Image_1_in_48 = NULL;
main_Image_1_in_49 = NULL;
Executive("product version 4 1 3");
$sync
main();
