Mercurial > repos > tyty > structurefold
comparison predict/ct_to_dot.py @ 102:43ea205aa49e draft
Uploaded
author | tyty |
---|---|
date | Thu, 19 Mar 2015 17:40:57 -0400 |
parents | |
children |
comparison
equal
deleted
inserted
replaced
101:a5c8b470bc9f | 102:43ea205aa49e |
---|---|
1 #!/usr/bin/env python | |
2 # -*- coding: utf-8 -*- | |
3 | |
4 import sys | |
5 import shlex | |
6 import os | |
7 import subprocess | |
8 from read_file import * | |
9 | |
10 ct_file = sys.argv[1] | |
11 path = sys.argv[2] | |
12 id_s = sys.argv[3] | |
13 result_file = sys.argv[4] | |
14 | |
15 h = file(result_file, 'w') | |
16 os.system('grep "'+id_s+'" '+ct_file+' |wc -l > '+path+'/count.txt') | |
17 count = read_t_file(path+'/count.txt') | |
18 comm = '' | |
19 for i in range(int(count[0][0])): | |
20 command = shlex.split('ct2dot %s %s %s' % (ct_file, str(i+1), os.path.join(path, 'db_file_%s.dbnn' % str(i+1)))) | |
21 subprocess.call(command) | |
22 comm = comm +' '+path+'/db_file_'+str(i+1)+'.dbnn' | |
23 | |
24 | |
25 | |
26 os.system('cat'+comm+' > '+result_file) | |
27 for i in range(int(count[0][0])): | |
28 command = shlex.split('rm %s' % (os.path.join(path, 'db_file_%s.dbnn' % str(i+1)))) | |
29 subprocess.call(command) | |
30 command = shlex.split('rm %s' % (os.path.join(path, 'count.txt'))) | |
31 subprocess.call(command) | |
32 | |
33 | |
34 | |
35 h.close() | |
36 |