Re: [GRASS-user] calculate the volume of a DEM with r.volume

Micha Silver Fri, 13 Jul 2018 04:19:20 -0700

Hi
I noticed you have a raster with a "minus" character in the name.
Sometimes that gives trouble. Not sure if it's relevant to the current problem, but it might be worth renaming to something like "dec_juin" (underscore instead of '-') and trying again.

On 07/13/2018 10:51 AM, marion-bru...@sfr.fr wrote:

‌Hello,

I am working on Windows and I tried using GRASS 7.2.2, GRASS 7.4.0 and GRASS 7.4.1 to see if there is any difference, but it is always the same.

I have a ratser map called dec-juin
First I did r.clump input=dec-juin@1217_Interpo0618 output=test
to get the clumps
Then I did r.volume input=dec-juin@1217_Interpo0618 clump=test@1217_Interpo0618
to have the volume.
The result I have is:

(Fri Jul 13 14:12:45 2018)
r.volume input=dec-juin@1217_Interpo0618 clump=test@1217_Interpo0618
Volume report on data from <dec-juin@1217_Interpo0618> using clumps on <test@1217_Interpo0618> raster map
Category   Average   Data   # Cells        Centroid             Total
Number     in clump Total in clump   Easting     Northing     Volume
-----------------------------------------------------------------------------
       1     -8.02       -72       9   191262.28 1222178.97           -22.59
       2     -9.02      -622      69   191265.64 1222177.29          -194.78
       3     -8.04      -394      49   191268.44 1222176.73          -123.30
       4     -7.08      -276      39   191271.23 1222177.85           -86.46
       5     -5.94      -427      72   191274.59 1222176.73          -133.75
       6     -4.30      -275      64   191285.22 1222178.97           -86.09
       7     -5.10     -2714     532   191288.58 1222174.50          -849.47
       8     -0.60        -7      12   191221.43 1222177.29            -2.27
       9     -8.95     -1030     115   191230.39 1222176.17          -322.31
      10    -10.06     -1217     121   191245.49 1222175.61          -380.88
      11    -10.93     -1956     179   191249.97 1222177.29          -612.34
      12    -10.24      -645      63   191259.48 1222175.61          -201.92
      13     -7.05      -621      88   191329.43 1222177.29          -194.28
      14     -7.51       -60       8   191336.14 1222177.29           -18.81
      15     -7.47       -30       4   191337.26 1222177.29            -9.35
      16     -8.04      -322      40   191340.06 1222177.29          -100.63
      17     -8.99      -108      12   191343.98 1222177.29           -33.78
      18      4.98        45       9   191182.26 1222175.05            14.03
      etc
      etc
   24458     -5.50       -22       4   191397.13 1221585.50            -6.88
   24459     -5.51       -33       6   191378.11 1221584.39           -10.34
   24460     -5.52       -11       2   191380.91 1221584.39            -3.45
   24461     -5.27    -30348    5764   191309.84 1221601.17         -9498.51
   24462     -5.49       -38       7   191358.53 1221582.15           -12.02
   24463     -5.49       -22       4   191363.56 1221583.83            -6.87
   24464     -5.48       -11       2   191385.38 1221583.83            -3.43
   24465     -5.49       -11       2   191387.06 1221583.83            -3.44
   24466     -5.50       -11       2   191388.74 1221583.83            -3.44
   24467     -5.50       -22       4   191360.20 1221582.71
           -6.88
   24468     -5.92   -523187   88358   191241.58 1221675.00       -163750.63
(Fri Jul 13 14:13:03 2018) La commande s'est terminée (18 sec)

As you can see, I have nothing at the end of the command... But I know it is finished, because when I created the clumps it said that it created 24468 clumps.
Do you need any more information?

Thank you very much

De : "Markus Neteler"
A : marion-bru...@sfr.fr,"grass-user"
Envoyé: jeudi 12 juillet 2018 16:43
Objet : Re: [GRASS-user] calculate the volume of a DEM with r.volume

Hello Marion,

On Wed, Jul 11, 2018 at 10:32 AM, <marion-bru...@sfr.fr> wrote:
> Hello,
>
> I tried using r.volume as you suggested. According to the documentation it
> is suitable for what I am trying to do.
>
> However, when I use it, I don't have the "Total volume" line at the end of
> the report, like it is shown in the documentation. I just have the total
> volume of each clump, but not of all the clumps... Has it ever happened to
> you when using this command?

That's weird! Which version/operating system do you use?

> in the command line
>
> I simply wroter.volume input=difference clump=clump in the command line

Here on command line:

GRASS 7.4.1svn (nc_spm_08_grass7):~ > g.region raster=elevation
GRASS 7.4.1svn (nc_spm_08_grass7):~ > r.volume input=elevation clump=geology_30m
100%

Volume report on data from <elevation> using clumps on <geology_30m> raster map

Category Average Data # Cells Centroid Total
Number in clump Total in clump Easting Northing Volume
-----------------------------------------------------------------------------
217 118.93 86288828 725562 635325.00 221535.00 8628882798.63
262 108.97 21650560 198684 638935.00 222495.00 2165056037.02
270 92.23 63578874 689373 642405.00 221485.00 6357887443.53
405 132.96 33732662 253710 631835.00 224095.00 3373266208.59
583 139.35 3011288 21609 630205.00 224665.00 301128821.55
720 124.30 599618 4824 634075.00 227995.00 59961816.06
766 132.43 936791 7074 631425.00 227845.00 93679120.08
862 118.31 7302317 61722 630505.00 218885.00 730231746.74
910 94.20 4235816 44964 639215.00 216365.00 423581613.11
921 135.22 1693985 12528 630755.00 215445.00 169398523.05
945 127.24 1145 9 630015.00 215015.00 114512.03
946 89.91 365748 4068 639085.00 215255.00 36574833.85
948 129.02 112632 873 630185.00 215115.00 11263181.57
-----------------------------------------------------------------------------
Total Volume = 22351026655.81

... which looks ok.

Can you say a bit more about your system?

Markus
_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

-- 
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918

_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] calculate the volume of a DEM with r.volume

Reply via email to