[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