MRIQC: number of observations cannot be determined on an empty distance matrix

Using MRIQC docker version: nipreps/mriqc:22.0.6. Encountered several edge cases where subjects have the following value error (error in 11 out of 15500 ran) for *.bold.nii.gz images. Looked at images and was not able to find anything that immediately stood out. Curious if anyone has encountered something similar?

ValueError: The number of observations cannot be determined on an empty distance matrix.

The full error is below:

nipype.pipeline.engine.nodes.NodeExecutionError: Exception raised while executing Node BigPlot.

Traceback (most recent call last):
File “/opt/conda/lib/python3.9/site-packages/nipype/interfaces/base/core.py”, line 398, in run
runtime = self._run_interface(runtime)
File “/opt/conda/lib/python3.9/site-packages/niworkflows/interfaces/plotting.py”, line 99, in _run_interface
fig = fMRIPlot(
File “/opt/conda/lib/python3.9/site-packages/niworkflows/viz/plots.py”, line 129, in plot
plot_carpet(
File “/opt/conda/lib/python3.9/site-packages/niworkflows/viz/plots.py”, line 236, in plot_carpet
children, _, n_leaves, _, distances = ward_tree(roi_data, return_distance=True)
File “/opt/conda/lib/python3.9/site-packages/sklearn/cluster/_agglomerative.py”, line 270, in ward_tree
out = hierarchy.ward(X)
File “/opt/conda/lib/python3.9/site-packages/scipy/cluster/hierarchy.py”, line 834, in ward
return linkage(y, method=‘ward’, metric=‘euclidean’)
File “/opt/conda/lib/python3.9/site-packages/scipy/cluster/hierarchy.py”, line 1068, in linkage
n = int(distance.num_obs_y(y))
File “/opt/conda/lib/python3.9/site-packages/scipy/spatial/distance.py”, line 2573, in num_obs_y
raise ValueError("The number of observations cannot be determined on "
ValueError: The number of observations cannot be determined on an empty distance matrix.

We’ve run across this here: ds002684-mriqc: "ValueError: The number of observations cannot be determined on an empty distance matrix." · Issue #39 · poldracklab/tacc-openneuro · GitHub. Should be resolved in FIX: Check for empty segments in plot_carpet by effigies · Pull Request #783 · nipreps/niworkflows · GitHub, so the next MRIQC release that includes niworkflows 1.7.4 or higher should be fixed.

1 Like