Osprey Error - Error using io_loadspec_twix

[daghbashian]

I was wondering if someone can please assist me with some processing errors I’m running into in Osprey.

I am unable to complete the load data function within Osprey and am getting the following errors:

Error using io_loadspec_twix
Too many output arguments.

Error in osp_LoadTwix (line 66)
[raw, raw_ref] = io_loadspec_twix(MRSCont.files{metab_ll,kk});

Error in OspreyLoad (line 197)
[MRSCont] = osp_LoadTwix(MRSCont);

Error in osp_onLoad (line 35)
MRSCont = OspreyLoad(MRSCont);

Error while evaluating UIControl Callback.

[Helge]

Hi @daghbashian,

Could you try downloading the most recent develop branch. I have pushed a fix for this last week (Update io_loadspec_twix · schorschinho/osprey@7d8b6f5 · GitHub).

Let me know if this solves the problem.

Best,
Helge

[khwood]

Good afternoon!

I am having this same error when trying to load data. I downloaded Osprey 2 days ago. Do I need to download this most recent develop branch?

Your guidance and input is greatly appreciated!
–Kimberly

[khwood]

Hi @Helge,

I checked my version of the io_loadspec_twix and it appears to be the most recent version (as seen on GitHub).

I’m still getting this Load data error… any thoughts on where else I need to look to resolve the issue?

Much thanks!!
–Kimberly

[admin]

What sequence is this exactly?

[khwood]

Sequence name: svs_se

(I’m new to MRS, but have the scanner protocol, what do you need to know exactly?)

[admin]

The protocol would be great, as well as the Siemens scanner software version, and ideally your Osprey job file!

[khwood]

Unfortunately, unless it’s in the scanner protocol, I do not know the scanner software version, but will see what I can find out! Unless, “XA30” is the software version…?

Scanner protocol and Osprey job file attached.

Much thanks!
–Kimberly
spss2_mrs_twix_1493_1_rtVS_khw.m (13.7 KB)
SPSSA XA30.pdf (155.1 KB)

[khwood]

Hi @admin,

The current operating system on the 3T Prisma is Syngo XA30.

Here’s a little extra info from our technologist: https://www.siemens-healthineers.com/en-us/magnetic-resonance-imaging/options-and-upgrades/upgrades/syngo-mr-xa30

[admin]

Ah, great - that is a use case we have not yet encountered, an XA30 sequence with an internal water reference scan (Sequence → Common → Ref. Scan Mode = Inline Correction)! This was not possible on older Siemens software versions. It means that the sequence itself acquires a water reference scan and stores it alongside the metabolite data (and you do not need to acquire a separate water reference scan).

I am not sure whether you are able to share these data institutionally (we have deidentification scripts in the Osprey folder), but they would be very helpful for us (and we need to modify the loader to fix the problem). If you can’t, phantom data acquired with the same protocol will do.

[khwood]

@admin,

So this means I do not need to load the non-water suppressed file (at all) in the analysis, and only use the water suppressed data?

The ACC (anterior cingulate cortex) data were acquired with a slightly different protocol (I’m working to get a copy from the scanner).

The rtVS (right ventral striatum) data were acquired with the protocol sent earlier.

The data can be found here and have been deidentified at acquisition.

Thanks so much for your prompt reply and guidance!
–Kimberly

[admin]

Yes, it’d mean you don’t need to use anything in files_ref.

Can’t access the data though:

[khwood]

Trying again a different way… found here

Alternatively, I could share with your email? (haven’t been able to find it online)

[khwood]

If the above doesn’t work, please try this Google drive

Thank you greatly!!
–Kimberly

[khwood]

@admin

Good afternoon!

Any luck with downloading the data?

Looking forward to hearing from you.

Thank you greatly!!
–Kimberly

[admin]

Luck - yes. Time - no It’s on my list.

[admin]

I just tried to reproduce the error you reported but the data (from the right ventral striatum) load and process totally fine. I’m noticing that in your job file you’re pointing to the exact same file (meas_MID00084_FID46817_RightVS_svs_se_30_TXRX_TR200.dat) for both files and files_ref, so I assume you made a simple copy/paste error. I used the water-unsuppressed file meas_MID00085_FID46818_RightVS_svs_se_30_TXRX_NWS_TR2000.dat in files_ref and everything is peachy.

Note the following:

• I have not investigated the inline water correction mode right now but the dataset will be useful to explore that once I have more time.
• I have not looked at the ACC data that you have acquired with a different protocol (please share that protocol).

I also have some feedback on your protocol in general:

• (IMPORTANT) Whenever you acquire a separate water reference, you need to set Sequence → Common → Delta Frequency to 0.0 ppm so that the water reference signal is co-localized to the water-suppressed acquisition (this counters the chemical shift displacement artefact).
• (SOMEWHAT LESS IMPORTANT) 12 transients for the water-unsuppressed scan is plenty. You can easily get away with 4 because you get so much signal.

Since these are protocol issues that I give feedback on very frequently, I should probably write a sticky post about them… maybe this weekend if I find time.

[khwood]

Thank you greatly for this feedback and guidance! Thank you for the feedback on our MRS protocol!

Could you share your job file? I have changed the file name for the “files_ref” input but continue to get an error when I try to load the data in the GUI.

[admin]

Relevant lines of my job file (where data_folder is the directory I’ve stored both files in):

files      = {fullfile(data_folder, 'meas_MID00084_FID46817_RightVS_svs_se_30_TXRX_TR200.dat')};
files_ref  = {fullfile(data_folder, 'meas_MID00085_FID46818_RightVS_svs_se_30_TXRX_NWS_TR2000.dat')};

What are yours?

[admin]

Just to check, what MATLAB version are you using?