[GRASS-user] v vol rst super slow

[Francois Chartier]

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)                             *

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
>>
>>
>>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20180523/0dec4e3d/attachment.html>