Node: qsiprep_1_0_wf.sub_FACT130_ses_post_pre_wf.dwi_preproc_ses_pre_wf.pre_hmc_wf.merge_and_denoise_wf.merge_dwis
Working directory: /work/qsiprep_1_0_wf/sub_FACT130_ses_post_pre_wf/dwi_preproc_ses_pre_wf/pre_hmc_wf/merge_and_denoise_wf/merge_dwis

Node inputs:

b0_refs = <undefined>
b0_threshold = 100
bids_dwi_files = ['/data/sub-FACT130/ses-pre/dwi/sub-FACT130_ses-pre_dwi.nii.gz']
bval_files = <undefined>
bvec_files = <undefined>
carpetplot_data = <undefined>
denoising_confounds = <undefined>
dwi_files = <undefined>
harmonize_b0_intensities = True
raw_concatenated_files = <undefined>
scan_metadata = {'/data/sub-FACT130/ses-pre/dwi/sub-FACT130_ses-pre_dwi.nii.gz': {'AcquisitionMatrixPE': 128, 'AcquisitionNumber': 1, 'AcquisitionTime': '12:21:16.752500', 'BandwidthPerPixelPhaseEncode': 19.531, 'BaseResolution': 128, 'ConversionSoftware': 'dcm2niix', 'ConversionSoftwareVersion': 'v1.0.20180622 GCC6.3.0', 'DerivedVendorReportedEchoSpacing': 0.00080001, 'DeviceSerialNumber': '35469', 'DwellTime': 2.8e-06, 'EchoTime': 0.057, 'EffectiveEchoSpacing': 0.000400005, 'FlipAngle': 90, 'ImageOrientationPatientDICOM': [0.996543, 0.0807755, -0.0194357, -0.0824693, 0.990091, -0.113657], 'ImageType': ['ORIGINAL', 'PRIMARY', 'M', 'ND', 'MOSAIC'], 'InPlanePhaseEncodingDirectionDICOM': 'COL', 'InstitutionAddress': 'Vassar_43_Cambridge_MA_US_02139', 'InstitutionName': 'McGovern_Brain_Institute', 'InstitutionalDepartmentName': 'Department', 'MRAcquisitionType': '2D', 'MagneticFieldStrength': 3, 'Manufacturer': 'Siemens', 'ManufacturersModelName': 'TrioTim', 'Modality': 'MR', 'ParallelReductionFactorInPlane': 2, 'PartialFourier': 0.75, 'PatientPosition': 'HFS', 'PercentPhaseFOV': 100, 'PhaseEncodingDirection': 'j-', 'PhaseEncodingSteps': 96, 'PhaseResolution': 1, 'PixelBandwidth': 1395, 'ProcedureStepDescription': 'INVESTIGATORS_Gabrieli', 'ProtocolName': 'TopUp_Distortion_Map_A-P-A_Diffusion', 'PulseSequenceDetails': '%CustomerSeq%_Andre_ep2d_se_dist_nav_map', 'ReceiveCoilName': '32Ch_Head', 'ReconMatrixPE': 128, 'RepetitionTime': 7.24, 'SAR': 0.212809, 'ScanOptions': 'PFP_FS', 'ScanningSequence': 'EP', 'SequenceName': 'epse2d1_128', 'SequenceVariant': 'SK_SP', 'SeriesDescription': 'TopUp_Distortion_Map_A-P-A_Diffusion', 'SeriesNumber': 17, 'ShimSetting': [6351, -28080, -24530, 440, -193, 187, 88, 113], 'SliceThickness': 2, 'SliceTiming': [3.635, 0, 3.7325, 0.0975, 3.83, 0.1975, 3.93, 0.295, 4.0275, 0.3925, 4.125, 0.49, 4.2225, 0.59, 4.3225, 0.6875, 4.42, 0.785, 4.5175, 0.885, 4.6175, 0.9825, 4.715, 1.08, 4.8125, 1.1775, 4.91, 1.2775, 5.01, 1.375, 5.1075, 1.4725, 5.205, 1.5725, 5.305, 1.67, 5.4025, 1.7675, 5.5, 1.865, 5.5975, 1.965, 5.6975, 2.0625, 5.795, 2.16, 5.8925, 2.26, 5.9925, 2.3575, 6.09, 2.455, 6.1875, 2.555, 6.285, 2.6525, 6.385, 2.75, 6.4825, 2.8475, 6.58, 2.9475, 6.68, 3.045, 6.7775, 3.1425, 6.875, 3.2425, 6.975, 3.34, 7.0725, 3.4375, 7.17, 3.535], 'SoftwareVersions': 'syngo_MR_B17', 'SpacingBetweenSlices': 2, 'StationName': 'MRC35469', 'TaskName': 'topupdwi', 'TotalReadoutTime': 0.0508007, 'TxRefAmp': 346.603}}

Traceback (most recent call last):
  File "/opt/conda/envs/qsiprep/lib/python3.10/site-packages/nipype/pipeline/plugins/multiproc.py", line 66, in run_node
    result["result"] = node.run(updatehash=updatehash)
  File "/opt/conda/envs/qsiprep/lib/python3.10/site-packages/nipype/pipeline/engine/nodes.py", line 525, in run
    result = self._run_interface(execute=True)
  File "/opt/conda/envs/qsiprep/lib/python3.10/site-packages/nipype/pipeline/engine/nodes.py", line 643, in _run_interface
    return self._run_command(execute)
  File "/opt/conda/envs/qsiprep/lib/python3.10/site-packages/nipype/pipeline/engine/nodes.py", line 769, in _run_command
    raise NodeExecutionError(msg)
nipype.pipeline.engine.nodes.NodeExecutionError: Exception raised while executing Node merge_dwis.

Traceback:
	Traceback (most recent call last):
	  File "/opt/conda/envs/qsiprep/lib/python3.10/site-packages/nipype/interfaces/base/core.py", line 401, in run
	    runtime = self._run_interface(runtime)
	  File "/opt/conda/envs/qsiprep/lib/python3.10/site-packages/qsiprep/interfaces/dwi_merge.py", line 76, in _run_interface
	    to_concat, b0_means, corrections = harmonize_b0s(
	  File "/opt/conda/envs/qsiprep/lib/python3.10/site-packages/qsiprep/interfaces/dwi_merge.py", line 700, in harmonize_b0s
	    b0_mean = index_img(dwi_nii, b0_indices).get_fdata().mean()
	  File "/opt/conda/envs/qsiprep/lib/python3.10/site-packages/nilearn/image/image.py", line 669, in index_img
	    return _index_img(imgs, index)
	  File "/opt/conda/envs/qsiprep/lib/python3.10/site-packages/nilearn/_utils/niimg_conversions.py", line 77, in _index_img
	    img, _get_data(img)[:, :, :, index], img.affine, copy_header=True
	IndexError: index 4 is out of bounds for axis 3 with size 4