forked from Kitware/VTK
-
Notifications
You must be signed in to change notification settings - Fork 26
/
vtkImageRFFT.h
83 lines (56 loc) · 2.79 KB
/
vtkImageRFFT.h
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
/*=========================================================================
Program: Visualization Toolkit
Module: vtkImageRFFT.h
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.
=========================================================================*/
// .NAME vtkImageRFFT - Reverse Fast Fourier Transform.
// .SECTION Description
// vtkImageRFFT implements the reverse fast Fourier transform. The input
// can have real or complex data in any components and data types, but
// the output is always complex doubles with real values in component0, and
// imaginary values in component1. The filter is fastest for images that
// have power of two sizes. The filter uses a butterfly fitlers for each
// prime factor of the dimension. This makes images with prime number dimensions
// (i.e. 17x17) much slower to compute. Multi dimensional (i.e volumes)
// FFT's are decomposed so that each axis executes in series.
// In most cases the RFFT will produce an image whose imaginary values are all
// zero's. In this case vtkImageExtractComponents can be used to remove
// this imaginary components leaving only the real image.
// .SECTION See Also
// vtkImageExtractComponenents
#ifndef __vtkImageRFFT_h
#define __vtkImageRFFT_h
#include "vtkImageFourierFilter.h"
class VTK_IMAGING_EXPORT vtkImageRFFT : public vtkImageFourierFilter
{
public:
static vtkImageRFFT *New();
vtkTypeMacro(vtkImageRFFT,vtkImageFourierFilter);
// Description:
// For streaming and threads. Splits output update extent into num pieces.
// This method needs to be called num times. Results must not overlap for
// consistent starting extent. Subclass can override this method. This
// method returns the number of pieces resulting from a successful split.
// This can be from 1 to "total". If 1 is returned, the extent cannot be
// split.
int SplitExtent(int splitExt[6], int startExt[6],
int num, int total);
protected:
vtkImageRFFT() {};
~vtkImageRFFT() {};
virtual int IterativeRequestInformation(vtkInformation* in,
vtkInformation* out);
virtual int IterativeRequestUpdateExtent(vtkInformation* in,
vtkInformation* out);
void ThreadedExecute(vtkImageData *inData, vtkImageData *outData,
int outExt[6], int threadId);
private:
vtkImageRFFT(const vtkImageRFFT&); // Not implemented.
void operator=(const vtkImageRFFT&); // Not implemented.
};
#endif