<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>