161
edits
Line 50: | Line 50: | ||
3dDFT -prefix {$ec}r{$r}vrdtblFFT -detrend junkr{$r}+orig | 3dDFT -prefix {$ec}r{$r}vrdtblFFT -detrend junkr{$r}+orig | ||
end | end | ||
The retinotopic stimulus has a periodicity at frequency 5, so calculate phase at only that frequency. Then, calculate amplitude squared at all frequencies | |||
foreach r (1 2 3 4 5 6 7 8) | |||
3dcalc -datum short -prefix {$ec}r{$r}vrdtblFFTf5 -cx2r PHASE -a {$ec}r{$r}vrdtblFFT+orig'[5]' -expr "a" | |||
3drefit -sublabel 0 FFTr{$r}f5 -subrepkey 0 FFTr{$r}f5 {$ec}r{$r}vrdtblFFTf5+orig | |||
3dcalc -datum short -prefix {$ec}r{$r}vrdtblFFTMAG -cx2r ABS -a {$ec}r{$r}vrdtblFFT+orig -expr "a*a" | 3dcalc -datum short -prefix {$ec}r{$r}vrdtblFFTMAG -cx2r ABS -a {$ec}r{$r}vrdtblFFT+orig -expr "a*a" | ||
end | end |
edits