W =3.9 E = 15 S= 8 and N = 11

#just to check the data:
r.in.xyz -s 10_20_11_Bank1_RR.txt out=Bank1_RR_10_20_11m fs=, --o

# set region, import the data to create rough DEM
g.region n=11.20 s=8.00 e=16.3 w=4.8 t=5 b=0 res=0.03 tbres=1 -ap3
r.in.xyz 10_20_11_Bank1_RR.txt out=Bank1_RR_10_20_11m fs=, --o

# set 0.01m res region, import the point data to create high res DEM
g.region res=0.01
v.in.ascii -zt 10_20_11_Bank1_RR.txt out=Bank1_10_20_11RR fs=, x=1 y=2 z=3
v.surf.rst -z Bank1_10_20_11RR elev=10_20_11_Bank1_RRelrst npmin=100

# do it for all and then create space-time cube
g.region b=1 t=7 tbres=1 res3=0.03 -p3
r.to.rast3 in=12_10_10_Bank1_RRelrst,01_29_11_Bank1_RRelrst,04_70_11_Bank1_RRelrst,08_24_11_Bank1_RRelrst,10_20_11_Bank1_RRelrst out=Bank1_time

# process filtered data
v.in.ascii -zt 10_20_11_DTF.txt out=bank110_20_11_DTF fs=space x=1 y=2 z=3
r.mapcalc bveg_10_20_11=10_20_11_Bank1_RRelrst-bank110_20_11_DTF
v.in.ascii -zt 8_24_11_DTF.txt out=bank108_24_11_DTF fs=space x=1 y=2 z=3
v.in.ascii -zt 4_7_11_DTF.txt out=bank104_07_11_DTF fs=space x=1 y=2 z=3
v.in.ascii -zt 1_29_11_DTF.txt out=bank101_29_11_DTF fs=space x=1 y=2 z=3
v.in.ascii -zt 12_10_10_DTF.txt out=bank112_10_10_DTF fs=space x=1 y=2 z=3
v.surf.rst -z bank110_20_11_DTF elev=bank110_20_11_DTF npmin=100 --o
v.surf.rst -z bank108_24_11_DTF elev=bank108_24_11_DTF npmin=100 --o
v.surf.rst -z bank104_07_11_DTF elev=bank104_07_11_DTF npmin=100 --o
v.surf.rst -z bank101_29_11_DTF elev=bank101_29_11_DTF npmin=100 --o
v.surf.rst -z bank112_10_10_DTF elev=bank112_10_10_DTF npmin=100 --o

r.series in=bank112_10_10_DTF,bank101_29_11_DTF,bank104_07_11_DTF,bank108_24_11_DTF,bank110_20_11_DTF out=bank1filt_min method=minimum
r.series in=bank112_10_10_DTF,bank101_29_11_DTF,bank104_07_11_DTF,bank108_24_11_DTF,bank110_20_11_DTF out=bank1filt_max method=maximum
r.series in=bank112_10_10_DTF,bank101_29_11_DTF,bank104_07_11_DTF,bank108_24_11_DTF,bank110_20_11_DTF out=bank1filt_range method=range
r.series in=bank112_10_10_DTF,bank101_29_11_DTF,bank104_07_11_DTF,bank108_24_11_DTF,bank110_20_11_DTF out=bank1filt_min_time method=min_raster
r.series in=bank112_10_10_DTF,bank101_29_11_DTF,bank104_07_11_DTF,bank108_24_11_DTF,bank110_20_11_DTF out=bank1filt_max_time method=max_raster

#map range of loss and range of gain
r.mapcalc b1f_loss="if(bank1filt_min_time>bank1filt_max_time,bank1filt_range,null())"
r.mapcalc b1f_gain="if(bank1filt_min_time<bank1filt_max_time,bank1filt_range,null())"
# volumes
r.sum b1f_loss
r.sum b1f_gain
# 110546 x 0.01 x 0.01 = 11m3, 6073 x 0.0001 = 0.6m3

d.grid size=1000 origin=0,8.185
d.grid size=1000 origin=0,8.94 co=brown
d.grid size=1000 origin=0,10.85 co=orange

r.colors bank112_10_10_DTF rast=Bank1raw_min
r.colors bank101_29_11_DTF rast=Bank1raw_min
r.colors bank104_07_11_DTF rast=Bank1raw_min
r.colors bank108_24_11_DTF rast=Bank1raw_min
r.colors bank110_20_11_DTF rast=Bank1raw_min
r.shaded.relief bank112_10_10_DTF shaded=bank112_10_10_DTFsh --o
r.shaded.relief bank101_29_11_DTF shaded=bank101_29_11_DTFsh --o
r.shaded.relief bank104_07_11_DTF shaded=bank104_07_11_DTFsh --o
r.shaded.relief bank108_24_11_DTF shaded=bank108_24_11_DTFsh --o
r.shaded.relief bank110_20_11_DTF shaded=bank110_20_11_DTFsh --o

