SpectrumProcess.cc
3.86 KB
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
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
/*
* To change this license header, choose License Headers in Project Properties.
* To change this template file, choose Tools | Templates
* and open the template in the editor.
*/
/*
* File: SpectrumProcess.cpp
* Author: hacene
*
* Created on December 31, 2020, 3:42 PM
*/
#include <stdlib.h>
#include <string>
#include <boost/lexical_cast.hpp>
#include "DicError.hh"
#include "AMDA_exception.hh"
#include "Operation.hh"
#include "ParameterManager.hh"
#include "ParamMgr.hh"
#include "ParamTable.hh"
#include "ParamInfo.hh"
#include "ParamMgr.hh"
#include "ParameterCreatorFromExpression.hh"
#include "SpectrumProcess.hh"
using namespace std;
using namespace boost;
using namespace log4cxx;
namespace AMDA {
namespace Parameters {
SpectrumProcess::SpectrumProcess(Parameter ¶meter) : SingleParamProcess_CRTP(parameter) {
_type = SpectrumCreator::SpectrumType::FOURIER;
_paramInfoId = parameter.getId();
parameter.setInfoId(_paramInfoId);
}
SpectrumProcess::SpectrumProcess(const SpectrumProcess& pProcess, Parameter ¶meter) : SingleParamProcess_CRTP(pProcess, parameter) {
_type = SpectrumCreator::SpectrumType::FOURIER;
_paramInfoId = parameter.getId();
parameter.setInfoId(_paramInfoId);
}
SpectrumProcess::~SpectrumProcess() {
}
TimeStamp SpectrumProcess::init() {
if (_attributList.size() < 1) {
BOOST_THROW_EXCEPTION(AMDA::AMDA_exception() << AMDA::errno_code(AMDA_PROCESS_ERR) << AMDA::ex_msg(std::string("ProcessSpectrum::parse require 1 attribute'")));
}
_windowtime = atof(_attributList[0].c_str());
_sampling = "center";
if (_attributList.size() > 1)
_sampling = _attributList[1];
/*****/
TimeStamp time = _parameterInput->init(this, _timeIntervalList);
_paramInput = _parameterInput->getParamData(this).get();
/** in case we need to sample
ParameterSPtr sampledParamInput = _parameterInput->getParameterManager().getSampledParameter(_parameterInput->getId(), "classic", _parameterInput->getGapThreshold(), false);
sampledParamInput->openConnection(this);
time = sampledParamInput->init(this, _timeIntervalList);
ParamData* sampledParamData = sampledParamInput->getParamData(this).get();
SpectrumCreator lCreator(*this, _timeIntervalList, *sampledParamData, _type, _windowtime, _sampling);
_operation = lCreator.getOperation();
*/
SpectrumCreator lCreator(*this, _timeIntervalList, *_paramInput, _type, _windowtime, _sampling);
_operation = lCreator.getOperation();
_paramData = ParamDataSPtr(_operation->getParamOutput());
_paramData->setMinSampling(_paramInput->getMinSampling());
int size_ = std::ceil(_windowtime/_paramInput->getMinSampling()) +1;
std::vector<double> frequencies(size_);
for(int i =0; i<size_; i++ )
frequencies[i] = (double) 1./_paramInput->getMinSampling()/(double) (i+1);
AMDA::Info::ParamInfoSPtr paramInfo = AMDA::Info::ParamMgr::getInstance()->getParamInfoFromId(_paramInfoId);
boost::shared_ptr<AMDA::Info::ParamTable> table(new AMDA::Info::EmbeddedCenterTable(
_paramInfoId.c_str(), 0, frequencies, false));
table->addTableParam(AMDA::Info::EmbeddedCenterTable::_paramKey, _paramInfoId );
table->setName("Frequencies");
table->setUnits("Hz");
table->setIsVariable(false);
table->setIsFullVariable(false);
paramInfo->addTable(0, table);
return time;
}
}
}