<div dir="ltr"><div>Thank you, </div><div><br></div><div>v.vol.rst input=BHMat1All_May12@Toronto wcolumn=cat dmin=5 elevation=may23trial<br></div><div><br></div><div>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. </div><div>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,</div><div><br></div><div><br></div><div><div>----------------------------------------------------------------------------+</div><div> | Layer: may23trial@Toronto Date: Wed May 23 21:46:30 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: may23trial |</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: 61 |</div><div> | Total Cells: 466650 |</div><div> | Total size: 853368 Bytes |</div><div> | Number of tiles: 64 |</div><div> | Mean tile size: 13333 Bytes |</div><div> | Tile size in memory: 32384 Bytes |</div><div> | Number of tiles in x, y and z: 4, 4, 4 |</div><div> | Dimension of a tile in x, y, z: 22, 23, 16 |</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: 5.03278689 |</div><div> | Range of data: min = -0.18226054 max = 32.72385025 |</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=5 wscale=1.0 zscale=1.\ |</div><div> | 0 elevation="may23trial" |</div><div> | |</div><div> +----------------------------------------------------------------------------+</div><div>(Wed May 23 21:54:13 2018) Command finished (0 sec) </div></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">2018-05-23 21:42 GMT-04:00 Anna Petrášová <span dir="ltr"><<a href="mailto:kratochanna@gmail.com" target="_blank">kratochanna@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">On Wed, May 23, 2018 at 9:02 PM, Francois Chartier<br>
<<a href="mailto:fra.chartier@gmail.com">fra.chartier@gmail.com</a>> wrote:<br>
> How do you change the resolution to something coarser?<br>
<br>
</span>With g.region, res= is horizontal and tbres= is vertical resolution, e.g:<br>
<br>
g.region res=200 tbres=200 -p3<br>
<div class="HOEnZb"><div class="h5"><br>
><br>
> 2018-05-23 15:30 GMT-04:00 Markus Metz <<a href="mailto:markus.metz.giswork@gmail.com">markus.metz.giswork@gmail.com</a><wbr>>:<br>
>><br>
>><br>
>><br>
>> On Wed, May 23, 2018 at 1:09 PM, Francois Chartier<br>
>> <<a href="mailto:fra.chartier@gmail.com">fra.chartier@gmail.com</a>> 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 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 disk space<br>
>> > for temp files<br>
>> > WARNING: Points are more dense than specified 'DMIN'--ignored 65 points<br>
>> > (remain 76822)<br>
>> > WARNING: Unable to create <(null)> raster map without cross_input raster<br>
>> > map being specified<br>
>> ><br>
>> ><br>
>> ><br>
>> > <a href="http://v.info" rel="noreferrer" target="_blank">v.info</a><br>
>> > (Wed May 23 07:06:24 2018)<br>
>> > <a href="http://v.info" rel="noreferrer" target="_blank">v.info</a> map=BHMat1All_May12@Toronto<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>
>> > |-----------------------------<wbr>------------------------------<wbr>-----------------|<br>
>> > | Map format: native<br>
>> > |<br>
>> ><br>
>> > |-----------------------------<wbr>------------------------------<wbr>-----------------|<br>
>> > | Type of map: vector (level: 2)<br>
>> > |<br>
>> > |<br>
>> > |<br>
>> > | Number of points: 76887 Number of centroids: 0<br>
>> > |<br>
>> > | Number of lines: 0 Number of boundaries: 0<br>
>> > |<br>
>> > | Number of areas: 0 Number of islands: 0<br>
>> > |<br>
>> > | Number of faces: 0 Number of 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>
>> > +-----------------------------<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 input<br>
>> points<br>
>><br>
>> According to the output of <a href="http://v.info" rel="noreferrer" target="_blank">v.info</a>, a point covers on average about 7643645<br>
>> cubic meters which corresponds to voxels with edge length of 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 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 <<a href="mailto:mlennert@club.worldonline.be">mlennert@club.worldonline.be</a>><wbr>:<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 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>
>> >> - the output of <a href="http://v.info" rel="noreferrer" target="_blank">v.info</a> on your point dataset<br>
>> >><br>
>> >> ?<br>
>> >><br>
>> >> Moritz<br>
>> >><br>
>> >>><br>
>> >>> 2018-05-22 7:58 GMT-04:00 Moritz Lennert <<a href="mailto:mlennert@club.worldonline.be">mlennert@club.worldonline.be</a><br>
>> >>> <mailto:<a href="mailto:mlennert@club.worldonline.be">mlennert@club.<wbr>worldonline.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">fra.chartier@gmail.com</a> <mailto:<a href="mailto:fra.chartier@gmail.com">fra.chartier@gmail.com</a><wbr>>>:<br>
>> >>> >Hi,<br>
>> >>> ><br>
>> >>> >I am interpolating 2D data using V Vol RST. After 90 min the<br>
>> >>> process<br>
>> >>> >progressed 15%, and I stopped it.<br>
>> >>> >There are approx 370 points, and it should not take this long.<br>
>> >>> I am<br>
>> >>> >not<br>
>> >>> >sure what I am doing incorrectly. I have set the 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>
>> > <a href="mailto:grass-user@lists.osgeo.org">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/<wbr>mailman/listinfo/grass-user</a><br>
><br>
><br>
><br>
> ______________________________<wbr>_________________<br>
> grass-user mailing list<br>
> <a href="mailto:grass-user@lists.osgeo.org">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/<wbr>mailman/listinfo/grass-user</a><br>
</div></div></blockquote></div><br></div>