[GRASSLIST:6443] d.profile: phenology curve plotting

Hello,

is there a function which plots the phenology of a time series?

d.profile is doing it spatially but I would like to look at the temporal curve
of various raster which are a time series.

e.g. d.profile input=raster.jan, raster.feb, .... raster.dec and then select a
pixel/area in the image and receive a time profile for each pixel which I
selected.

thanks in advance, Martin

--
Martin Wegmann

DLR - German Aerospace Center
German Remote Sensing Data Center
@
Dept.of Geography
Remote Sensing and Biodiversity Unit
University of Wuerzburg
Am Hubland
97074 Würzburg

phone: +49-(0)931 - 888 4797
fax: +49-(0)931 - 888 4961
http://www.biota-africa.org
http://www.biogis.de

On Wed, Apr 13, 2005 at 01:58:30PM -0400, Martin Wegmann wrote:

Hello,

is there a function which plots the phenology of a time series?

d.profile is doing it spatially but I would like to look at the temporal curve
of various raster which are a time series.

e.g. d.profile input=raster.jan, raster.feb, .... raster.dec and then select a
pixel/area in the image and receive a time profile for each pixel which I
selected.

hi Martin,

you may use

i.spectral help

Description:
displays spectral response at user specified locations in images

Modifying this script a bit it you could use g.mlist to insert
your time series and visualize. Needs gnuplot.

Example:
http://mpa.itc.it/markus/grass60/i_spectral_lsat7_spearfish.jpg
(LANDSAT7 for Spearfish)

Markus

Hello Markus,

thanks, that is exactly what I am looking for.

However I looked at the script but cannot figure out how/where to insert
g.mlist.

I assume I have to substitute RASTERMAPS="$GIS_OPT_band1 ... " somehow with
the g.mlist output, right?

Martin

On Thursday 14 April 2005 06:19, Markus Neteler wrote:

On Wed, Apr 13, 2005 at 01:58:30PM -0400, Martin Wegmann wrote:
> Hello,
>
> is there a function which plots the phenology of a time series?
>
> d.profile is doing it spatially but I would like to look at the temporal
> curve of various raster which are a time series.
>
> e.g. d.profile input=raster.jan, raster.feb, .... raster.dec and then
> select a pixel/area in the image and receive a time profile for each
> pixel which I selected.

hi Martin,

you may use

i.spectral help

Description:
displays spectral response at user specified locations in images

Modifying this script a bit it you could use g.mlist to insert
your time series and visualize. Needs gnuplot.

Example:
http://mpa.itc.it/markus/grass60/i_spectral_lsat7_spearfish.jpg
(LANDSAT7 for Spearfish)

Markus

--
Martin Wegmann

DLR - German Aerospace Center
German Remote Sensing Data Center
@
Dept.of Geography
Remote Sensing and Biodiversity Unit
University of Wuerzburg
Am Hubland
97074 Würzburg

phone: +49-(0)931 - 888 4797
fax: +49-(0)931 - 888 4961
http://www.biota-africa.org
http://www.biogis.de

However I looked at the script but cannot figure out how/where to
insert g.mlist.

I assume I have to substitute RASTERMAPS="$GIS_OPT_band1 ... " somehow
with the g.mlist output, right?

replace the two RASTERMAPS lines with:

RASTERMAPS="`g.mlist sep=, pattern='*'`"

Backwards `quotes` fills that part of the string with the output from
the program within them, run the inner g.mlist part first to see what
will happen.

Double "quotes" around the string are just precautionary.

If you are using less than 8 maps, just use i.spectral without
modification.

Hamish