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 <[email protected]> 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@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 <[email protected] > > <mailto:[email protected]>> a écrit : > > > > does anybody know how to resolve this? > > > > On Thu, Nov 8, 2018, 06:23 Francois Chartier <[email protected] > > <mailto:[email protected]> 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^��� > > > > > <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>
_______________________________________________ grass-user mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/grass-user
