[GRASS5] [bug #3246] (grass) d.vect.thematic fails when the map is in another mapset

GRASS GIS GRASS GIS developer
1

sorry for this doubled post Michael - I forgot to CC the grass5 list the
previous time

Michael

Unfortunatelly the bug remains - it is impossible to use d.vect thematic on a
map which is not in the same mapset where the session was started. There is
always " Table '(...)' doesn't exist." information. Could you please try to
fix it?

The GRASSRGB switch works like a charm. Thanks!

But there is incostistency in the legend provided by d.vect.thematic -l:

Value range: 57.920000 - 69.400000
Mapped by 1 intervals of 11.5948

Color(R:G:B) Value
============ ==========
250:0:0 >57.3408 - 69.5148

As you can see the calculated "Value range" is different than in the legend's
"Value" column. No matter what interval is used (here nint=1) there is always
a difference. Could it be fixed too?

Also I think that the ">" sign between the color and the the value in the
legend should be removed - it doesn't help in anything and worsens the clarity
of the information. How do you think?

Maciek

-------------------------------------------- Managed by Request Tracker

2

1. This math interval difference only happens for the first interval
(min-value) for invterval mapping. This is to make sure that ALL values are
included and each value is included only once. If I specify [>lower value
and <= higher value] for each interval, I miss the minimum value of the
total range. If I specify [>=lower value and <=higher value] for each range,
some values get counted twice. So the actual range of values mapped is
correct, but the lower value shown in the legend is slightly below the
actual minimum. The "11.5948" instead of "11.60" is some kind of internal
rounding error that I have no idea how to correct.

2. I use the notation >lower value - higher value because of several
requests by users to do so. This makes it clearer that the interval is
greater (i.e., not inclusive of) the minimum and less than or equal to
(i.e., inclusive of) the maximum.

Michael
______________________________
Michael Barton, Professor of Anthropology
School of Human Evolution and Social Change
Arizona State University
Tempe, AZ 85287-2402
USA

voice: 480-965-6262; fax: 480-965-7671
www: http://www.public.asu.edu/~cmbarton

From: Maciek Sieczka via RT <grass-bugs@intevation.de>
Reply-To: Maciek Sieczka via RT <grass-bugs@intevation.de>
Date: Tue, 17 May 2005 17:59:21 +0200 (CEST)
Cc: <grass5@grass.itc.it>, <michael.barton@asu.edu>
Subject: [bug #3246] (grass) d.vect.thematic fails when the map is in another
mapset

sorry for this doubled post Michael - I forgot to CC the grass5 list the

previous time

Michael

Unfortunatelly the bug remains - it is impossible to use d.vect thematic on a

map which is not in the same mapset where the session was started. There is

always " Table '(...)' doesn't exist." information. Could you please try to

fix it?

The GRASSRGB switch works like a charm. Thanks!

But there is incostistency in the legend provided by d.vect.thematic -l:

Value range: 57.920000 - 69.400000

Mapped by 1 intervals of 11.5948

Color(R:G:B) Value

============ ==========

250:0:0 >57.3408 - 69.5148

As you can see the calculated "Value range" is different than in the legend's

"Value" column. No matter what interval is used (here nint=1) there is always

a difference. Could it be fixed too?

Also I think that the ">" sign between the color and the the value in the

legend should be removed - it doesn't help in anything and worsens the clarity

of the information. How do you think?

Maciek

-------------------------------------------- Managed by Request Tracker

3

From: "Michael Barton" <michael.barton@asu.edu>

1. This math interval difference only happens for the first interval
(min-value) for invterval mapping.

Actually it also happens for the last interval too - max in the legend and
the actual max are never the same. Is it also by design?

This is to make sure that ALL values are
included and each value is included only once. If I specify [>lower value
and <= higher value] for each interval, I miss the minimum value of the
total range. If I specify [>=lower value and <=higher value] for each
range,
some values get counted twice. So the actual range of values mapped is
correct, but the lower value shown in the legend is slightly below the
actual minimum. The "11.5948" instead of "11.60" is some kind of internal
rounding error that I have no idea how to correct.

2. I use the notation >lower value - higher value because of several
requests by users to do so. This makes it clearer that the interval is
greater (i.e., not inclusive of) the minimum and less than or equal to
(i.e., inclusive of) the maximum.

Michael

Understood. But anyway if I use the d.vect.thematic -l output for
publication it may look supsicious for the eventual reader. I'm not
complaining here - I rerally apprecaite your effort and just trying to be helpfull. This legend looks suspicious and if somebody could fix it would great.

And what about the bug remaining (as in the title)? Any chances it will be
fixed?

Maciek

______________________________
Michael Barton, Professor of Anthropology
School of Human Evolution and Social Change
Arizona State University
Tempe, AZ 85287-2402
USA

voice: 480-965-6262; fax: 480-965-7671
www: http://www.public.asu.edu/~cmbarton

From: Maciek Sieczka via RT <grass-bugs@intevation.de>
Reply-To: Maciek Sieczka via RT <grass-bugs@intevation.de>
Date: Tue, 17 May 2005 17:59:21 +0200 (CEST)
Cc: <grass5@grass.itc.it>, <michael.barton@asu.edu>
Subject: [bug #3246] (grass) d.vect.thematic fails when the map is in
another
mapset

sorry for this doubled post Michael - I forgot to CC the grass5 list the

previous time

Michael

Unfortunatelly the bug remains - it is impossible to use d.vect thematic
on a

map which is not in the same mapset where the session was started. There
is

always " Table '(...)' doesn't exist." information. Could you please try
to

fix it?

The GRASSRGB switch works like a charm. Thanks!

But there is incostistency in the legend provided by d.vect.thematic -l:

Value range: 57.920000 - 69.400000

Mapped by 1 intervals of 11.5948

Color(R:G:B) Value

============ ==========

250:0:0 >57.3408 - 69.5148

As you can see the calculated "Value range" is different than in the
legend's

"Value" column. No matter what interval is used (here nint=1) there is
always

a difference. Could it be fixed too?

Also I think that the ">" sign between the color and the the value in the

legend should be removed - it doesn't help in anything and worsens the
clarity

of the information. How do you think?

Maciek

-------------------------------------------- Managed by Request Tracker

_______________________________________________
grass5 mailing list
grass5@grass.itc.it
http://grass.itc.it/mailman/listinfo/grass5

4

Maciek,

I fixed the problem with d.vect.thematic not working with maps from a
different mapset. I think I've also improved the interval thematic mapping.
Please check out the newest version I just committed to the cvs and posted
to my website.

I tend to agree with you that

min - max

looks better than

min - max

although the second one is more accurate.

What does everyone else think?

Michael
______________________________
Michael Barton, Professor of Anthropology
School of Human Evolution and Social Change
Arizona State University
Tempe, AZ 85287-2402
USA

voice: 480-965-6262; fax: 480-965-7671
www: http://www.public.asu.edu/~cmbarton

From: Maciek Sieczka <werchowyna@epf.pl>
Date: Thu, 19 May 2005 19:39:31 +0200
To: Michael Barton <michael.barton@asu.edu>, Paolo Cavallini via RT
<grass-bugs@intevation.de>
Cc: grass devel <grass5@grass.itc.it>
Subject: Re: [GRASS5] Re: [bug #3246] (grass) d.vect.thematic fails when the
map is in another mapset

From: "Michael Barton" <michael.barton@asu.edu>

1. This math interval difference only happens for the first interval
(min-value) for invterval mapping.

Actually it also happens for the last interval too - max in the legend and
the actual max are never the same. Is it also by design?

This is to make sure that ALL values are
included and each value is included only once. If I specify [>lower value
and <= higher value] for each interval, I miss the minimum value of the
total range. If I specify [>=lower value and <=higher value] for each
range,
some values get counted twice. So the actual range of values mapped is
correct, but the lower value shown in the legend is slightly below the
actual minimum. The "11.5948" instead of "11.60" is some kind of internal
rounding error that I have no idea how to correct.

2. I use the notation >lower value - higher value because of several
requests by users to do so. This makes it clearer that the interval is
greater (i.e., not inclusive of) the minimum and less than or equal to
(i.e., inclusive of) the maximum.

Michael

Understood. But anyway if I use the d.vect.thematic -l output for
publication it may look supsicious for the eventual reader. I'm not
complaining here - I rerally apprecaite your effort and just trying to be
helpfull. This legend looks suspicious and if somebody could fix it would
great.

And what about the bug remaining (as in the title)? Any chances it will be
fixed?

Maciek

______________________________
Michael Barton, Professor of Anthropology
School of Human Evolution and Social Change
Arizona State University
Tempe, AZ 85287-2402
USA

voice: 480-965-6262; fax: 480-965-7671
www: http://www.public.asu.edu/~cmbarton

From: Maciek Sieczka via RT <grass-bugs@intevation.de>
Reply-To: Maciek Sieczka via RT <grass-bugs@intevation.de>
Date: Tue, 17 May 2005 17:59:21 +0200 (CEST)
Cc: <grass5@grass.itc.it>, <michael.barton@asu.edu>
Subject: [bug #3246] (grass) d.vect.thematic fails when the map is in
another
mapset

sorry for this doubled post Michael - I forgot to CC the grass5 list the

previous time

Michael

Unfortunatelly the bug remains - it is impossible to use d.vect thematic
on a

map which is not in the same mapset where the session was started. There
is

always " Table '(...)' doesn't exist." information. Could you please try
to

fix it?

The GRASSRGB switch works like a charm. Thanks!

But there is incostistency in the legend provided by d.vect.thematic -l:

Value range: 57.920000 - 69.400000

Mapped by 1 intervals of 11.5948

Color(R:G:B) Value

============ ==========

250:0:0 >57.3408 - 69.5148

As you can see the calculated "Value range" is different than in the
legend's

"Value" column. No matter what interval is used (here nint=1) there is
always

a difference. Could it be fixed too?

Also I think that the ">" sign between the color and the the value in the

legend should be removed - it doesn't help in anything and worsens the
clarity

of the information. How do you think?

Maciek

-------------------------------------------- Managed by Request Tracker

_______________________________________________
grass5 mailing list
grass5@grass.itc.it
http://grass.itc.it/mailman/listinfo/grass5

5

From: Michael Barton <michael.barton@asu.edu>
Sent: Fri, 20 May 2005 10:12:41 -0700
[...]

I tend to agree with you that

min - max

looks better than

>min - max

although the second one is more accurate.

What does everyone else think?

I'd like the math standard :

[min ; cut1]
]cut1 ; cut2]
...
]cutn ; max]

Note the *included* min. I'll see if I can put that into your current script,
Michael.

Daniel.

6

If you want to change your local version, please go right ahead and see how
it looks.

Let's not change the cvs version until we have a consensus (if possible).
For work in my field (archaeology) and from what I know of geography, the
notation

]min ; max]

