Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

FLIRT 6.0 (FSL 6.0.2) linear tests fail #52

Closed
smoia opened this issue Nov 15, 2019 · 4 comments
Closed

FLIRT 6.0 (FSL 6.0.2) linear tests fail #52

smoia opened this issue Nov 15, 2019 · 4 comments
Milestone

Comments

@smoia
Copy link
Contributor

smoia commented Nov 15, 2019

In test_linear.py, test_apply_linear_transform fails for all the parametrisations sw_tool='fsl', with a fresh install of FSL 6.0.2 (FLIRT 6.0).
It's an assertion error due to the difference being higher than the tolerance.
This doesn't happen in the travis CI tests, so I assume it is due to the different version of FSL.

@oesteban
Copy link
Collaborator

oesteban commented Feb 14, 2022

This could be related, although we are using 3D images - Washington-University/HCPpipelines#118

I'm investigating this right now (cc @effigies, for any quick insights I might want to check).

EDIT: BTW, it would be nice to maybe consider some of these new behaviors (bugs?) in FSL 6 within nipype (pinging @satra to learn if he knows more about these issues, his opinion about how nipype should evolve considering these issues, and any other enlightening comment he may have).

@oesteban
Copy link
Collaborator

Okay, I'm chatting with Matthew Webster and Fidel Alfaro-Almagro about this issue and will report back ASAP. It seems there is a bug in FSL 6 resampling. If feels like the option -setbackground 0 has stopped being honored or maybe has changed its behavior. The outcome artifact is that there's a rotation frame around the resampled images that is filled with zeros for nitransforms and FSL 5, but it's a constant nonzero number for FSL 6.

For the moment, and to make all these tests less brittle, I'm going to look for intensity differences only within the brain mask.

@satra
Copy link
Member

satra commented Feb 14, 2022

there is another issue we came across recently using the applyisoxfm in flirt, which gave us weird issues that we could only uncover in the tsnr maps. they affected the std dev maps and looked like a checkerboard pattern. i replaced the script with a nitransforms version that did not have the effect. i was going to write this up and send to the FSL folks, but haven't had the time. this wasn't a direct case of before and after, but it's possible that something fundamental changed in the resampling behavior in flirt.

@oesteban
Copy link
Collaborator

The plot thickens - when I looked at the differences within the brain mask, there are still differences between FSL 5 and 6. These were masked by the out-of-bounds constant value I first checked.

oesteban added a commit that referenced this issue Feb 15, 2022
Makes the tests less brittle so that they don't fail with FSL 6.0.

Closes #52.
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

3 participants