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")
# bag_file = "C:/Users/gmasetti/Google Drive/QC Tools/test data/H12730/H12730_MB_1m_MLLW.bag"
bag_file = "C:/Users/gmasetti/Google Drive/QC Tools/test data/H12657/H12657_MB_2m_MLLW.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)
print("has density? %s" % bag.has_density())
bag_density = bag.density(mask_nan=True)
print(type(bag.density(mask_nan=True)), bag.density(mask_nan=True).shape, bag.density(mask_nan=True).dtype)
print("min: %s" % np.nanmin(bag_density))
print("max: %s" % np.nanmax(bag_density))
# from hydroffice.bag.density import Density2Gdal
# Density2Gdal(bag_density=bag_density, bag_meta=bag_meta, fmt="ascii")
# Density2Gdal(bag_density=bag_density, bag_meta=bag_meta, fmt="geotiff")
# Density2Gdal(bag_density=bag_density, bag_meta=bag_meta, fmt="xyz")