would be meaningless--although it might be perfect for your field

Michael
______________________________
Michael Barton, Professor of Anthropology
School of Human Evolution and Social Change
Arizona State University
Tempe, AZ 85287-2402
USA

voice: 480-965-6262; fax: 480-965-7671
www: http://www.public.asu.edu/~cmbarton

From: Daniel Calvelo Aros <dcalvelo@minag.gob.pe>
Reply-To: <daniel.calvelo@minag.gob.pe>
Date: Fri, 20 May 2005 12:22:08 -0500
To: Michael Barton <michael.barton@asu.edu>, Maciek Sieczka
<werchowyna@epf.pl>, Paolo Cavallini via RT <grass-bugs@intevation.de>
Cc: grass devel <grass5@grass.itc.it>
Subject: Re: [GRASS5] Re: [bug #3246] (grass) d.vect.thematic fails when the
map is in another mapset

From: Michael Barton <michael.barton@asu.edu>
Sent: Fri, 20 May 2005 10:12:41 -0700
[...]

I tend to agree with you that

min - max

looks better than

min - max

although the second one is more accurate.

What does everyone else think?

I'd like the math standard :

[min ; cut1]
]cut1 ; cut2]
...
]cutn ; max]

Note the *included* min. I'll see if I can put that into your current script,
Michael.

