Questions regarding three fMRI datasets with different TR

Dear neuroimaging experts,

I have three fMRI datasets preprocessed with fMRIprep having different volumes due to different TR: 1) 91x109x91x192 (TR=2510, slice thickness=3); 2) 91x109x91x570 (TR= 720, slice thickness=2); 3) 91X109X91X240 (TR=2000, slice thickness=2.5). I would like to use neurocombat to harmonise data which requires converting the 4D image to a vector. Given neuroCombat requires the same vector length, I wonder how can I resclice/resample the data to have the same vector length. Any suggestions are appriciated!

Best wishes
Joyce