test.t.merge.sh 2.0 KB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152
  1. #!/bin/sh
  2. # Tests the merging of space time datasets
  3. # We need to set a specific region in the
  4. # @preprocess step of this test.
  5. # The region setting should work for UTM and LL test locations
  6. g.region s=0 n=80 w=0 e=120 b=0 t=50 res=10 res3=10 -p3
  7. r.mapcalc --o expr="prec_1 = rand(10, 150)"
  8. r.mapcalc --o expr="prec_2 = rand(20, 250)"
  9. r.mapcalc --o expr="prec_3 = rand(30, 350)"
  10. r.mapcalc --o expr="prec_4 = rand(40, 450)"
  11. r.mapcalc --o expr="prec_5 = rand(50, 550)"
  12. r.mapcalc --o expr="prec_6 = rand(60, 650)"
  13. # Register maps in temporal database
  14. t.register -i --o maps=prec_1,prec_2,prec_3,prec_4,prec_5,prec_6 \
  15. start="2001-01-01" increment="1 month"
  16. # We need to create three space time dataset and registe the maps
  17. # in several space time datasets
  18. t.create --o type=strds temporaltype=absolute output=precip_abs1 \
  19. title="Test" descr="This is the 1 test strds" semantictype=sum
  20. t.register -i --o input=precip_abs1 maps=prec_1,prec_2,prec_3
  21. t.create --o type=strds temporaltype=absolute output=precip_abs2 \
  22. title="Test" descr="This is the 2 test strds" semantictype=sum
  23. t.register --o input=precip_abs2 maps=prec_3,prec_4,prec_5
  24. t.create --o type=strds temporaltype=absolute output=precip_abs3 \
  25. title="Test" descr="This is the 3 test strds" semantictype=sum
  26. t.register --o input=precip_abs3 maps=prec_4,prec_5,prec_6
  27. # @test to merge two and three space time datasets
  28. t.merge inputs=precip_abs1 output=precip_abs4
  29. t.info precip_abs4
  30. t.rast.list precip_abs4
  31. t.unregister type=rast maps=prec_1,prec_2,prec_3 input=precip_abs4
  32. t.merge --o inputs=precip_abs1,precip_abs2,precip_abs3 output=precip_abs4
  33. t.info precip_abs4
  34. t.rast.list precip_abs4
  35. t.merge --o inputs=precip_abs1,precip_abs2,precip_abs3,precip_abs4 output=precip_abs4
  36. t.info precip_abs4
  37. t.rast.list precip_abs4
  38. t.remove type=strds input=precip_abs1,precip_abs2,precip_abs3,precip_abs4
  39. t.unregister type=rast maps=prec_1,prec_2,prec_3,prec_4,prec_5,prec_6
  40. g.remove rast=prec_1,prec_2,prec_3,prec_4,prec_5,prec_6