[GRASS-user] v.vol.rst - DMIN and Segment Size

Francois Chartier fra.chartier at gmail.com
Wed Nov 14 18:06:32 PST 2018


So i rerun v.vol.rst with new g.region -p3

C:\>g.region -p3
projection: 1 (UTM)
zone:       17
datum:      nad83
ellipsoid:  grs80
north:      4862000
south:      4825700
west:       609200
east:       651700
top:        250.00000000
bottom:     50.00000000
nsres:      10
nsres3:     504.16666667
ewres:      10
ewres3:     500
tbres:      1
rows:       3630
rows3:      72
cols:       4250
cols3:      85
depths:     200
cells:      15427500
cells3:     1224000

V.VOL.RST run (below) with increasing DMIN, from 0.5 to 200. finally worked
with 200 m of min spacing. as you can see the majority of the points are
excluded.
I am using 180k points over a 40 km 40 km domain.
I have not set any smoothing parameters. also using the default tension
parameters and as data points are very close vertically (every m or less)
Thanks

(Wed Nov 14 20:50:41 2018)

v.vol.rst input=geotechNov8 at Toronto wcolumn=dbl_2 segmax=10 npmin=20
npmax=100 dmin=0.5 elevation=geotechnov8nov14
188692 records selected from table
WARNING: Some points outside of region -- will ignore...
ERROR: Points are too concentrated -- please increase DMIN
(Wed Nov 14 20:50:50 2018) Command finished (8 sec)

(Wed Nov 14 20:51:08 2018)

v.vol.rst input=geotechNov8 at Toronto wcolumn=dbl_2 segmax=10 npmin=20
npmax=100 dmin=5 elevation=geotechnov8nov14
188692 records selected from table
WARNING: Some points outside of region -- will ignore...
ERROR: Points are too concentrated -- please increase DMIN
(Wed Nov 14 20:51:17 2018) Command finished (8 sec)

(Wed Nov 14 20:51:25 2018)

v.vol.rst input=geotechNov8 at Toronto wcolumn=dbl_2 segmax=10 npmin=20
npmax=100 dmin=20 elevation=geotechnov8nov14
188692 records selected from table
WARNING: Some points outside of region -- will ignore...
ERROR: Points are too concentrated -- please increase DMIN
(Wed Nov 14 20:51:34 2018) Command finished (8 sec)

(Wed Nov 14 20:51:46 2018)

v.vol.rst input=geotechNov8 at Toronto wcolumn=dbl_2 segmax=10 npmin=20
npmax=100 dmin=50 elevation=geotechnov8nov14
188692 records selected from table
WARNING: Some points outside of region -- will ignore...
ERROR: Points are too concentrated -- please increase DMIN
(Wed Nov 14 20:51:59 2018) Command finished (12 sec)

(Wed Nov 14 20:52:11 2018)

v.vol.rst input=geotechNov8 at Toronto wcolumn=dbl_2 segmax=10 npmin=20
npmax=100 dmin=70 elevation=geotechnov8nov14
188692 records selected from table
WARNING: Some points outside of region -- will ignore...
ERROR: Points are too concentrated -- please increase DMIN
(Wed Nov 14 20:52:20 2018) Command finished (8 sec)

(Wed Nov 14 20:55:19 2018)

v.vol.rst input=geotechNov8 at Toronto wcolumn=dbl_2 segmax=10 npmin=20
npmax=100 dmin=100 elevation=geotechnov8nov14
188692 records selected from table
WARNING: Some points outside of region -- will ignore...
ERROR: Points are too concentrated -- please increase DMIN
(Wed Nov 14 20:55:28 2018) Command finished (8 sec)

(Wed Nov 14 20:55:46 2018)

