From a4959e8f177a807558464727c257e149096a7af5 Mon Sep 17 00:00:00 2001
From: Anson Liu <ansonl@users.noreply.github.com>
Date: Tue, 7 Nov 2023 18:56:14 -0800
Subject: [PATCH] Add bounds script

---
 scripts/getTIFBounds.py | 15 +++++++++++++++
 1 file changed, 15 insertions(+)
 create mode 100644 scripts/getTIFBounds.py

diff --git a/scripts/getTIFBounds.py b/scripts/getTIFBounds.py
new file mode 100644
index 0000000..a384587
--- /dev/null
+++ b/scripts/getTIFBounds.py
@@ -0,0 +1,15 @@
+from osgeo import gdal
+from osgeo import gdalconst
+from osgeo.gdalconst import GA_ReadOnly
+
+import sys
+
+# Get TIF extents to 3d print
+data = gdal.Open(sys.argv[1], GA_ReadOnly)
+geoTransform = data.GetGeoTransform()
+minx = geoTransform[0]
+maxy = geoTransform[3]
+maxx = minx + geoTransform[1] * data.RasterXSize
+miny = maxy + geoTransform[5] * data.RasterYSize
+print([minx, miny, maxx, maxy])
+print([f"bllat {miny}", f"bllon {minx}", f"trlat {maxy}", f"trlon {maxx}"])
\ No newline at end of file