These codes are used in this project: The spatial-temporal distribution and variation of GOCI-derived suspended sediment along coastal water induced by a typhoon.
Build_array, lmfit, regressionn are used to do exponential regression, and time2original finds out how long it takes to recover to the (assumed) initial state.
(exponential form: TSM_concentration = aexp(-btime)+c)
Also, I wrote another scrips to plot the result of each pixel to make sure the results are reasonable.