Hello HCP users,
We're doing our own preprocessing of a variety of datasets and are planning
on running HCP data through our pipeline. I'm trying to take one of the raw
images from HCP and do gradient distortion correction (using the HCP
gradunwarp package at https://github.com/Washington-University/gradunwarp).
I'm using the command-line tool, specifying a gradient coefficient file,
but am running into a problem.
Here's the error:
Traceback (most recent call last):
File "/usr/local/bin/gradient_unwarp.py", line 124, in <module>
grad_unwarp.run()
File "/usr/local/bin/gradient_unwarp.py", line 113, in run
self.unwarper.run()
File
"/usr/local/lib/python2.7/dist-packages/gradunwarp/core/unwarp_resample.py",
line 152, in run
m_rcs2lai, m_rcs2lai_nohalf, g_xyz2rcs)
File
"/usr/local/lib/python2.7/dist-packages/gradunwarp/core/unwarp_resample.py",
line 271, in non_linear_unwarp_siemens
order=self.order)
File "/usr/lib/python2.7/dist-packages/scipy/ndimage/interpolation.py",
line 299, in map_coordinates
raise RuntimeError('invalid shape for coordinate array')
RuntimeError: invalid shape for coordinate array
Digging into the code in some of those files, the problem seems to be
occurring with the 4-dimensional volume that is being passed to the NumPy
interpolation.
I'm running Ubuntu 14.04 with Python version 2.7.6, NumPy version 1.8.2 and
Scipy version 0.13.3.
Any help would be greatly appreciated!
Michael
_______________________________________________
HCP-Users mailing list
[email protected]
http://lists.humanconnectome.org/mailman/listinfo/hcp-users