Daniel.

7

Michael,

I would like to play with d.vect.thematic. I have installed Lorenzo Moretti's GRASS distribution for Mac OS X for both 5.4 & 6.0 and use both. What do I need to do to try it out? Thanks for your help.

Regards,
Tom

Michael Barton wrote:

If you want to change your local version, please go right ahead and see how
it looks.

Let's not change the cvs version until we have a consensus (if possible).
For work in my field (archaeology) and from what I know of geography, the
notation

]min ; max]

would be meaningless--although it might be perfect for your field

Michael
______________________________
Michael Barton, Professor of Anthropology
School of Human Evolution and Social Change
Arizona State University
Tempe, AZ 85287-2402
USA

voice: 480-965-6262; fax: 480-965-7671
www: http://www.public.asu.edu/~cmbarton

From: Daniel Calvelo Aros <dcalvelo@minag.gob.pe>
Reply-To: <daniel.calvelo@minag.gob.pe>
Date: Fri, 20 May 2005 12:22:08 -0500
To: Michael Barton <michael.barton@asu.edu>, Maciek Sieczka
<werchowyna@epf.pl>, Paolo Cavallini via RT <grass-bugs@intevation.de>
Cc: grass devel <grass5@grass.itc.it>
Subject: Re: [GRASS5] Re: [bug #3246] (grass) d.vect.thematic fails when the
map is in another mapset

From: Michael Barton <michael.barton@asu.edu>
Sent: Fri, 20 May 2005 10:12:41 -0700
[...]
   

I tend to agree with you that

min - max

looks better than

min - max
       

although the second one is more accurate.

What does everyone else think?
     

I'd like the math standard :

[min ; cut1]
]cut1 ; cut2]
...
]cutn ; max]

