Hot answers tagged

5

The module r.univar can take a comma separated list of rasters as input, and calculates stats on the whole list. Does that help? MAPLIST=`g.list rast pattern="rast*" separator=comma` g.region -p $MAPLIST r.univar input=${MAPLIST} percentile=90


3

UnicodeDecodeError: 'charmap' codec can't decode byte 0x8c in position 113: character maps to I have had problems with GRASS with a non ASCII charachter (like ö) in the path. Try changing it.


2

Unless your terrain is very/extremely(?) mountainous you are not going to have many slopes > 60 degrees, so maybe almost all your slopes are 0-60? Try: Slope Raster calculator on your slope. I calculate slopes <=30 to 1. Steeper slopes gets value 0. Style the output using "Paletted/Unique values and add only the value 1:


2

Micha's answer is right. If instead you do a for loop with r.univar, you get the same result than with t.rast.univar which is a wrapper for r.univar for all the maps within a time series.


1

What the problem was After inspecting your data, I could identify the problem: you have a lot of duplicate points as well as a few points that are "almost" duplicates (extremely close to each other). From originally 954 features in the point layer you provided, only 229 (ca. 24%) are really unique points (thus not duplicates or quasi-duplicates). ...


1

For r.in.lidar, you need to have the Liblas library installed (before creating the Mac binary package of GRASS GIS). But Liblas is rather outdated. As an alternative, you may install PDAL and then the addon r.in.pdal (https://grass.osgeo.org/grass7/manuals/addons/r.in.pdal.html) which relies on it.


1

The problem was identified in this thread at the GRASS mail list. The manual page for the v.surf.idw module states the following: **column**=name Name of attribute column with values to interpolate If not given and input is 2D vector map then category values are used. If input is 3D vector map then z-coordinates are used. The points ...


1

Eventually I found the simple and effective way. The expression array_contains( overlay_nearest( @layer , max_distance:=1000 , expression:="name" , limit:= 40 ) , "name") will return true if in a radius of 1000 layer units there is at least ...


1

You can also calculate the length of a line taking into account the topography the following way: Use the "Profiles from line" (SAGA) tool, which adds points along the line, giving the length from the line start to each point ("DIST" as the length on a flat surface, "DIST_SURF" with heights included). Select the point with the ...


Only top voted, non community-wiki answers of a minimum length are eligible