SpectralNotch.cxx
Go to the documentation of this file.00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022 #include "SpectralNotch.hxx"
00023 #include "ProcessingFactory.hxx"
00024
00025 namespace CLAM
00026 {
00027
00028 namespace Hidden
00029 {
00030 static const char * metadata[] = {
00031 "key", "SpectralNotch",
00032
00033
00034 0
00035 };
00036 static FactoryRegistrator<ProcessingFactory, SpectralNotch> reg = metadata;
00037 }
00038
00039 bool SpectralNotch::Do(const Spectrum& in, Spectrum& out)
00040 {
00041 out = in;
00042
00043 DataArray& outMag = out.GetMagBuffer();
00044
00045 int spectrumSize = in.GetSize();
00046
00047 TData spectralResolution = spectrumSize/in.GetSpectralRange();
00048 int centerFreq = Round(mFreq.GetLastValue()*spectralResolution);
00049 int band = Round(mBandwidthCtl.GetLastValue()*spectralResolution);
00050 TData gain = log2lin(mGainCtl.GetLastValue());
00051 int n;
00052 int harmonic = 1;
00053 int leftLimit = centerFreq-band;
00054 if(leftLimit<0) leftLimit = 0;
00055 int rightLimit = centerFreq+band+1;
00056 if(rightLimit>spectrumSize) rightLimit=spectrumSize;
00057 for (n=leftLimit; n<rightLimit;n++)
00058 {
00059 outMag[n] *= gain ;
00060 }
00061
00062 outMag[leftLimit-1] *= (gain*0.5);
00063 outMag[rightLimit] *= (gain*0.5);
00064 return true;
00065 }
00066
00067
00068 }
00069