How do you change the resolution to something coarser? 2018-05-23 15:30 GMT-04:00 Markus Metz <[email protected]>:
> > > On Wed, May 23, 2018 at 1:09 PM, Francois Chartier <[email protected]> > wrote: > > > > 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) > > > The current target resolutions of > > nsres3: 5 > ewres3: 5 > tbres: 1 > > are much too high for the given extents and the coverage of the input > points > > According to the output of v.info, 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. > > Markus M > > > > > > 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
