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@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@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@gmail.com <mailto:fra.chartier@gmail.com>> a écrit :
does anybody know how to resolve this?
On Thu, Nov 8, 2018, 06:23 Francois Chartier <fra.chartier@gmail.com
<mailto:fra.chartier@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@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@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^���
