1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
|
/*
* vp_renderC.c
*
* Function to render classified volumes.
*
* Copyright (c) 1994 The Board of Trustees of The Leland Stanford
* Junior University. All rights reserved.
*
* Permission to use, copy, modify and distribute this software and its
* documentation for any purpose is hereby granted without fee, provided
* that the above copyright notice and this permission notice appear in
* all copies of this software and that you do not sell the software.
* Commercial licensing is available by contacting the author.
*
* THE SOFTWARE IS PROVIDED "AS IS" AND WITHOUT WARRANTY OF ANY KIND,
* EXPRESS, IMPLIED OR OTHERWISE, INCLUDING WITHOUT LIMITATION, ANY
* WARRANTY OF MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE.
*
* Author:
* Phil Lacroute
* Computer Systems Laboratory
* Electrical Engineering Dept.
* Stanford University
*/
/*
* $Date: 1994/12/30 23:52:38 $
* $Revision: 1.27 $
*/
#include "vp_global.h"
#define COMP_AC1PB_FUNC VPCompAC1PB
extern void VPCompAC1PB();
#define COMP_AC3PB_FUNC VPCompAC3PB
extern void VPCompAC3PB();
#ifdef COMP_AC11B
#define COMP_AC11B_FUNC VPCompAC11B
extern void VPCompAC11B();
#else
#define COMP_AC11B_FUNC VPCompAC1NB
#endif
#ifdef COMP_AC31B
#define COMP_AC31B_FUNC VPCompAC31B
extern void VPCompAC31B();
#else
#define COMP_AC31B_FUNC VPCompAC3NB
#endif
#ifdef COMP_AC12B
#define COMP_AC12B_FUNC VPCompAC12B
extern void VPCompAC12B();
#else
#define COMP_AC12B_FUNC VPCompAC1NB
#endif
#ifdef COMP_AC32B
#define COMP_AC32B_FUNC VPCompAC32B
extern void VPCompAC32B();
#else
#define COMP_AC32B_FUNC VPCompAC3NB
#endif
#define COMP_AC1NB_FUNC VPCompAC1NB
extern void VPCompAC1NB();
#define COMP_AC3NB_FUNC VPCompAC3NB
extern void VPCompAC3NB();
#define COMP_AC1PS_FUNC VPCompAC1PB
#define COMP_AC3PS_FUNC VPCompAC3PB
#ifdef COMP_AC11S
#define COMP_AC11S_FUNC VPCompAC11S
extern void VPCompAC11S();
#else
#define COMP_AC11S_FUNC VPCompAC1NS
#endif
#ifdef COMP_AC31S
#define COMP_AC31S_FUNC VPCompAC31S
extern void VPCompAC31S();
#else
#define COMP_AC31S_FUNC VPCompAC3NS
#endif
#ifdef COMP_AC12S
#define COMP_AC12S_FUNC VPCompAC12S
extern void VPCompAC12S();
#else
#define COMP_AC12S_FUNC VPCompAC1NS
#endif
#ifdef COMP_AC32S
#define COMP_AC32S_FUNC VPCompAC32S
extern void VPCompAC32S();
#else
#define COMP_AC32S_FUNC VPCompAC3NS
#endif
#define COMP_AC1NS_FUNC VPCompAC1NS
extern void VPCompAC1NS();
#define COMP_AC3NS_FUNC VPCompAC3NS
extern void VPCompAC3NS();
#ifdef INDEX_VOLUME
extern void VPCompAI11B();
#endif
#define SHADOWS_OFF 0
#define SHADOWS_ON 1
#define SHADOW_OPTS 2
#define MATERIAL_CALLBACK 0
#define MATERIAL_ONE 1
#define MATERIAL_TWO 2
#define MATERIAL_MORE 3
#define MATERIAL_OPTS 4
#define COLOR_GRAY 0
#define COLOR_RGB 1
#define COLOR_OPTS 2
static void (*AffineProcTable[SHADOW_OPTS][MATERIAL_OPTS][COLOR_OPTS])() = {
{
{ COMP_AC1PB_FUNC, COMP_AC3PB_FUNC },
{ COMP_AC11B_FUNC, COMP_AC31B_FUNC },
{ COMP_AC12B_FUNC, COMP_AC32B_FUNC },
{ COMP_AC1NB_FUNC, COMP_AC3NB_FUNC }
},
{
{ COMP_AC1PS_FUNC, COMP_AC3PS_FUNC },
{ COMP_AC11S_FUNC, COMP_AC31S_FUNC },
{ COMP_AC12S_FUNC, COMP_AC32S_FUNC },
{ COMP_AC1NS_FUNC, COMP_AC3NS_FUNC }
}
};
/*
* vpRenderClassifiedVolume
*
* Render a classified, run-length encoded volume using the shear-warp
* algorithm.
*/
vpResult
vpRenderClassifiedVolume(vpc)
vpContext *vpc;
{
int retcode;
void (*composite_func)();
int shadow_option, material_option, color_option;
#ifdef INDEX_VOLUME
int have_index;
#endif
/* check for errors and initialize */
if ((retcode = VPCheckShader(vpc)) != VP_OK)
return(retcode);
if ((retcode = VPFactorView(vpc)) != VP_OK)
return(retcode);
if ((retcode = VPCheckClassifiedVolume(vpc, vpc->best_view_axis)) != VP_OK)
return(retcode);
if ((retcode = VPCheckImage(vpc)) != VP_OK)
return(retcode);
if ((retcode = VPCheckShadows(vpc)) != VP_OK)
return(retcode);
#ifdef INDEX_VOLUME
switch (vpc->best_view_axis) {
case VP_X_AXIS:
if (vpc->rle_x->voxel_index != NULL)
have_index = 1;
else
have_index = 0;
break;
case VP_Y_AXIS:
if (vpc->rle_y->voxel_index != NULL)
have_index = 1;
else
have_index = 0;
break;
case VP_Z_AXIS:
if (vpc->rle_z->voxel_index != NULL)
have_index = 1;
else
have_index = 0;
break;
default:
VPBug("invalid viewing axis in vpRenderClassifiedVolume");
}
#endif /* INDEX_VOLUME */
Debug((vpc, VPDEBUG_RENDER, "Algorithm: affine RLEvolume\n"));
/* determine which options are enabled */
if (vpc->enable_shadows)
shadow_option = SHADOWS_ON;
else
shadow_option = SHADOWS_OFF;
if (vpc->shading_mode == CALLBACK_SHADER)
material_option = MATERIAL_CALLBACK;
else if (vpc->num_materials == 1)
material_option = MATERIAL_ONE;
else if (vpc->num_materials == 2)
material_option = MATERIAL_TWO;
else
material_option = MATERIAL_MORE;
if (vpc->color_channels == 1)
color_option = COLOR_GRAY;
else
color_option = COLOR_RGB;
/* render */
if (vpc->affine_view) {
/* choose a compositing function */
composite_func = AffineProcTable[shadow_option][material_option]
[color_option];
#ifdef INDEX_VOLUME
if (have_index && shadow_option == SHADOWS_OFF &&
material_option == MATERIAL_ONE && color_option == COLOR_GRAY) {
composite_func = VPCompAI11B;
}
#endif
VPRenderAffine(vpc, USE_RLEVOLUME, composite_func);
} else {
/* XXX perspective rendering not available yet */
return(VPSetError(vpc, VPERROR_BAD_OPTION));
}
return(VP_OK);
}
|
