1 | #!/usr/bin/env python |
---|
2 | |
---|
3 | from __future__ import division |
---|
4 | |
---|
5 | import sys |
---|
6 | from galaxy import eggs |
---|
7 | import pkg_resources; pkg_resources.require( "bx-python" ) |
---|
8 | from bx.arrays.array_tree import * |
---|
9 | from bx.arrays.wiggle import WiggleReader |
---|
10 | |
---|
11 | BLOCK_SIZE = 100 |
---|
12 | |
---|
13 | def main(): |
---|
14 | |
---|
15 | input_fname = sys.argv[1] |
---|
16 | out_fname = sys.argv[2] |
---|
17 | |
---|
18 | reader = WiggleReader( open( input_fname ) ) |
---|
19 | |
---|
20 | # Fill array from reader |
---|
21 | d = array_tree_dict_from_reader( reader, {}, block_size = BLOCK_SIZE ) |
---|
22 | |
---|
23 | for array_tree in d.itervalues(): |
---|
24 | array_tree.root.build_summary() |
---|
25 | |
---|
26 | FileArrayTreeDict.dict_to_file( d, open( out_fname, "w" ) ) |
---|
27 | |
---|
28 | if __name__ == "__main__": |
---|
29 | main() |
---|