[GRASS-user] v.vol.rst - DMIN and Segment Size
Moritz Lennert
mlennert at club.worldonline.be
Wed Nov 14 05:14:49 PST 2018
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^���
>
More information about the grass-user
mailing list