Note the *included* min. I'll see if I can put that into your current script,
Michael.

Daniel.
   
_______________________________________________
grass5 mailing list
grass5@grass.itc.it
http://grass.itc.it/mailman/listinfo/grass5

--
Thomas E Adams
National Weather Service
Ohio River Forecast Center
1901 South State Route 134
Wilmington, OH 45177

EMAIL: thomas.adams@noaa.gov

VOICE: 937-383-0528
FAX: 937-383-0033

8

I'm on MacOSX right now. Quick and dirty, since "." is in the PATH:

1. Unpack Michael's thematic.tgz;
2. start grass (I use 6);
3. within grass cd to the thematic/d.vect.thematic directory;
4. issue a 'cp ../v.univar.sh/v.univar.sh .'
5. presto.

You won't see colors in the graphical legend, though (you need CVS for that).

Daniel

-- Daniel Calvelo Aros

---------- Original Message -----------
From: "Thomas Adams" <Thomas.Adams@noaa.gov>
To: Michael Barton <michael.barton@asu.edu>
Cc: grass5@grass.itc.it
Sent: Fri, 20 May 2005 14:07:20 -0400
Subject: [GRASS5] d.vect.thematic on Mac OS X

Michael,

I would like to play with d.vect.thematic. I have installed
LorenzoMoretti's GRASS distribution for Mac OS X for both 5.4 & 6.0
anduse both. What do I need to do to try it out? Thanks for your help.

Regards,
Tom

Michael Barton wrote:
If you want to change your local version, please go right ahead and
see how it looks.

Let's not change the cvs version until we have a consensus (if
possible). For work in my field (archaeology) and from what I know
of geography, the notation

]min ; max]

would be meaningless--although it might be perfect for your field

Michael
______________________________
Michael Barton, Professor of Anthropology
School of Human Evolution and Social Change
Arizona State University
Tempe, AZ 85287-2402
USA

voice: 480-965-6262; fax: 480-965-7671
www: http://www.public.asu.edu/~cmbarton

From: Daniel Calvelo Aros <dcalvelo@minag.gob.pe>
Reply-To: <daniel.calvelo@minag.gob.pe>
Date: Fri, 20 May 2005 12:22:08 -0500
To: Michael Barton <michael.barton@asu.edu>, Maciek Sieczka
<werchowyna@epf.pl>, Paolo Cavallini via RT <grass-
bugs@intevation.de> Cc: grass devel <grass5@grass.itc.it> Subject:
Re: [GRASS5] Re: [bug #3246] (grass) d.vect.thematic fails when the
map is in another mapset

From: Michael Barton <michael.barton@asu.edu>
Sent: Fri, 20 May 2005 10:12:41 -0700
[...]
I tend to agree with you that

min - max

looks better than

min - max
although the second one is more accurate.

What does everyone else think?
I'd like the math standard :

[min ; cut1]
]cut1 ; cut2]
...
]cutn ; max]

Note the *included* min. I'll see if I can put that into your
current script, Michael.

Daniel.

_______________________________________________
grass5 mailing list
grass5@grass.itc.it
http://grass.itc.it/mailman/listinfo/grass5

--
Thomas E Adams
National Weather Service
Ohio River Forecast Center
1901 South State Route 134
Wilmington, OH 45177

EMAIL: thomas.adams@noaa.gov

VOICE: 937-383-0528
FAX: 937-383-0033

------- End of Original Message -------

9

From: Michael Barton <michael.barton.asu.edu>
Sent: Fri, 20 May 2005 10:39:53 -0700> If you want to change your local
version, please go right ahead and

see how it looks.

Let's not change the cvs version until we have a consensus (if
possible). For work in my field (archaeology) and from what I know
of geography, the notation

]min ; max]

