This document describes the HydrOffice BAG library and tools (0.2). For the source code, go here.

Source code for hydroffice.bag.tools.bag_elevation

from __future__ import absolute_import, division, print_function, unicode_literals

import os
import logging


[docs]def main(): logger = logging.getLogger() logger.setLevel(logging.NOTSET) import argparse from hydroffice.bag import BAGFile, is_bag, __version__ app_name = "bag_elevation" app_info = "Extraction of elevation layer from an OpenNS BAG file, using hydroffice.bag r%s" % __version__ formats = ['ascii', 'geotiff', 'xyz'] parser = argparse.ArgumentParser(prog=app_name, description=app_info) parser.add_argument("bag_file", type=str, help="a valid BAG file from which to extract metadata") parser.add_argument("-f", "--format", help="one of the available file format: " + ", ".join(formats), choices=formats, default="geotiff", metavar='') parser.add_argument("-o", "--output", help="the output file", type=str) parser.add_argument("-v", "--verbose", help="increase output verbosity", action="store_true") args = parser.parse_args() if args.verbose: print("> verbosity: ON") ch = logging.StreamHandler() ch.setLevel(logging.DEBUG) # change to WARNING to reduce verbosity, DEBUG for high verbosity ch_formatter = logging.Formatter('%(levelname)-9s %(name)s.%(funcName)s:%(lineno)d > %(message)s') ch.setFormatter(ch_formatter) logger.addHandler(ch) if args.verbose: print("> input: %s" % args.bag_file) if args.output: args.output = os.path.abspath(args.output) print("> output: %s" % args.output) else: args.output = None print("> output: [default]") print("> format: %s" % args.format) if not os.path.exists(args.bag_file): parser.exit(1, "ERROR: the input valid does not exist: %s" % args.bag_file) if not is_bag(args.bag_file): parser.exit(1, "ERROR: the input valid does not seem a BAG file: %s" % args.bag_file) bf = BAGFile(args.bag_file) bag_meta = None try: bag_meta = bf.populate_metadata() except Exception as e: parser.exit(1, "ERROR: issue in metadata population: %s" % e) bag_elevation = None try: bag_elevation = bf.elevation(mask_nan=False) except Exception as e: parser.exit(1, "ERROR: issue in elevation population: %s" % e) try: from hydroffice.bag.elevation import Elevation2Gdal Elevation2Gdal(bag_elevation=bag_elevation, bag_meta=bag_meta, fmt=args.format, out_file=args.output) except Exception as e: parser.exit(1, "ERROR: issue in output creation: %s" % e) if args.verbose: print("> DONE")
if __name__ == "__main__": main()