Newer
Older
/* MELODIC - Multivariate exploratory linear optimized decomposition into
independent components
melgmix.h - class for Gaussian/Gamma Mixture Model
Christian F. Beckmann, FMRIB Image Analysis Group
Copyright (C) 1999-2002 University of Oxford */
/* CCOPYRIGHT */
#ifndef __MELGMIX_h
#define __MELGMIX_h
#include "newimage/newimageall.h"
#include "utils/log.h"
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
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
#include "meloptions.h"
//#include "melreport.h"
using namespace Utilities;
using namespace NEWIMAGE;
namespace Melodic{
class MelGMix
{
public:
//constructor
/* MelGMix(MelodicOptions &popts, Log &plogger): */
/* opts(popts), */
/* logger(plogger), */
/* mainhtml() */
/* { */
/* } */
MelGMix(MelodicOptions &popts, Log &plogger):
opts(popts),
logger(plogger)
{
}
~MelGMix() {
//mainhtml << endl << "<hr></CENTER></BODY></HTML>" << endl;
}
void save();
void setup(const RowVector& dat, volumeinfo inf, const string dirname,
int here, volume<float> themask,
volume<float> themean, int num_mix = 3,
float eps = 0.0, bool fixdim = false);
void gmmfit();
void ggmfit();
inline void fit(string mtype = string("GGM"))
{
mmtype = mtype;
if(mmtype==string("GGM"))
this->ggmfit();
else
this->gmmfit();
}
inline Matrix threshold(string levels)
{return this->threshold(data, levels);}
inline Matrix threshold(RowVector& levels)
{return this->threshold(data, levels);}
Matrix threshold(const RowVector& dat, Matrix& levels);
Matrix threshold(const RowVector& dat, string levels);
void status(const string &txt);
inline RowVector& get_means() {return means;}
inline void set_means(RowVector& Arg) {means = Arg;}
inline RowVector& get_vars() {return vars;}
inline void set_vars(RowVector& Arg) {vars = Arg;}
inline RowVector& get_pi() {return props;}
inline void set_pi(RowVector& Arg) {props = Arg;}
inline RowVector& get_data() {return data;}
inline void set_data(RowVector& Arg) {data = Arg;}
inline RowVector& get_prob() {return probmap;}
inline float get_eps() {return epsilon;}
inline void set_eps(float Arg) {epsilon = Arg;}
inline Matrix& get_threshmaps() {return threshmaps;}
inline void set_threshmaps(Matrix& Arg) {threshmaps = Arg;}
inline bool isfitted(){return fitted;}
inline int mixtures(){return nummix;}
inline string get_type() { return mmtype;}
inline void set_type(string Arg) { mmtype = Arg;}
inline string get_prefix() { return prefix;}
inline void set_prefix(string Arg) { prefix = Arg;}
inline RowVector get_probmap() {return probmap;}
inline float get_offset() {return offset;}
inline void set_offset(float Arg) {offset = Arg;}
inline void flipres(int num){
means = -means;
data = -data;
threshmaps = -threshmaps;
if(mmtype=="GGM"){
float tmp;
tmp= means(2);means(2)= means(3);means(3)=tmp;
tmp=vars(2);vars(2)=vars(3);vars(3)=tmp;
tmp=props(2);props(2)=props(3);props(3)=tmp;
}
}
void create_rep();
inline void add_infstr(string what){
threshinfo.push_back(what);
}
inline string get_infstr(int num){
if((threshinfo.size()<(unsigned int)(num-1))||(num<1))
return string("");
else
return threshinfo[num-1];
}
inline int size_infstr(){
return threshinfo.size();
}
inline void clear_infstr(){
threshinfo.clear();
}
private:
MelodicOptions &opts;
Log &logger; //global log file
//Log mainhtml;
void gmmupdate();
float gmmevidence();
void gmmreducemm();
void add_params(Matrix& mu, Matrix& sig, Matrix& pi,
float logLH, float MDL, float Evi, bool advance = false);
void get_params(int index, Matrix& mu, Matrix& sig, Matrix& pi,
float logLH, float MDL, float Evi);
Matrix Params;
Matrix threshmaps;
RowVector means;
RowVector vars;
RowVector props;
RowVector data;
RowVector probmap;
volume<float> Mean;
volume<float> Mask;
float epsilon;
float logprobY;
float MDL;
float Evi;
float offset;
int nummix;
int numdata;
int cnumber;
bool fitted;
bool fixdim;
string prefix;
string mmtype;
string dirname;
volumeinfo bginfo;
vector<string> threshinfo;
};
}
#endif