@@ -338,10 +338,10 @@ function create_timecourses(self,varargin)
338
338
data = zscore(data );
339
339
mag_d = sqrt(sum(data .^ 2 ));
340
340
341
- results.R = zeros(self .n_total ,1 );
342
- results.X = zeros(self .n_total ,1 );
343
- results.Y = zeros(self .n_total ,1 );
344
- results.Sigma = zeros(self .n_total ,1 );
341
+ results.corr_fit = zeros(self .n_total ,1 );
342
+ results.mu_x = zeros(self .n_total ,1 );
343
+ results.mu_y = zeros(self .n_total ,1 );
344
+ results.sigma = zeros(self .n_total ,1 );
345
345
346
346
if size(self .hrf ,2 )==1
347
347
hrf_fft = fft(repmat([self .hrf ;...
@@ -355,12 +355,12 @@ function create_timecourses(self,varargin)
355
355
(mag_tc * mag_d(v ));
356
356
id = isinf(CS ) | isnan(CS );
357
357
CS(id ) = 0 ;
358
- [results .R (v ),j ] = max(CS );
359
- results .X (v ) = cos(self .pa(self .idx(j ,1 ))) * ...
358
+ [results .corr_fit (v ),j ] = max(CS );
359
+ results .mu_x (v ) = cos(self .pa(self .idx(j ,1 ))) * ...
360
360
self .ecc(self .idx(j ,2 ));
361
- results .Y (v ) = sin(self .pa(self .idx(j ,1 ))) * ...
361
+ results .mu_y (v ) = sin(self .pa(self .idx(j ,1 ))) * ...
362
362
self .ecc(self .idx(j ,2 ));
363
- results .Sigma (v ) = self .ecc(self .idx(j ,2 )) * ...
363
+ results .sigma (v ) = self .ecc(self .idx(j ,2 )) * ...
364
364
self .slope(self .idx(j ,3 ));
365
365
end
366
366
waitbar(v / self .n_total ,wb )
@@ -379,23 +379,23 @@ function create_timecourses(self,varargin)
379
379
(mag_tc * mag_d(v ));
380
380
id = isinf(CS ) | isnan(CS );
381
381
CS(id ) = 0 ;
382
- [results .R (v ),j ] = max(CS );
383
- results .X (v ) = cos(self .pa(self .idx(j ,1 ))) * ...
382
+ [results .corr_fit (v ),j ] = max(CS );
383
+ results .mu_x (v ) = cos(self .pa(self .idx(j ,1 ))) * ...
384
384
self .ecc(self .idx(j ,2 ));
385
- results .Y (v ) = sin(self .pa(self .idx(j ,1 ))) * ...
385
+ results .mu_y (v ) = sin(self .pa(self .idx(j ,1 ))) * ...
386
386
self .ecc(self .idx(j ,2 ));
387
- results .Sigma (v ) = self .ecc(self .idx(j ,2 )) * ...
387
+ results .sigma (v ) = self .ecc(self .idx(j ,2 )) * ...
388
388
self .slope(self .idx(j ,3 ));
389
389
end
390
390
waitbar(v / self .n_total ,wb )
391
391
end
392
392
end
393
- results.R = reshape(results .R ,self .n_rows ,self .n_cols ,self .n_slices );
394
- results.X = reshape(results .X ,self .n_rows ,self .n_cols ,self .n_slices );
395
- results.Y = reshape(results .Y ,self .n_rows ,self .n_cols ,self .n_slices );
396
- results.Sigma = reshape(results .Sigma ,self .n_rows ,self .n_cols ,self .n_slices );
397
- results.Eccentricity = abs(results .X + results .Y * 1i );
398
- results.Polar_Angle = angle(results .X + results .Y * 1i );
393
+ results.corr_fit = reshape(results .corr_fit ,self .n_rows ,self .n_cols ,self .n_slices );
394
+ results.mu_x = reshape(results .mu_x ,self .n_rows ,self .n_cols ,self .n_slices );
395
+ results.mu_y = reshape(results .mu_y ,self .n_rows ,self .n_cols ,self .n_slices );
396
+ results.sigma = reshape(results .sigma ,self .n_rows ,self .n_cols ,self .n_slices );
397
+ results.Eccentricity = abs(results .mu_x + results .mu_y * 1i );
398
+ results.Polar_Angle = angle(results .mu_x + results .mu_y * 1i );
399
399
close(wb )
400
400
fprintf(' done\n ' );
401
401
end
0 commit comments