RE: [GRASS5] rinascii & routascii

Hello,

Here is the r.in.ascii displacement problem and the failure
of the dp paramater in routascii. There were no compile problems
for either module on my machine (intel with Mandrake 7.2). I
wonder why it fails on my machine but different on others...?

GRASS:/usr/local/grassedit > more rinascii.txt
north: 4504755
south: 4404275
east: 351345
west: 350985
rows: 4
cols: 3
0.0 1.0 0.0
1.0 1.0 1.0
0.0 1.0 1.0
0.0 0.0 1.0

GRASS:/usr/local/grassedit > r.in.ascii -f input=rinascii.txt output=test
CREATING SUPPORT FILES FOR test

GRASS:/usr/local/grassedit > r.out.ascii map=test dp=7 |more
north: 4504755
south: 4504275
east: 351345
west: 350985
rows: 4
cols: 3
0 1 0
0 1 0
0 1 0
0 1 0
GRASS:/usr/local/grassedit > r.out.ascii map=test dp=2 |more
north: 4504755
south: 4504275
east: 351345
west: 350985
rows: 4
cols: 3
0 1 0
0 1 0
0 1 0
0 1 0
GRASS:/usr/local/grassedit > r.out.ascii map=test dp=0 |more
north: 4504755
south: 4504275
east: 351345
west: 350985
rows: 4
cols: 3
0 1 0
0 1 0
0 1 0
0 1 0
GRASS:/usr/local/grassedit > r.report map=test units=c
r.stats: 100%
+---------------------------------------------------------------------------
--+
| RASTER MAP CATEGORY REPORT
|
|LOCATION: blgrassdata Sun Dec 17 12:52:51
2000|
|---------------------------------------------------------------------------
--|
| north: 4504755 east: 351345
|
|REGION south: 4504275 west: 350985
|
| res: 120 res: 120
|
|---------------------------------------------------------------------------
--|
|MASK:none
|
|---------------------------------------------------------------------------
--|
|MAP: (untitled) (test in leproj)
|
|---------------------------------------------------------------------------
--|
| Category Information |
cell|
| #|description
|count|
|---------------------------------------------------------------------------
--|
|0-0.003922|from to . . . . . . . . . . . . . . . . . . . . . . . . . .|
8|
|0.996078-1|from to . . . . . . . . . . . . . . . . . . . . . . . . . .|
4|
|---------------------------------------------------------------------------
--|
|TOTAL |
12|
+---------------------------------------------------------------------------
--+
GRASS:/usr/local/grassedit > r.info test

+---------------------------------------------------------------------------
-+
| Layer: test Date: Sun Dec 17 12:48:29 2000
|
| Mapset: leproj Login of Creator: root
|
| Location: blgrassdata
|
| DataBase: /usr/local/grass5/data
|
| Title: ( test )
|

|---------------------------------------------------------------------------
-|
|
|
| Type of Map: raster Number of Categories: 255
|
| Data Type: FCELL
|
| Rows: 4
|
| Columns: 3
|
| Total Cells: 12
|
| Projection: UTM (zone 13)
|
| N: 4504755 S: 4404275 Res: 25120
|
| E: 351345 W: 350985 Res: 120
|
|
|
| Data Source:
|
|
|
|
|
|
|
| Data Description:
|
| generated by r.in.ascii
|
|
|
|
|

+---------------------------------------------------------------------------
-+

-----Original Message-----
From: Markus Neteler [mailto:neteler@geog.uni-hannover.de]
Sent: Sunday, December 17, 2000 12:53 PM
To: grass5 developers list
Subject: Re: [GRASS5] nviz

Helena,

thanks for the hint. I have forwarded it to my student.
I feel he has a library violation somewhere.

On Sun, Dec 17, 2000 at 12:59:30PM -0600, Helena Mitasova - staff wrote:

>nviz2.2 -q el=hoehen_test

I am not sure whether that is a problem, but -q does not
expect anything after it - it is supposed to open
nviz without any data (and you can then load everything
from e.g. your state file saved from your previous work.
-q did not work in pre-beta10 versions, it was supposed
to be fixed, but maybe the fix is not 100% done.

so it should be either nviz2.2 -q
or nviz2.2 hoehen_test

>module: r.out.ascii
>
>bugreport: dp paramater doesn't seem to work. All output ascii files
>are integer regardless whether the source map is integer,
>float, or double.

I have used r.out.ascii with dp without problems in beta8 - has
anything changed that cold have messed it up in beta10?
(I haven't dowloaded beta10 yet, I hope to get to it next week.)

Well, it's somehow now working:

r.out.ascii m=dgm12.5 dp=7 |more
north: 5770690
south: 5768830
east: 3567625
west: 3565765
rows: 372
cols: 372
287.24 287.24 289.22 289.22 290.75 290.75 290.75 291.86 291.86 292.32 292.32
292.32 292.62 292.62 292.51 292.51 292.51 291.03

r.out.ascii m=dgm12.5 dp=1 |more
north: 5770690
south: 5768830
east: 3567625
west: 3565765
rows: 372
cols: 372
287.2 287.2 289.2 289.2 290.8 290.8

r.out.ascii m=dgm12.5 dp=3 |more
north: 5770690
south: 5768830
east: 3567625
west: 3565765
rows: 372
cols: 372
287.24 287.24 289.22 289.22

dgm12.5 is a DCELL map. However, the module wasn't changed. Very strange.

r.info dgm12.5

+---------------------------------------------------------------------------
-+
| Layer: dgm12.5 Date: Thu Oct 19 13:20:50 2000
| Type of Map: raster Number of Categories: 373
|
| Data Type: DCELL
|

Oh, I would be glad to have a fully working GRASS once...

Markus

----------------------------------------
If you want to unsubscribe from GRASS Development Team mailing list write
to:
minordomo@geog.uni-hannover.de with
subject 'unsubscribe grass5'

----------------------------------------
If you want to unsubscribe from GRASS Development Team mailing list write to:
minordomo@geog.uni-hannover.de with
subject 'unsubscribe grass5'

On Sun, Dec 17, 2000 at 01:04:29PM -0700, William L. Baker wrote:

Hello,

Here is the r.in.ascii displacement problem and the failure
of the dp paramater in routascii. There were no compile problems
for either module on my machine (intel with Mandrake 7.2). I
wonder why it fails on my machine but different on others...?

GRASS:/usr/local/grassedit > more rinascii.txt
north: 4504755
south: 4404275
east: 351345
west: 350985
rows: 4
cols: 3
0.0 1.0 0.0
1.0 1.0 1.0
0.0 1.0 1.0
0.0 0.0 1.0

GRASS:/usr/local/grassedit > r.in.ascii -f input=rinascii.txt output=test
CREATING SUPPORT FILES FOR test

Put some non-zero digits after those decimal points.

GRASS:/usr/local/grassedit > r.out.ascii map=test dp=7 |more
north: 4504755
south: 4504275
east: 351345
west: 350985
rows: 4
cols: 3
0 1 0
0 1 0
0 1 0
0 1 0

If all of the digits past the decimal point are zero, they get removed
in the ascii output (along with the decimal point). This is a
"feature" I guess (there's a lib function that does this stripping).
It looks like your output is not the same as the input though (was there
something else done?). Note that the header for rinascii.txt should have
a different "south" -- the north south res there is 25120 instead of 120
(like E-W).

--
Eric G. Miller <egm2@jps.net>

----------------------------------------
If you want to unsubscribe from GRASS Development Team mailing list write to:
minordomo@geog.uni-hannover.de with
subject 'unsubscribe grass5'