Skip to content

Conversation

@jan-petr
Copy link
Contributor

@jan-petr jan-petr commented Dec 9, 2021

Linked issue

#978

@jan-petr jan-petr self-assigned this Dec 9, 2021
@jan-petr jan-petr linked an issue Dec 9, 2021 that may be closed by this pull request
9 tasks
@jan-petr jan-petr changed the title Bug #978 fix multi pl dquant Bug #978 fix multi PLD quant Dec 9, 2021
Copy link
Contributor

@MichaelStritt MichaelStritt left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

You nicely cleaned up the code. Especially the switch statements and the io of the functions seems a lot better now. All headers look fine. I added two minor comments. Would be good if someone with an FSL/BASIL installation (@BeatrizPadrela or @MDijsselhof) could test this though (either now or during the release testing).

@MichaelStritt MichaelStritt added the bug Something isn't working label Dec 10, 2021
@jan-petr
Copy link
Contributor Author

jan-petr commented Dec 10, 2021

We have several flavors to test this with. I've tried them, but a re-test by others would be welcome:
These work nicely and give CBF and TT maps

  • GE_PCASL_3Dspiral_25.0LX_eASLPLD3_1
  • GE_PCASL_3Dspiral_25.0LX_eASLPLD7_1

This one runs and outputs everything. But for a very weird reason, ATT map is messed up for upper slices - needs a bit more investigation

  • Philips_PCASL_2DEPI_5.6.1.2_Hadamard_1

These two crashes in FABBER as multi-TE is not yet supported - we will do that later

  • Siemens_PCASL_3DGRASE_VD13A_Hadamard8FME_1
  • Siemens_PCASL_3DGRASE_VD13A_Hadamard4FME_1

This one crashes badly, but the reason is that it doesn't have M0 and has BSup. And our BSup correction for M0 doesn't yet work for multi-PLD. I already created a future issue for that in 2.1.0

  • Siemens_PCASL_3DGRASE_VB17A_multiPLD_1

This one used to crash because there are way too many volumes for motion correction. However, I have disabled motion correction for Look-Locker, so it should now run. And we will adapt the motion correction for Look-Locker in the future - an issue for that has been already created:

  • Philips_PASL_2DEPI_5.4.0.2-LookLocker_1

@MichaelStritt Can you install FSL in WSL?
@BeatrizPadrela Can you test these yourself as well?
@MDijsselhof Can you test your testing set if it works?
@everybody: Any other ideas or testing data?

Copy link
Member

@HenkMutsaerts HenkMutsaerts left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

In general, this issue contains multiple issues, making it harder to review this.

xASL_wrp_Resample already saved both PWI and PWI4D, and it should still do this, irregardless of the number of PLDs. At the end of the quantification wrapper we delete PWI4D as it can take up quite some space. And for BASIL we always provide PWI4D, regardless the number of PLDs. (at least, this is how it should be and how I programmed this)

@MichaelStritt MichaelStritt changed the title Bug #978 fix multi PLD quant Closes #978 Multi PLD quantification Dec 13, 2021
@ExploreASL ExploreASL deleted a comment from HenkMutsaerts Dec 13, 2021
@MichaelStritt
Copy link
Contributor

MichaelStritt commented Dec 14, 2021

Copy link
Member

@HenkMutsaerts HenkMutsaerts left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Few minor fixes left

Copy link
Contributor

@MDijsselhof MDijsselhof left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Couple of questions!


xASL_spm_deformations(x, {x.P.Path_CBF}, {x.P.Pop_Path_CBF}, [], [], AffineTransfPath, x.P.Path_y_ASL);
if xASL_exist(x.P.Path_TT,'file')
xASL_spm_deformations(x, {x.P.Path_TT}, {x.P.Pop_Path_TT}, [], [], AffineTransfPath, x.P.Path_y_ASL);
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I apologise for opening this again, but what is the difference? Both produce ATT maps?

@MichaelStritt
Copy link
Contributor

@jan-petr, @BeatrizPadrela, @MDijsselhof & @HenkMutsaerts: Let's merge as is and test in release branch. Easier and better IMHO 👍

Copy link
Member

@HenkMutsaerts HenkMutsaerts left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I added 2 commits, for clarity and to avoid running through xASL_wrp_Realign when bMotionCorrection

@jan-petr jan-petr force-pushed the bug-#978_FixMultiPLDquant branch from 28eb183 to 3314efe Compare December 17, 2021 09:15
@jan-petr jan-petr merged commit 3314efe into develop Dec 17, 2021
@jan-petr jan-petr deleted the bug-#978_FixMultiPLDquant branch December 17, 2021 09:16
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

bug Something isn't working

Projects

None yet

Development

Successfully merging this pull request may close these issues.

Fix multiPLD quantification

6 participants