I am not sure what you mean about getting multiple samples.
Indeed, “The algorithm relates the change in backscatter intensity to the change in moisture” by computing the minimum, maximum and sum of
VV values over the time range. Then it uses the latest acquisition in your time-range, as well as the min and sensitivity to compute soil moisture:
Mv = ((samples.VV) - min)/(sensitivity);
If you pull out
Mv for earlier acquisitions in the script you will be computing your values with different time-ranges, and I am not sure that would be a robust method.
If you want to do this nevertheless, you can modify
Mv = ((samples.VV) - min)/(sensitivity); for it to read different acquisitions, for example
samples.VV and return the response to a different variable. Then you would return these variables as bands.