65#ifndef vtkFixedPointVolumeRayCastMapper_h
66#define vtkFixedPointVolumeRayCastMapper_h
68#include "vtkRenderingVolumeModule.h"
69#include "vtkThreads.h"
72#define VTKKW_FP_SHIFT 15
73#define VTKKW_FPMM_SHIFT 17
74#define VTKKW_FP_MASK 0x7fff
75#define VTKKW_FP_SCALE 32767.0
241 void LookupColorUC(
unsigned short* colorTable,
unsigned short* scalarOpacityTable,
242 unsigned short index,
unsigned char color[4]);
244 unsigned short index[4],
int components,
unsigned char color[4]);
246 unsigned short* scalarOpacityTable[4],
unsigned short index[4],
float weights[4],
247 int components,
unsigned char color[4]);
275 int x,
int y,
unsigned int pos[3],
unsigned int dir[3],
unsigned int* numSteps);
299 double viewDirection[3],
double viewUp[3]);
370 void ComputeMatrices(
double inputOrigin[3],
double inputSpacing[3],
int inputExtent[6],
466 double rayStart[3],
double rayEnd[3],
float rayDirection[3],
double bounds[6]);
476 double rayStart[3],
double rayEnd[3],
int numClippingPlanes,
float* clippingPlanes);
542 return ((dir < 0.0) ? (
static_cast<unsigned int>(-dir *
VTKKW_FP_SCALE + 0.5))
543 : (0x80000000 +
static_cast<unsigned int>(dir *
VTKKW_FP_SCALE + 0.5)));
547 float in[3],
unsigned int out[3])
549 out[0] = ((in[0] < 0.0) ? (
static_cast<unsigned int>(-in[0] *
VTKKW_FP_SCALE + 0.5))
550 : (0x80000000 +
static_cast<unsigned int>(in[0] *
VTKKW_FP_SCALE + 0.5)));
551 out[1] = ((in[1] < 0.0) ? (
static_cast<unsigned int>(-in[1] *
VTKKW_FP_SCALE + 0.5))
552 : (0x80000000 +
static_cast<unsigned int>(in[1] *
VTKKW_FP_SCALE + 0.5)));
553 out[2] = ((in[2] < 0.0) ? (
static_cast<unsigned int>(-in[2] *
VTKKW_FP_SCALE + 0.5))
554 : (0x80000000 +
static_cast<unsigned int>(in[2] *
VTKKW_FP_SCALE + 0.5)));
558 unsigned int position[3],
unsigned int increment[3])
560 if (increment[0] & 0x80000000)
562 position[0] += (increment[0] & 0x7fffffff);
566 position[0] -= increment[0];
568 if (increment[1] & 0x80000000)
570 position[1] += (increment[1] & 0x7fffffff);
574 position[1] -= increment[1];
576 if (increment[2] & 0x80000000)
578 position[2] += (increment[2] & 0x7fffffff);
582 position[2] -= increment[2];
594 unsigned int v[3],
unsigned int* ptr)
602 unsigned int in[3],
unsigned int out[3])
613 mmpos[1] *
static_cast<vtkIdType>(this->MinMaxVolumeSize[0]) + mmpos[0]) +
616 return ((*(this->
MinMaxVolume + 3 * offset + 2)) & 0x00ff);
620 unsigned int mmpos[3],
int c,
unsigned short maxIdx,
int flip)
624 mmpos[1] *
static_cast<vtkIdType>(this->MinMaxVolumeSize[0]) + mmpos[0]) +
635 return (*(this->
MinMaxVolume + 3 * offset + 1) > maxIdx);
645 unsigned short* scalarOpacityTable,
unsigned short index,
unsigned char color[4])
647 unsigned short alpha = scalarOpacityTable[index];
648 color[0] =
static_cast<unsigned char>(
649 (colorTable[3 * index] * alpha + 0x7fff) >> (2 *
VTKKW_FP_SHIFT - 8));
650 color[1] =
static_cast<unsigned char>(
651 (colorTable[3 * index + 1] * alpha + 0x7fff) >> (2 *
VTKKW_FP_SHIFT - 8));
652 color[2] =
static_cast<unsigned char>(
653 (colorTable[3 * index + 2] * alpha + 0x7fff) >> (2 *
VTKKW_FP_SHIFT - 8));
654 color[3] =
static_cast<unsigned char>(alpha >> (
VTKKW_FP_SHIFT - 8));
658 unsigned short* scalarOpacityTable,
unsigned short index[4],
int components,
659 unsigned char color[4])
661 unsigned short alpha;
665 alpha = scalarOpacityTable[index[1]];
666 color[0] =
static_cast<unsigned char>(
667 (colorTable[3 * index[0]] * alpha + 0x7fff) >> (2 *
VTKKW_FP_SHIFT - 8));
668 color[1] =
static_cast<unsigned char>(
669 (colorTable[3 * index[0] + 1] * alpha + 0x7fff) >> (2 *
VTKKW_FP_SHIFT - 8));
670 color[2] =
static_cast<unsigned char>(
671 (colorTable[3 * index[0] + 2] * alpha + 0x7fff) >> (2 *
VTKKW_FP_SHIFT - 8));
672 color[3] =
static_cast<unsigned char>(alpha >> (
VTKKW_FP_SHIFT - 8));
675 alpha = scalarOpacityTable[index[3]];
676 color[0] =
static_cast<unsigned char>((index[0] * alpha + 0x7fff) >>
VTKKW_FP_SHIFT);
677 color[1] =
static_cast<unsigned char>((index[1] * alpha + 0x7fff) >>
VTKKW_FP_SHIFT);
678 color[2] =
static_cast<unsigned char>((index[2] * alpha + 0x7fff) >>
VTKKW_FP_SHIFT);
679 color[3] =
static_cast<unsigned char>(alpha >> (
VTKKW_FP_SHIFT - 8));
685 unsigned short* colorTable[4],
unsigned short* scalarOpacityTable[4],
unsigned short index[4],
686 float weights[4],
int components,
unsigned char color[4])
688 unsigned int tmp[4] = { 0, 0, 0, 0 };
690 for (
int i = 0; i < components; i++)
692 unsigned short alpha =
693 static_cast<unsigned short>(
static_cast<float>(scalarOpacityTable[i][index[i]]) * weights[i]);
694 tmp[0] +=
static_cast<unsigned char>(
695 ((colorTable[i][3 * index[i]]) * alpha + 0x7fff) >> (2 *
VTKKW_FP_SHIFT - 8));
696 tmp[1] +=
static_cast<unsigned char>(
697 ((colorTable[i][3 * index[i] + 1]) * alpha + 0x7fff) >> (2 *
VTKKW_FP_SHIFT - 8));
698 tmp[2] +=
static_cast<unsigned char>(
699 ((colorTable[i][3 * index[i] + 2]) * alpha + 0x7fff) >> (2 *
VTKKW_FP_SHIFT - 8));
700 tmp[3] +=
static_cast<unsigned char>(alpha >> (
VTKKW_FP_SHIFT - 8));
703 color[0] =
static_cast<unsigned char>((tmp[0] > 255) ? (255) : (tmp[0]));
704 color[1] =
static_cast<unsigned char>((tmp[1] > 255) ? (255) : (tmp[1]));
705 color[2] =
static_cast<unsigned char>((tmp[2] > 255) ? (255) : (tmp[2]));
706 color[3] =
static_cast<unsigned char>((tmp[3] > 255) ? (255) : (tmp[3]));
Defines a transfer function for mapping a property to an RGB color value.
encode a direction into a one or two byte value
Compute shading tables for encoded normals.
Use finite differences to estimate gradient.
helper class for a ray cast image
A helper that generates composite images for the volume ray cast mapper.
A helper that generates composite images for the volume ray cast mapper.
A helper that generates composite images for the volume ray cast mapper.
A helper that generates composite images for the volume ray cast mapper.
A helper that generates MIP images for the volume ray cast mapper.
void SetRayCastImage(vtkFixedPointRayCastImage *)
Set / Get the underlying image object.
vtkRayCastImageDisplayHelper * ImageDisplayHelper
vtkMatrix4x4 * ViewToWorldMatrix
vtkPiecewiseFunction * SavedGradientOpacityFunction[4]
vtkMatrix4x4 * ViewToVoxelsMatrix
vtkEncodedGradientShader * GradientShader
vtkFixedPointVolumeRayCastMIPHelper * MIPHelper
vtkFixedPointVolumeRayCastCompositeGOHelper * CompositeGOHelper
float OldImageSampleDistance
int ShouldUseNearestNeighborInterpolation(vtkVolume *vol)
void ReleaseGraphicsResources(vtkWindow *) override
WARNING: INTERNAL METHOD - NOT INTENDED FOR GENERAL USE Release any graphics resources that are being...
vtkPiecewiseFunction * SavedScalarOpacityFunction[4]
vtkMultiThreader * Threader
unsigned short ScalarOpacityTable[4][32768]
float RetrieveRenderTime(vtkRenderer *ren)
float MaximumImageSampleDistance
unsigned short * GetColorTable(int c)
unsigned short * GetSpecularShadingTable(int c)
float ImageSampleDistance
int SavedColorChannels[4]
vtkPiecewiseFunction * SavedGrayFunction[4]
vtkTransform * VoxelsTransform
unsigned short * ContiguousGradientNormal
float MinimumImageSampleDistance
vtkFixedPointVolumeRayCastMapper()
vtkImageData * SavedMinMaxInput
unsigned short DiffuseShadingTable[4][65536 *3]
void UpdateCroppingRegions()
void GetFloatTripleFromPointer(float v[3], float *ptr)
unsigned short * MinMaxVolume
unsigned short GradientOpacityTable[4][256]
unsigned char ** GradientMagnitude
vtkColorTransferFunction * SavedRGBFunction[4]
void CaptureZBuffer(vtkRenderer *ren)
void ApplyFinalColorWindowLevel()
float ComputeRequiredImageSampleDistance(float desiredTime, vtkRenderer *ren)
What is the image sample distance required to achieve the desired time?
void ComputeRayInfo(int x, int y, unsigned int pos[3], unsigned int dir[3], unsigned int *numSteps)
vtkFixedPointVolumeRayCastCompositeShadeHelper * CompositeShadeHelper
unsigned short * GetScalarOpacityTable(int c)
vtkMatrix4x4 * VoxelsToWorldMatrix
unsigned short ColorTable[4][32768 *3]
vtkTimeStamp SavedGradientsMTime
unsigned short * GetGradientOpacityTable(int c)
vtkRenderWindow * RenderWindow
float SavedScalarOpacityDistance[4]
float MinimumViewDistance
void StoreRenderTime(vtkRenderer *ren, vtkVolume *vol, float t)
vtkFiniteDifferenceGradientEstimator * GradientEstimator
vtkFixedPointVolumeRayCastCompositeGOShadeHelper * CompositeGOShadeHelper
vtkTypeBool IntermixIntersectingGeometry
int ComputeRowBounds(vtkRenderer *ren, int imageFlag, int rowBoundsFlag, int inputExtent[6])
vtkTypeBool AutoAdjustSampleDistances
vtkImageData * MinMaxVolumeCache
float GetZBufferValue(int x, int y)
void DisplayRenderedImage(vtkRenderer *, vtkVolume *)
void FillInMaxGradientMagnitudes(int fullDim[3], int smallDim[3])
unsigned short SpecularShadingTable[4][65536 *3]
unsigned int ToFixedPointDirection(float dir)
float WorldToVoxelsArray[16]
vtkMatrix4x4 * VolumeMatrix
void InitializeRayInfo(vtkVolume *vol)
void SetNumberOfThreads(int num)
Set/Get the number of threads to use.
vtkImageData * SavedParametersInput
void Render(vtkRenderer *, vtkVolume *) override
WARNING: INTERNAL METHOD - NOT INTENDED FOR GENERAL USE Initialize rendering for this volume.
void ComputeGradients(vtkVolume *vol)
int UpdateShadingTable(vtkRenderer *ren, vtkVolume *vol)
void GetUIntTripleFromPointer(unsigned int v[3], unsigned int *ptr)
float ComputeRequiredImageSampleDistance(float desiredTime, vtkRenderer *ren, vtkVolume *vol)
What is the image sample distance required to achieve the desired time?
~vtkFixedPointVolumeRayCastMapper() override
float SavedSampleDistance
static vtkFixedPointVolumeRayCastMapper * New()
float GetEstimatedRenderTime(vtkRenderer *ren, vtkVolume *vol)
Get an estimate of the rendering time for a given volume / renderer.
float RetrieveRenderTime(vtkRenderer *ren, vtkVolume *vol)
int GradientOpacityRequired
int CheckMinMaxVolumeFlag(unsigned int pos[3], int c)
vtkFixedPointVolumeRayCastCompositeHelper * CompositeHelper
vtkTransform * PerspectiveTransform
unsigned int ToFixedPointPosition(float val)
int CheckMIPMinMaxVolumeFlag(unsigned int pos[3], int c, unsigned short maxIdx, int flip)
float GetEstimatedRenderTime(vtkRenderer *ren)
unsigned short * GetDiffuseShadingTable(int c)
vtkDirectionEncoder * DirectionEncoder
vtkMatrix4x4 * PerspectiveMatrix
unsigned char ** GetGradientMagnitude()
float InteractiveSampleDistance
float VoxelsToWorldArray[16]
friend VTK_THREAD_RETURN_TYPE vtkFPVRCMSwitchOnDataType(void *arg)
void PerSubVolumeInitialization(vtkRenderer *, vtkVolume *, int)
vtkImageData * SavedGradientsInput
vtkTimeStamp SavedParametersMTime
unsigned short ** GradientNormal
void ComputeMatrices(double inputOrigin[3], double inputSpacing[3], int inputExtent[6], vtkRenderer *ren, vtkVolume *vol)
float ViewToVoxelsArray[16]
void PrintSelf(ostream &os, vtkIndent indent) override
Methods invoked by print to print information about the object including superclasses.
void LookupColorUC(unsigned short *colorTable, unsigned short *scalarOpacityTable, unsigned short index, unsigned char color[4])
vtkTransform * VoxelsToViewTransform
void FixedPointIncrement(unsigned int position[3], unsigned int increment[3])
float GradientMagnitudeScale[4]
void LookupAndCombineIndependentColorsUC(unsigned short *colorTable[4], unsigned short *scalarOpacityTable[4], unsigned short index[4], float weights[4], int components, unsigned char color[4])
vtkMatrix4x4 * WorldToVoxelsMatrix
vtkVolumeRayCastSpaceLeapingImageFilter * SpaceLeapFilter
void PerVolumeInitialization(vtkRenderer *, vtkVolume *)
int ClipRayAgainstVolume(double rayStart[3], double rayEnd[3], float rayDirection[3], double bounds[6])
void CreateCanonicalView(vtkVolume *volume, vtkImageData *image, int blend_mode, double viewDirection[3], double viewUp[3])
int PerImageInitialization(vtkRenderer *, vtkVolume *, int, double *, double *, int *)
unsigned int FixedPointCroppingRegionPlanes[6]
vtkFixedPointRayCastImage * RayCastImage
float * TransformedClippingPlanes
vtkVolume ** RenderVolumeTable
void UpdateMinMaxVolume(vtkVolume *vol)
friend VTK_THREAD_RETURN_TYPE FixedPointVolumeRayCastMapper_CastRays(void *arg)
vtkDataArray * PreviousScalars
vtkRenderer ** RenderRendererTable
void ShiftVectorDown(unsigned int in[3], unsigned int out[3])
int UpdateColorTable(vtkVolume *vol)
int NumTransformedClippingPlanes
unsigned char * ContiguousGradientMagnitude
int ClipRayAgainstClippingPlanes(double rayStart[3], double rayEnd[3], int numClippingPlanes, float *clippingPlanes)
int GetNumberOfThreads()
Set/Get the number of threads to use.
int UpdateGradients(vtkVolume *vol)
vtkDataArray * CurrentScalars
vtkTypeBool LockSampleDistanceToInputSpacing
unsigned int CroppingRegionMask[27]
unsigned short ** GetGradientNormal()
void LookupDependentColorUC(unsigned short *colorTable, unsigned short *scalarOpacityTable, unsigned short index[4], int components, unsigned char color[4])
float GradientMagnitudeShift[4]
int CheckIfCropped(unsigned int pos[3])
int NumberOfGradientSlices
vtkMatrix4x4 * VoxelsToViewMatrix
topologically and geometrically regular array of data
a simple class to control print indentation
represent and manipulate 4x4 transformation matrices
A class for performing multithreaded execution.
Defines a 1D piecewise function.
maintain a list of planes
helper class that draws the image to the screen
create a window for renderers to draw into
abstract specification for renderers
record modification and/or execution time
Timer support and logging.
int CroppingRegionFlags
Cropping variables, and a method for converting the world coordinate cropping region planes to voxel ...
Builds the space leaping data structure.
represents a volume (data & properties) in a rendered scene
window superclass for vtkRenderWindow
VTK_THREAD_RETURN_TYPE vtkFPVRCMSwitchOnDataType(void *arg)
VTK_THREAD_RETURN_TYPE FixedPointVolumeRayCastMapper_CastRays(void *arg)