X-Git-Url: http://git.kpe.io/?a=blobdiff_plain;f=include%2Ffilter.h;h=5527b3eaddabf57cd094809de8db781cf042d3e1;hb=04c6595c9b3a36151671fd8f697f8a9a38dce735;hp=9ee7390a4acd83c96bd9a72b7f9990fbfe7bede9;hpb=dfa390de2efc04d85b03718a6480f735516df0e8;p=ctsim.git diff --git a/include/filter.h b/include/filter.h index 9ee7390..5527b3e 100644 --- a/include/filter.h +++ b/include/filter.h @@ -9,7 +9,7 @@ ** This is part of the CTSim program ** Copyright (C) 1983-2000 Kevin Rosenberg ** -** $Id: filter.h,v 1.12 2000/07/07 15:30:59 kevin Exp $ +** $Id: filter.h,v 1.15 2000/07/20 11:17:31 kevin Exp $ ** ** This program is free software; you can redistribute it and/or modify ** it under the terms of the GNU General Public License (version 2) as @@ -29,6 +29,17 @@ #define FILTER_H +#ifdef HAVE_CONFIG_H +#include "config.h" +#endif +#ifdef HAVE_FFTW +#include +#include +#endif + +#include + + class SignalFilter { public: @@ -64,35 +75,61 @@ class SignalFilter { DOMAIN_SPATIAL } DomainID; - static const char FILTER_ABS_BANDLIMIT_STR[]= "abs_bandlimit"; - static const char FILTER_ABS_SINC_STR[]= "abs_sinc"; - static const char FILTER_ABS_COS_STR[]= "abs_cos"; - static const char FILTER_ABS_HAMMING_STR[]= "abs_hamming"; - static const char FILTER_SHEPP_STR[]= "shepp"; - static const char FILTER_BANDLIMIT_STR[]= "bandlimit"; - static const char FILTER_SINC_STR[]= "sinc"; - static const char FILTER_COS_STR[]= "cos"; - static const char FILTER_HAMMING_STR[]= "hamming"; - static const char FILTER_TRIANGLE_STR[]= "triangle"; + // Filters + static const char FILTER_ABS_BANDLIMIT_STR[]; + static const char FILTER_ABS_SINC_STR[]; + static const char FILTER_ABS_COS_STR[]; + static const char FILTER_ABS_HAMMING_STR[]; + static const char FILTER_SHEPP_STR[]; + static const char FILTER_BANDLIMIT_STR[]; + static const char FILTER_SINC_STR[]; + static const char FILTER_COS_STR[]; + static const char FILTER_HAMMING_STR[]; + static const char FILTER_TRIANGLE_STR[]; + + static const char FILTER_ABS_BANDLIMIT_TITLE_STR[]; + static const char FILTER_ABS_SINC_TITLE_STR[]; + static const char FILTER_ABS_COS_TITLE_STR[]; + static const char FILTER_ABS_HAMMING_TITLE_STR[]; + static const char FILTER_SHEPP_TITLE_STR[]; + static const char FILTER_BANDLIMIT_TITLE_STR[]; + static const char FILTER_SINC_TITLE_STR[]; + static const char FILTER_COS_TITLE_STR[]; + static const char FILTER_HAMMING_TITLE_STR[]; + static const char FILTER_TRIANGLE_TITLE_STR[]; - static const char FILTER_METHOD_CONVOLUTION_STR[]= "convolution"; - static const char FILTER_METHOD_FOURIER_STR[]= "fourier"; - static const char FILTER_METHOD_FOURIER_TABLE_STR[]="fourier_table"; - static const char FILTER_METHOD_FFT_STR[]= "fft"; + // Filter Methods + static const char FILTER_METHOD_CONVOLUTION_STR[]; + static const char FILTER_METHOD_FOURIER_STR[]; + static const char FILTER_METHOD_FOURIER_TABLE_STR[]; + static const char FILTER_METHOD_FFT_STR[]; #if HAVE_FFTW - static const char FILTER_METHOD_FFTW_STR[]= "fftw"; - static const char FILTER_METHOD_RFFTW_STR[]= "rfftw"; + static const char FILTER_METHOD_FFTW_STR[]; + static const char FILTER_METHOD_RFFTW_STR[]; #endif - static const char DOMAIN_FREQUENCY_STR[]="frequency"; - static const char DOMAIN_SPATIAL_STR[]="spatial"; + static const char FILTER_METHOD_CONVOLUTION_TITLE_STR[]; + static const char FILTER_METHOD_FOURIER_TITLE_STR[]; + static const char FILTER_METHOD_FOURIER_TABLE_TITLE_STR[]; + static const char FILTER_METHOD_FFT_TITLE_STR[]; +#if HAVE_FFTW + static const char FILTER_METHOD_FFTW_TITLESTR[]; + static const char FILTER_METHOD_RFFTW_TITLE_STR[]; +#endif + // Domains + static const char DOMAIN_FREQUENCY_STR[]; + static const char DOMAIN_SPATIAL_STR[]; - SignalFilter (const char* filterName, const char* filterMethodName,double bw, double signalIncrement, int n, double param, const char* domainName, const int zeropad = 0, const int numIntegral = 0); + static const char DOMAIN_FREQUENCY_TITLE_STR[]; + static const char DOMAIN_SPATIAL_TITLE_STR[]; - SignalFilter (const FilterID filt_type, FilterMethodID filterMethodID, double bw, double signalIncrement, int n, double param, const DomainID domain, const int zeropad = 0, const int numIntegral = 0); - SignalFilter (const char* filterName, const char* domainName, double bw, double param, int numIntegral = 0); + SignalFilter (const char* filterName, const char* filterMethodName,double bw, double signalIncrement, int n, double param, const char* domainName, const int zeropad = 0, const int preinterpolationFactor = 1); + + SignalFilter (const FilterID filt_type, FilterMethodID filterMethodID, double bw, double signalIncrement, int n, double param, const DomainID domain, const int zeropad = 0, const int preinterpolationFactor = 1); + + SignalFilter (const char* filterName, const char* domainName, double bw, double param); ~SignalFilter (void); @@ -132,14 +169,16 @@ class SignalFilter { double response (double x); - static double spatialResponse (FilterID fType, double bw, double x, double param, int nIntegral = 0); + static double spatialResponse (FilterID fType, double bw, double x, double param); static double frequencyResponse (FilterID fType, double bw, double u, double param); - static double spatialResponseCalc (FilterID fType, double bw, double x, double param, int n); - static double spatialResponseAnalytic (FilterID fType, double bw, double x, double param); + static double spatialResponseCalc (FilterID fType, double bw, double x, double param, int nIntegral); + + static void setNumIntegral(int nIntegral) {N_INTEGRAL = nIntegral;} + private: double m_bw; int m_nFilterPoints; @@ -153,9 +192,9 @@ class SignalFilter { double* m_vecFourierSinTable; complex* m_complexVecFilter; #ifdef HAVE_FFTW - fftw_real* m_vecRealFftInput; + fftw_real* m_vecRealFftInput, *m_vecRealFftSignal; rfftw_plan m_realPlanForward, m_realPlanBackward; - fftw_complex* m_vecComplexFftInput; + fftw_complex* m_vecComplexFftInput, *m_vecComplexFftSignal; fftw_plan m_complexPlanForward, m_complexPlanBackward; #endif @@ -168,10 +207,14 @@ class SignalFilter { FilterMethodID m_idFilterMethod; DomainID m_idDomain; double m_filterParam; - int m_numIntegral; int m_traceLevel; int m_zeropad; + int m_nOutputPoints; + int m_preinterpolationFactor; + + static int N_INTEGRAL; + static const bool haveAnalyticSpatial (const FilterID filterID); static const FilterID convertFilterNameToID (const char* filterName); static const char* convertFilterIDToName (const FilterID filterID); static const FilterMethodID convertFilterMethodNameToID (const char* filterMethodName); @@ -179,9 +222,9 @@ class SignalFilter { static const DomainID convertDomainNameToID (const char* domainName); static const char* convertDomainIDToName (const DomainID domainID); - void init (const FilterID filt_type, const FilterMethodID filterMethod, double bw, double signalIncrement, int n, double param, const DomainID domain, const int zeropad, const int numInt); + void init (const FilterID filt_type, const FilterMethodID filterMethod, double bw, double signalIncrement, int n, double param, const DomainID domain, const int zeropad, const int preInterpScale); -double spatialResponseCalc (double x, double param, int n) const; + double spatialResponseCalc (double x, double param) const; double spatialResponseAnalytic (double x, double param) const; @@ -194,5 +237,4 @@ double spatialResponseCalc (double x, double param, int n) const; }; - #endif