r.mapcalc b1f_diff012911_121010=bank101_29_11_DTF-bank112_10_10_DTF
r.mapcalc b1f_diff047011_012911=bank104_07_11_DTF-bank101_29_11_DTF
r.mapcalc b1f_diff082411_047011=bank108_24_11_DTF-bank104_07_11_DTF
r.mapcalc b1f_diff102011_082411=bank110_20_11_DTF-bank108_24_11_DTF

r.to.rast3 in=bank112_10_10_DTF,bank101_29_11_DTF,bank104_07_11_DTF,bank108_24_11_DTF,bank110_20_11_DTF out=bank1filt_time
xganim view1=bank112_10_10_DTF,bank101_29_11_DTF,bank104_07_11_DTF,bank108_24_11_DTF,bank110_20_11_DTF
xganim view1=bank112_10_10_DTFsh,bank101_29_11_DTFsh,bank104_07_11_DTFsh,bank108_24_11_DTFsh,bank110_20_11_DTFsh


#time series analysis raw
r.mapcalc b1diff012911_121010=12_10_10_Bank1_RRelrst-01_29_11_Bank1_RRelrst
r.mapcalc b1diff047011_012911=01_29_11_Bank1_RRelrst-04_70_11_Bank1_RRelrst

r.series in=12_10_10_Bank1_RRelrst,01_29_11_Bank1_RRelrst,04_70_11_Bank1_RRelrst,08_24_11_Bank1_RRelrst,10_20_11_Bank1_RRelrst out=Bank1raw_min method=minimum
r.series in=12_10_10_Bank1_RRelrst,01_29_11_Bank1_RRelrst,04_70_11_Bank1_RRelrst,08_24_11_Bank1_RRelrst,10_20_11_Bank1_RRelrst out=Bank1raw_max method=maximum
r.series in=12_10_10_Bank1_RRelrst,01_29_11_Bank1_RRelrst,04_70_11_Bank1_RRelrst,08_24_11_Bank1_RRelrst,10_20_11_Bank1_RRelrst out=Bank1raw_range method=range
r.series in=12_10_10_Bank1_RRelrst,01_29_11_Bank1_RRelrst,04_70_11_Bank1_RRelrst,08_24_11_Bank1_RRelrst,10_20_11_Bank1_RRelrst out=Bank1raw_timemax method=max_raster
r.series in=12_10_10_Bank1_RRelrst,01_29_11_Bank1_RRelrst,04_70_11_Bank1_RRelrst,08_24_11_Bank1_RRelrst,10_20_11_Bank1_RRelrst out=Bank1raw_timemin method=min_raster

r.colors bank108_24_11_DTF rast=
r.mapcalc diff

r.colors 12_10_10_Bank1_RRelrst rast=Bank1raw_min
r.colors 01_29_11_Bank1_RRelrst rast=Bank1raw_min
r.colors 04_70_11_Bank1_RRelrst rast=Bank1raw_min
r.colors 08_24_11_Bank1_RRelrst rast=Bank1raw_min
r.colors 10_20_11_Bank1_RRelrst rast=Bank1raw_min

r.shaded.relief 12_10_10_Bank1_RRelrst shaded=12_10_10_Bank1_RRelrstsh
r.shaded.relief 01_29_11_Bank1_RRelrst shaded=01_29_11_Bank1_RRelrstsh
r.shaded.relief 04_70_11_Bank1_RRelrst shaded=04_70_11_Bank1_RRelrstsh
r.shaded.relief 08_24_11_Bank1_RRelrst shaded=08_24_11_Bank1_RRelrstsh
r.shaded.relief 10_20_11_Bank1_RRelrst shaded=10_20_11_Bank1_RRelrstsh

d.his h=12_10_10_Bank1_RRelrst i=12_10_10_Bank1_RRelrstsh
d.his h=01_29_11_Bank1_RRelrst i=01_29_11_Bank1_RRelrstsh
d.his h=04_70_11_Bank1_RRelrst i=04_70_11_Bank1_RRelrstsh
d.his h=08_24_11_Bank1_RRelrst i=08_24_11_Bank1_RRelrstsh
d.his h=10_20_11_Bank1_RRelrst i=10_20_11_Bank1_RRelrstsh