[GRASS-user] v.vol.rst - DMIN and Segment Size
Francois Chartier
fra.chartier at gmail.com
Tue Nov 13 16:28:56 PST 2018
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
+----------------------------------------------------------------------------+
|
|----------------------------------------------------------------------------|
|
|
| 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> a
écrit :
> does anybody know how to resolve this?
>
> On Thu, Nov 8, 2018, 06:23 Francois Chartier <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 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)
>>
>>
>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20181113/7a5604b9/attachment-0001.html>
More information about the grass-user
mailing list