[GRASS-dev] Re: [GRASS GIS] #802: v.vol.rst reverts mask coordinates

[GRASS GIS]

#802: v.vol.rst reverts mask coordinates
------------------------------------+---------------------------------------
 Reporter:  helena                  |       Owner:  grass-dev@…              
     Type:  defect                  |      Status:  new                      
 Priority:  normal                  |   Milestone:  6.5.0                    
Component:  Default                 |     Version:  svn-develbranch6         
 Keywords:  mask, 3D interpolation  |    Platform:  All                      
      Cpu:  x86-64                  |  
------------------------------------+---------------------------------------
Changes (by helena):

  * platform:  Linux => All


Comment:

 Both 2D and 3D masks have problems.[[BR]]

 v.vol.rst has option maskmap to define a 2D raster as a mask for volume.
 It appears that when output is a 2D raster (e.g. example from grassbook
 for precipitation with influence of topography) the result is correct.
 However, for the 3D raster output the result is not written correctly - it
 looks like wrong order of rows or columns discussed by Soeren recently.
 So if v.vol.rst is run with the mask shown in pink in the linked image,
 the result when converted to 2D series is like the yellow strip and the
 isosurfaces are shown in the second image.[[BR]]

 [http://skagit.meas.ncsu.edu/~helena/grasswork/volrstmask_bug2dX.jpg Mask
 and result in 2D][[BR]]
 [http://skagit.meas.ncsu.edu/~helena/grasswork/volrstmask_bug.tif Masked
 isosurfaces][[BR]]

 if a 2d MASK or 3d G3D_MASK is set, v.vol.rst finishes with segfault
 before writing the output. It finishes without segfault in GRASS7 but the
 mask is not applied.

 Helena

-- 
Ticket URL: <http://trac.osgeo.org/grass/ticket/802#comment:2>
GRASS GIS <http://grass.osgeo.org>