search

nextgis / dhi

Collection of scripts to create DHI data
MIT License
1 stars 2 forks source link

'Good pixels' calculation script #9

Open simgislab opened 8 years ago

simgislab commented 8 years ago

Use this as a base: https://github.com/nextgis/dhi/blob/master/general/create_combined_dhi.py

  1. For each pixel at each time step calculate the number of years with valid (not-nodata) pixels (value 1-12). This part is done, see https://github.com/nextgis/dhi/blob/master/general/create_combined_dhi.py#L96.
  2. Result will be as many rasters as time steps, i.e. 46 rasters for FPAR8.
  3. Reclassify each one of them according to the rule, if value > 3: then assign 1, else: assign 0.
  4. At each pixel, sum them all up! This will tell for each pixel how many time steps are passing the threshold of > 3 (value 1-46). I.e. how many steps on the phenology curve are kept because of the threshold.
  5. Export these counts to a TIF file (just one)
simgislab commented 8 years ago

check calculation for NDVI/EVI/GPP, too many good pixels