Hi, So now i am able to run the v.vol.rst (thank you). the 3d raster exists (see below) and shows up in the layer manager. However, i cannot show in the Map Display Window. I have assigned a color table. any ideas how to show the 3D raster with color coding based on the properties?
----------------------------------------------------------------------------+ | Layer: BHALL_May24v2@Toronto Date: Sun May 27 16:43:46 2018 | | Mapset: Toronto Login of Creator: Francois Charti | | Location: Toronto | | DataBase: C:\Users\Francois Chartier\Documents\grassdata | | Title: BHALL_May24v2 | | Units: none | | Vertical unit: units | | Timestamp: none | |----------------------------------------------------------------------------| | | | Type of Map: raster_3d Number of Categories: 0 | | Data Type: FCELL | | Rows: 90 | | Columns: 85 | | Depths: 307 | | Total Cells: 2348550 | | Total size: 3180868 Bytes | | Number of tiles: 336 | | Mean tile size: 9466 Bytes | | Tile size in memory: 30420 Bytes | | Number of tiles in x, y and z: 6, 7, 8 | | Dimension of a tile in x, y, z: 15, 13, 39 | | | | Projection: UTM (zone 17) | | N: 4870795 S: 4825715 Res: 500.88888889 | | E: 651669 W: 609204 Res: 499.58823529 | | T: 352 B: 45 Res: 1 | | Range of data: min = -8.99099445 max = 33.89262009 | | | | 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=2 wscale=1.0 zscale=1.\ | | 0 elevation="BHALL_May24v2" | | | +----------------------------------------------------------------------------+ (Sun May 27 16:57:19 2018) Command finished (0 sec) 2018-05-24 4:52 GMT-04:00 Moritz Lennert <[email protected]>: > On 24/05/18 04:05, Francois Chartier 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. >> > > Watch out: according to the info below, your horizontal resolution is > around 500, not 200. > > FYI, if you want to define your region's extension with a vector, but also > want a fixed resolution (here 200 and 5), you can use the -a flag of > g.region to force g.region to use this precise resolution instead of > defining a resolution that is an exact divisor of the vector's extension. > > Moritz > > >> >> ------------------------------------------------------------ >> ----------------+ >> | 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] <mailto: >> [email protected]>>: >> >> On Wed, May 23, 2018 at 9:02 PM, Francois Chartier >> <[email protected] <mailto:[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] <mailto:[email protected] >> >>: >> >> >> >> >> >> >> >> On Wed, May 23, 2018 at 1:09 PM, Francois Chartier >> >> <[email protected] <mailto:[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 <http://v.info> >> >> > (Wed May 23 07:06:24 2018) >> >> > v.info <http://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 <http://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] <mailto:[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 <http://v.info> on your point dataset >> >> >> >> >> >> ? >> >> >> >> >> >> Moritz >> >> >> >> >> >>> >> >> >>> 2018-05-22 7:58 GMT-04:00 Moritz Lennert >> <[email protected] <mailto:[email protected]> >> >> >>> <mailto:[email protected] >> <mailto:[email protected]>>>: >> >> >>> >> >> >>> >> >> >>> >> >> >>> Am 22. Mai 2018 04:35:13 MESZ schrieb Francois Chartier >> >> >>> <[email protected] <mailto:[email protected]> >> <mailto:[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] <mailto:[email protected]> >> >> > https://lists.osgeo.org/mailman/listinfo/grass-user >> <https://lists.osgeo.org/mailman/listinfo/grass-user> >> > >> > >> > >> > _______________________________________________ >> > grass-user mailing list >> > [email protected] <mailto:[email protected]> >> > https://lists.osgeo.org/mailman/listinfo/grass-user >> <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
