import logging from matplotlib import pyplot as plt logger = logging.getLogger() logger.setLevel(logging.NOTSET) 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) from hydroffice.bag import BAGFile from hydroffice.bag import BAGError from hydroffice.bag.helper import Helper bag_file = os.path.join(Helper.samples_folder(), "bdb_00.bag") if os.path.exists(bag_file): print("- file_bag_0: %s" % bag_file) bag = BAGFile(bag_file) bag_meta = bag.populate_metadata() print(bag_meta) bag_uncertainty = bag.uncertainty(mask_nan=True) print(type(bag.elevation(mask_nan=True)), bag.elevation(mask_nan=True).shape, bag.elevation(mask_nan=True).dtype) from hydroffice.bag.uncertainty import Uncertainty2Gdal Uncertainty2Gdal(bag_uncertainty=bag_uncertainty, bag_meta=bag_meta, fmt="ascii", epsg=4326) Uncertainty2Gdal(bag_uncertainty=bag_uncertainty, bag_meta=bag_meta, fmt="geotiff", epsg=4326) Uncertainty2Gdal(bag_uncertainty=bag_uncertainty, bag_meta=bag_meta, fmt="xyz", epsg=4326)
import subprocess import os from hydroffice.bag.helper import Helper from hydroffice.bag import tools tools_folder = os.path.abspath(os.path.dirname(tools.__file__)) tool_path = os.path.join(tools_folder, 'bag_validate.py') print("tool: %s" % tool_path) # help print("\n\n# -h") subprocess.call("python %s -h" % tool_path) # verbose + test file file_bag_0 = os.path.join(Helper.samples_folder(), "bdb_00.bag") print("\n\n# -v %s" % file_bag_0) subprocess.call("python %s -v %s" % (tool_path, file_bag_0)) # verbose + fake file file_bag_0 = os.path.join(Helper.samples_folder(), "not_present_00.bag") print("\n\n# -v %s" % file_bag_0) subprocess.call("python %s -v %s" % (tool_path, file_bag_0)) # test file file_bag_0 = os.path.join(Helper.samples_folder(), "bdb_00.bag") print("\n\n# %s" % file_bag_0) subprocess.call("python %s %s" % (tool_path, file_bag_0))
import os import logging logger = logging.getLogger() logger.setLevel(logging.NOTSET) 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) from hydroffice.bag.helper import BAGError, Helper try: raise BAGError("test") except BAGError as e: print(e) data_folder = Helper.samples_folder() if os.path.exists(data_folder): print("data folder: %s" % data_folder) iso_folder = Helper.iso19139_folder() if os.path.exists(iso_folder): print("iso folder: %s" % iso_folder) iso_folder = Helper.iso19757_3_folder() if os.path.exists(iso_folder): print("iso folder: %s" % iso_folder)
def test_bag_samples_folder(self): assert os.path.exists(Helper.samples_folder())
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) from hydroffice.bag import BAGFile from hydroffice.bag import BAGError from hydroffice.bag.helper import Helper # Example that use bag.BAGFile to: # - open a BAG file # - read the whole elevation and uncertainty layers # - read a selected range of rows for the elevation and uncertainty layers file_bag_0 = os.path.join(Helper.samples_folder(), "bdb_00.bag") if os.path.exists(file_bag_0): logger.debug("- file_bag_0: %s" % file_bag_0) # - open a file bag_0 = BAGFile(file_bag_0) logger.debug("\n%s\n" % bag_0) # - get elevation shape logger.info("elevation shape: %s" % (bag_0.elevation_shape(), )) # - read the full elevation full_elevation = bag_0.elevation(mask_nan=True) logger.info("elevation array:\n type: %s\n shape: %s\n dtype: %s" % (type(full_elevation), full_elevation.shape, full_elevation.dtype)) ax = plt.contourf(full_elevation) plt.colorbar(ax)
def setUp(self): self.file_bag_0 = os.path.join(Helper.samples_folder(), "bdb_00.bag") self.file_bag_1 = os.path.join(Helper.samples_folder(), "bdb_01.bag") self.file_fake_0 = os.path.join(Helper.samples_folder(), "fake_00.bag")