from qsiprep.workflows.anatomical import init_dwi_reference_wf
wf = init_dwi_reference_wf(omp_nthreads=1,
                           gen_report=True,
                           source_file="sub-1_dwi.nii.gz",
                           register_t1=True)