From d6cf1a980363d23e81f992acf93f20ee1103fc67 Mon Sep 17 00:00:00 2001
From: giumas <giuseppemasetti@gmail.com>
Date: Sat, 21 Sep 2024 22:51:13 +0200
Subject: [PATCH] updated example

---
 examples/bag_metadata_vertical_datum.py | 21 ++++++++++++++++++++-
 1 file changed, 20 insertions(+), 1 deletion(-)

diff --git a/examples/bag_metadata_vertical_datum.py b/examples/bag_metadata_vertical_datum.py
index 880af6e..63b2ca4 100644
--- a/examples/bag_metadata_vertical_datum.py
+++ b/examples/bag_metadata_vertical_datum.py
@@ -1,5 +1,6 @@
 import os
 import logging
+from osgeo import osr
 
 from hyo2.abc2.lib.logging import set_logging
 from hyo2.bag.bag import BAGFile
@@ -9,7 +10,8 @@
 set_logging(ns_list=['hyo2.bag'])
 logger = logging.getLogger(__name__)
 
-file_bag_0 = os.path.join(Helper.samples_folder(), "bdb_01.bag")
+# file_bag_0 = os.path.join(Helper.samples_folder(), "bdb_01.bag")
+file_bag_0 = r"C:\Users\gmasetti\Downloads\H13745_MB_4m_MLLW_1of1 (1).bag"
 if os.path.exists(file_bag_0):
     logger.debug("- file_bag_0: %s" % file_bag_0)
 
@@ -19,3 +21,20 @@
 meta = Meta(bag_0.metadata())
 
 logger.debug("vertical datum: %s (%s)" % (meta.wkt_vertical_datum, type(meta.wkt_vertical_datum)))
+
+vrs = osr.SpatialReference()
+ret = vrs.ImportFromWkt(meta.wkt_vertical_datum)
+logger.debug("wkt import: %s" % ret)
+
+is_vertical = vrs.IsVertical() == 1
+logger.debug('wkt is vertical: %s' % is_vertical)
+
+is_epsg = vrs.GetAttrValue("AUTHORITY", 0).lower() == 'epsg'
+if is_epsg:
+    epsg = vrs.GetAttrValue("AUTHORITY", 1)
+    logger.debug("EPSG:%s" % epsg)
+
+is_depth = vrs.GetAttrValue("AXIS", 0).lower() == 'depth'
+if is_depth:
+    direction = vrs.GetAttrValue("AXIS", 1)
+    logger.debug("depth is %s" % direction)