would be meaningless--although it might be perfect for your field

Ok, let's wait for things to settle. I attach a heavily overhauled version of
d.vect.thematic that does the following:

- takes care of the min included (it *is* in the first interval)

- is prepared for any bracketing needed, but uses the "a - b" notation currently

- solves the problem of using stdevs and having eg. min>mean-2stdev

The code, though, is much shorter since I've changed the way the script is
architectured. That's why I send it in toto instead of patch:

- I added an (internal, e.g. not accesible to users) "annotation" facility,
which miught be important for stdev-type classifications.

- I've factored all the query draw and print sections to work from a set of
breakpoint definitions, which are calculated differently for each of the
themecalc methods. It complicates thing a little but it allows for easier
classifier additions later, and is less cumbersome for further enhancements
(only one spot for changes).

- I've left the points part alone, since I don't have point data at hand to
test against, and I'd like your comments before refactoring further.

Please test and comment.

Daniel.

PS. this is not standalone as it still requires v.univar.sh

Michael
______________________________
Michael Barton, Professor of Anthropology
School of Human Evolution and Social Change
Arizona State University
Tempe, AZ 85287-2402
USA

voice: 480-965-6262; fax: 480-965-7671
www: http://www.public.asu.edu/~cmbarton

> From: Daniel Calvelo Aros <dcalvelo@minag.gob.pe>
> Reply-To: <daniel.calvelo@minag.gob.pe>
> Date: Fri, 20 May 2005 12:22:08 -0500
> To: Michael Barton <michael.barton@asu.edu>, Maciek Sieczka
> <werchowyna@epf.pl>, Paolo Cavallini via RT <grass-bugs@intevation.de>
> Cc: grass devel <grass5@grass.itc.it>
> Subject: Re: [GRASS5] Re: [bug #3246] (grass) d.vect.thematic fails when the
> map is in another mapset
>
> From: Michael Barton <michael.barton@asu.edu>
> Sent: Fri, 20 May 2005 10:12:41 -0700
> [...]
>>
>> I tend to agree with you that
>>
>> min - max
>>
>> looks better than
>>
>>> min - max
>>
>> although the second one is more accurate.
>>
>> What does everyone else think?
>
> I'd like the math standard :
>
> [min ; cut1]
> ]cut1 ; cut2]
> ...
> ]cutn ; max]
>
> Note the *included* min. I'll see if I can put that into your current script,
> Michael.
>
> Daniel.

_______________________________________________
grass5 mailing list
grass5@grass.itc.it
http://grass.itc.it/mailman/listinfo/grass5

------- End of Original Message -------

(attachments)

d.vect.thematic.dca (29.8 KB)

10

Thanks Daniel.

Here is a more "Mac" way to install it.

This assumes that you are using Lorenzo Moretti's Mac binaries.

1. Go to your grass folder (in Applications).

2. Select the GRASS program file. I did this for grass61cvs, but it should
work for grass60 also.

3. Hold down the control key and click with the mouse (or right click on a
3-button mouse).

4. Select "Show package contents"

5. Navigate to Contents/Resources/grass-6.1cvs/scripts (or
.../grass-6.0/scripts)

6. Put d.vect.thematic and v.univar.sh in the scripts folder

7. Voila, you can run them from the GRASS command line (I'll add them to the
GIS Manager menu soon).

Note that sometimes the execute permission gets turned off when these files
travel across the ether. You can not turn this on via the Mac "get info"
dialog for some weird reason. You can either use the chmod command from the
terminal to do this or download a nifty utility called BatChmod from
www.versiontracker.com that will allow you to set the execute permissions.

Michael
______________________________
Michael Barton, Professor of Anthropology
School of Human Evolution and Social Change
Arizona State University
Tempe, AZ 85287-2402
USA

voice: 480-965-6262; fax: 480-965-7671
www: http://www.public.asu.edu/~cmbarton

From: Daniel Calvelo Aros <dcalvelo@minag.gob.pe>
Reply-To: <daniel.calvelo@minag.gob.pe>
Date: Fri, 20 May 2005 14:22:50 -0500
To: Thomas Adams <Thomas.Adams@noaa.gov>, Michael Barton
<michael.barton@asu.edu>
Cc: <grass5@grass.itc.it>
Subject: Re: [GRASS5] d.vect.thematic on Mac OS X

I'm on MacOSX right now. Quick and dirty, since "." is in the PATH:

1. Unpack Michael's thematic.tgz;
2. start grass (I use 6);
3. within grass cd to the thematic/d.vect.thematic directory;
4. issue a 'cp ../v.univar.sh/v.univar.sh .'
5. presto.

You won't see colors in the graphical legend, though (you need CVS for that).

Daniel

-- Daniel Calvelo Aros

---------- Original Message -----------
From: "Thomas Adams" <Thomas.Adams@noaa.gov>
To: Michael Barton <michael.barton@asu.edu>
Cc: grass5@grass.itc.it
Sent: Fri, 20 May 2005 14:07:20 -0400
Subject: [GRASS5] d.vect.thematic on Mac OS X

Michael,

I would like to play with d.vect.thematic. I have installed
LorenzoMoretti's GRASS distribution for Mac OS X for both 5.4 & 6.0
anduse both. What do I need to do to try it out? Thanks for your help.

Regards,
Tom

Michael Barton wrote:
If you want to change your local version, please go right ahead and
see how it looks.

Let's not change the cvs version until we have a consensus (if
possible). For work in my field (archaeology) and from what I know
of geography, the notation

]min ; max]

would be meaningless--although it might be perfect for your field

Michael
______________________________
Michael Barton, Professor of Anthropology
School of Human Evolution and Social Change
Arizona State University
Tempe, AZ 85287-2402
USA

voice: 480-965-6262; fax: 480-965-7671
www: http://www.public.asu.edu/~cmbarton

From: Daniel Calvelo Aros <dcalvelo@minag.gob.pe>
Reply-To: <daniel.calvelo@minag.gob.pe>
Date: Fri, 20 May 2005 12:22:08 -0500
To: Michael Barton <michael.barton@asu.edu>, Maciek Sieczka
<werchowyna@epf.pl>, Paolo Cavallini via RT <grass-
bugs@intevation.de> Cc: grass devel <grass5@grass.itc.it> Subject:
Re: [GRASS5] Re: [bug #3246] (grass) d.vect.thematic fails when the
map is in another mapset

From: Michael Barton <michael.barton@asu.edu>
Sent: Fri, 20 May 2005 10:12:41 -0700
[...]
I tend to agree with you that

min - max

looks better than

min - max
although the second one is more accurate.

What does everyone else think?
I'd like the math standard :

[min ; cut1]
]cut1 ; cut2]
...
]cutn ; max]

Note the *included* min. I'll see if I can put that into your
current script, Michael.

Daniel.

_______________________________________________
grass5 mailing list
grass5@grass.itc.it
http://grass.itc.it/mailman/listinfo/grass5

--
Thomas E Adams
National Weather Service
Ohio River Forecast Center
1901 South State Route 134
Wilmington, OH 45177

EMAIL: thomas.adams@noaa.gov

VOICE: 937-383-0528
FAX: 937-383-0033

------- End of Original Message -------

11

Daniel,

I tried out your code changes to d.vect.thematic.

There were a few odd display glitches, but overall it worked well.

I've gone ahead and rewritten the graduated points routine along the lines
of your graduated colors routine. It revealed a problem I didn't notice
previously with the graduated points, so I fixed it. I made a couple of
other improvements to the points legend.

I also added a flag option for math notation brackets, so everyone should be
happy (or at least happier).

Could you try it out and see what you think? If you want to test the
graduated point routine, just use areas and select centroids as the vector
type.

Lorenzo Moretti suggested an option to write out the thematic mapping as a
group file for the GIS Manager. This is intriguing and I might try to add
this in the next update if there is some interest.

Michael

On 5/20/05 12:39 PM, "Daniel Calvelo Aros" <dcalvelo@minag.gob.pe> wrote:

From: Michael Barton <michael.barton.asu.edu>
Sent: Fri, 20 May 2005 10:39:53 -0700> If you want to change your local
version, please go right ahead and

see how it looks.

Let's not change the cvs version until we have a consensus (if
possible). For work in my field (archaeology) and from what I know
of geography, the notation

]min ; max]

would be meaningless--although it might be perfect for your field

Ok, let's wait for things to settle. I attach a heavily overhauled version of
d.vect.thematic that does the following:

- takes care of the min included (it *is* in the first interval)

- is prepared for any bracketing needed, but uses the "a - b" notation
currently

- solves the problem of using stdevs and having eg. min>mean-2stdev

The code, though, is much shorter since I've changed the way the script is
architectured. That's why I send it in toto instead of patch:

- I added an (internal, e.g. not accesible to users) "annotation" facility,
which miught be important for stdev-type classifications.

- I've factored all the query draw and print sections to work from a set of
breakpoint definitions, which are calculated differently for each of the
themecalc methods. It complicates thing a little but it allows for easier
classifier additions later, and is less cumbersome for further enhancements
(only one spot for changes).

- I've left the points part alone, since I don't have point data at hand to
test against, and I'd like your comments before refactoring further.

Please test and comment.

Daniel.

PS. this is not standalone as it still requires v.univar.sh

Michael
______________________________
Michael Barton, Professor of Anthropology
School of Human Evolution and Social Change
Arizona State University
Tempe, AZ 85287-2402
USA

voice: 480-965-6262; fax: 480-965-7671
www: http://www.public.asu.edu/~cmbarton

From: Daniel Calvelo Aros <dcalvelo@minag.gob.pe>
Reply-To: <daniel.calvelo@minag.gob.pe>
Date: Fri, 20 May 2005 12:22:08 -0500
To: Michael Barton <michael.barton@asu.edu>, Maciek Sieczka
<werchowyna@epf.pl>, Paolo Cavallini via RT <grass-bugs@intevation.de>
Cc: grass devel <grass5@grass.itc.it>
Subject: Re: [GRASS5] Re: [bug #3246] (grass) d.vect.thematic fails when the
map is in another mapset

From: Michael Barton <michael.barton@asu.edu>
Sent: Fri, 20 May 2005 10:12:41 -0700
[...]

I tend to agree with you that

min - max

looks better than

min - max

although the second one is more accurate.

What does everyone else think?

I'd like the math standard :

[min ; cut1]
]cut1 ; cut2]
...
]cutn ; max]

Note the *included* min. I'll see if I can put that into your current
script,
Michael.

Daniel.

_______________________________________________
grass5 mailing list
grass5@grass.itc.it
http://grass.itc.it/mailman/listinfo/grass5

------- End of Original Message -------

____________________
C. Michael Barton, Professor of Anthropology
School of Human Evolution and Social Change
PO Box 872402
Arizona State University
Tempe, AZ 85287-2402
USA

Phone: 480-965-6262
Fax: 480-965-7671
www: <www.public.asu.edu/~cmbarton>

(attachments)

d.vect.thematic (18.9 KB)

12

Michael & Daniel,

Thanks for your help!

Regards,
Tom

Michael Barton wrote:

Thanks Daniel.

Here is a more "Mac" way to install it.

This assumes that you are using Lorenzo Moretti's Mac binaries.

1. Go to your grass folder (in Applications).

2. Select the GRASS program file. I did this for grass61cvs, but it should
work for grass60 also.

3. Hold down the control key and click with the mouse (or right click on a
3-button mouse).

4. Select "Show package contents"

5. Navigate to Contents/Resources/grass-6.1cvs/scripts (or
.../grass-6.0/scripts)

6. Put d.vect.thematic and v.univar.sh in the scripts folder

7. Voila, you can run them from the GRASS command line (I'll add them to the
GIS Manager menu soon).

Note that sometimes the execute permission gets turned off when these files
travel across the ether. You can not turn this on via the Mac "get info"
dialog for some weird reason. You can either use the chmod command from the
terminal to do this or download a nifty utility called BatChmod from
www.versiontracker.com that will allow you to set the execute permissions.

Michael
______________________________
Michael Barton, Professor of Anthropology
School of Human Evolution and Social Change
Arizona State University
Tempe, AZ 85287-2402
USA

voice: 480-965-6262; fax: 480-965-7671
www: http://www.public.asu.edu/~cmbarton

From: Daniel Calvelo Aros <dcalvelo@minag.gob.pe>
Reply-To: <daniel.calvelo@minag.gob.pe>
Date: Fri, 20 May 2005 14:22:50 -0500
To: Thomas Adams <Thomas.Adams@noaa.gov>, Michael Barton
<michael.barton@asu.edu>
Cc: <grass5@grass.itc.it>
Subject: Re: [GRASS5] d.vect.thematic on Mac OS X

I'm on MacOSX right now. Quick and dirty, since "." is in the PATH:

1. Unpack Michael's thematic.tgz;
2. start grass (I use 6);
3. within grass cd to the thematic/d.vect.thematic directory;
4. issue a 'cp ../v.univar.sh/v.univar.sh .'
5. presto.

You won't see colors in the graphical legend, though (you need CVS for that).

Daniel

