<div dir="ltr"><div>Hi, </div><div>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?</div><div><br></div><div>----------------------------------------------------------------------------+</div><div> | Layer:    BHALL_May24v2@Toronto          Date: Sun May 27 16:43:46 2018    |</div><div> | Mapset:   Toronto                        Login of Creator: Francois Charti |</div><div> | Location: Toronto                                                          |</div><div> | DataBase: C:\Users\Francois Chartier\Documents\grassdata                   |</div><div> | Title:    BHALL_May24v2                                                    |</div><div> | Units:    none                                                             |</div><div> | Vertical unit: units                                                       |</div><div> | Timestamp: none                                                            |</div><div> |----------------------------------------------------------------------------|</div><div> |                                                                            |</div><div> |   Type of Map:  raster_3d            Number of Categories: 0               |</div><div> |   Data Type:    FCELL                                                      |</div><div> |   Rows:         90                                                         |</div><div> |   Columns:      85                                                         |</div><div> |   Depths:       307                                                        |</div><div> |   Total Cells:  2348550                                                    |</div><div> |   Total size:           3180868 Bytes                                      |</div><div> |   Number of tiles:      336                                                |</div><div> |   Mean tile size:       9466 Bytes                                         |</div><div> |   Tile size in memory:  30420 Bytes                                        |</div><div> |   Number of tiles in x, y and  z:   6, 7, 8                                |</div><div> |   Dimension of a tile in x, y, z:   15, 13, 39                             |</div><div> |                                                                            |</div><div> |        Projection: UTM (zone 17)                                           |</div><div> |            N:    4870795    S:    4825715   Res: 500.88888889              |</div><div> |            E:     651669    W:     609204   Res: 499.58823529              |</div><div> |            T:        352    B:         45   Res:     1                     |</div><div> |   Range of data:   min = -8.99099445 max = 33.89262009                     |</div><div> |                                                                            |</div><div> |   Data Source:                                                             |</div><div> |                                                                            |</div><div> |                                                                            |</div><div> |                                                                            |</div><div> |   Data Description:                                                        |</div><div> |    generated by v.vol.rst                                                  |</div><div> |                                                                            |</div><div> |   Comments:                                                                |</div><div> |    v.vol.rst input="BHMat1All_May12@Toronto" wcolumn="cat" tension=40. \   |</div><div> |    smooth=0.1 segmax=50 npmin=200 npmax=700 dmin=2 wscale=1.0 zscale=1.\   |</div><div> |    0 elevation="BHALL_May24v2"                                             |</div><div> |                                                                            |</div><div> +----------------------------------------------------------------------------+</div><div>(Sun May 27 16:57:19 2018) Command finished (0 sec)   </div></div><div class="gmail_extra"><br><div class="gmail_quote">2018-05-24 4:52 GMT-04:00 Moritz Lennert <span dir="ltr"><<a href="mailto:mlennert@club.worldonline.be" target="_blank">mlennert@club.worldonline.be</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">On 24/05/18 04:05, Francois Chartier wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Thank you,<br>
<br>
v.vol.rst input=BHMat1All_May12@Toronto wcolumn=cat dmin=5 elevation=may23trial<br>
<br>
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.<br>
</blockquote>
<br></span>
Watch out: according to the info below, your horizontal resolution is around 500, not 200.<br>
<br>
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.<br>
<br>
Moritz<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5">
<br>
<br>
------------------------------<wbr>------------------------------<wbr>----------------+<br>
  | Layer:    may23trial@Toronto             Date: Wed May 23 21:46:30 2018    |<br>
  | Mapset:   Toronto                        Login of Creator: Francois Charti |<br>
  | Location: Toronto                                                           |<br>
  | DataBase: C:\Users\Francois Chartier\Documents\grassdata                    |<br>
  | Title:    may23trial                                                        |<br>
  | Units:    none                                                              |<br>
  | Vertical unit: units                                                        |<br>
  | Timestamp: none                                                             |<br>
  |---------------------------<wbr>------------------------------<wbr>-------------------|<br>
  |                                                                             |<br>
  |   Type of Map:  raster_3d            Number of Categories: 0                |<br>
  |   Data Type:    FCELL                                                       |<br>
  |   Rows:         90                                                          |<br>
  |   Columns:      85                                                          |<br>
  |   Depths:       61                                                          |<br>
  |   Total Cells:  466650                                                      |<br>
  |   Total size:           853368 Bytes                                        |<br>
  |   Number of tiles:      64                                                  |<br>
  |   Mean tile size:       13333 Bytes                                         |<br>
  |   Tile size in memory:  32384 Bytes                                         |<br>
  |   Number of tiles in x, y and  z:   4, 4, 4                                 |<br>
  |   Dimension of a tile in x, y, z:   22, 23, 16                              |<br>
  |                                                                             |<br>
  |        Projection: UTM (zone 17)                                            |<br>
  |            N:    4870795    S:    4825715   Res: 500.88888889               |<br>
  |            E:     651669    W:     609204   Res: 499.58823529               |<br>
  |            T:        352    B:         45   Res: 5.03278689                 |<br>
  |   Range of data:   min = -0.18226054 max = 32.72385025                      |<br>
  |                                                                             |<br>
  |   Data Source:                                                              |<br>
  |                                                                             |<br>
  |                                                                             |<br>
  |                                                                             |<br>
  |   Data Description:                                                         |<br>
  |    generated by v.vol.rst                                                   |<br>
  |                                                                             |<br>
  |   Comments:                                                                 |<br>
  |    v.vol.rst input="BHMat1All_May12@Toronto<wbr>" wcolumn="cat" tension=40. \   |<br>
  |    smooth=0.1 segmax=50 npmin=200 npmax=700 dmin=5 wscale=1.0 zscale=1.\   |<br>
  |    0 elevation="may23trial"                                                 |<br>
  |                                                                             |<br>
  +---------------------------<wbr>------------------------------<wbr>-------------------+<br>
