/usr/include/mia-2.4/mia/core/ica.hh is in libmia-2.4-dev 2.4.3-5.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
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 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 | /* -*- mia-c++ -*-
*
* This file is part of MIA - a toolbox for medical image analysis
* Copyright (c) Leipzig, Madrid 1999-2016 Gert Wollny
*
* MIA is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation; either version 3 of the License, or
* (at your option) any later version.
*
* This program is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with MIA; if not, see <http://www.gnu.org/licenses/>.
*
*/
#ifndef mia_core_ica_hh
#define mia_core_ica_hh
#include <set>
#include <vector>
#include <itpp/signal/fastica.h>
#include <mia/core/defines.hh>
#include <mia/core/slopevector.hh>
#include <mia/core/icaanalysisbase.hh>
#include <boost/concept/requires.hpp>
#include <boost/concept_check.hpp>
NS_MIA_BEGIN
/**
\ingroup perf
\brief class for ICA analysis and use of such data.
This class implements basic operations for of ICA. It makes use of the ITPP implementation of FastICA.
*/
class EXPORT_CORE CICAAnalysisITPP : public CICAAnalysis {
public:
typedef CICAAnalysis::IndexSet IndexSet;
/// The type of a vector as used by IT++
typedef itpp::Vec<itpp::mat::value_type> itppvector;
/**
Initialize an ICA based of predefined data - this is only used for test cases.
*/
CICAAnalysisITPP(const itpp::mat& ic, const itpp::mat& mix, const std::vector<double>& mean );
CICAAnalysisITPP();
~CICAAnalysisITPP();
/**
Initialize the ICA
\param series_length number of data sets that will be provided
\param slice_size number of elements each set containes
*/
void initialize(unsigned int series_length, unsigned int slice_size);
/**
Run the independed component analysis using the given numbers of components
\param nica number of indentepended components
\param guess initial guess for the ICA, pass an empty vector of you
don't want to use this feature
*/
bool run(unsigned int nica, std::vector<std::vector<float> > guess);
/// \returns the feature vector of \a row
std::vector<float> get_feature_row(unsigned int row)const;
/// \returns the mixing signal curve of the feature \a row
std::vector<float> get_mix_series(unsigned int row)const;
/// \returns the complete mixed signal at series index \a idx
std::vector<float> get_mix(unsigned int idx)const;
/** Evaluate an incomplete mixed signal. Here the features are given that are \a not to be used.
\sa get_partial_mix
\param idx series index
\param skip a set of feature indices that will be skipped when evaluating the mix
\returns the mixed signal
*/
std::vector<float> get_incomplete_mix(unsigned int idx, const IndexSet& skip)const;
/** Evaluate an incomplete mixed signal. Here the features are given that are \a used to create the mix.
\sa get_incomplete_mix
\param idx series index
\param use the set of feature indices that will be used to evaluate the mix
\returns an incolmplete mixed signal.
*/
std::vector<float> get_partial_mix(unsigned int idx, const IndexSet& use)const;
/** Evaluate a mix of the feature signals by adding and subtractig individual features.
\param plus features o be added
\param minus features to be subtracted
\returns the feature mix
*/
std::vector<float> get_delta_feature(const IndexSet& plus, const IndexSet& minus)const;
/**
Replace a mixing signal curve
\param index of the curve to be replaced
\param series new data for mixing curve
*/
void set_mixing_series(unsigned int index, const std::vector<float>& series);
/// \returns a vector containing all mixing curves
CSlopeColumns get_mixing_curves() const;
/**
Normalize the ICs in the following manner:
* Scale and shift the range of the ICs to [-1, 1]
* Scale the mixing curved to compensate for the required scaling
* move the means of the time points to compensate for the shifting.
*/
void normalize_ICs();
/**
Normalize the mixing curves to have a zero mean. As a result a mean image is created that
containes the sum of the ICs weighted by the required mean shift.
*/
std::vector<float> normalize_Mix();
/// \returns the number of actual ICs
unsigned int get_ncomponents() const;
/**
sets the number of iterations in the ICA
\param n
*/
void set_max_iterations(int n);
/**
Set the ICA approach to either FICA_APPROACH_DEFL(default) or FICA_APPROACH_SYMM.
\param approach
*/
void set_approach(EApproach approach);
private:
void set_row_internal(unsigned row, const std::vector<double>& buffer, double mean);
struct CICAAnalysisITPPImpl *impl;
};
class EXPORT_CORE CICAAnalysisITPPFactory: public CICAAnalysisFactory {
public:
PICAAnalysis create() const;
};
typedef std::shared_ptr<CICAAnalysisITPPFactory> PICAAnalysisITPPFactory;
NS_MIA_END
#endif
|