[GRASS-dev] [GRASS GIS] #1392: r3.to.rast converts FCELL map data to 0

GRASS GRASS developer
1

#1392: r3.to.rast converts FCELL map data to 0
------------------------+---------------------------------------------------
Reporter: hamish | Owner: grass-dev@…
     Type: defect | Status: new
Priority: normal | Milestone: 6.4.2
Component: Raster | Version: svn-develbranch6
Keywords: r3.to.rast | Platform: Linux
      Cpu: x86-64 |
------------------------+---------------------------------------------------
tested in 6.5svn.

man page test case:
{{{
cat << EOF > test.g3d
north: 4925010.000000
south: 4924890.000000
east: 596760.000000
west: 596610.000000
top: 1.000000
bottom: 0.000000
rows: 4
cols: 5
levels: 1
1204.74 1204.48 1204.19 1203.81 1203.39
1203.89 1203.67 1203.34 1202.98 1202.43
1203.05 1202.80 1202.51 1202.11 1201.48
1202.10 1201.92 1201.62 1201.27 1200.68
EOF

r3.in.ascii in=test.g3d out=test3d.f type=float

r3.info test3d.f
  | Type of Map: 3d cell Number of Categories: 0
|
  | Data Type: float
|
  | Rows: 4
|
  | Columns: 5
|
  | Depths: 1
|
  | Total Cells: 20
|
  | Projection: UTM (zone 13)
|
  | N: 4925010 S: 4924890 Res: 30
|
  | E: 596760 W: 596610 Res: 30
|
  | T: 1 B: 0 Res: 1
|
  | Range of data: min = 1200.68 max = 1204.74
|

g.region rast3d=test3d.f
r3.to.rast test3d.f out=test3d.f
r.info -r test3d.f_00001
min=0
max=0
}}}

but double works:
{{{
r3.in.ascii in=test.g3d out=test3d.d type=double
r3.to.rast test3d.d out=test3d.d
      #try it twice: r3.to.rast has --overwrite set by default?!

r.info -r test3d.d_00001
min=1200.68
max=1204.74
}}}

Hamish

--
Ticket URL: <https://trac.osgeo.org/grass/ticket/1392&gt;
GRASS GIS <http://grass.osgeo.org>

2

#1392: r3.to.rast converts FCELL map data to 0
------------------------+---------------------------------------------------
Reporter: hamish | Owner: grass-dev@…
     Type: defect | Status: new
Priority: normal | Milestone: 6.4.2
Component: Raster3D | Version: svn-develbranch6
Keywords: r3.to.rast | Platform: Linux
      Cpu: x86-64 |
------------------------+---------------------------------------------------
Changes (by hamish):

  * component: Raster => Raster3D

--
Ticket URL: <https://trac.osgeo.org/grass/ticket/1392#comment:1&gt;
GRASS GIS <http://grass.osgeo.org>

3

#1392: r3.to.rast converts FCELL map data to 0
------------------------+---------------------------------------------------
Reporter: hamish | Owner: grass-dev@…
     Type: defect | Status: new
Priority: normal | Milestone: 6.4.2
Component: Raster3D | Version: svn-develbranch6
Keywords: r3.to.rast | Platform: Linux
      Cpu: x86-64 |
------------------------+---------------------------------------------------

Comment(by hamish):

note

  * both FCELL and DCELL 3D rasters have "Number of Categories: 0" in
r3.info.

  * r.info for both F,D results of r3.to.rast report 255 cats.

??,
H

--
Ticket URL: <https://trac.osgeo.org/grass/ticket/1392#comment:2&gt;
GRASS GIS <http://grass.osgeo.org>

4

#1392: r3.to.rast converts FCELL map data to 0
-------------------------+--------------------------------------------------
  Reporter: hamish | Owner: grass-dev@…
      Type: defect | Status: closed
  Priority: normal | Milestone: 6.4.2
Component: Raster3D | Version: svn-develbranch6
Resolution: worksforme | Keywords: r3.to.rast
  Platform: Linux | Cpu: x86-64
-------------------------+--------------------------------------------------
Changes (by huhabla):

* cc: huhabla (added)
  * status: new => closed
  * resolution: => worksforme

Comment:

The FCELL bug should be fixed in trunk r46745. Additional float test is
implemented. In case this works on other machine it should backported to
6.4 and 6.5.

The number of categories is in r46745 equal between raster and raster3d ->
0.

Best
Soeren

--
Ticket URL: <https://trac.osgeo.org/grass/ticket/1392#comment:3&gt;
GRASS GIS <http://grass.osgeo.org>