Note: This document is for an older version of GRASS GIS that will be discontinued soon. You should upgrade, and read the current manual page.
Curretnly, only STRDS with absolute temporal type are supported.
# Check time MODIS LST raster time series to sample t.rast.list LST_Day_monthly@modis_lst # Set computational region to first raster map of the time series g.region -up raster=MOD11B3.A2015001.h11v05.single_LST_Day_6km \ align=MOD11B3.A2015001.h11v05.single_LST_Day_6km # Create sampling points map with timestamp column for sampling time echo "cat,sampling_time,x,y 1,2015-01-12,-125870.664090128,591821.149756026 2,2015-02-01,77667.7974483332,571718.585653462 3,2015-02-27,351565.233345769,481257.047191923 4,2015-03-10,565154.976935513,375718.585653462 5,2015-04-07,366642.156422692,277718.585653462 6,2015-05-22,178180.617961154,363154.483089359 7,2015-06-12,12334.464115,395821.149756026 8,2016-03-04,112847.28462782,204846.790781667 9,2016-05-12,311360.105140641,26436.5343714103 10,2017-01-03,467154.976935513,8846.79078166676 11,2018-01-01,600334.464115,129462.175397051 12,2016-06-07,690796.002576538,262641.662576539 13,2016-08-16,685770.361550897,441051.918986795 14,2016-12-11,625462.669243205,594333.970268846" > sampling_points.txt # Import sampling points v.in.ascii in=sampling_points.txt out=sampling_points x=3 y=4 cat=1 \ columns='cat int, sampling_time date, x double precision, y double precision' \ separator="," skip=1 # Assign LST values to sampling points at time of sampling v.what.strds.timestamp input=sampling_points \ timestamp_column=sampling_time column=lst_at_sampling_time \ strds=LST_Day_monthly@modis_lst # Check result v.db.select sampling_points cat|sampling_time|x|y|lst_at_sampling_time 1|2015-01-12|-125870.664090128|591821.149756026|13811 2|2015-02-01|77667.7974483332|571718.585653462|13486 3|2015-02-27|351565.233345769|481257.047191923|13793 4|2015-03-10|565154.976935513|375718.585653462|14355 5|2015-04-07|366642.156422692|277718.585653462|14565 6|2015-05-22|178180.617961154|363154.483089359|14864 7|2015-06-12|12334.464115|395821.149756026|15091 8|2016-03-04|112847.28462782|204846.790781667|14641 9|2016-05-12|311360.105140641|26436.5343714103|14982 10|2017-01-03|467154.976935513|8846.79078166676| 11|2018-01-01|600334.464115|129462.175397051| 12|2016-06-07|690796.002576538|262641.662576539|15052 13|2016-08-16|685770.361550897|441051.918986795|15052 14|2016-12-11|625462.669243205|594333.970268846|13920 # Check result for point with category 5 r.what -n map=MOD11B3.A2015091.h11v05.single_LST_Day_6km@modis_lst \ coordinates=366642.156422692,277718.585653462 # add extra column for LST in degree Celsius v.db.addcolumn map=sampling_points column="lst_celsius double precision" # rescale original MODIS LST to degree Celsius v.db.update map=sampling_points column="lst_celsius" query_column="lst_at_sampling_time * 0.02 - 273.15" v.db.select sampling_points
Available at: v.what.strds.timestamp source code (history)
Latest change: Monday Jun 28 07:54:09 2021 in commit: 1cfc0af029a35a5d6c7dae5ca7204d0eb85dbc55
Main index | Vector index | Topics index | Keywords index | Graphical index | Full index
© 2003-2023 GRASS Development Team, GRASS GIS 7.8.9dev Reference Manual