intital commit

Author: ssmahto

InfeRes_package/MASK.py

@@ -67,6 +67,7 @@ def mask(res_name, yearOFcommission, max_wl, point, boundary, dem_file_path, res
 
 
     # [1] =============================  CLIP LANDSAT IMAGES BY THE UTM BOUNDING BOX 
+    print('============ [1] CLIP LANDSAT IMAGES BY THE UTM BOUNDING BOX ===============')
     print("Clipping Landsat images by the bounding box ...")
     clip_count = 0 
     os.chdir(res_directory + "/LandsatData")
@@ -102,6 +103,7 @@ def mask(res_name, yearOFcommission, max_wl, point, boundary, dem_file_path, res
 
 
     # [2] =============================== DELETE >80% cloudy (over the reservoir) images
+    print('============ [2] DELETE >80% cloudy (over the reservoir) images ===============')
     print("Estimating cloud fraction...")
     class_count = 0 
     cloud_threshold = 80
@@ -182,6 +184,7 @@ def mask(res_name, yearOFcommission, max_wl, point, boundary, dem_file_path, res
 
 
     # [3] ===================================  NDWI CALCULATION (adding cloud mask)
+    print('============ [3] NDWI CALCULATION (adding cloud mask) ===============')
     print("Adding cloud mask to NDWI images...")
     count = 0 
     os.chdir(res_directory + "/LandsatData_Clip")
@@ -241,6 +244,7 @@ def mask(res_name, yearOFcommission, max_wl, point, boundary, dem_file_path, res
 
     # [4] ==============================  CREATE DEM-BASED MAX WATER EXTENT MASK    
     # DEM is preprocessed to have the same cell size and alignment with Landsat images 
+    print('============ [4] CREATE DEM-BASED MAX WATER EXTENT MASK ===============')
     print("Creating DEM-based max water extent mask ...") 
     os.chdir(res_directory +  "/Outputs") 
     res_dem_file = res_name + "DEM.tif"
@@ -268,6 +272,7 @@ def mask(res_name, yearOFcommission, max_wl, point, boundary, dem_file_path, res
 
 
     # [5] ============================  CREATE LANDSAT-BASED MAX WATER EXTENT MASK
+    print('============ [5] CREATE LANDSAT-BASED MAX WATER EXTENT MASK ===============')
     print("Creating Landsat-based max water extent mask ...")
     os.chdir(res_directory +  "/Outputs") 
     res_dem_file = res_name + "DEM.tif"
@@ -292,7 +297,7 @@ def mask(res_name, yearOFcommission, max_wl, point, boundary, dem_file_path, res
                 cloud_percentage = round(np.nansum(ndwi)/np.nansum(res_iso)*100,2)
                 print(str(cloud_percentage) + '% cloud')
                 if cloud_percentage < 20:
-                    print('Image qualified for creating Landsat-based MASK')
+                    print('-------------------------------------------------')
                     ndwi = gdal_array.LoadFile(filename).astype(np.float32)
                     water = ndwi
                     water[np.where(ndwi >= 0)] = 1
@@ -340,6 +345,7 @@ def mask(res_name, yearOFcommission, max_wl, point, boundary, dem_file_path, res
 
 
     # [6] ======================  CREATE EXPANDED MASK (by 3 pixels surrounding each of water pixels)
+    print('============ [6] CREATE EXPANDED MASK ===============')
     print("Creating expanded mask ...")
     os.chdir(res_directory +  "/Outputs")
     mask_1 = gdal_array.LoadFile("Landsat_Mask.tif").astype(np.float32)
@@ -366,6 +372,7 @@ def mask(res_name, yearOFcommission, max_wl, point, boundary, dem_file_path, res
 
 
     # [7] =================================  CREATE 50-ZONE MAP (FREQUENCE MAP)
+    print('============ [7] CREATE 50-ZONE MAP (FREQUENCE MAP) ===============')
     print("Creating 50-zone map (frequence map) ...")
     os.chdir(res_directory +  "/Outputs")
     count = gdal_array.LoadFile("Count.tif").astype(np.float32)

InfeRes_package/data_download.py

@@ -52,7 +52,7 @@ def clip_image_to_geometry(image):
 
     slno = l8images_clip.toList(l8images_clip.size())
     count = 1
-    for i in range(num_images):  # range(num_images)
+    for i in range(2):  # range(num_images)
         print(count)
         selected_image = ee.Image(slno.get(i))
 
@@ -142,7 +142,7 @@ def clip_image_to_geometry(image):
     collection_id = "LANDSAT/LC08/C02/T1_TOA"      
     start_data = '2001-01-01'
     end_date = '2022-12-31'
-    #get_landsat_images(Satellite, row, path, point_coordinates, boundary_coordinates, collection_id, start_data, end_date)
+    get_landsat_images(Satellite, row, path, point_coordinates, boundary_coordinates, collection_id, start_data, end_date)
     print("Congratulations...all Landsat-8 files have successfully downloaded!")
 
     #====================================>> USER INPUT PARAMETERS (Landsat-7 Image Specifications)

InfeRes_package/data_processing.py

@@ -24,12 +24,12 @@
     os.chdir(parent_directory)
     res_name = "Xiaowan"                        # Name of the reservoir 
     # A point within the reservoir [longitude, latitude]
-    point = [99.95, 24.745]
+    point = [100.420, 24.66]
     # Upper-Left and Lower-right coordinates. Example coordinates [longitude, latitude]
-    boundary = [99.20, 25.60, 100.25, 24.65]
-    max_wl = 1236 
-    dead_wl = 1162
-    yearOFcommission = 2010
+    boundary = [100.115, 24.79, 100.51, 24.61]
+    max_wl = 1020 
+    dead_wl = 980
+    yearOFcommission = 2000
     Number_of_tiles = 1                             
     os.makedirs(res_name, exist_ok=True)                  
     os.chdir(parent_directory + res_name)