Hi, your region settings are interesting. Do you really need a resolution of 5m in x,y plane and 1m in z direction? As a result you have 23,507,876,616 voxel to compute, but only 76887 points for interpolation.
I would suggest you reduce the resolution of your region significantly, at least to test the interpolation settings. Best regards Soeren 2018-05-23 13:09 GMT+02:00 Francois Chartier <[email protected]>: > 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 > (Wed May 23 07:06:24 2018) > 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) > > > 2018-05-23 6:07 GMT-04:00 Moritz Lennert <[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 on your point dataset >> >> ? >> >> Moritz >> >>> >>> 2018-05-22 7:58 GMT-04:00 Moritz Lennert <[email protected] >>> <mailto:[email protected]>>: >>> >>> >>> >>> Am 22. Mai 2018 04:35:13 MESZ schrieb Francois Chartier >>> <[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] > https://lists.osgeo.org/mailman/listinfo/grass-user _______________________________________________ grass-user mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/grass-user
