[GRASS-user] v.vol.rst interpolation results in rectangular isosurfaces

[Markus Metz]

On Sat, Jan 26, 2019 at 7:43 PM Francois Chartier <fra.chartier at gmail.com>
wrote:
>
> Hi,
>
> Using v.vol.rst, (see command below, and g.region), the interpolation
results in rectangular areas in certain areas, see attached photo.
> Does this have to do with the segmentation process or something else?

Most probably yes, try to increase npmin. Maybe you also need to increase
segmax, but segmax should not be increased without increasing npmin,
otherwise you might get more of these block artefacts.

Markus M
>
>  v.vol.rst --overwrite input="Dec16CleanedJjan12 at Toronto" wcolumn="db\   |
>  |    l_2" tension=10 smooth=0.1 segmax=50 npmin=200 npmax=700 dmin=0.1
ws\   |
>  |    cale=1.0 zscale=1.0 elevation="jan25"
>
> (Sat Jan 26 13:37:21 2019) Command finished (0 sec)

> (Sat Jan 26 13:37:44 2019)

> r3.info map=Jan25 at Toronto

>
 +----------------------------------------------------------------------------+
>  | Layer:    Jan25 at Toronto                  Date: Fri Jan 25 22:12:59
2019    |
>  | Mapset:   Toronto                        Login of Creator: Francois
Charti |
>  | Location: Toronto
     |
>  | DataBase: C:\Users\Francois Chartier\Documents\grassdata
    |
>  | Title:    jan25
     |
>  | Units:    none
    |
>  | Vertical unit: units
    |
>  | Timestamp: none
     |
>
 |----------------------------------------------------------------------------|
>  |
     |
>  |   Type of Map:  raster_3d            Number of Categories: 0
    |
>  |   Data Type:    FCELL
     |
>  |   Rows:         400
     |
>  |   Columns:      450
     |
>  |   Depths:       65
    |
>  |   Total Cells:  11700000
    |
>  |   Total size:           27646661 Bytes
    |
>  |   Number of tiles:      1575
    |
>  |   Mean tile size:       17553 Bytes
     |
>  |   Tile size in memory:  32400 Bytes
     |
>  |   Number of tiles in x, y and  z:   15, 15, 7
     |
>  |   Dimension of a tile in x, y, z:   30, 27, 10
    |
>  |
     |
>  |        Projection: UTM (zone 17)
    |
>  |            N:    4837000    S:    4833000   Res:    10
    |
>  |            E:     632500    W:     628000   Res:    10
    |
>  |            T:        115    B:         50   Res:     1
    |
>  |   Range of data:   min = -25.42733002 max = 121.86579132
    |
>  |
     |
>  |   Data Source:
    |
>  |
     |
>  |
     |
>  |
     |
>  |   Data Description:
     |
>  |    generated by v.vol.rst
     |
>  |
     |
>  |   Comments:
     |
>  |    v.vol.rst --overwrite input="Dec16CleanedJjan12 at Toronto"
wcolumn="db\   |
>  |    l_2" tension=10 smooth=0.1 segmax=50 npmin=200 npmax=700 dmin=0.1
ws\   |
>  |    cale=1.0 zscale=1.0 elevation="jan25"
    |
>  |
     |
>
 +----------------------------------------------------------------------------+
> (Sat Jan 26 13:37:45 2019) Command finished (0 sec)

>
> C:\Users\Francois Chartier>g.region -p3
> projection: 1 (UTM)
> zone:       17
> datum:      nad83
> ellipsoid:  grs80
> north:      4837000
> south:      4833000
> west:       628000
> east:       632500
> top:        115.00000000
> bottom:     50.00000000
> nsres:      10
> nsres3:     10
> ewres:      10
> ewres3:     10
> tbres:      1
> rows:       400
> rows3:      400
> cols:       450
> cols3:      450
> depths:     65
> cells:      180000
> cells3:     11700000
>
> C:\Users\Francois Chartier>
>
> Virus-free. www.avast.com
> _______________________________________________
> grass-user mailing list
> grass-user at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/grass-user
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20190131/2e4a3929/attachment.html>