-
Notifications
You must be signed in to change notification settings - Fork 2
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
- Loading branch information
1 parent
2f0c4a7
commit 8ec7431
Showing
4 changed files
with
231 additions
and
27 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,145 @@ | ||
| """ | ||
| Anna Tips 4 MODIS | ||
| - all bands are in a single file | ||
| - Format: hdf5 file | ||
| - resolution dependent [0.5km, 0.25km, 1km] | ||
| - downloader has the correct Level! | ||
| - time will be painful | ||
| * multiple files 4 multiple SWATHS | ||
| * very little revisit time during the day | ||
| - Filtering Locations / Bounding Boxes / Tiles | ||
| - Day & Night Flag: file sizes, filler value | ||
| - Level 1B - SWATH Product | ||
| potentially useful packages: | ||
| - modis-tools: https://github.com/fraymio/modis-tools | ||
| - pyMODIS: http://www.pymodis.org and https://github.com/lucadelu/pyModis | ||
| ------------------------------------------------------------------------- | ||
| Summary: The download script to interact directly with the goes2go package. | ||
| We only want to specify what is necessary and compatible with the goes2go package. | ||
| In general, we want to satellite, the spatial domain, and the period. | ||
| Args: | ||
| satellite_number: | ||
| **Input Parameters** | ||
| Downloading: | ||
| - satellite number: int --> (16,17,18) | ||
| - spatial extent: str --> full disk (F), CONUS (C), Mesoscale domains (M, M1, M2) | ||
| - goes instrument: str --> e.g. ABI radiance (or SUVI for helio) | ||
| - preprocessing level: str --> e.g. level-1b | ||
| - directory: str | ||
| - return xarray dataset or list of files --> return as file list works better? | ||
| - band specifications: list[int] --> download all or subset only | ||
| - start time (of range of times to be downloaded): str | ||
| - end time: str | ||
| - timesteps/number of files: str | ||
| - day vs. night mode: --> e.g. for only downloading day mode images | ||
| ---------- | ||
| --- | ||
| Basic Processing: | ||
| - resolution: --> downscale all bands to common resolution (e.g. 2 km) | ||
| - coordinate system transformations | ||
| - etc. | ||
| ================= | ||
| INPUT PARAMETERS | ||
| ================= | ||
| # LIST OF DATES | START DATE, END DATE, STEP | ||
| np.arange, np.linspace | ||
| t0, t1, dt | num_files | ||
| timestamps = [t0, t1, t2] | ||
| # create list of dates | ||
| list_of_dates: list = ["2020-10-19 12:00:00", "2020-10-12 12:00:00", ...] | ||
| # SATELLITE INFORMATION | ||
| satellite_number: int = (16, 17, 18)i | ||
| instrument: str = (ABI, ...) | ||
| processing_level: str = (level-1b,): str = (level-1b, ...) | ||
| data_product: str = (radiances, ...) | ||
| # LIST OF BANDS | ||
| list_of_bands: list = [1, 2, ..., 15, 16] | ||
| # TARGET-GRID | ||
| target_grid: xr.Dataset = ... | ||
| % =============== | ||
| HOW DO WE CHECK DAYTIME HOURS? | ||
| * Get Centroid for SATELLITE FOV - FIXED | ||
| * Get Radius points for SATELLITE FOV - FIXED | ||
| * Check if centroid and/or radius points for FOV is within daytime hours | ||
| * Add warning if chosen date is before GOES orbit was stabilized | ||
| * True: | ||
| download that date | ||
| False: | ||
| Skippppp / Continue | ||
| @dataclass | ||
| class SatelliteFOV: | ||
| lon_min: float | ||
| lon_max: float | ||
| lat_min: float | ||
| lat_max: float | ||
| viewing_angle: ... # [0.15, -0.15] [rad] | ||
| @property | ||
| def get_radius(self): | ||
| ... | ||
| class GOES16(SatelliteFOV): | ||
| ... | ||
| ================= | ||
| ALGORITHMS | ||
| ================= | ||
| for itime in timestamps: | ||
| for iband in list_of_bands: | ||
| # ------------------------------------------------- | ||
| # download to folder - use GOES2GO loader | ||
| # ------------------------------------------------- | ||
| # open data in folder | ||
| # ------------------------------------------------- | ||
| # quality check 1 - did it download | ||
| # ------------------------------------------------- | ||
| if download_criteria: | ||
| continue if allow_missing else break | ||
| # quality check 2 - day and/or night specification | ||
| if day_night_criteria: | ||
| continue if allow_missing else break | ||
| # ------------------------------------------------- | ||
| # CRS Transformation (Optional, preferred) | ||
| # ------------------------------------------------- | ||
| # load dataset | ||
| ds: xr.Dataset = xr.load_dataset(...) | ||
| # coordinate transformation | ||
| ds: xr.Dataset = crs_transform(ds, target_crs, *args, **kwargs) | ||
| # resave | ||
| ds.to_netcdf(...) | ||
| # ------------------------------------------------- | ||
| # downsample/upscale/lower-res (optional, preferred) | ||
| # ------------------------------------------------- | ||
| # load dataset | ||
| ds: xr.Dataset = xr.load_dataset(...) | ||
| # resample | ||
| ds: xr.Dataset = downsample(ds, target_grid, *args, **kwargs) | ||
| ds: xr.Dataset = transform_coords(ds, target_coords) | ||
| # resave | ||
| ds.to_netcdf(...) | ||
| """ |