-- Daniel Calvelo Aros

---------- Original Message -----------
From: "Thomas Adams" <Thomas.Adams@noaa.gov>
To: Michael Barton <michael.barton@asu.edu>
Cc: grass5@grass.itc.it
Sent: Fri, 20 May 2005 14:07:20 -0400
Subject: [GRASS5] d.vect.thematic on Mac OS X

Michael,

I would like to play with d.vect.thematic. I have installed
LorenzoMoretti's GRASS distribution for Mac OS X for both 5.4 & 6.0
anduse both. What do I need to do to try it out? Thanks for your help.

Regards,
Tom

Michael Barton wrote:
If you want to change your local version, please go right ahead and
see how it looks.

Let's not change the cvs version until we have a consensus (if
possible). For work in my field (archaeology) and from what I know
of geography, the notation

]min ; max]

would be meaningless--although it might be perfect for your field

Michael
______________________________
Michael Barton, Professor of Anthropology
School of Human Evolution and Social Change
Arizona State University
Tempe, AZ 85287-2402
USA

voice: 480-965-6262; fax: 480-965-7671
www: http://www.public.asu.edu/~cmbarton

From: Daniel Calvelo Aros <dcalvelo@minag.gob.pe>
Reply-To: <daniel.calvelo@minag.gob.pe>
Date: Fri, 20 May 2005 12:22:08 -0500
To: Michael Barton <michael.barton@asu.edu>, Maciek Sieczka
<werchowyna@epf.pl>, Paolo Cavallini via RT <grass-
bugs@intevation.de> Cc: grass devel <grass5@grass.itc.it> Subject:
Re: [GRASS5] Re: [bug #3246] (grass) d.vect.thematic fails when the
map is in another mapset

From: Michael Barton <michael.barton@asu.edu>
Sent: Fri, 20 May 2005 10:12:41 -0700
[...]
I tend to agree with you that

min - max

looks better than

min - max
although the second one is more accurate.

What does everyone else think?
I'd like the math standard :

[min ; cut1]
]cut1 ; cut2]
...
]cutn ; max]

Note the *included* min. I'll see if I can put that into your
current script, Michael.

Daniel.

_______________________________________________
grass5 mailing list
grass5@grass.itc.it
http://grass.itc.it/mailman/listinfo/grass5

--
Thomas E Adams
National Weather Service
Ohio River Forecast Center
1901 South State Route 134
Wilmington, OH 45177

EMAIL: thomas.adams@noaa.gov

VOICE: 937-383-0528
FAX: 937-383-0033
     

------- End of Original Message -------

--
Thomas E Adams
National Weather Service
Ohio River Forecast Center
1901 South State Route 134
Wilmington, OH 45177

EMAIL: thomas.adams@noaa.gov

VOICE: 937-383-0528
FAX: 937-383-0033

13

On Sun, May 22, 2005 18:49, Michael Barton said:

Daniel,

I tried out your code changes to d.vect.thematic.

There were a few odd display glitches, but overall it worked well.

Here both versions (your's and Daniel's) work fine.

I've gone ahead and rewritten the graduated points routine along the lines
of your graduated colors routine. It revealed a problem I didn't notice
previously with the graduated points, so I fixed it. I made a couple of
other improvements to the points legend.

The graduated points legend only displays boxes, even when I chose other
symbols. Is this normal ?

I also added a flag option for math notation brackets, so everyone should be
happy (or at least happier).

Yes, this seems like the best solution to me.

Lorenzo Moretti suggested an option to write out the thematic mapping as a
group file for the GIS Manager. This is intriguing and I might try to add
this in the next update if there is some interest.

Sounds good. This would allow using the GIS Manager print command to print the
results to a ps file.

Thanks a lot to both of you for the great work !

Moritz

14

The graduated points legend only displays boxes, even when I chose other
symbols. Is this normal ?

Yes. I'm using d.graph to create the legend. Currently, it only allows the
old GRASS 5 default (with o displaying as an open box). With Hamish's update
to d.graph, it will display the box in the selected color. Perhaps, it is
possible to update d.graph to show the new GRASS 6 icons?

Michael

Thanks a lot to both of you for the great work !

Moritz

____________________
C. Michael Barton, Professor of Anthropology
School of Human Evolution and Social Change
PO Box 872402
Arizona State University
Tempe, AZ 85287-2402
USA

Phone: 480-965-6262
Fax: 480-965-7671
www: <www.public.asu.edu/~cmbarton>