[GRASS-user] v.vol.rst - DMIN and Segment Size
Francois Chartier
fra.chartier at gmail.com
Wed Nov 14 17:40:00 PST 2018
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>
<#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>
<#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20181114/beccb98f/attachment-0001.html>
More information about the grass-user
mailing list