Merge branch 'main' of https://github.com/spaceml-org/rs_tools

rs_tools/_src/data/modis/bands.py

@@ -8,4 +8,5 @@
 
 MODIS_CHANNEL_NUMBERS = [
     str(i) for i in range(1, 39)
-] + ['13lo', '14lo', '13hi', '14hi']
+] + ['13lo', '14lo', '13hi', '14hi']
+}

scripts/goes-download.py

@@ -389,7 +389,7 @@ def _check_save_dir(save_dir: str) -> bool:
         return True
     else:
         try:
-            os.mkdir(save_dir)
+            os.makedirs(save_dir)
             return True
         except:
             msg = "Save directory does not exist"

scripts/modis-download.py

@@ -327,7 +327,7 @@ def _check_save_dir(save_dir: str) -> bool:
         return True
     else:
         try:
-            os.mkdir(save_dir)
+            os.makedirs(save_dir)
             return True
         except:
             msg = "Save directory does not exist"

scripts/pipeline/preprocess_modis.py

@@ -34,23 +34,28 @@ class MODISGeoProcessing:
     save_path: str = "./"
     channels: list[str] = MODIS_CHANNEL_NUMBERS
 
-    def geoprocess_files(self, params: GeoProcessingParams, files: List[str]):
-        save_path: str = params.save_path.join(self.save_path)
-        target_grid: np.ndarray = create_latlon_grid(params.region, self.resolution) # create using np.linspace?
+    # def geoprocess_files(self, params: GeoProcessingParams, files: List[str]):
+    # resolution: float = 0.25 # in degrees?
+    # save_path: str = "./"
+
+    # def geoprocess(self, params: GeoProcessingParams, files: List[str]):
+    #     save_path: str = params.save_path.join(self.save_path)
+    #     target_grid: np.ndarray = create_latlon_grid(params.region, self.resolution) # create using np.linspace?
 
-        # loop through files
-        # open dataset
-        # stack variables to channels
-        # resample
-        # convert units (before or after resampling???)
-        # save as netcdf
+    #     # loop through files
+    #     # open dataset
+    #     # stack variables to channels
+    #     # resample
+    #     # convert units (before or after resampling???)
+    #     # save as netcdf
 
     def geoprocess_file(self, file: str, save_name: str="test") -> None:
 
         # load MODIS SCENE
         ds, time_stamp = self.load_modis_scene(file)
         # convert units
         ds = self.convert_units(file)
+        # TODO: resample
         # save to raster
         time_stamp = time_stamp.strftime("%Y%m%d%H%M")
         save_name =f"modis_{time_stamp}.tif"