From 9a4e5cfb48755f588e93c3056b8b944b4cb6ec9d Mon Sep 17 00:00:00 2001 From: Amir Souri Date: Thu, 7 Oct 2021 08:53:09 -0400 Subject: [PATCH] adding crs --- GEOAkaze/GEOAkaze_mod.py | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/GEOAkaze/GEOAkaze_mod.py b/GEOAkaze/GEOAkaze_mod.py index 6e0358d..bc64221 100644 --- a/GEOAkaze/GEOAkaze_mod.py +++ b/GEOAkaze/GEOAkaze_mod.py @@ -677,7 +677,7 @@ def read_gee_tiff(self): for fname in geefname: try: - src = rasterio.open(fname) + src = rasterio.open(fname,crs='EPSG:3857') except: continue out_trans = src.transform @@ -715,14 +715,14 @@ def read_gee_tiff(self): if (not within_box) and (intersect_box): src_appended = [] for int_box in range(len(intersect_box)): - src = rasterio.open(intersect_box[int_box]) + src = rasterio.open(intersect_box[int_box],crs='EPSG:3857') src_appended.append(src) msi_img, out_trans = rasterio.merge(src_appended) print('Two tiles are chosen from the clim ' + intersect_box) # if there is at least one master to fully enclose the slave elif within_box: print('The chosen clim MSI is ' + within_box[0]) - src = rasterio.open(within_box[0]) + src = rasterio.open(within_box[0],crs='EPSG:3857') out_trans = src.transform msi_img = src.read(1)