v.vol.rst input=geotechNov8 at Toronto wcolumn=dbl_2 segmax=10 npmin=20
npmax=100 dmin=200 elevation=geotechnov8nov14
188692 records selected from table
WARNING: Some points outside of region -- will ignore...
Processing all selected output files will require 24480 bytes of disk space
for temp files
WARNING: There are points outside specified 2D/3D region--ignored 3533
points (total points: 188692)
WARNING: Points are more dense than specified 'DMIN'--ignored 175580 points
(remain 13112)
WARNING: Unable to create <(null)> raster map without cross_input raster
map being specified
Finished interpolating
3D raster map <geotechnov8nov14> created
The number of points in vector map is 188692
The number of points outside of 2D/3D region 3533
The number of points used (after reduction) is 9579
(Wed Nov 14 20:55:58 2018) Command finished (12 sec)

<https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail>
Virus-free.
www.avast.com
<https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail>
<#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>

Le mer. 14 nov. 2018 à 20:40, Francois Chartier <fra.chartier at gmail.com> a
écrit :

> good point.  i created the 3d raster before. see g.region -p3 below.
> I will rerun v.vol.rst with the higher 3d g.region resolution, and keep
> increasing the dmin until it runs.
>
> C:\>g.region -p3
> projection: 1 (UTM)
> zone:       17
> datum:      nad83
> ellipsoid:  grs80
> north:      4862000
> south:      4825700
> west:       609200
> east:       651700
> top:        250.00000000
> bottom:     50.00000000
> nsres:      10
> nsres3:     504.16666667
> ewres:      10
> ewres3:     500
> tbres:      1
> rows:       3630
> rows3:      72
> cols:       4250
> cols3:      85
> depths:     200
> cells:      15427500
> cells3:     1224000
>
>
> <https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail> Virus-free.
> www.avast.com
> <https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail>
> <#m_3960299417046487252_DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
>
> Le mer. 14 nov. 2018 à 09:38, Moritz Lennert <mlennert at club.worldonline.be>
> a écrit :
>
>> On 14/11/18 01:28, Francois Chartier wrote:
>> > Hi,
>> >
>> > v.vol.rst starts working when npmin is 100 m or more.  I dont believe
>> > that this algorythm only works with data points spaced by 100 m that
>> > does not make sense.  is there a max total number for a dataset? is it
>> > because i have too few rows and columns?
>> > see metadata below.  i also showed a lower dmin of 0.9.  g.region is
>> > printed at the end of the email
>>
>>
>> Your g.region output does not correspond to the region that was set when
>> you ran v.vol.rst as the map created has a resolution of 500 in x and y
>> direction, but your region says 10.
>>
>> Could you provide the output of g.region -p3 to also get the 3D settings ?
>>
>> I don't have an exact idea, but the map you created has a resolution of
>> 500 x 500 x 176, so maybe if you have a point per meter this just means
>> that the number of points you have per target pixel is much too high.
>>
>> Maybe you need to think about the appropriate 3D region first and make
>> sure to set that before running v.vol.rst ?
>>
>> Moritz
>>
>> >
>> >
>> >
>>  +----------------------------------------------------------------------------+
>> >                                              |
>> >
>>  |----------------------------------------------------------------------------|
>> >   |
>>
>> >        |
>> >   |   Type of Map:  raster_3d            Number of Categories: 0
>> >         |
>> >   |   Data Type:    FCELL
>>
>> >        |
>> >   |   Rows:         90
>> >         |
>> >   |   Columns:      85
>> >         |
>> >   |   Depths:       2
>>
>> >        |
>> >   |   Total Cells:  15300
>>
>> >        |
>> >   |   Total size:           48229 Bytes
>>
>> >        |
>> >   |   Number of tiles:      4
>>
>> >        |
>> >   |   Mean tile size:       12057 Bytes
>>
>> >        |
>> >   |   Tile size in memory:  15480 Bytes
>>
>> >        |
>> >   |   Number of tiles in x, y and  z:   2, 2, 1
>>
>> >        |
>> >   |   Dimension of a tile in x, y, z:   43, 45, 2
>>
>> >        |
>> >   |
>>
>> >        |
>> >   |        Projection: UTM (zone 17)
>> >         |
>> >   |            N:    4870782    S:    4825713   Res: 500.76666667
>>
>> >        |
>> >   |            E:     651697    W:     609201   Res: 499.95294118
>>
>> >        |
>> >   |            T:        352    B:          0   Res:   176
>> >         |
>> >   |   Range of data:   min = -1.17721057 max = 89.22615051
>> >         |
>> >   |
>>
>> >        |
>> >   |   Data Source:
>> >         |
>> >   |
>>
>> >        |
>> >   |
>>
>> >        |
>> >   |
>>
>> >        |
>> >   |   Data Description:
>>
>> >        |
>> >   |    generated by v.vol.rst
>>
>> >        |
>> >   |
>>
>> >        |
>> >   |   Comments:
>>
>> >        |
>> >   |    v.vol.rst --overwrite input="geotechNov8 at Toronto"
>> wcolumn="dbl_2"
>> > te\   |
>> >   |    nsion=40. smooth=0.1 segmax=50 npmin=200 npmax=700 dmin=100
>> > wscale=1\   |
>> >   |    .0 zscale=1.0 elevation="geotechnov8"
>> >         |
>> >   |
>>
>> >        |
>> >
>>  +----------------------------------------------------------------------------+
>> > (Tue Nov 13 19:18:12 2018) Command finished (0 sec)
>> >
>> > *i have also tried with low dmin of 0.9 m:*
>> >
>> >
>> > (Tue Nov 13 19:22:49 2018)
>> > v.vol.rst --overwrite input=geotechNov8 at Toronto wcolumn=dbl_2 npmin=60
>> > npmax=100 dmin=0.9 elevation=geotechnov13
>> > 188692 records selected from table
>> > WARNING: Some points outside of region -- will ignore...
>> > ERROR: Points are too concentrated -- please increase DMIN
>> > (Tue Nov 13 19:23:00 2018) Command finished (10 sec)
>> >
>> > METADATA of Vector points
>> >
>> > Thu Nov 08 06:06:56 2018                                  |
>> >   | Timestamp (first layer): none
>>
>> >        |
>> >
>>  |----------------------------------------------------------------------------|
>> >   | Map format:      native
>>
>> >        |
>> >
>>  |----------------------------------------------------------------------------|
>> >   |   Type of map: vector (level: 2)
>> >         |
>> >   |
>>
>> >        |
>> >   |   Number of points:       188692          Number of centroids:  0
>>
>> >        |
>> >   |   Number of lines:        0               Number of boundaries: 0
>>
>> >        |
>> >   |   Number of areas:        0               Number of islands:    0
>>
>> >        |
>> >   |   Number of faces:        0               Number of kernels:    0
>>
>> >        |
>> >   |   Number of volumes:      0               Number of holes:      0
>>
>> >        |
>> >   |
>>
>> >        |
>> >   |   Map is 3D:              Yes
>>
>> >        |
>> >   |   Number of dblinks:      1
>>
>> >        |
>> >   |
>>
>> >        |
>> >   |   Projection: UTM (zone 17)
>>
>> >        |
>> >   |
>>
>> >        |
>> >   |               N:           4870782    S:           4825713
>> >         |
>> >   |               E:            651697    W:            609201
>> >         |
>> >   |               B:                 0    T:               352
>> >         |
>> >   |
>>
>> >        |
>> >   |   Digitization threshold: 0
>>
>> >        |
>> >   |   Comment:
>> >         |
>> >   |
>>
>> >        |
>> >
>>  +----------------------------------------------------------------------------+
>> > (Tue Nov 13 19:25:30 2018) Command finished (0 sec)
>> >
>> > g.region - p
>> >
>> > projection: 1 (UTM)
>> > zone:       17
>> > datum:      nad83
>> > ellipsoid:  grs80
>> > north:      4862000
>> > south:      4825700
>> > west:       609200
>> > east:       651700
>> > nsres:      10
>> > ewres:      10
>> > rows:       3630
>> > cols:       4250
>> > cells:      15427500
>> >
>> > Le ven. 9 nov. 2018 à 07:16, Francois Chartier <fra.chartier at gmail.com
>> > <mailto:fra.chartier at gmail.com>> a écrit :
>> >
>> >     does anybody know how to resolve this?
>> >
>> >     On Thu, Nov 8, 2018, 06:23 Francois Chartier <
>> fra.chartier at gmail.com
>> >     <mailto:fra.chartier at gmail.com> wrote:
>> >
>> >         Hi,
>> >
>> >         I am trying to use V.VOL.RST to interpolate particle size for
>> >         soil samples.   i obtain the error message to increase DMin.
>> >         unfortunately i have a data point every meter vertically and a
>> >         total of 190K.  if grass is asking to increase DMIN, which i
>> >         have increased to 30 m to try to see if this was the pb, can I
>> >         resolve this by reducing the Segment size?
>> >
>> >         * COMMAND OUTPUT*
>> >         (Thu Nov 08 06:17:32 2018) Command finished (10 sec)
>> >         (Thu Nov 08 06:20:40 2018)
>> >         v.vol.rst input=geotechNov8 at Toronto wcolumn=dbl_2 npmin=50
>> >         npmax=200 dmin=30 elevation=Geotech3DNov8
>> >         188692 records selected from table
>> >         WARNING: Some points outside of region -- will ignore...
>> >         ERROR: Points are too concentrated -- please increase DMIN
>> >         (Thu Nov 08 06:20:51 2018) Command finished (10 sec)
>> >
>> >         *METADATA*
>> >         Thu Nov 08 06:11:32 2018) Command finished (0 sec)
>> >         (Thu Nov 08 06:19:17 2018)
>> >         v.info <http://v.info> map=geotechNov8 at Toronto
>> >
>>  +----------------------------------------------------------------------------+
>> >           | Name:            geotechNov8
>> >                         |
>> >           | Mapset:          Toronto
>> >                         |
>> >           | Location:        Toronto
>> >                         |
>> >           | Database:        C:\Users\Francois
>> >         Chartier\Documents\grassdata            |
>> >           | Title:
>> >                         |
>> >           | Map scale:       1:1
>> >                         |
>> >           | Name of creator: Francois Chartier
>> >                         |
>> >           | Organization:
>> >                        |
>> >           | Source date:     Thu Nov 08 06:06:56 2018
>> >                        |
>> >           | Timestamp (first layer): none
>> >                        |
>> >
>>  |----------------------------------------------------------------------------|
>> >           | Map format:      native
>> >                        |
>> >
>>  |----------------------------------------------------------------------------|
>> >           |   Type of map: vector (level: 2)
>> >                         |
>> >           |
>> >                        |
>> >           |   Number of points:       188692          Number of
>> >         centroids:  0          |
>> >           |   Number of lines:        0               Number of
>> >         boundaries: 0          |
>> >           |   Number of areas:        0               Number of
>> >         islands:    0          |
>> >           |   Number of faces:        0               Number of
>> >         kernels:    0          |
>> >           |   Number of volumes:      0               Number of holes:
>> >              0          |
>> >           |
>> >                        |
>> >           |   Map is 3D:              Yes
>> >                        |
>> >           |   Number of dblinks:      1
>> >                        |
>> >           |
>> >                        |
>> >           |   Projection: UTM (zone 17)
>> >                        |
>> >           |
>> >                        |
>> >           |               N:           4870782    S:           4825713
>> >                         |
>> >           |               E:            651697    W:            609201
>> >                         |
>> >           |               B:                 0    T:               352
>> >                         |
>> >           |
>> >                        |
>> >           |   Digitization threshold: 0
>> >                        |
>> >           |   Comment:
>> >                         |
>> >           |
>> >                        |
>> >
>>  +----------------------------------------------------------------------------+
>> >         (Thu Nov 08 06:19:18 2018) Command finished (0 sec)
>> >
>> >
>> > N �n�r����)em�h�yhiם�w^���
>> >
>>
>>
>>
>
> <https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail> Virus-free.
> www.avast.com
> <https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail>
> <#m_3960299417046487252_DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20181114/49b563cd/attachment-0001.html>


More information about the grass-user mailing list