00001 /* 00002 * Copyright (c) 2001-2004 MUSIC TECHNOLOGY GROUP (MTG) 00003 * UNIVERSITAT POMPEU FABRA 00004 * 00005 * 00006 * This program is free software; you can redistribute it and/or modify 00007 * it under the terms of the GNU General Public License as published by 00008 * the Free Software Foundation; either version 2 of the License, or 00009 * (at your option) any later version. 00010 * 00011 * This program is distributed in the hope that it will be useful, 00012 * but WITHOUT ANY WARRANTY; without even the implied warranty of 00013 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 00014 * GNU General Public License for more details. 00015 * 00016 * You should have received a copy of the GNU General Public License 00017 * along with this program; if not, write to the Free Software 00018 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA 00019 * 00020 */ 00021 00022 #include "SpectralCombDelta.hxx" 00023 #include "ProcessingFactory.hxx" 00024 00025 namespace CLAM 00026 { 00027 00028 namespace Hidden 00029 { 00030 static const char * metadata[] = { 00031 "key", "SpectralCombDelta", 00032 // "category", "Spectral Transformations", 00033 // "description", "SpectralCombDelta", 00034 0 00035 }; 00036 static FactoryRegistrator<ProcessingFactory, SpectralCombDelta> reg = metadata; 00037 } 00038 00039 bool SpectralCombDelta::Do(const Spectrum& in, Spectrum& out) 00040 { 00041 if (!mConfig.GetPreserveOuts()) 00042 { 00043 out = in; //TODO big cludge for streaming 00044 } 00045 DataArray& inMag = in.GetMagBuffer(); 00046 DataArray& outMag = out.GetMagBuffer(); 00047 00048 int spectrumSize = in.GetSize(); 00049 mMag.Resize(spectrumSize); 00050 mMag.SetSize(spectrumSize); 00051 00052 TData spectralResolution = spectrumSize/in.GetSpectralRange(); 00053 int fundamental = Round(mFundamentalCtl.GetLastValue()* spectralResolution); 00054 TData overallFactor = log2lin(mLogOverallFactor.GetLastValue()); 00055 int i; 00056 int n; 00057 //first we set everything to 0 (it could be done more efficiently) 00058 for(i = 0; i<spectrumSize; i++) 00059 { 00060 mMag[i] = 0.; 00061 } 00062 n = 1; 00063 for(i = 0; i<spectrumSize; i++) 00064 { 00065 if(i == int (TData(fundamental) * n)) 00066 { 00067 mMag[i] = inMag[i]*overallFactor; 00068 n++; 00069 } 00070 } 00071 out.SetMagBuffer(mMag); 00072 return true; 00073 } 00074 00075 /*bool SpectralCombDelta::Do(const Spectrum& in, Spectrum& out) 00076 { 00077 if (!mConfig.GetPreserveOuts()) 00078 { 00079 out = in; //TODO big cludge for streaming 00080 } 00081 DataArray& inMag = in.GetMagBuffer(); 00082 DataArray& outMag = out.GetMagBuffer(); 00083 00084 int spectrumSize = in.GetSize(); 00085 mMag.Resize(spectrumSize); 00086 mMag.SetSize(spectrumSize); 00087 00088 TData spectralResolution = spectrumSize/in.GetSpectralRange(); 00089 TData freqFactor = mFreqSpaceFactorCtl.GetLastValue(); 00090 int fundamental = Round(mFundamentalCtl.GetLastValue()* spectralResolution); 00091 TData overallFactor = log2lin(mLogOverallFactor.GetLastValue()); 00092 int i; 00093 TData n; 00094 //first we set everything to 0 (it could be done more efficiently) 00095 for(i = 0; i<spectrumSize; i++) 00096 { 00097 mMag[i] = 0.; 00098 } 00099 n = 1.; 00100 for(i = 0; i<spectrumSize; i++) 00101 { 00102 if(i == int (TData(fundamental) * n)) 00103 { 00104 mMag[i] = inMag[i]*overallFactor; 00105 n*=freqFactor; 00106 } 00107 } 00108 out.SetMagBuffer(mMag); 00109 return true; 00110 } 00111 */ 00112 00113 00114 } 00115