For me personally dcm2niix still crashes (even though the up-to-date flavor library was used):

==============================================================================================
Importing subject = Sub1_1:   
Subject = Sub1, visit = 1, session = ASL_1, scan = ASL4D
executing: ["M:\SoftwareDevelopment\MATLAB\m.stritt\Server_xASL\ExploreASL\External\MRIcron\20190902\dcm2nii.exe" 
-f "%f_%p_%t_%s_%r" -o "M:\SoftwareDevelopment\MATLAB\m.stritt\Server_xASL\FlavorDatabase\Siemens_PCASL_3DGRASE_VD13A_Hadamard4FME_1\derivatives\ExploreASL\temp\Sub1\ASL_1\dcm2nii_temp_ASL4D"
"M:\SoftwareDevelopment\MATLAB\m.stritt\Server_xASL\FlavorDatabase\Siemens_PCASL_3DGRASE_VD13A_Hadamard4FME_1\sourcedata\Sub1\ASL\IM0007.0001.4914180.IMA"]
==============================================================================================

==============================================================================================
status: -1073740791
Warning: No NIfTIs found! Either illegal characters in filenames or other code defect... 
> In xASL_io_dcm2nii (line 251)

Okay guys, I had to put an if loop here

because with Amnah's data (FME), we only have 2 repetitions, so this line makes sense and works for the next for loop of the averaging across repetitions.
However, with Philips data, since we have 12 repetitions, the size of ASL_im, divided by the 12 repetitions, gives the PLDs per repetition (7). If I used the same NumberofPLDs line that I used for FME data, in the foor loop, in line 251, in the last iteration, we would have 12 (number of rep) * 12 (84 / Had8-1 = 84/7 = 12), which would overfit the matrix size

Tell me what you think @MichaelStritt @jan-petr

@jan-petr & @BeatrizPadrela: Possibly a stupid question from my side, but wouldn't it be easier to determine that based on the matrix size? Then your for loop can't crash.

By the way, there are no if loops, we call them if statements.

@MichaelStritt that's a good idea!
And yes indeed, if statements and for loops, yesterday was already late :)