Thanks for all the contributions!
I was running kspaceFirstOrder3DG
using a sensor mask defined by cuboid corners, and I defined the sensors' frequency response.
Error message popped up when it ran into gaussianFilter
, saying that there are "Too many input arguments".
After some inspection, I found that in kspaceFirstOrder3D
, gaussianFilter
was applied to sensor data (where sensor_data.p
is a 2-D matrix) before re-ordering, so everything was fine.
However, in kspaceFirstOrder3DG
, the sensor data read back from h5 file is already in 4-D format, such that gaussianFilter
can not work properly.
I think I can remove the sensors' frequency response for now and filter the signal after kspaceFirstOrder3DG
finishes. Or is there any other alternative modifications that I can do?
Thanks again!
P.S. I have modified the h5writeatt
part according to http://www.k-wave.org/forum/topic/error-reading-h5-files-when-using-binaries.