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__()