Merge pull request #1416 from GeoscienceAustralia/deaintertidal_graphs_bugfixing

Add debug fields to DEA Intertidal to investigate x and y coordinates

Author: erialC-P

dev/services/wms/ows_refactored/sea_ocean_coast/intertidal_c3/style_intertidal_cfg.py

@@ -20,6 +20,7 @@
     "abstract": "Intertidal elevation in metres above Mean Sea Level",
     "custom_includes": {
         "tide_graph_path": "ows_refactored.sea_ocean_coast.intertidal_c3.utils_intertidal.tide_graph_path",  # add custom metadata field
+        "featureinfo_coords": "ows_refactored.sea_ocean_coast.intertidal_c3.utils_intertidal.featureinfo_coords",  # add custom metadata field
     },
     "index_function": {
         "function": "datacube_ows.band_utils.single_band",

dev/services/wms/ows_refactored/sea_ocean_coast/intertidal_c3/utils_intertidal.py

@@ -94,3 +94,16 @@ def tide_graph_path(data, ds):
 
     # Combine into a string
     return f"{base_dir}/{product}/{version}/{region_split}/{year}--P1Y/{product}_{region}_{year}--P1Y_final_tide_graph.png"
+
+
+def featureinfo_coords(data, ds):
+    """
+    Temporary function to test what x and y coordinates are returned on `data`
+    """
+    # Get point from GetFeatureInfo data. On DEA Maps, we
+    # can assume that `data` coordinates are always in "EPSG:3857"
+    epsg3857_y, epsg3857_x = data.y.item(), data.x.item()
+    point_albers = point(y=epsg3857_y, x=epsg3857_x, crs="EPSG:3857").to_crs("EPSG:3577").geom
+    epsg3577_y, epsg3577_x = point_albers.y, point_albers.x
+
+    return f"EPSG:3857 Y: {epsg3857_y}, EPSG:3857 X: {epsg3857_x}, EPSG:3577 Y: {epsg3577_y}, EPSG:3577 X: {epsg3577_x}"