On Wed, May 23, 2018 at 10:05 PM, Francois Chartier <[email protected]> wrote: > 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,
You can view it in 3D view, see https://youtu.be/CmK-fEyf2SY or use r3.to.rast (https://grass.osgeo.org/grass74/manuals/r3.to.rast.html) to slice it into a series of 2D rasters and view them. > > > ----------------------------------------------------------------------------+ > | 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 > > _______________________________________________ grass-user mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/grass-user
