r3.to.rast converts FCELL map data to 0

[hamish]

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

[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

[huhabla]

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