Back to the main page.
Bug 2906 - ft_volumereslice fails for mris not in millimeters
| Status | CLOSED FIXED |
| Reported | 2015-06-04 16:18:00 +0200 |
| Modified | 2019-08-10 12:33:13 +0200 |
| Product: | FieldTrip |
| Component: | core |
| Version: | unspecified |
| Hardware: | All |
| Operating System: | All |
| Importance: | P5 major |
| Assigned to: | Tzvetan Popov |
| URL: | |
| Tags: | |
| Depends on: | |
| Blocks: | |
| See also: |
Thomas Hartmann - 2015-06-04 16:18:02 +0200
as described in the title, ft_volumereslice only handles mris in mm correctly: %% clear... clear all global close all %% load mri.... load standard_mri %% do unit conversion to mm.... mri_m = ft_convert_units(mri, 'm'); mri_mm = ft_convert_units(mri_m, 'mm'); mri_cm = ft_convert_units(mri_m, 'cm'); %% try reslice.... mri_resliced = ft_volumereslice([], mri_m); % fails mri_resliced = ft_volumereslice([], mri_cm); % fails without error message producing garbage output mri_resliced = ft_volumereslice([], mri_mm); % works as expected
Tzvetan Popov - 2015-06-23 17:41:11 +0200
Hi Thomas, I've been able to reproduce it and created a test script. I guess the problem arises in ft_sourceinterpolate during the mldivide converting the anatomical voxel positions into voxel indices into the functional volume. I'll check cheerzen tzvetan
Tzvetan Popov - 2015-06-24 10:43:35 +0200
Sending ft_volumereslice.m Transmitting file data. Committed revision 10477
Tzvetan Popov - 2016-02-24 09:22:33 +0100
Sending ft_volumereslice.m Transmitting file data. Committed revision 10477