Skip to content
Snippets Groups Projects
Commit 034e7a79 authored by Christian Beckmann's avatar Christian Beckmann
Browse files

*** empty log message ***

parent 6579e3ff
No related branches found
No related tags found
No related merge requests found
...@@ -2,7 +2,7 @@ ...@@ -2,7 +2,7 @@
include ${FSLCONFDIR}/default.mk include ${FSLCONFDIR}/default.mk
OPTFLAGS = -O3 OPTFLAGS = -O3
PROJNAME = melodic2 PROJNAME = melodic2
......
...@@ -133,17 +133,19 @@ namespace Melodic{ ...@@ -133,17 +133,19 @@ namespace Melodic{
{//plot frequency {//plot frequency
miscplot newplot; miscplot newplot;
RowVector empty(1);
empty = 0.0;
int fact = int(std::pow(10.0,int(std::log10(float(melodat.get_mix().Nrows()))))); int fact = int(std::pow(10.0,int(std::log10(float(melodat.get_mix().Nrows())))));
if(opts.tr.value()>0.0) if(opts.tr.value()>0.0)
newplot.timeseries(melodat.get_fmix().Column(cnum).t(), newplot.timeseries(empty | melodat.get_fmix().Column(cnum).t(),
report.appendDir(string("f")+ report.appendDir(string("f")+
num2str(cnum)+".png"), num2str(cnum)+".png"),
string("FFT of timecourse (in Hz / ") +num2str(fact)+")", string("FFT of timecourse (in Hz / ") +num2str(fact)+")",
fact/(opts.tr.value()*melodat.get_mix().Nrows()), 150, fact/(opts.tr.value()*melodat.get_mix().Nrows()), 150,
0,2); 0,2);
else else
newplot.timeseries(melodat.get_fmix().Column(cnum).t(), newplot.timeseries(empty | melodat.get_fmix().Column(cnum).t(),
report.appendDir(string("f")+ report.appendDir(string("f")+
num2str(cnum)+".png"), num2str(cnum)+".png"),
string(string("FFT of timecourse (in cycles); ") string(string("FFT of timecourse (in cycles); ")
...@@ -154,7 +156,7 @@ namespace Melodic{ ...@@ -154,7 +156,7 @@ namespace Melodic{
+"* TR)/cycles")); +"* TR)/cycles"));
write_ascii_matrix(report.appendDir(string("f") write_ascii_matrix(report.appendDir(string("f")
+num2str(cnum)+".txt"), +num2str(cnum)+".txt"),
melodat.get_mix().Column(cnum)); melodat.get_fmix().Column(cnum));
IChtml << "<A HREF=\"" << string("f") IChtml << "<A HREF=\"" << string("f")
+num2str(cnum)+".txt" << "\"> "; +num2str(cnum)+".txt" << "\"> ";
IChtml << "<img BORDER=0 SRC=\"" IChtml << "<img BORDER=0 SRC=\""
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment