% Example batch script for Philips MEGA-PRESS data with accompanying water
% references and structural images Remember to check GannetPreInitialise.m
% has the appropriate settings for your data
data_dir = '/Users/username/Documents/my_project/data';
metab = {'S01_GABA_68_act.sdat'
'S02_GABA_68_act.sdat'
'S03_GABA_68_act.sdat'};
metab = fullfile(data_dir, metab); % It is good practice to include the full path of files
water = {'S01_GABA_68_ref.sdat'
'S02_GABA_68_ref.sdat'
'S03_GABA_68_ref.sdat'};
water = fullfile(data_dir, water);
anat = {'S01_T1_anat.nii'
'S02_T1_anat.nii'
'S03_T1_anat.nii'};
anat = fullfile(data_dir, anat);
MRS = GannetLoad(metab, water);
MRS = GannetFit(MRS);
MRS = GannetCoRegister(MRS, anat);
MRS = GannetSegment(MRS);
MRS = GannetQuantify(MRS);
LS0tCnRpdGxlOiAiQmF0Y2ggcHJvY2Vzc2luZyBzY3JpcHQiCmRhdGU6ICJMYXN0IHVwZGF0ZWQ6IGByIGZvcm1hdChTeXMudGltZSgpLCAnJUIgJWQsICVZJylgIgpvdXRwdXQ6IGh0bWxfZG9jdW1lbnQKLS0tCgpgYGB7ciBzZXR1cCwgaW5jbHVkZSA9IEZBTFNFfQprbml0cjo6b3B0c19jaHVuayRzZXQoZWNobyA9IFRSVUUpCmBgYAoKYGBge3IsIGNoaWxkID0gYygianMvYmFjay10by10b3AuaHRtbCIsICJqcy9jb3B5LXRvLWNsaXBib2FyZC5odG1sIil9CmBgYAoKYGBge2NzcywgZWNobyA9IEZBTFNFfQpib2R5IC5tYWluLWNvbnRhaW5lciB7CiAgbWF4LXdpZHRoOiA5MDBweDsKfQpgYGAKCjxicj4KCmBgYHtvY3RhdmUsIGV2YWwgPSBGQUxTRX0KJSBFeGFtcGxlIGJhdGNoIHNjcmlwdCBmb3IgUGhpbGlwcyBNRUdBLVBSRVNTIGRhdGEgd2l0aCBhY2NvbXBhbnlpbmcgd2F0ZXIKJSByZWZlcmVuY2VzIGFuZCBzdHJ1Y3R1cmFsIGltYWdlcyBSZW1lbWJlciB0byBjaGVjayBHYW5uZXRQcmVJbml0aWFsaXNlLm0KJSBoYXMgdGhlIGFwcHJvcHJpYXRlIHNldHRpbmdzIGZvciB5b3VyIGRhdGEKCmRhdGFfZGlyID0gJy9Vc2Vycy91c2VybmFtZS9Eb2N1bWVudHMvbXlfcHJvamVjdC9kYXRhJzsKCm1ldGFiID0geydTMDFfR0FCQV82OF9hY3Quc2RhdCcKICAgICAgICAgJ1MwMl9HQUJBXzY4X2FjdC5zZGF0JwogICAgICAgICAnUzAzX0dBQkFfNjhfYWN0LnNkYXQnfTsKbWV0YWIgPSBmdWxsZmlsZShkYXRhX2RpciwgbWV0YWIpOyAlIEl0IGlzIGdvb2QgcHJhY3RpY2UgdG8gaW5jbHVkZSB0aGUgZnVsbCBwYXRoIG9mIGZpbGVzCiAgICAgICAgIAp3YXRlciA9IHsnUzAxX0dBQkFfNjhfcmVmLnNkYXQnCiAgICAgICAgICdTMDJfR0FCQV82OF9yZWYuc2RhdCcKICAgICAgICAgJ1MwM19HQUJBXzY4X3JlZi5zZGF0J307CndhdGVyID0gZnVsbGZpbGUoZGF0YV9kaXIsIHdhdGVyKTsKCmFuYXQgPSB7J1MwMV9UMV9hbmF0Lm5paScKICAgICAgICAnUzAyX1QxX2FuYXQubmlpJwogICAgICAgICdTMDNfVDFfYW5hdC5uaWknfTsKYW5hdCA9IGZ1bGxmaWxlKGRhdGFfZGlyLCBhbmF0KTsKCk1SUyA9IEdhbm5ldExvYWQobWV0YWIsIHdhdGVyKTsKTVJTID0gR2FubmV0Rml0KE1SUyk7Ck1SUyA9IEdhbm5ldENvUmVnaXN0ZXIoTVJTLCBhbmF0KTsKTVJTID0gR2FubmV0U2VnbWVudChNUlMpOwpNUlMgPSBHYW5uZXRRdWFudGlmeShNUlMpOwpgYGAK
Built with
R Markdown in
RStudio
Copyright © 2020–2026, Mark Mikkelsen