Mercurial > repos > rv43 > tomo
view tomo_reduce.py @ 73:adb84c450789 draft
planemo upload for repository https://github.com/rolfverberg/galaxytools commit 366e516aef0735af2998c6ff3af037181c8d5213-dirty
author | rv43 |
---|---|
date | Mon, 20 Mar 2023 17:30:24 +0000 |
parents | 1cf15b61cd83 |
children | 4f4ee8db5f67 |
line wrap: on
line source
#!/usr/bin/env python3 import logging import argparse import pathlib import sys #import tracemalloc from workflow.run_tomo import Tomo #from memory_profiler import profile #@profile def __main__(): # Parse command line arguments parser = argparse.ArgumentParser( description='Reduce tomography data') # parser.add_argument('-i', '--input_file', # required=True, # type=pathlib.Path, # help='''Full or relative path to the input file (in yaml or nxs format).''') # parser.add_argument('-o', '--output_file', # required=True, # type=pathlib.Path, # help='''Full or relative path to the output file (in Nexus format).''') # parser.add_argument('--galaxy_flag', # action='store_true', # help='''Use this flag to run the scripts as a galaxy tool.''') # parser.add_argument('--img_x_bounds', # required=False, # nargs=2, # type=int, # help='Vertical data reduction image range') parser.add_argument('-l', '--log', # type=argparse.FileType('w'), default=sys.stdout, help='Logging stream or filename') # parser.add_argument('--log_level', # choices=logging._nameToLevel.keys(), # default='INFO', # help='''Specify a preferred logging level.''') args = parser.parse_args() # Set log configuration # When logging to file, the stdout log level defaults to WARNING logging_format = '%(asctime)s : %(levelname)s - %(module)s : %(funcName)s - %(message)s' level = logging.getLevelName(args.log_level) if args.log is sys.stdout: logging.basicConfig(format=logging_format, level=level, force=True, handlers=[logging.StreamHandler()]) else: if isinstance(args.log, str): logging.basicConfig(filename=f'{args.log}', filemode='w', format=logging_format, level=level, force=True) elif isinstance(args.log, io.TextIOWrapper): logging.basicConfig(filemode='w', format=logging_format, level=level, stream=args.log, force=True) else: raise(ValueError(f'Invalid argument --log: {args.log}')) stream_handler = logging.StreamHandler() logging.getLogger().addHandler(stream_handler) stream_handler.setLevel(logging.WARNING) stream_handler.setFormatter(logging.Formatter(logging_format)) # Start memory monitoring # tracemalloc.start() # Log command line arguments # logging.info(f'input_file = {args.input_file}') # logging.info(f'output_file = {args.output_file}') # logging.info(f'galaxy_flag = {args.galaxy_flag}') # logging.info(f'img_x_bounds = {args.img_x_bounds}') # logging.debug(f'log = {args.log}') # logging.debug(f'is log stdout? {args.log is sys.stdout}') # logging.debug(f'log_level = {args.log_level}') logging.info(f'log = {args.log}') logging.info(f'is log stdout? {args.log is sys.stdout}') print('Done') # Instantiate Tomo object # tomo = Tomo(galaxy_flag=args.galaxy_flag) # Read input file # data = tomo.read(args.input_file) # Generate reduced tomography images # data = tomo.gen_reduced_data(data, img_x_bounds=args.img_x_bounds) # Write output file # data = tomo.write(data, args.output_file) # Displaying memory usage # logging.info(f'Memory usage: {tracemalloc.get_traced_memory()}') # Stop memory monitoring # tracemalloc.stop() if __name__ == "__main__": __main__()