fMRIPrep: how to resample data from T1w space into MNI space

In this situation I would use the outputs already on MNI space, and for the warp, an all-zeros XxYxZx3 NIfTI with the same shape and affine of the input.

import numpy as np
import nibabel as nb

im = nb.load('<bids_entities>_space-MNI152NLin2009cAsym_desc-preproc_bold.nii.gz')
hdr = im.header.copy()
hdr.set_intent(2006)  # Set NIFTI intent code to "FSL_FNIRT_DISPLACEMENT_FIELD"
nb.Nifti1Image(np.zeros(list(im.shape[:3]) + [3], dtype=float), im.affine, hdr).to_filename('zero_warpfield.nii.gz')