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,
)