[GRASS5] r.to.sites: change proposal for -a
Markus Neteler
neteler at geog.uni-hannover.de
Tue Oct 24 10:30:33 EDT 2000
On Tue, Oct 24, 2000 at 08:38:04AM -0500, Bill Brown - staff wrote:
> >
> >Hi all,
> >
> >a proposal for r.to.sites:
> >Currently when converting raster data to sites format you
> >get a list like this:
> >
> >3571556.25|5775981.25|#106
> >3571731.25|5775981.25|#123.8
> >3572243.75|5775981.25|#109.5
> >3570556.25|5775968.75|#132.5
> >
> >Here we have two problems:
> > - #cat may be INT only
> > - usually you want %doubles format in sites list
> >
>
> I didn't read all the earlier messages, but last I knew,
> cat can be any RASTER_MAP_TYPE. Did you change this?
Hi Bill,
I didn't change that, was just guessing (I should read the programmer's
manual). But I am expecting problems with s.surf.rst & friends if it is not
stored as %double. Or am I wrong?
The s.info only accepts INT cats:
s.info ktb12.5
[...]
Number of DIMENSIONS: 2
--------------------
- - MIN - - - - MAX - -
dim 1 3558006.250000 3574281.250000 Easting
dim 2 5756006.250000 5775993.750000 Northing
Type of CATEGORY information: CELL_TYPE
----------------------------
- - MIN - - - - MAX - -
80 370 No Label
[...]
although there are double cats. At least this is inconsistent
currently.
And if I want to reconvert the list back to raster with s.to.rast,
this module will make an INT raster map:
s.to.rast in=ktb12.5 out=test
Using size option: 0
Finding and opening site list...
ERROR: selected decimal field column no. 1 not present in sites list.
If interpolating, the s.surf.rst is also expecting doubles:
s.surf.rst help
[...]
input Name of the site file with input x|y|%z [%sm]
So it might be easier to modify the r.to.sites than all the other site
modules.
My 0.02 Euro...
Regards
Markus
----------------------------------------
If you want to unsubscribe from GRASS Development Team mailing list write to:
minordomo at geog.uni-hannover.de with
subject 'unsubscribe grass5'
More information about the grass-dev
mailing list