We are having trouble getting the latest version (2.5.0) of osprey to process Siemens MEGA-PRESS data. When running OspreyLoad, we get
| Reader version: 1496740353 (UTC: 06-Jun-2017 09:12:33)
| Scan 1/2, read all mdhs:
| 35.4 MB read in 1 s
| Scan 2/2, read all mdhs:
| 532.2 MB read in 1 s
|
| No water reference data found. Returning empty field
and then when running OspreyProcess
| Unrecognized field name “names”.
|
| Error in osp_XReferencing (line 42)
| ind = find(strcmp(dataToFit.names,‘sum’));
|
| Error in op_preref (line 178)
| [refShift, ~] =
| osp_XReferencing(temp_spec,frequencies,polarity,lim,0);%
| determine frequency shift
|
| Error in OspreyProcess (line 248)
| [refShift_ind_ini]=op_preref(raw,seq);
|
These errors occur even with data that was successfully processed with earlier versions of Osprey. Any help with this would be much appreciated.
So the first message is normal. The no water reference data found. Returning empty field message means that your sequence did not have an integrated water reference in the TWIX .dat file. This is a feature of, for example, the CMRR sequence. In this case, the TWIX .dat file contains water reference data which is then separated by Osprey.
For the second message. Could you make a screenshot of MRSCont.raw{1, 1}.flags and MRSCont.raw{1, 1}.seq. It seems like the sequence is not correctly identified.
I coincidentally was using Osprey for the first time and ran into the exact same error messages as Ben using Siemens MEGA-PRESS data. I will send you a direct message with a DropBox link to download an example data file.
Error in osp_plotAllPDF (line 76)
osp_plotModule(MRSCont, ‘OspreyFit’, kk,[bb ss], Names{mm});
Error in OspreyFit (line 228)
osp_plotAllPDF(MRSCont, ‘OspreyFit’);
which was the result of water fit results being stored as MRSCont.fit.results.w.fitParams. rather than MRSCont.fit.results.water.fitParams. Changing line 606 in osp_plotModule.m to:
The corrections you made are correct ( fixed them on GitHub now). One comment about this, I would recommend running a unsuppressed water reference with the same settings as the metabolite data (TE etc) to perform the eddy current correction.
The second error is a bit tricky because the test data you shared did not have a T1 scan. Could you open the GUI after you run OsrpeyCoreg and post the output? Also, could you check the MRSCont flags didCoreg, too?