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
|
/*=========================================================================
Program: Visualization Toolkit
Module: ImageWeightedSum.cxx
Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen
All rights reserved.
See Copyright.txt or http://www.kitware.com/Copyright.htm for details.
This software is distributed WITHOUT ANY WARRANTY; without even
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
PURPOSE. See the above copyright notice for more information.
=========================================================================*/
#include <vtkColorTransferFunction.h>
#include <vtkGlyph3D.h>
#include <vtkImageData.h>
#include <vtkImageReslice.h>
#include <vtkImageMapToColors.h>
#include <vtkImageAppendComponents.h>
#include <vtkInformation.h>
#include <vtkNew.h>
#include <vtkSmartPointer.h>
#include <vtkSphereSource.h>
#include <vtkStreamingDemandDrivenPipeline.h>
int TestUpdateExtentReset(int vtkNotUsed(argc), char * vtkNotUsed(argv) [] )
{
vtkSmartPointer<vtkImageData> img = vtkSmartPointer<vtkImageData>::New();
img->SetDimensions(100, 100, 100);
img->AllocateScalars(VTK_FLOAT, 1);
float *scalars = static_cast<float *>(img->GetScalarPointer());
vtkIdType n = 100*100*100;
for (vtkIdType i = 0; i < n; i++)
{
scalars[i] = 0.0;
}
vtkSmartPointer<vtkImageReslice> reslicer = vtkSmartPointer<vtkImageReslice>::New();
reslicer->SetInputData(img);
reslicer->SetOutputExtent(0, 100, 0, 100, 0, 0);
vtkSmartPointer<vtkImageMapToColors> colors =
vtkSmartPointer<vtkImageMapToColors>::New();
colors->SetInputConnection(reslicer->GetOutputPort());
vtkSmartPointer<vtkColorTransferFunction> ctf =
vtkSmartPointer<vtkColorTransferFunction>::New();
ctf->AddRGBPoint(0, 1., 0., 0.);
colors->SetLookupTable(ctf);
vtkSmartPointer<vtkImageAppendComponents> append =
vtkSmartPointer<vtkImageAppendComponents>::New();
append->SetInputConnection(0 ,colors->GetOutputPort());
colors->Update();
append->Update();
colors->Update();
// At this point, the COMBINED_UPDATE_EXTENT of the output of reslicer must be
// reset to {0, -1, 0, -1, 0, -1}, otherwise, the following will fail because
// when computing the output extent it will take into account the old (not reset)
// COMBINED_UPDATE_EXTENT value.
reslicer->SetOutputExtent(0, 100, 0, 80, 0, 0);
colors->Update();
vtkNew<vtkGlyph3D> polyDataFilter;
polyDataFilter->SetInputConnection(0, colors->GetOutputPort());
vtkNew<vtkSphereSource> sphere;
polyDataFilter->SetSourceConnection(sphere->GetOutputPort());
polyDataFilter->Update();
// After Update(), the COMBINED_UPDATE_EXTENT of the output of reslicer must be
// reset to {0, -1, 0, -1, 0, -1}, otherwise, the following will fail because
// when computing the output extent it will take into account the old (not reset)
// COMBINED_UPDATE_EXTENT value.
int combinedExtent[6];
reslicer->GetExecutive()->GetOutputInformation(0)->Get(
vtkStreamingDemandDrivenPipeline::COMBINED_UPDATE_EXTENT(), combinedExtent);
if (combinedExtent[0] <= combinedExtent[1] ||
combinedExtent[2] <= combinedExtent[3] ||
combinedExtent[4] <= combinedExtent[5])
{
return EXIT_FAILURE;
}
reslicer->SetOutputExtent(0, 100, 0, 50, 0, 0);
// For some reason there was no error reported when the combined extent was
// was wrong, however you could check in vtkImageReslice::ThreadRequestData
// and you'll see that the output extent is still {0, 100, 0, 80, 0, 0}
append->Update();
return EXIT_SUCCESS;
}
|
