Steps followed to retrieve SM from CyGNSS data over a 36 km grid (TxSON) using Artificial Neural Networks are given below. Our CyGNSS-derived SM product is also compared with TxSON in-situ data, SMAP and NASA's L3 SM product.
- Download CyGNSS data -> get_cygnss_raw_or.sh
- Extract CyGNSS data over a 36 km region (This can be changed within the file) -> within_36km.py (conda activate py3.8_cygnsslib)
- Download NDVI -> MOD13A1 from us.earthdata.nasa (MODIS/TERRA Vegetation Indices 16 day L3 Global 500m SIN Grid V006
- Convert NDVI .hdf file to .tiff file in ArcGIS
- Reproject NDVI .tiff file to EPSG:4326 -> reproject_ndvi.ipynb (conda activate geo_env_2)
- Upscale TxSON (in-situ) data to a 9km grid at an hourly scale using Voronoi Method -> upscaled_TxSON_9km.mat
- Make the dataset for ML models (Contains code to aggregate txson insitu data to a daily scale) -> make_ml_dataset.ipynb
- Train ANN on data for all months of the year -> ann_all_months_9km.ipynb
- Comparison with TxSON, SMAP and NASA's CyGNSS L3 product (This comparison is at a daily 36 km grid resolution) -> l3_smap_comp_all_months.ipynb