Add TODO.

preprocessing/calculate_ndvi.py

@@ -99,6 +99,7 @@ def calculate_ndvi(extracted_paths: Union[List[str], str], config_file: str, chu
     nir = xr.open_mfdataset(nir_paths, combine = 'nested', concat_dim = 'time', chunks = chunk_size, parallel = True).squeeze(dim = ['band'], drop = True).rename({'band_data': 'nir'}).sortby('y', ascending = False)
     mask = xr.open_mfdataset(mask_paths, combine = 'nested', concat_dim = 'time', chunks = chunk_size, parallel = True).squeeze(dim = ['band'], drop = True).rename({'band_data': 'mask'}).sortby('y', ascending = False)
     
+    # TODO: replace interp by coarsen for 20m resolution: use 20m band as grid reference
     # Get masking condition and resolution management based on provider, select only data inside the bounds of the input shapefile
     if preferred_provider == 'copernicus':
         red = red.sel(x = slice(bounds[0], bounds[2]), y = slice(bounds[3], bounds[1]))