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
|
/*=========================================================================
*
* Copyright Insight Software Consortium
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0.txt
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*
*=========================================================================*/
#include <iostream>
#include "itkImageFilterToVideoFilterWrapper.h"
#include "itkImage.h"
#include "itkRecursiveGaussianImageFilter.h"
#include "itkVideoFileReader.h"
#include "itkVideoFileWriter.h"
#include "itkImageFileReader.h"
#include "itkFileListVideoIO.h"
#include "itkFileListVideoIOFactory.h"
#include "itkTestingComparisonImageFilter.h"
/**
* Main test
*/
int itkImageFilterToVideoFilterWrapperTest( int argc, char* argv[] )
{
// Check parameters
if (argc < 3)
{
std::cerr << "Usage: " << argv[0] << " input_video output_video" << std::endl;
return EXIT_FAILURE;
}
// Get the lists of input and output files
std::vector<std::string> inputFiles = itk::FileListVideoIO::SplitFileNames(argv[1]);
std::vector<std::string> outputFiles = itk::FileListVideoIO::SplitFileNames(argv[2]);
if (inputFiles.size() != outputFiles.size())
{
std::cerr << "Must specify the same number of input and output frames" << std::endl;
return EXIT_FAILURE;
}
// Typedefs
typedef unsigned char PixelType;
typedef itk::Image<PixelType, 2> FrameType;
typedef itk::VideoStream< FrameType > VideoType;
typedef itk::RecursiveGaussianImageFilter< FrameType, FrameType > GaussianImageFilterType;
typedef itk::ImageFilterToVideoFilterWrapper< GaussianImageFilterType > GaussianVideoFilterType;
typedef itk::VideoFileReader< VideoType > VideoReaderType;
typedef itk::VideoFileWriter< VideoType > VideoWriterType;
// Register FileListIO with the factory -- shouldn't have to do this. Needs fixing
itk::ObjectFactoryBase::RegisterFactory( itk::FileListVideoIOFactory::New() );
// Set up reader and writer
VideoReaderType::Pointer reader = VideoReaderType::New();
VideoWriterType::Pointer writer = VideoWriterType::New();
reader->SetFileName(argv[1]);
writer->SetFileName(argv[2]);
// Instantiate a new video filter and an image filter
GaussianImageFilterType::Pointer imgGauss = GaussianImageFilterType::New();
GaussianVideoFilterType::Pointer vidGauss = GaussianVideoFilterType::New();
// Set the parameters on the image filter and plug it into the video filter
imgGauss->SetSigma(3);
vidGauss->SetImageFilter(imgGauss);
// String the pipeline together
vidGauss->SetInput(reader->GetOutput());
writer->SetInput(vidGauss->GetOutput());
// Run the pipeline
writer->Update();
//
// Check output
//
typedef itk::ImageFileReader< FrameType > ImageReaderType;
typedef itk::Testing::ComparisonImageFilter< FrameType, FrameType > DifferenceFilterType;
ImageReaderType::Pointer imReader1 = ImageReaderType::New();
ImageReaderType::Pointer imReader2 = ImageReaderType::New();
DifferenceFilterType::Pointer differ = DifferenceFilterType::New();
imgGauss->SetInput(imReader1->GetOutput());
differ->SetValidInput(imgGauss->GetOutput());
differ->SetTestInput(imReader2->GetOutput());
for (unsigned int i = 0; i < inputFiles.size(); ++i)
{
imReader1->SetFileName(inputFiles[i]);
imReader2->SetFileName(outputFiles[i]);
differ->Update();
if (itk::Math::NotAlmostEquals( differ->GetTotalDifference(), 0) )
{
std::cerr << "Frame " << i << " didn't produce the correct output. Difference = "
<< differ->GetTotalDifference() << std::endl;
return EXIT_FAILURE;
}
}
//////
// Return successfully
//////
return EXIT_SUCCESS;
}
|
