1 /*=========================================================================
2
3 Program: Visualization Toolkit
4 Module: TestGPURayCastVolumePolyData.cxx
5
6 Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen
7 All rights reserved.
8 See Copyright.txt or http://www.kitware.com/Copyright.htm for details.
9
10 This software is distributed WITHOUT ANY WARRANTY; without even
11 the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
12 PURPOSE. See the above copyright notice for more information.
13
14 =========================================================================*/
15 // This test covers additive method.
16 // This test volume renders a synthetic dataset with unsigned char values,
17 // with the additive method.
18
19 #include <vtkCamera.h>
20 #include <vtkColorTransferFunction.h>
21 #include <vtkDataArray.h>
22 #include <vtkGPUVolumeRayCastMapper.h>
23 #include <vtkImageData.h>
24 #include <vtkImageReader.h>
25 #include <vtkImageShiftScale.h>
26 #include <vtkNew.h>
27 #include <vtkOutlineFilter.h>
28 #include <vtkPiecewiseFunction.h>
29 #include <vtkPointData.h>
30 #include <vtkPolyDataMapper.h>
31 #include <vtkRegressionTestImage.h>
32 #include <vtkRenderWindow.h>
33 #include <vtkRenderWindowInteractor.h>
34 #include <vtkRenderer.h>
35 #include <vtkSmartPointer.h>
36 #include <vtkSphereSource.h>
37 #include <vtkTestUtilities.h>
38 #include <vtkTimerLog.h>
39 #include <vtkVolumeProperty.h>
40 #include <vtkXMLImageDataReader.h>
41
TestGPURayCastVolumePolyData(int argc,char * argv[])42 int TestGPURayCastVolumePolyData(int argc, char *argv[])
43 {
44 double scalarRange[2];
45
46 vtkNew<vtkActor> outlineActor;
47 vtkNew<vtkPolyDataMapper> outlineMapper;
48 vtkNew<vtkGPUVolumeRayCastMapper> volumeMapper;
49
50 vtkNew<vtkXMLImageDataReader> reader;
51 const char* volumeFile = vtkTestUtilities::ExpandDataFileName(
52 argc, argv, "Data/vase_1comp.vti");
53 reader->SetFileName(volumeFile);
54 volumeMapper->SetInputConnection(reader->GetOutputPort());
55
56 // Add outline filter
57 vtkNew<vtkOutlineFilter> outlineFilter;
58 outlineFilter->SetInputConnection(reader->GetOutputPort());
59 outlineMapper->SetInputConnection(outlineFilter->GetOutputPort());
60 outlineActor->SetMapper(outlineMapper.GetPointer());
61
62 volumeMapper->GetInput()->GetScalarRange(scalarRange);
63 volumeMapper->SetSampleDistance(0.1);
64 volumeMapper->SetAutoAdjustSampleDistances(0);
65 volumeMapper->SetBlendModeToComposite();
66
67 vtkNew<vtkRenderWindow> renWin;
68 renWin->SetMultiSamples(0);
69 vtkNew<vtkRenderer> ren;
70 renWin->AddRenderer(ren.GetPointer());
71 renWin->SetSize(400, 400);
72 ren->SetBackground(0.2, 0.2, 0.5);
73
74 vtkNew<vtkRenderWindowInteractor> iren;
75 iren->SetRenderWindow(renWin.GetPointer());
76
77 vtkNew<vtkPiecewiseFunction> scalarOpacity;
78 scalarOpacity->AddPoint(50, 0.0);
79 scalarOpacity->AddPoint(75, 1.0);
80
81 vtkNew<vtkVolumeProperty> volumeProperty;
82 volumeProperty->ShadeOn();
83 volumeProperty->SetInterpolationType(VTK_LINEAR_INTERPOLATION);
84 volumeProperty->SetScalarOpacity(scalarOpacity.GetPointer());
85
86 vtkSmartPointer<vtkColorTransferFunction> colorTransferFunction =
87 volumeProperty->GetRGBTransferFunction(0);
88 colorTransferFunction->RemoveAllPoints();
89 colorTransferFunction->AddRGBPoint(scalarRange[0], 0.6, 0.4, 0.1);
90
91 vtkSmartPointer<vtkVolume> volume = vtkSmartPointer<vtkVolume>::New();
92 volume->SetMapper(volumeMapper.GetPointer());
93 volume->SetProperty(volumeProperty.GetPointer());
94
95 /// Add sphere in the center of volume
96 int dims[3];
97 double spacing[3], center[3], origin[3];
98 reader->Update();
99 vtkSmartPointer<vtkImageData> im = reader->GetOutput();
100 im->GetDimensions(dims);
101 im->GetOrigin(origin);
102 im->GetSpacing(spacing);
103
104 center[0] = origin[0] + spacing[0]*dims[0]/2.0;
105 center[1] = origin[1] + spacing[1]*dims[1]/2.0;
106 center[2] = origin[2] + spacing[2]*dims[2]/2.0;
107
108 vtkNew<vtkSphereSource> sphereSource;
109 sphereSource->SetCenter(center);
110 sphereSource->SetRadius(dims[1]/3.0);
111 vtkNew<vtkPolyDataMapper> sphereMapper;
112 vtkNew<vtkActor> sphereActor;
113 sphereMapper->SetInputConnection(sphereSource->GetOutputPort());
114 sphereActor->SetMapper(sphereMapper.GetPointer());
115
116 ren->AddViewProp(volume.GetPointer());
117 ren->AddActor(outlineActor.GetPointer());
118 ren->AddActor(sphereActor.GetPointer());
119
120 renWin->Render();
121 ren->ResetCamera();
122
123 iren->Initialize();
124
125 int retVal = vtkRegressionTestImage( renWin.GetPointer() );
126 if( retVal == vtkRegressionTester::DO_INTERACTOR)
127 {
128 iren->Start();
129 }
130
131 return !retVal;
132 }
133