Thank you, v.vol.rst input=BHMat1All_May12@Toronto wcolumn=cat dmin=5 elevation=may23trial
it would appear to have worked. I changed the resolution to 200 m horizontally and 5 m for z using g.region as a trial. now the big question is how do you see the data. I am able to view a 2D raster by changing the color scale but 3d raster is blank, ----------------------------------------------------------------------------+ | Layer: may23trial@Toronto Date: Wed May 23 21:46:30 2018 | | Mapset: Toronto Login of Creator: Francois Charti | | Location: Toronto | | DataBase: C:\Users\Francois Chartier\Documents\grassdata | | Title: may23trial | | Units: none | | Vertical unit: units | | Timestamp: none | |----------------------------------------------------------------------------| | | | Type of Map: raster_3d Number of Categories: 0 | | Data Type: FCELL | | Rows: 90 | | Columns: 85 | | Depths: 61 | | Total Cells: 466650 | | Total size: 853368 Bytes | | Number of tiles: 64 | | Mean tile size: 13333 Bytes | | Tile size in memory: 32384 Bytes | | Number of tiles in x, y and z: 4, 4, 4 | | Dimension of a tile in x, y, z: 22, 23, 16 | | | | Projection: UTM (zone 17) | | N: 4870795 S: 4825715 Res: 500.88888889 | | E: 651669 W: 609204 Res: 499.58823529 | | T: 352 B: 45 Res: 5.03278689 | | Range of data: min = -0.18226054 max = 32.72385025 | | | | Data Source: | | | | | | | | Data Description: | | generated by v.vol.rst | | | | Comments: | | v.vol.rst input="BHMat1All_May12@Toronto" wcolumn="cat" tension=40. \ | | smooth=0.1 segmax=50 npmin=200 npmax=700 dmin=5 wscale=1.0 zscale=1.\ | | 0 elevation="may23trial" | | | +----------------------------------------------------------------------------+ (Wed May 23 21:54:13 2018) Command finished (0 sec) 2018-05-23 21:42 GMT-04:00 Anna Petrášová <[email protected]>: > On Wed, May 23, 2018 at 9:02 PM, Francois Chartier > <[email protected]> wrote: > > How do you change the resolution to something coarser? > > With g.region, res= is horizontal and tbres= is vertical resolution, e.g: > > g.region res=200 tbres=200 -p3 > > > > > 2018-05-23 15:30 GMT-04:00 Markus Metz <[email protected]>: > >> > >> > >> > >> On Wed, May 23, 2018 at 1:09 PM, Francois Chartier > >> <[email protected]> wrote: > >> > > >> > see below > >> > g.region -3p > >> > > >> > projection: 1 (UTM) > >> > zone: 17 > >> > datum: nad83 > >> > ellipsoid: grs80 > >> > north: 4870795 > >> > south: 4825715 > >> > west: 609204 > >> > east: 651669 > >> > top: 352.00000000 > >> > bottom: 45.00000000 > >> > nsres: 5 > >> > nsres3: 5 > >> > ewres: 5 > >> > ewres3: 5 > >> > tbres: 1 > >> > rows: 9016 > >> > rows3: 9016 > >> > cols: 8493 > >> > cols3: 8493 > >> > depths: 307 > >> > cells: 76572888 > >> > cells3: 23507876616 > >> > > >> > for v.vol.rst command line - i am not sure where you get this so i > just > >> > restarted the simulation. > >> > Wed May 23 07:07:59 2018) > >> > v.vol.rst input=BHMat1All_May12@Toronto wcolumn=cat dmin=0.1 > >> > 25 records selected from table > >> > Processing all selected output files will require 0 bytes of disk > space > >> > for temp files > >> > WARNING: Points are more dense than specified 'DMIN'--ignored 65 > points > >> > (remain 76822) > >> > WARNING: Unable to create <(null)> raster map without cross_input > raster > >> > map being specified > >> > > >> > > >> > > >> > v.info > >> > (Wed May 23 07:06:24 2018) > >> > v.info map=BHMat1All_May12@Toronto > >> > > >> > +----------------------------------------------------------- > -----------------+ > >> > | Name: BHMat1All_May12 > >> > | > >> > | Mapset: Toronto > >> > | > >> > | Location: Toronto > >> > | > >> > | Database: C:\Users\Francois Chartier\Documents\grassdata > >> > | > >> > | Title: > >> > | > >> > | Map scale: 1:1 > >> > | > >> > | Name of creator: Francois Chartier > >> > | > >> > | Organization: > >> > | > >> > | Source date: Sat May 12 15:45:44 2018 > >> > | > >> > | Timestamp (first layer): none > >> > | > >> > > >> > |----------------------------------------------------------- > -----------------| > >> > | Map format: native > >> > | > >> > > >> > |----------------------------------------------------------- > -----------------| > >> > | Type of map: vector (level: 2) > >> > | > >> > | > >> > | > >> > | Number of points: 76887 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: 4870795 S: 4825715 > >> > | > >> > | E: 651669 W: 609204 > >> > | > >> > | B: 45 T: 352 > >> > | > >> > | > >> > | > >> > | Digitization threshold: 0 > >> > | > >> > | Comment: > >> > | > >> > | > >> > | > >> > > >> > +----------------------------------------------------------- > -----------------+ > >> > (Wed May 23 07:06:24 2018) Command finished (0 sec) > >> > >> The current target resolutions of > >> > >> nsres3: 5 > >> ewres3: 5 > >> tbres: 1 > >> > >> are much too high for the given extents and the coverage of the input > >> points > >> > >> According to the output of v.info, a point covers on average about > 7643645 > >> cubic meters which corresponds to voxels with edge length of about 200 > >> meters. Therefore a good starting point would be a 3D region with > >> resolutions in all 3 dimensions of 200 meter. If the results make > sense, the > >> resolutions can be slowly increased. > >> > >> Markus M > >> > > >> > > >> > 2018-05-23 6:07 GMT-04:00 Moritz Lennert < > [email protected]>: > >> >> > >> >> Hi François, > >> >> > >> >> Please keep all discussions on the list. > >> >> > >> >> On 23/05/18 02:26, Francois Chartier wrote: > >> >>> > >> >>> Hi Moritz, > >> >>> > >> >>> see below. seems like i have 6M cells > >> >>> > >> >>> projection: 1 (UTM) > >> >>> zone: 17 > >> >>> datum: nad83 > >> >>> ellipsoid: grs80 > >> >>> north: 4868980 > >> >>> south: 4829680 > >> >>> west: 609620 > >> >>> east: 651628 > >> >>> nsres: 5 > >> >>> ewres: 4.99976196 > >> >>> rows: 7860 > >> >>> cols: 8402 > >> >>> cells: 66039720 > >> >> > >> >> > >> >> > >> >> Sorry, my fault: as you are interpolating in 3D we need the 3D > region, > >> >> so the output of 'g.region -3p'. > >> >> > >> >> And could you also provide us with: > >> >> > >> >> - the exact v.vol.rst command line you used > >> >> - the output of v.info on your point dataset > >> >> > >> >> ? > >> >> > >> >> Moritz > >> >> > >> >>> > >> >>> 2018-05-22 7:58 GMT-04:00 Moritz Lennert < > [email protected] > >> >>> <mailto:[email protected]>>: > >> >>> > >> >>> > >> >>> > >> >>> Am 22. Mai 2018 04:35:13 MESZ schrieb Francois Chartier > >> >>> <[email protected] <mailto:[email protected]>>: > >> >>> >Hi, > >> >>> > > >> >>> >I am interpolating 2D data using V Vol RST. After 90 min the > >> >>> process > >> >>> >progressed 15%, and I stopped it. > >> >>> >There are approx 370 points, and it should not take this long. > >> >>> I am > >> >>> >not > >> >>> >sure what I am doing incorrectly. I have set the g.region to > >> >>> match the > >> >>> >vector, and joined the vector attribute to the database. > >> >>> > >> >>> What is the resolution of the region ? > >> >>> > >> >>> Please provide the output of g.region -p. > >> >>> > >> >>> Moritz > >> >>> > >> >>> > >> >> > >> >> > >> > > >> > > >> > _______________________________________________ > >> > grass-user mailing list > >> > [email protected] > >> > https://lists.osgeo.org/mailman/listinfo/grass-user > > > > > > > > _______________________________________________ > > grass-user mailing list > > [email protected] > > https://lists.osgeo.org/mailman/listinfo/grass-user >
(Wed May 23 21:46:14 2018) v.vol.rst input=BHMat1All_May12@Toronto wcolumn=cat dmin=0.5 elevation=may23trial 25 records selected from table ERROR: Points are too concentrated -- please increase DMIN (Wed May 23 21:46:15 2018) Command finished (0 sec) (Wed May 23 21:46:29 2018) v.vol.rst input=BHMat1All_May12@Toronto wcolumn=cat dmin=5 elevation=may23trial 25 records selected from table WARNING: Some points outside of region -- will ignore... Processing all selected output files will require 30600 bytes of disk space for temp files WARNING: There are points outside specified 2D/3D region--ignored 19 points (total points: 76887) WARNING: Points are more dense than specified 'DMIN'--ignored 60377 points (remain 16510) WARNING: Unable to create <(null)> raster map without cross_input raster map being specified Finished interpolating 3D raster map <may23trial> created The number of points in vector map is 76887 The number of points outside of 2D/3D region 19 The number of points used (after reduction) is 16491 (Wed May 23 21:47:38 2018) Command finished (1 min 8 sec)
_______________________________________________ grass-user mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/grass-user
