from qsiprep.workflows.dwi.hmc_sdc import init_qsiprep_hmcsdc_wf
wf = init_qsiprep_hmcsdc_wf({'dwi_series':['dwi1.nii', 'dwi2.nii'],
                            'fieldmap_info': {'suffix': None},
                            'dwi_series_pedir': 'j'},
                            source_file='/data/sub-1/dwi/sub-1_dwi.nii.gz',
                            b0_threshold=100,
                            hmc_transform='Affine',
                            hmc_model='3dSHORE',
                            hmc_align_to='iterative',
                            template='MNI152NLin2009cAsym',
                            shoreline_iters=1,
                            impute_slice_threshold=0,
                            omp_nthreads=1,
                            fmap_bspline=False,
                            fmap_demean=False,
                            use_syn=True,
                            force_syn=False,
                            name='qsiprep_hmcsdc_wf',
                            dwi_metadata={})