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

GRASS GIS trac at osgeo.org
Mon Jun 20 06:54:12 EDT 2011


#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>
GRASS GIS <http://grass.osgeo.org>



More information about the grass-dev mailing list