55 #include <fiff/fiff_evoked.h>
65 #include <QtCore/QCoreApplication>
92 int main(
int argc,
char *argv[])
94 QCoreApplication a(argc, argv);
109 QFile t_fileEvoked(
"E:/Data/sl_data/MEG/mind006/mind006_051209_auditory01_raw-ave.fif");
110 QFile t_fileInv(
"E:/Data/sl_data/MEG/mind006/mind006_051209_auditory01_raw-oct-6p-meg-inv.fif");
113 QString method(
"dSPM");
114 QString t_sFileNameStc(
"");
117 for(qint32 i = 0; i < argc; ++i)
119 if(strcmp(argv[i],
"-snr") == 0 || strcmp(argv[i],
"--snr") == 0)
122 snr = atof(argv[i+1]);
124 else if(strcmp(argv[i],
"-method") == 0 || strcmp(argv[i],
"--method") == 0)
127 method = QString::fromUtf8(argv[i+1]);
129 else if(strcmp(argv[i],
"-stc") == 0 || strcmp(argv[i],
"--stc") == 0)
132 t_sFileNameStc = QString::fromUtf8(argv[i+1]);
136 double lambda2 = 1.0 / pow(snr, 2);
137 qDebug() <<
"Start calculation with: SNR" << snr <<
"; Lambda" << lambda2 <<
"; Method" << method <<
"; stc:" << t_sFileNameStc;
142 fiff_int_t setno = 0;
143 QPair<QVariant, QVariant> baseline(QVariant(), 0);
144 FiffEvoked evoked(t_fileEvoked, setno, baseline);
156 MinimumNorm minimumNorm(inverse_operator, lambda2, method);
162 std::cout <<
"\npart ( block( 0, 0, 10, 10) ) of the inverse solution:\n" << sourceEstimate.
data.block(0,0,10,10) << std::endl;
163 printf(
"tmin = %f s\n", sourceEstimate.
tmin);
164 printf(
"tstep = %f s\n", sourceEstimate.
tstep);
167 if(!t_sFileNameStc.isEmpty())
169 QFile t_fileStc(t_sFileNameStc);
170 sourceEstimate.
write(t_fileStc);
175 std::cout <<
"\npart ( block( 0, 0, 10, 10) ) of the inverse solution:\n" << readSourceEstimate.data.block(0,0,10,10) << std::endl;
176 printf(
"tmin = %f s\n", readSourceEstimate.tmin);
177 printf(
"tstep = %f s\n", readSourceEstimate.tstep);
Minimum norm class declaration.
bool write(QIODevice &p_IODevice)
MNEInverseOperator class declaration.
FiffEvokedSet class declaration.
MNESourceEstimate class declaration.