[GRASS-user] v vol rst super slow

Sören Gebbert soerengebbert at googlemail.com
Wed May 23 04:18:25 PDT 2018


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 <fra.chartier at gmail.com>:
> 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
> (Wed May 23 07:06:24 2018)
> 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)
>
>
> 2018-05-23 6:07 GMT-04:00 Moritz Lennert <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 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>>:
>>>
>>>
>>>
>>>     Am 22. Mai 2018 04:35:13 MESZ schrieb Francois Chartier
>>>     <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
> https://lists.osgeo.org/mailman/listinfo/grass-user


More information about the grass-user mailing list