def create_polynomial_drift(bold_shape, tr, drift_order, drift_var): p = buildPolyMat(drift_order, bold_shape[0], tr) nvox = bold_shape[1] coeff = drift_var**.5 * randn(p.shape[1], nvox) drift = np.dot(p, coeff) logger.info('Drift shape: %s', str(drift.shape)) return drift
def create_polynomial_drift(bold_shape, tr, drift_order, drift_var): p = buildPolyMat(drift_order, bold_shape[0], tr) nvox = bold_shape[1] coeff = drift_var ** 0.5 * randn(p.shape[1], nvox) drift = np.dot(p, coeff) logger.info("Drift shape: %s", str(drift.shape)) return drift
def create_polynomial_drift(bold_shape, tr, drift_order, drift_var): p = buildPolyMat(drift_order, bold_shape[0], tr) nvox = bold_shape[1] coeff = drift_var**.5 * randn(p.shape[1],nvox) drift = np.dot(p, coeff) pyhrf.verbose(3, 'Drift shape: %s' %str(drift.shape)) return drift
def create_polynomial_drift_from_coeffs(bold_shape, tr, drift_order, drift_coeffs, drift_mean=0.0, drift_amplitude=1.0): p = buildPolyMat(drift_order, bold_shape[0], tr) if 0: print "%%%%%%%%%%%%%%%%%%%%%%%%%" print "p:", p print "%%%%%%%%%%%%%%%%%%%%%%%%%" print "drift_coeffs:", drift_coeffs.shape drift = np.dot(p, drift_coeffs) * drift_amplitude + drift_mean logger.info("Drift shape: %s", str(drift.shape)) return drift
def create_polynomial_drift_from_coeffs(bold_shape, tr, drift_order, drift_coeffs, drift_mean=0., drift_amplitude=1.): p = buildPolyMat(drift_order, bold_shape[0], tr) if 0: print '%%%%%%%%%%%%%%%%%%%%%%%%%' print 'p:', p print '%%%%%%%%%%%%%%%%%%%%%%%%%' print 'drift_coeffs:', drift_coeffs.shape drift = np.dot(p, drift_coeffs) * drift_amplitude + drift_mean logger.info('Drift shape: %s', str(drift.shape)) return drift