Workshop at OpenGeoHub Summer School 2019 about the analysis of space-time satellite data for disease ecology applications with GRASS GIS and R stats.

Recently, the National Aeronautics and Space Administration (NASA) has re-processed all the MODIS archive and made available the sixth version for all land products. In this post, I share a script that combines different FOSS4G tools to download, mosaic, re-project, subset, import, apply quality assessment flags and build a time series of MODIS LST daily products. First, we use the pymodis tool modis_download to automatically download 4 tiles of the MOD11A1 product for the period 2010-2016.

Several MODIS products come in “8-day” compositions. Suppose we have daily data and we want to aggregate it with this MODIS-like granularity. How can we achieve that? First idea would be to use the great t.rast.aggregate module in GRASS GIS with granularity="8 days”. But this presents a problem. It does not re-start every year, but aggregates the whole series with 8-day granularity, not considering the year. Meanwhile, MODIS 8-day products re-start the aggregation every January 1st.