[GRASS-user] More maps with same legend (for visually compare)

GRASS GRASS user
1

Dear GRASS-User List,

I have a CSV file with X,Y,Z coords of 150 point in a time series (50 measuring).
I imported it in a for cycle of $i (1…50):

v.in.ascii input=su1_data.csv output=su1_v_$i format=point fs=‘;’ skip=2 x=2
y=3 z=$iplus4 cat=0 --overwrite

and made surface with:

v.surf.rst input=su1_v_$i layer=1 dmax=25 dmin=5 zmult=1 tension=40 segmax=200
npmin=100 zcolumn=dbl_$i elev=su1_v_surf_rst_$i --overwrite

I displayed the maps with:

d.erase
d.rast su1_v_surf_rst_$i
d.legend map=su1_v_surf_rst_$i

and would like compare it visually.

The problem is, that the scaling
and coloring of surfaces is different by each map.
Is it possible to set with a parameter of v.surf.rst that
all maps legend should have use the legend
with the same highest and deepest point?
I don’t find so a parameter.
Is there another solution to visualize the time series?

If somebody have the time to run it, the full script is here:
http://born.try.hu/grass/su1.sh

And the data file is here:
http://born.try.hu/grass/su1_data.csv

Thanks for all help:
Imre Bornemisza

Bornemisza Imre . born@freemail.hu . +36 30 288 24 64

2

On 07/11/08 15:06, Bornemisza Imre wrote:

Dear GRASS-User List,

I have a CSV file with X,Y,Z coords of 150 point in a time series (50 measuring).
I imported it in a for cycle of $i (1..50):

v.in.ascii input=su1_data.csv output=su1_v_$i format=point fs=';' skip=2 x=2
y=3 z=$iplus4 cat=0 --overwrite

and made surface with:

v.surf.rst input=su1_v_$i layer=1 dmax=25 dmin=5 zmult=1 tension=40 segmax=200
npmin=100 zcolumn=dbl_$i elev=su1_v_surf_rst_$i --overwrite

I displayed the maps with:

d.erase
d.rast su1_v_surf_rst_$i
d.legend map=su1_v_surf_rst_$i

and would like compare it visually.

The problem is, that the scaling
and coloring of surfaces is different by each map.
Is it possible to set with a parameter of v.surf.rst that
all maps legend should have use the legend
with the same highest and deepest point?
I don't find so a parameter.
Is there another solut ion to visualize the time series?

This is not a v.surf.rst issue, but a r.colors issues. You can either create a customized color table in a file and use r.colors rules= to apply these rules to all your maps, or, if you have one raster map which is the reference map, you can also use r.colors raster=

Moritz

3

On Fri, Nov 7, 2008 at 3:18 PM, Moritz Lennert
<mlennert@club.worldonline.be> wrote:

On 07/11/08 15:06, Bornemisza Imre wrote:

Dear GRASS-User List,

I have a CSV file with X,Y,Z coords of 150 point in a time series (50
measuring).
I imported it in a for cycle of $i (1..50):

v.in.ascii input=su1_data.csv output=su1_v_$i format=point fs=';' skip=2
x=2
y=3 z=$iplus4 cat=0 --overwrite

and made surface with:

v.surf.rst input=su1_v_$i layer=1 dmax=25 dmin=5 zmult=1 tension=40
segmax=200
npmin=100 zcolumn=dbl_$i elev=su1_v_surf_rst_$i --overwrite

I displayed the maps with:

d.erase
d.rast su1_v_surf_rst_$i
d.legend map=su1_v_surf_rst_$i

and would like compare it visually.

The problem is, that the scaling
and coloring of surfaces is different by each map.
Is it possible to set with a parameter of v.surf.rst that
all maps legend should have use the legend
with the same highest and deepest point?
I don't find so a parameter.
Is there another solut ion to visualize the time series?

This is not a v.surf.rst issue, but a r.colors issues. You can either create
a customized color table in a file and use r.colors rules= to apply these
rules to all your maps, or, if you have one raster map which is the
reference map, you can also use r.colors raster=

Just as a note, I have recently added a rast parameter to d.rast.leg:

Parameters:
....
          rast Name of raster map to generate legend from

Of course this reference raster map should have a range covering all
maps in question.

Markus

4

Bornemisza Imre wrote:

> The problem is, that the scaling
> and coloring of surfaces is different by each map.
> Is it possible to set with a parameter of v.surf.rst that
> all maps legend should have use the legend
> with the same highest and deepest point?
> I don't find so a parameter.
> Is there another solut ion to visualize the time series?

Moritz Lennert:

This is not a v.surf.rst issue, but a r.colors issues. You
can either create a customized color table in a file and use
r.colors rules= to apply these rules to all your maps, or,
if you have one raster map which is the reference map, you
can also use r.colors raster=

added as a wiki FAQ:
http://grass.osgeo.org/wiki/Time_series#Common_legends_for_many_raster_maps

Hamish