| 1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950 |
- #!/bin/sh
- # Test the temporal and spatial sampling/observation of raster maps by vector point maps
- # We need to set a specific region in the
- # @preprocess step of this test.
- # The region setting should work for UTM and LL test locations
- g.region s=0 n=80 w=0 e=120 b=0 t=50 res=10 res3=10 -p3
- export GRASS_OVERWRITE=1
- r.mapcalc expr="prec_1 = 100.0"
- r.mapcalc expr="prec_2 = 200.0"
- r.mapcalc expr="prec_3 = 300"
- r.mapcalc expr="prec_4 = 400"
- r.mapcalc expr="prec_5 = 500.0"
- r.mapcalc expr="prec_6 = 600.0"
- r.mapcalc expr="prec_7 = 400"
- r.mapcalc expr="prec_8 = 500.0"
- r.mapcalc expr="prec_9 = 600.0"
- v.random output=prec n=5 seed=1
- v.random -z output=test_1 column=test n=5 seed=1
- t.create type=strds temporaltype=absolute output=precip_abs1 title="A test" descr="A test"
- t.register -i input=precip_abs1 maps=prec_1,prec_2,prec_3,prec_4,prec_5,prec_6 start="2001-03-01 00:00:00" increment="1 months"
- t.create type=strds temporaltype=absolute output=precip_abs2 title="A test" descr="A test"
- t.register -i input=precip_abs2 maps=prec_7,prec_8,prec_9 start="2001-05-01 00:00:00" increment="1 months"
- # The @test
- t.vect.observe.strds input=prec strds=precip_abs1 output=prec_observer vector=prec_observer column="test_val"
- v.info prec_observer
- t.info type=stvds input=prec_observer
- t.vect.list input=prec_observer
- t.vect.db.select input=prec_observer
- t.vect.observe.strds columns=test1,test2,test3 input=test_1 \
- strds=precip_abs1,precip_abs1,precip_abs2 output=test_1_observer \
- vector=test_1_observer
- v.info test_1_observer
- t.info type=stvds input=test_1_observer
- t.vect.list input=test_1_observer
- t.vect.db.select input=test_1_observer columns=cat,test1,test2,test3
- # @postprocess
- t.remove -rf type=strds input=precip_abs1,precip_abs2
- t.remove -rf type=stvds input=prec_observer,test_1_observer
|
