[GRASS-user] v vol rst super slow
Francois Chartier
fra.chartier at gmail.com
Sun May 27 14:01:04 PDT 2018
Hi,
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?
----------------------------------------------------------------------------+
| Layer: BHALL_May24v2 at Toronto Date: Sun May 27 16:43:46 2018
|
| Mapset: Toronto Login of Creator: Francois
Charti |
| Location: Toronto
|
| DataBase: C:\Users\Francois Chartier\Documents\grassdata
|
| Title: BHALL_May24v2
|
| Units: none
|
| Vertical unit: units
|
| Timestamp: none
|
|----------------------------------------------------------------------------|
|
|
| Type of Map: raster_3d Number of Categories: 0
|
| Data Type: FCELL
|
| Rows: 90
|
| Columns: 85
|
| Depths: 307
|
| Total Cells: 2348550
|
| Total size: 3180868 Bytes
|
| Number of tiles: 336
|
| Mean tile size: 9466 Bytes
|
| Tile size in memory: 30420 Bytes
|
| Number of tiles in x, y and z: 6, 7, 8
|
| Dimension of a tile in x, y, z: 15, 13, 39
|
|
|
| Projection: UTM (zone 17)
|
| N: 4870795 S: 4825715 Res: 500.88888889
|
| E: 651669 W: 609204 Res: 499.58823529
|
| T: 352 B: 45 Res: 1
|
| Range of data: min = -8.99099445 max = 33.89262009
|
|
|
| Data Source:
|
|
|
|
|
|
|
| Data Description:
|
| generated by v.vol.rst
|
|
|
| Comments:
|
| v.vol.rst input="BHMat1All_May12 at Toronto" wcolumn="cat" tension=40.
\ |
| smooth=0.1 segmax=50 npmin=200 npmax=700 dmin=2 wscale=1.0
zscale=1.\ |
| 0 elevation="BHALL_May24v2"
|
|
|
+----------------------------------------------------------------------------+
(Sun May 27 16:57:19 2018) Command finished (0 sec)
2018-05-24 4:52 GMT-04:00 Moritz Lennert <mlennert at club.worldonline.be>:
> On 24/05/18 04:05, Francois Chartier wrote:
>
>> Thank you,
>>
>> v.vol.rst input=BHMat1All_May12 at Toronto wcolumn=cat dmin=5
>> elevation=may23trial
>>
>> 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.
>>
>
> Watch out: according to the info below, your horizontal resolution is
> around 500, not 200.
>
> 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.
>
> Moritz
>
>
>>
>> ------------------------------------------------------------
>> ----------------+
>> | Layer: may23trial at Toronto Date: Wed May 23 21:46:30
>> 2018 |
>> | Mapset: Toronto Login of Creator: Francois
>> Charti |
>> | Location: Toronto
>> |
>> | DataBase: C:\Users\Francois Chartier\Documents\grassdata
>> |
>> | Title: may23trial
>> |
>> | Units: none
>> |
>> | Vertical unit: units
>> |
>> | Timestamp: none
>> |
>> |---------------------------------------------------------
>> -------------------|
>> |
>> |
>> | Type of Map: raster_3d Number of Categories: 0
>> |
>> | Data Type: FCELL
>> |
>> | Rows: 90
>> |
>> | Columns: 85
>> |
>> | Depths: 61
>> |
>> | Total Cells: 466650
>> |
>> | Total size: 853368 Bytes
>> |
>> | Number of tiles: 64
>> |
>> | Mean tile size: 13333 Bytes
>> |
>> | Tile size in memory: 32384 Bytes
>> |
>> | Number of tiles in x, y and z: 4, 4, 4
>> |
>> | Dimension of a tile in x, y, z: 22, 23, 16
>> |
>> |
>> |
>> | Projection: UTM (zone 17)
>> |
>> | N: 4870795 S: 4825715 Res: 500.88888889
>> |
>> | E: 651669 W: 609204 Res: 499.58823529
>> |
>> | T: 352 B: 45 Res: 5.03278689
>> |
>> | Range of data: min = -0.18226054 max = 32.72385025
>> |
>> |
>> |
>> | Data Source:
>> |
>> |
>> |
>> |
>> |
>> |
>> |
>> | Data Description:
>> |
>> | generated by v.vol.rst
>> |
>> |
>> |
>> | Comments:
>> |
>> | v.vol.rst input="BHMat1All_May12 at Toronto" wcolumn="cat"
>> tension=40. \ |
>> | smooth=0.1 segmax=50 npmin=200 npmax=700 dmin=5 wscale=1.0
>> zscale=1.\ |
>> | 0 elevation="may23trial"
>> |
>> |
>> |
>> +---------------------------------------------------------
>> -------------------+
>> (Wed May 23 21:54:13 2018) Command finished (0 sec)
>>
>>
>> 2018-05-23 21:42 GMT-04:00 Anna Petrášová <kratochanna at gmail.com <mailto:
>> kratochanna at gmail.com>>:
>>
>> On Wed, May 23, 2018 at 9:02 PM, Francois Chartier
>> <fra.chartier at gmail.com <mailto:fra.chartier at gmail.com>> wrote:
>> > How do you change the resolution to something coarser?
>>
>> With g.region, res= is horizontal and tbres= is vertical resolution,
>> e.g:
>>
>> g.region res=200 tbres=200 -p3
>>
>> >
>> > 2018-05-23 15:30 GMT-04:00 Markus Metz
>> <markus.metz.giswork at gmail.com <mailto:markus.metz.giswork at gmail.com
>> >>:
>> >>
>> >>
>> >>
>> >> On Wed, May 23, 2018 at 1:09 PM, Francois Chartier
>> >> <fra.chartier at gmail.com <mailto:fra.chartier at gmail.com>> 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 at 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 <http://v.info>
>> >> > (Wed May 23 07:06:24 2018)
>> >> > v.info <http://v.info> map=BHMat1All_May12 at 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 <http://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
>> <mlennert at club.worldonline.be <mailto:mlennert at club.worldonline.be>>:
>> >> >>
>> >> >> 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 <http://v.info> on your point dataset
>> >> >>
>> >> >> ?
>> >> >>
>> >> >> Moritz
>> >> >>
>> >> >>>
>> >> >>> 2018-05-22 7:58 GMT-04:00 Moritz Lennert
>> <mlennert at club.worldonline.be <mailto:mlennert at club.worldonline.be>
>> >> >>> <mailto:mlennert at club.worldonline.be
>> <mailto:mlennert at club.worldonline.be>>>:
>> >> >>>
>> >> >>>
>> >> >>>
>> >> >>> Am 22. Mai 2018 04:35:13 MESZ schrieb Francois Chartier
>> >> >>> <fra.chartier at gmail.com <mailto:fra.chartier at gmail.com>
>> <mailto:fra.chartier at gmail.com <mailto:fra.chartier at gmail.com>>>:
>> >> >>> >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
>> >> > grass-user at lists.osgeo.org <mailto:grass-user at lists.osgeo.org>
>> >> > https://lists.osgeo.org/mailman/listinfo/grass-user
>> <https://lists.osgeo.org/mailman/listinfo/grass-user>
>> >
>> >
>> >
>> > _______________________________________________
>> > grass-user mailing list
>> > grass-user at lists.osgeo.org <mailto:grass-user at lists.osgeo.org>
>> > https://lists.osgeo.org/mailman/listinfo/grass-user
>> <https://lists.osgeo.org/mailman/listinfo/grass-user>
>>
>>
>>
>>
>> _______________________________________________
>> grass-user mailing list
>> grass-user at lists.osgeo.org
>> https://lists.osgeo.org/mailman/listinfo/grass-user
>>
>>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20180527/e0cd1264/attachment-0001.html>
More information about the grass-user
mailing list