annotate galaxy-tools/biobank/tools/get_studies_details.py @ 0:ba6cf6ede027 draft default tip

Uploaded
author ric
date Wed, 28 Sep 2016 06:03:30 -0400
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
1 '''
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
2 '''
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
3
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
4 import argparse, csv, sys
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
5
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
6 from bl.vl.kb import KnowledgeBase as KB
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
7 import bl.vl.utils.ome_utils as vlu
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
8 from bl.vl.utils import get_logger, LOG_LEVELS
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
9
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
10
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
11 def make_parser():
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
12 parser = argparse.ArgumentParser(description='retrieve studies details')
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
13 parser.add_argument('--logfile', type=str, help='log file (default=stderr)')
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
14 parser.add_argument('--loglevel', type=str, choices=LOG_LEVELS,
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
15 help='logging level (default=INFO)', default='INFO')
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
16 parser.add_argument('-H', '--host', type=str, help='omero hostname')
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
17 parser.add_argument('-U', '--user', type=str, help='omero user')
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
18 parser.add_argument('-P', '--passwd', type=str, help='omero passwd')
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
19 parser.add_argument('--out_file', type=str, required=True,
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
20 help='output file with studies details')
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
21 return parser
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
22
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
23
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
24 def main(argv):
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
25 parser = make_parser()
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
26 args = parser.parse_args(argv)
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
27
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
28 logger = get_logger('get_studies', level=args.loglevel,
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
29 filename=args.logfile)
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
30
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
31 try:
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
32 host = args.host or vlu.ome_host()
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
33 user = args.user or vlu.ome_user()
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
34 passwd = args.passwd or vlu.ome_passwd()
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
35 except ValueError, ve:
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
36 logger.critical(ve)
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
37 sys.exit(ve)
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
38
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
39 kb = KB(driver='omero')(host, user, passwd)
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
40
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
41 logger.info('Retrieving studies')
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
42 studies = kb.get_objects(kb.Study)
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
43 logger.info('Retrieved {} studies'.format(len(studies)))
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
44 records = []
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
45 for s in studies:
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
46 enr = kb.get_enrolled(s)
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
47 rec = {'label': s.label,
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
48 'description': s.description,
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
49 'enrolled': len(enr)}
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
50 records.append(rec)
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
51
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
52 with open(args.out_file, 'w') as outfile:
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
53 writer = csv.DictWriter(outfile, ['label', 'description', 'enrolled'],
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
54 delimiter='\t')
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
55 # writer.writeheader()
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
56 writer.writerows(sorted(records, key=lambda key: key['label']))
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
57
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
58 if __name__ == '__main__':
ba6cf6ede027 Uploaded
ric
parents:
diff changeset
59 main(sys.argv[1:])