<div dir="ltr">How do you change the resolution to something coarser?</div><div class="gmail_extra"><br><div class="gmail_quote">2018-05-23 15:30 GMT-04:00 Markus Metz <span dir="ltr"><<a href="mailto:markus.metz.giswork@gmail.com" target="_blank">markus.metz.giswork@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div class="h5"><br><br>On Wed, May 23, 2018 at 1:09 PM, Francois Chartier <<a href="mailto:fra.chartier@gmail.com" target="_blank">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 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 for temp files<br>> WARNING: Points are more dense than specified 'DMIN'--ignored 65 points (remain 76822)<br>> WARNING: Unable to create <(null)> raster map without cross_input raster map being specified<br>><br>><br>><br>> <a href="http://v.info" target="_blank">v.info</a><br>> (Wed May 23 07:06:24 2018) <br>> <a href="http://v.info" target="_blank">v.info</a> map=BHMat1All_May12@Toronto <br>> +----------------------------<wbr>------------------------------<wbr>------------------+<br>> | Name: BHMat1All_May12 |<br>> | Mapset: Toronto |<br>> | Location: Toronto |<br>> | Database: C:\Users\Francois Chartier\Documents\grassdata |<br>> | Title: |<br>> | Map scale: 1:1 |<br>> | Name of creator: Francois Chartier |<br>> | Organization: |<br>> | Source date: Sat May 12 15:45:44 2018 |<br>> | Timestamp (first layer): none |<br>> |----------------------------<wbr>------------------------------<wbr>------------------|<br>> | Map format: native |<br>> |----------------------------<wbr>------------------------------<wbr>------------------|<br>> | Type of map: vector (level: 2) |<br>> | |<br>> | Number of points: 76887 Number of centroids: 0 |<br>> | Number of lines: 0 Number of boundaries: 0 |<br>> | Number of areas: 0 Number of islands: 0 |<br>> | Number of faces: 0 Number of kernels: 0 |<br>> | Number of volumes: 0 Number of holes: 0 |<br>> | |<br>> | Map is 3D: Yes |<br>> | Number of dblinks: 1 |<br>> | |<br>> | Projection: UTM (zone 17) |<br>> | |<br>> | N: 4870795 S: 4825715 |<br>> | E: 651669 W: 609204 |<br>> | B: 45 T: 352 |<br>> | |<br>> | Digitization threshold: 0 |<br>> | Comment: |<br>> | |<br>> +----------------------------<wbr>------------------------------<wbr>------------------+<br>> (Wed May 23 07:06:24 2018) Command finished (0 sec) <br><br></div></div><div>The current target resolutions of<br><br>nsres3: 5<br>ewres3: 5<br>tbres: 1<br><br></div>are much too high for the given extents and the coverage of the input points<br><div></div><br>According to the output of <a href="http://v.info" target="_blank">v.info</a>, 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.<br><br></div><div>Markus M<br></div><div><div><div class="h5">><br>><br>> 2018-05-23 6:07 GMT-04:00 Moritz Lennert <<a href="mailto:mlennert@club.worldonline.be" target="_blank">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, 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" 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" target="_blank">mlennert@club.worldonline.be</a> <mailto:<a href="mailto:mlennert@club.worldonline.be" target="_blank">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" target="_blank">fra.chartier@gmail.com</a> <mailto:<a href="mailto:fra.chartier@gmail.com" target="_blank">fra.chartier@gmail.com</a><wbr>>>:<br>>>> >Hi,<br>>>> ><br>>>> >I am interpolating 2D data using V Vol RST. After 90 min the process<br>>>> >progressed 15%, and I stopped it.<br>>>> >There are approx 370 points, and it should not take this long. 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></div></div><span class="">> ______________________________<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" target="_blank">https://lists.osgeo.org/<wbr>mailman/listinfo/grass-user</a><br></span></div></div>
</blockquote></div><br></div>