SpectralReverse.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 "SpectralReverse.hxx"
00023 #include "ProcessingFactory.hxx"
00024
00025 namespace CLAM
00026 {
00027
00028 namespace Hidden
00029 {
00030 static const char * metadata[] = {
00031 "key", "SpectralReverse",
00032
00033
00034 0
00035 };
00036 static FactoryRegistrator<ProcessingFactory, SpectralReverse> reg = metadata;
00037 }
00038
00039 bool SpectralReverse::Do(const Spectrum& in, Spectrum& out)
00040 {
00041 if (!mConfig.GetPreserveOuts())
00042 {
00043 out = in;
00044 }
00045 DataArray& inPhase = in.GetPhaseBuffer();
00046 DataArray& inMag = in.GetMagBuffer();
00047
00048 int spectrumSize = in.GetSize();
00049
00050 TData spectralResolution = spectrumSize/in.GetSpectralRange();
00051
00052 int mirrorPoint = Round(mAmount.GetLastValue()* spectralResolution);
00053
00054 mMag.Resize(spectrumSize);
00055 mMag.SetSize(spectrumSize);
00056 mPhase.Resize(spectrumSize);
00057 mPhase.SetSize(spectrumSize);
00058 int i;
00059 for(i = 0; i<mirrorPoint; i++)
00060 {
00061 mMag[i] = inMag[mirrorPoint-i];
00062 mPhase[i] = inPhase[mirrorPoint-i];
00063 }
00064 for(i = mirrorPoint; i<spectrumSize; i++)
00065 {
00066 mMag[i] = inMag[spectrumSize-i+mirrorPoint];
00067 mPhase[i] = inPhase[spectrumSize-i+mirrorPoint];
00068 }
00069 out.SetMagBuffer(mMag);
00070 out.SetPhaseBuffer(mPhase);
00071 return true;
00072 }
00073
00074
00075 }
00076