(Wed May 23 21:54:13 2018) Command finished (0 sec)<br>
<br>
<br></div></div>
2018-05-23 21:42 GMT-04:00 Anna Petrášová <<a href="mailto:kratochanna@gmail.com" target="_blank">kratochanna@gmail.com</a> <mailto:<a href="mailto:kratochanna@gmail.com" target="_blank">kratochanna@gmail.com</a>><wbr>>:<span class=""><br>
<br>
    On Wed, May 23, 2018 at 9:02 PM, Francois Chartier<br></span><span class="">
    <<a href="mailto:fra.chartier@gmail.com" target="_blank">fra.chartier@gmail.com</a> <mailto:<a href="mailto:fra.chartier@gmail.com" target="_blank">fra.chartier@gmail.com</a><wbr>>> wrote:<br>
    > How do you change the resolution to something coarser?<br>
<br>
    With g.region, res= is horizontal and tbres= is vertical resolution,<br>
    e.g:<br>
<br>
    g.region res=200 tbres=200 -p3<br>
<br>
     ><br>
     > 2018-05-23 15:30 GMT-04:00 Markus Metz<br></span>
    <<a href="mailto:markus.metz.giswork@gmail.com" target="_blank">markus.metz.giswork@gmail.com</a> <mailto:<a href="mailto:markus.metz.giswork@gmail.com" target="_blank">markus.metz.giswork@gm<wbr>ail.com</a>>>:<span class=""><br>
     >><br>
     >><br>
     >><br>
     >> On Wed, May 23, 2018 at 1:09 PM, Francois Chartier<br></span><div><div class="h5">
     >> <<a href="mailto:fra.chartier@gmail.com" target="_blank">fra.chartier@gmail.com</a> <mailto:<a href="mailto:fra.chartier@gmail.com" target="_blank">fra.chartier@gmail.com</a><wbr>>> wrote:<br>
     >> ><br>
     >> > see below<br>
     >> > g.region -3p<br>
     >> ><br>
     >> > projection: 1 (UTM)<br>
     >> > zone:       17<br>
     >> > datum:      nad83<br>
     >> > ellipsoid:  grs80<br>
     >> > north:      4870795<br>
     >> > south:      4825715<br>
     >> > west:       609204<br>
     >> > east:       651669<br>
     >> > top:        352.00000000<br>
     >> > bottom:     45.00000000<br>
     >> > nsres:      5<br>
     >> > nsres3:     5<br>
     >> > ewres:      5<br>
     >> > ewres3:     5<br>
     >> > tbres:      1<br>
     >> > rows:       9016<br>
     >> > rows3:      9016<br>
     >> > cols:       8493<br>
     >> > cols3:      8493<br>
     >> > depths:     307<br>
     >> > cells:      76572888<br>
     >> > cells3:     23507876616<br>
     >> ><br>
     >> > for v.vol.rst command line - i am not sure where you get this<br>
    so i just<br>
     >> > restarted the simulation.<br>
     >> > Wed May 23 07:07:59 2018)<br>
     >> > v.vol.rst input=BHMat1All_May12@Toronto wcolumn=cat dmin=0.1<br>
     >> > 25 records selected from table<br>
     >> > Processing all selected output files will require 0 bytes of<br>
    disk space<br>
     >> > for temp files<br>
     >> > WARNING: Points are more dense than specified 'DMIN'--ignored<br>
    65 points<br>
     >> > (remain 76822)<br>
     >> > WARNING: Unable to create <(null)> raster map without<br>
    cross_input raster<br>
     >> > map being specified<br>
     >> ><br>
     >> ><br>
     >> ><br></div></div>
     >> > <a href="http://v.info" rel="noreferrer" target="_blank">v.info</a> <<a href="http://v.info" rel="noreferrer" target="_blank">http://v.info</a>><span class=""><br>
     >> > (Wed May 23 07:06:24 2018)<br></span>
     >> > <a href="http://v.info" rel="noreferrer" target="_blank">v.info</a> <<a href="http://v.info" rel="noreferrer" target="_blank">http://v.info</a>> map=BHMat1All_May12@Toronto<div><div class="h5"><br>
     >> ><br>
     >> ><br>
    +-----------------------------<wbr>------------------------------<wbr>-----------------+<br>
     >> >  | Name:            BHMat1All_May12<br>
     >> > |<br>
     >> >  | Mapset:          Toronto<br>
     >> > |<br>
     >> >  | Location:        Toronto<br>
     >> > |<br>
     >> >  | Database:        C:\Users\Francois Chartier\Documents\grassdata<br>
     >> > |<br>
     >> >  | Title:<br>
     >> > |<br>
     >> >  | Map scale:       1:1<br>
     >> > |<br>
     >> >  | Name of creator: Francois Chartier<br>
     >> > |<br>
     >> >  | Organization:<br>
     >> > |<br>
     >> >  | Source date:     Sat May 12 15:45:44 2018<br>
     >> > |<br>
     >> >  | Timestamp (first layer): none<br>
     >> > |<br>
     >> ><br>
     >> ><br>
    |-----------------------------<wbr>------------------------------<wbr>-----------------|<br>
     >> >  | Map format:      native<br>
     >> > |<br>
     >> ><br>
     >> ><br>
    |-----------------------------<wbr>------------------------------<wbr>-----------------|<br>
     >> >  |   Type of map: vector (level: 2)<br>
     >> > |<br>
     >> >  |<br>
     >> > |<br>
     >> >  |   Number of points:       76887           Number of<br>
    centroids:  0<br>
     >> > |<br>
     >> >  |   Number of lines:        0               Number of<br>
    boundaries: 0<br>
     >> > |<br>
     >> >  |   Number of areas:        0               Number of<br>
    islands:    0<br>
     >> > |<br>
     >> >  |   Number of faces:        0               Number of<br>
    kernels:    0<br>
     >> > |<br>
     >> >  |   Number of volumes:      0               Number of holes:          0<br>
     >> > |<br>
     >> >  |<br>
     >> > |<br>
     >> >  |   Map is 3D:              Yes<br>
     >> > |<br>
     >> >  |   Number of dblinks:      1<br>
     >> > |<br>
     >> >  |<br>
     >> > |<br>
     >> >  |   Projection: UTM (zone 17)<br>
     >> > |<br>
     >> >  |<br>
     >> > |<br>
     >> >  |               N:           4870795    S:           4825715<br>
     >> > |<br>
     >> >  |               E:            651669    W:            609204<br>
     >> > |<br>
     >> >  |               B:                45    T:               352<br>
     >> > |<br>
     >> >  |<br>
     >> > |<br>
     >> >  |   Digitization threshold: 0<br>
     >> > |<br>
     >> >  |   Comment:<br>
     >> > |<br>
     >> >  |<br>
     >> > |<br>
     >> ><br>
     >> ><br>
    +-----------------------------<wbr>------------------------------<wbr>-----------------+<br>
     >> > (Wed May 23 07:06:24 2018) Command finished (0 sec)<br>
     >><br>
     >> The current target resolutions of<br>
     >><br>
     >> nsres3:     5<br>
     >> ewres3:     5<br>
     >> tbres:      1<br>
     >><br>
     >> are much too high for the given extents and the coverage of the<br>
    input<br>
     >> points<br>
     >><br></div></div>
     >> According to the output of <a href="http://v.info" rel="noreferrer" target="_blank">v.info</a> <<a href="http://v.info" rel="noreferrer" target="_blank">http://v.info</a>>, a point<span class=""><br>
    covers on average about 7643645<br>
     >> cubic meters which corresponds to voxels with edge length of<br>
    about 200<br>
     >> meters. Therefore a good starting point would be a 3D region with<br>
     >> resolutions in all 3 dimensions of 200 meter. If the results<br>
    make sense, the<br>
     >> resolutions can be slowly increased.<br>
     >><br>
     >> Markus M<br>
     >> ><br>
     >> ><br>
     >> > 2018-05-23 6:07 GMT-04:00 Moritz Lennert<br></span>
    <<a href="mailto:mlennert@club.worldonline.be" target="_blank">mlennert@club.worldonline.be</a> <mailto:<a href="mailto:mlennert@club.worldonline.be" target="_blank">mlennert@club.worldonl<wbr>ine.be</a>>>:<span class=""><br>
     >> >><br>
     >> >> Hi François,<br>
     >> >><br>
     >> >> Please keep all discussions on the list.<br>
     >> >><br>
     >> >> On 23/05/18 02:26, Francois Chartier wrote:<br>
     >> >>><br>
     >> >>> Hi Moritz,<br>
     >> >>><br>
     >> >>> see below.  seems like i have 6M cells<br>
     >> >>><br>
     >> >>> projection: 1 (UTM)<br>
     >> >>> zone:       17<br>
     >> >>> datum:      nad83<br>
     >> >>> ellipsoid:  grs80<br>
     >> >>> north:      4868980<br>
     >> >>> south:      4829680<br>
     >> >>> west:       609620<br>
     >> >>> east:       651628<br>
     >> >>> nsres:      5<br>
     >> >>> ewres:      4.99976196<br>
     >> >>> rows:       7860<br>
     >> >>> cols:       8402<br>
     >> >>> cells:      66039720<br>
     >> >><br>
     >> >><br>
     >> >><br>
     >> >> Sorry, my fault: as you are interpolating in 3D we need the<br>
    3D region,<br>
     >> >> so the output of 'g.region -3p'.<br>
     >> >><br>
     >> >> And could you also provide us with:<br>
     >> >><br>
     >> >> - the exact v.vol.rst command line you used<br></span>
     >> >> - the output of <a href="http://v.info" rel="noreferrer" target="_blank">v.info</a> <<a href="http://v.info" rel="noreferrer" target="_blank">http://v.info</a>> on your point dataset<span class=""><br>
     >> >><br>
     >> >> ?<br>
     >> >><br>
     >> >> Moritz<br>
     >> >><br>
     >> >>><br>
     >> >>> 2018-05-22 7:58 GMT-04:00 Moritz Lennert<br>
    <<a href="mailto:mlennert@club.worldonline.be" target="_blank">mlennert@club.worldonline.be</a> <mailto:<a href="mailto:mlennert@club.worldonline.be" target="_blank">mlennert@club.worldonl<wbr>ine.be</a>><br></span>
     >> >>> <mailto:<a href="mailto:mlennert@club.worldonline.be" target="_blank">mlennert@club.worldonl<wbr>ine.be</a><span class=""><br>
    <mailto:<a href="mailto:mlennert@club.worldonline.be" target="_blank">mlennert@club.worldonl<wbr>ine.be</a>>>>:<br>
     >> >>><br>
     >> >>><br>
     >> >>><br>
     >> >>>     Am 22. Mai 2018 04:35:13 MESZ schrieb Francois Chartier<br>
     >> >>>     <<a href="mailto:fra.chartier@gmail.com" target="_blank">fra.chartier@gmail.com</a> <mailto:<a href="mailto:fra.chartier@gmail.com" target="_blank">fra.chartier@gmail.com</a><wbr>><br></span>
    <mailto:<a href="mailto:fra.chartier@gmail.com" target="_blank">fra.chartier@gmail.com</a> <mailto:<a href="mailto:fra.chartier@gmail.com" target="_blank">fra.chartier@gmail.com</a><wbr>>>>:<span class=""><br>
     >> >>>      >Hi,<br>
     >> >>>      ><br>
     >> >>>      >I am interpolating 2D data using V Vol RST.  After 90<br>
    min the<br>
     >> >>> process<br>
     >> >>>      >progressed 15%, and I stopped it.<br>
     >> >>>      >There are approx 370 points, and it should not take<br>
    this long.<br>
     >> >>> I am<br>
     >> >>>      >not<br>
     >> >>>      >sure what I am doing incorrectly.  I have set the<br>
    g.region to<br>
     >> >>>     match the<br>
     >> >>>      >vector, and joined the vector attribute to the database.<br>
     >> >>><br>
     >> >>>     What is the resolution of the region ?<br>
     >> >>><br>
     >> >>>     Please provide the output of g.region -p.<br>
     >> >>><br>
     >> >>>     Moritz<br>
     >> >>><br>
     >> >>><br>
     >> >><br>
     >> >><br>
     >> ><br>
     >> ><br>
     >> > ______________________________<wbr>_________________<br>
     >> > grass-user mailing list<br></span>
     >> > <a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a> <mailto:<a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo<wbr>.org</a>><br>
     >> > <a href="https://lists.osgeo.org/mailman/listinfo/grass-user" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailma<wbr>n/listinfo/grass-user</a><span class=""><br>
    <<a href="https://lists.osgeo.org/mailman/listinfo/grass-user" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailm<wbr>an/listinfo/grass-user</a>><br>
     ><br>
     ><br>
     ><br>
     > ______________________________<wbr>_________________<br>
     > grass-user mailing list<br></span>
     > <a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a> <mailto:<a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo<wbr>.org</a>><br>
     > <a href="https://lists.osgeo.org/mailman/listinfo/grass-user" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailma<wbr>n/listinfo/grass-user</a><span class=""><br>
    <<a href="https://lists.osgeo.org/mailman/listinfo/grass-user" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailm<wbr>an/listinfo/grass-user</a>><br>
<br>
<br>
<br>
<br>
______________________________<wbr>_________________<br>
grass-user mailing list<br>
<a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/grass-user" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailma<wbr>n/listinfo/grass-user</a><br>
<br>
</span></blockquote>
<br>
<br>
</blockquote></div><br></div>