예제 #1
0
def demo_fsl_feeds(output_dir):
    """Demo for FSL Feeds data.

    Parameters
    ----------
    output_dir: string
        where output will be written to

    """
    output_dir = os.path.join(output_dir, "fsl_feeds_mrimc_output")
    fsl_feeds = fetch_fsl_feeds()
    subject_id = "sub001"
    subjects = [SubjectData(subject_id=subject_id,
                            func=fsl_feeds.func,
                            output_dir=os.path.join(output_dir, subject_id))]
    _demo_runner(subjects, "FSL FEEDS")
예제 #2
0
def demo_fsl_feeds(output_dir):
    """Demo for FSL Feeds data.

    Parameters
    ----------
    output_dir: string
        where output will be written to

    """
    output_dir = os.path.join(output_dir, "fsl_feeds_mrimc_output")
    fsl_feeds = fetch_fsl_feeds()
    subject_id = "sub001"
    subjects = [
        SubjectData(subject_id=subject_id,
                    func=fsl_feeds.func,
                    output_dir=os.path.join(output_dir, subject_id))
    ]
    _demo_runner(subjects, "FSL FEEDS")
예제 #3
0
def demo_fsl_feeds(output_dir="/tmp/fsl_feeds_mrimc_output"):
    """Demo for FSL Feeds data.

    Parameters
    ----------
    data_dir: string, optional
        where the data is located on your disk, where it will be
        downloaded to
    output_dir: string, optional
        where output will be written to

    """
    fsl_feeds = fetch_fsl_feeds()
    subject_id = "sub001"
    subjects = [SubjectData(subject_id=subject_id,
                            func=fsl_feeds.func,
                            output_dir=os.path.join(output_dir, subject_id))]
    _demo_runner(subjects, "FSL FEEDS")
예제 #4
0
def demo_fsl_feeds(output_dir="/tmp/fsl_feeds_mrimc_output"):
    """Demo for FSL Feeds data.

    Parameters
    ----------
    data_dir: string, optional
        where the data is located on your disk, where it will be
        downloaded to
    output_dir: string, optional
        where output will be written to

    """
    fsl_feeds = fetch_fsl_feeds()
    subject_id = "sub001"
    subjects = [
        SubjectData(subject_id=subject_id,
                    func=fsl_feeds.func,
                    output_dir=os.path.join(output_dir, subject_id))
    ]
    _demo_runner(subjects, "FSL FEEDS")
                         'duration': duration})
frametimes = np.linspace(0, (n_scans - 1) * TR, n_scans)
maximum_epoch_duration = max(EV1_epoch_duration, EV2_epoch_duration)
hfcut = 1.5 * maximum_epoch_duration  # why ?

"""construct design matrix"""
drift_model = 'Cosine'
hrf_model = 'Canonical With Derivative'
design_matrix = make_design_matrix(frame_times=frametimes,
                                   paradigm=paradigm,
                                   hrf_model=hrf_model,
                                   drift_model=drift_model,
                                   period_cut=hfcut)

"""fetch input data"""
_subject_data = fetch_fsl_feeds()
subject_data = SubjectData()
subject_data.subject_id = "sub001"
subject_data.func = _subject_data.func
subject_data.anat = _subject_data.anat

output_dir = os.path.join(_subject_data.data_dir, "pypreprocess_output")
if not os.path.exists(output_dir):
    os.makedirs(output_dir)
subject_data.output_dir = os.path.join(
    output_dir, subject_data.subject_id)



"""preprocess the data"""
results = do_subjects_preproc(
예제 #6
0
                         'duration': duration})
frametimes = np.linspace(0, (n_scans - 1) * TR, n_scans)
maximum_epoch_duration = max(EV1_epoch_duration, EV2_epoch_duration)
hfcut = 1.5 * maximum_epoch_duration  # why ?

"""construct design matrix"""
drift_model = 'Cosine'
hrf_model = 'spm + derivative'
design_matrix = make_design_matrix(frame_times=frametimes,
                                   paradigm=paradigm,
                                   hrf_model=hrf_model,
                                   drift_model=drift_model,
                                   period_cut=hfcut)

"""fetch input data"""
_subject_data = fetch_fsl_feeds()
subject_data = SubjectData()
subject_data.subject_id = "sub001"
subject_data.func = _subject_data.func
subject_data.anat = _subject_data.anat

output_dir = os.path.join(_subject_data.data_dir, "pypreprocess_output")
if not os.path.exists(output_dir):
    os.makedirs(output_dir)
subject_data.output_dir = os.path.join(
    output_dir, subject_data.subject_id)



"""preprocess the data"""
results = do_subjects_preproc(