Source code for galaxy.datatypes.converters.wiggle_to_array_tree_converter

#!/usr/bin/env python

from __future__ import division

import sys
from galaxy import eggs
import pkg_resources; pkg_resources.require( "bx-python" )
from bx.arrays.array_tree import *
from bx.arrays.wiggle import WiggleReader

BLOCK_SIZE = 100

[docs]def main():

    input_fname = sys.argv[1]
    out_fname = sys.argv[2]

    reader = WiggleReader( open( input_fname ) )

    # Fill array from reader
    d = array_tree_dict_from_reader( reader, {}, block_size = BLOCK_SIZE )

    for array_tree in d.itervalues():
        array_tree.root.build_summary()

    FileArrayTreeDict.dict_to_file( d, open( out_fname, "w" ) )

if __name__ == "__main__":
    main()