annotate data_manager/vsnp_genbank_fetcher.py @ 0:36d828c2685b draft default tip

Uploaded
author greg
date Mon, 10 Feb 2020 14:35:22 +0000
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
36d828c2685b Uploaded
greg
parents:
diff changeset
1 import argparse
36d828c2685b Uploaded
greg
parents:
diff changeset
2 import json
36d828c2685b Uploaded
greg
parents:
diff changeset
3 import os
36d828c2685b Uploaded
greg
parents:
diff changeset
4 import sys
36d828c2685b Uploaded
greg
parents:
diff changeset
5 try:
36d828c2685b Uploaded
greg
parents:
diff changeset
6 # For Python 3.0 and later
36d828c2685b Uploaded
greg
parents:
diff changeset
7 from urllib.request import Request, urlopen
36d828c2685b Uploaded
greg
parents:
diff changeset
8 except ImportError:
36d828c2685b Uploaded
greg
parents:
diff changeset
9 # Fall back to Python 2 imports
36d828c2685b Uploaded
greg
parents:
diff changeset
10 from urllib2 import Request, urlopen
36d828c2685b Uploaded
greg
parents:
diff changeset
11
36d828c2685b Uploaded
greg
parents:
diff changeset
12
36d828c2685b Uploaded
greg
parents:
diff changeset
13 def url_download(url, workdir):
36d828c2685b Uploaded
greg
parents:
diff changeset
14 file_path = os.path.abspath(os.path.join(workdir, os.path.basename(url)))
36d828c2685b Uploaded
greg
parents:
diff changeset
15 src = None
36d828c2685b Uploaded
greg
parents:
diff changeset
16 dst = None
36d828c2685b Uploaded
greg
parents:
diff changeset
17 try:
36d828c2685b Uploaded
greg
parents:
diff changeset
18 req = Request(url)
36d828c2685b Uploaded
greg
parents:
diff changeset
19 src = urlopen(req)
36d828c2685b Uploaded
greg
parents:
diff changeset
20 with open(file_path, 'w') as dst:
36d828c2685b Uploaded
greg
parents:
diff changeset
21 while True:
36d828c2685b Uploaded
greg
parents:
diff changeset
22 chunk = src.read(2**10)
36d828c2685b Uploaded
greg
parents:
diff changeset
23 if chunk:
36d828c2685b Uploaded
greg
parents:
diff changeset
24 dst.write(chunk)
36d828c2685b Uploaded
greg
parents:
diff changeset
25 else:
36d828c2685b Uploaded
greg
parents:
diff changeset
26 break
36d828c2685b Uploaded
greg
parents:
diff changeset
27 except Exception as e:
36d828c2685b Uploaded
greg
parents:
diff changeset
28 sys.exit(str(e))
36d828c2685b Uploaded
greg
parents:
diff changeset
29 finally:
36d828c2685b Uploaded
greg
parents:
diff changeset
30 if src:
36d828c2685b Uploaded
greg
parents:
diff changeset
31 src.close()
36d828c2685b Uploaded
greg
parents:
diff changeset
32 return file_path
36d828c2685b Uploaded
greg
parents:
diff changeset
33
36d828c2685b Uploaded
greg
parents:
diff changeset
34
36d828c2685b Uploaded
greg
parents:
diff changeset
35 def download(dbkey, name, url, out_file):
36d828c2685b Uploaded
greg
parents:
diff changeset
36
36d828c2685b Uploaded
greg
parents:
diff changeset
37 with open(out_file) as fh:
36d828c2685b Uploaded
greg
parents:
diff changeset
38 params = json.loads(fh.read())
36d828c2685b Uploaded
greg
parents:
diff changeset
39
36d828c2685b Uploaded
greg
parents:
diff changeset
40 workdir = params['output_data'][0]['extra_files_path']
36d828c2685b Uploaded
greg
parents:
diff changeset
41 os.makedirs(workdir)
36d828c2685b Uploaded
greg
parents:
diff changeset
42 file_path = url_download(url, workdir)
36d828c2685b Uploaded
greg
parents:
diff changeset
43 entry_name = os.path.basename(file_path)
36d828c2685b Uploaded
greg
parents:
diff changeset
44
36d828c2685b Uploaded
greg
parents:
diff changeset
45 data_manager_json = {"data_tables": {}}
36d828c2685b Uploaded
greg
parents:
diff changeset
46 data_manager_entry = {}
36d828c2685b Uploaded
greg
parents:
diff changeset
47 data_manager_entry['value'] = dbkey
36d828c2685b Uploaded
greg
parents:
diff changeset
48 data_manager_entry['name'] = entry_name
36d828c2685b Uploaded
greg
parents:
diff changeset
49 data_manager_entry['path'] = file_path
36d828c2685b Uploaded
greg
parents:
diff changeset
50 data_manager_entry['description'] = "Genbank file for %s" % name
36d828c2685b Uploaded
greg
parents:
diff changeset
51 data_manager_json["data_tables"]["vsnp_genbank"] = data_manager_entry
36d828c2685b Uploaded
greg
parents:
diff changeset
52
36d828c2685b Uploaded
greg
parents:
diff changeset
53 with open(out_file, 'w') as fh:
36d828c2685b Uploaded
greg
parents:
diff changeset
54 fh.write(json.dumps(data_manager_json, sort_keys=True))
36d828c2685b Uploaded
greg
parents:
diff changeset
55
36d828c2685b Uploaded
greg
parents:
diff changeset
56
36d828c2685b Uploaded
greg
parents:
diff changeset
57 parser = argparse.ArgumentParser()
36d828c2685b Uploaded
greg
parents:
diff changeset
58
36d828c2685b Uploaded
greg
parents:
diff changeset
59 parser.add_argument('--dbkey', dest='dbkey', help='Genome reference dbkey')
36d828c2685b Uploaded
greg
parents:
diff changeset
60 parser.add_argument('--name', dest='name', help='Reference display name')
36d828c2685b Uploaded
greg
parents:
diff changeset
61 parser.add_argument('--url', dest='url', help='URL to download Genbank file')
36d828c2685b Uploaded
greg
parents:
diff changeset
62 parser.add_argument('--out_file', dest='out_file', help='JSON output file')
36d828c2685b Uploaded
greg
parents:
diff changeset
63
36d828c2685b Uploaded
greg
parents:
diff changeset
64 args = parser.parse_args()
36d828c2685b Uploaded
greg
parents:
diff changeset
65
36d828c2685b Uploaded
greg
parents:
diff changeset
66 download(args.dbkey, args.name, args.url, args.out_file)