/usr/include/mia-2.2/mia/2d/fftkernel.hh is in libmia-2.2-dev 2.2.7-3.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
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 | /* -*- mia-c++ -*-
*
* This file is part of MIA - a toolbox for medical image analysis
* Copyright (c) Leipzig, Madrid 1999-2015 Gert Wollny
*
* MIA is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation; either version 3 of the License, or
* (at your option) any later version.
*
* This program is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with MIA; if not, see <http://www.gnu.org/licenses/>.
*
*/
#ifndef mia_2d_fftkernel_hh
#define mia_2d_fftkernel_hh
#include <complex>
#include <fftw3.h>
#include <mia/core/defines.hh>
#include <mia/core/factory.hh>
#include <mia/core/spacial_kernel.hh>
#include <mia/2d/defines2d.hh>
#include <mia/2d/vector.hh>
NS_MIA_BEGIN
struct EXPORT_2D fft2d_kernel_data {
static const char *data_descr;
};
/**
\cond NEEDS_REHAUL
@ingroup filtering
\brief Base class for 2D FFT filters.
This class provides tha basic interface for filters that work within the
frequency domain by means of a FFT.
Filter applications are applied like this:
first a real-to-complex transformation is run, then run the filter on
the half-complex transform and run the back-transform.
*/
class EXPORT_2D CFFT2DKernel :public CProductBase {
public:
/// plugin search path helper type
typedef fft2d_kernel_data plugin_data;
/// plugin search path helper type
typedef kernel_plugin_type plugin_type;
CFFT2DKernel();
virtual ~CFFT2DKernel();
/**
Run the FFT by first calling the forward plan, then
call do_apply() and then running the backward plan.
Normalization is left to the caller,
*/
void apply() const;
/**
Prepare the FFT structures and return the buffer where the data needs to be
put when it should be processed.
*/
float *prepare(const C2DBounds& size);
private:
/*
free all the FFT structures
*/
void tear_down();
virtual void do_apply(const C2DBounds& m_size, size_t m_realsize_x,
fftwf_complex *m_cbuffer) const = 0;
C2DBounds m_size;
fftwf_complex *m_cbuffer;
float *m_fbuffer;
float m_scale;
fftwf_plan m_forward_plan;
fftwf_plan m_backward_plan;
size_t m_realsize_x;
};
typedef std::shared_ptr<CFFT2DKernel > PFFT2DKernel;
typedef TFactory<CFFT2DKernel> CFFT2DKernelPlugin;
typedef THandlerSingleton<TFactoryPluginHandler<CFFT2DKernelPlugin> > CFFT2DKernelPluginHandler;
/// \endcond
NS_MIA_END
#endif
|