#!/usr/bin/env python from hcpxnat.interface import HcpInterface idb = HcpInterface(config="/home/NRG/mhilem01/.hcpxnat_intradb.cfg") idb.project = "NKI" sessions = idb.getSessions(project="NKI") for s in sessions: idb.subject_label = idb.session_label = s.get("label") scans = idb.getSessionScans() for s in scans: idb.scan_id = s.get("ID") sd = s.get("series_description") if "MUSIC" in sd: print s
#'MPRAGE_SIEMENS': 'T1w', #'REST_1400': 'rfMRI', #'REST_CAP': 'rfMRI', #'CHECKERBOARD_645': 'tfMRI', #'CHECKERBOARD_1400': 'tfMRI', #'BREATH_HOLD_1400': 'tfMRI', #'DIFF_137_AP': 'dMRI' #'MPRAGE_SIEMENS_DEFACED': 'T1w', #'MPRAGE_SIEMENS_MUSIC ONLY': 'T1w', #'FLAIR - ch32 iPAT2 -CLINICAL': 'T2w', 'MPRAGE': 'T1w' } for s in sessions: idb.subject_label = idb.session_label = s.get('label') scans = idb.getSessionScans() print "\nModifying series_desc for session " + idb.session_label for s in scans: idb.scan_id = s.get('ID') sd = s.get('series_description') if not sd in mapping: print "+++ " + sd + " not in series desc mapping. SKIPPING..." continue print " %s %s - Scan type %s --> %s" % \ (s.get('ID'), sd, s.get('type'), mapping[sd]) idb.setScanElement('xnat:mrScanData', 'type', mapping[sd])