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

Hi, On Tue, Jul 10, 2018 at 6:04 AM, wrote: ... > I had two DEMS and I used r.mapcalc to get > the difference between them. I want to calculate the total volume of this > last raster. I wonder if https://grass.osgeo.org/grass74/manuals/r.volume.html could be of any use here? Markus

Re: [GRASS-user] script option

* frank [2018-07-10 08:05:29 +0200]: Hi dear grass users, I'm trying to get two field name of an attribute table (ctx_map) by the options of the gui. But I can access only at the first one and I don't find why...??. See my code below : #%option G_OPT_V_MAP #% key: ctx_map #% description:

[GRASS-user] script option

Hi dear grass users, I'm trying to get two field name of an attribute table (ctx_map) by the options of the gui. But I can access only at the first one and I don't find why...??. See my code below : #%option G_OPT_V_MAP #% key: ctx_map #% description: vector data #% required : yes #%

Re: [GRASS-user] r3.in.xyz

Hi François, On 10/07/18 04:40, Francois Chartier wrote: Hi Moritz, The idea is a bit novel.  I am working with over 2 Boreholes (BH) with soil information (silt, sand, gravel, Tills, and everything in between etc.).  These are not distinct geological contacts with recognizable

Re: [GRASS-user] Start in text mode: define new location with proj.4 string [RESOLVED]

On Mon, 9 Jul 2018, Rich Shepard wrote: When I initially imported the DEMs I pointed the GUI location wizard at the appropriate prj.adf file and applied v.in.ogr to the hdr.adf file. These were then re-projected to the project's location, patched together into a mosaic. Since these cover only a

Re: [GRASS-user] script option

Le 10/07/2018 à 12:25, Nikos Alexandris a écrit : * frank [2018-07-10 08:05:29 +0200]: Hi dear grass users, I'm trying to get two field name of an attribute table (ctx_map) by the options of the gui. But I can access only at the first one and I don't find why...??. See my code below :

[GRASS-user] g.remove failure

g.remove -f type=vect name=huc12 WARNING: No data base element files found and the map is not removed. But, I can display the map, use the 'query' button in the display monitor, and get the attribute information about individual hydrological units. Why might g.remove fail this request?

[GRASS-user] r.mask vs g.region

I have DEMs covering a large area and want to limit analyses to a single drainage basin. I can use 'g.region vect=basin_name' or 'r.mask vect=basin_name' and want to learn of differences by applying one in preference to tha other. Are there situations where one is preferred over the other? Does

[GRASS-user] question on r.out.gdal in a loop

Dear grass users, I'm using grass 7.4.1 trying to export all raster layers in the current mapset out of grass to a folder as GeoTiff files though a for loop from the python shell as follows and I'm getting the error in red... What could be the problem? Thanks a lot in advance. Regards, Gabriel

Re: [GRASS-user] g.remove failure

On Tue, 10 Jul 2018, Veronica Andreo wrote: Check if you are in the mapset where the map is, otherwise you cannot delete maps; i.e. maps can only be removed from the mapset they belong to Vero, Yes, that's probably the case. Thanks, Rich ___

[GRASS-user] Vector map point deleted but still displays

I removed a point from a vector map but it still shows up on the display monitor (see attached screenshot, point #2). The attribute table no longer contains that row: cat|y|x|label 1|45.299844|-122.353975|USGS gauge 3|45.302761|-122.358851|Building I've tried to remove this point from the

Re: [GRASS-user] Start in text mode: define new location with proj.4 string

  Reading the Startup page I do not see a way to apply the proj.4 string to create a new location on the command line. I know this is an option with the GUI location wizard but do not see the appropriate '-c' modifier to do this

Re: [GRASS-user] Start in text mode: define new location with proj.4 string

On Tue, 10 Jul 2018, Micha Silver wrote: What you've suggested in past posts seems to work fine (point to a shapefile with accompanying *.prj): Unless I'm not understanding the problem. Micha, My reason for writing was having had raster files imported incorrectly. Using 'gdalinfo

Re: [GRASS-user] g.remove failure

Check if you are in the mapset where the map is, otherwise you cannot delete maps; i.e. maps can only be removed from the mapset they belong to HTH Vero El mar., 10 jul. 2018 a las 20:44, Rich Shepard () escribió: > g.remove -f type=vect name=huc12 > WARNING: No data base element files found >

Re: [GRASS-user] limit number of confluent streams in r.steam.extract or r.watershed

Hello, Not sure if I understood your question. Are you trying to extract major streams only? r.stream.extract already returns a valid stream network with confluences snapped to the main stem. Also, limiting the number of confluent streams and moving (?) streams at the confluences are two

Re: [GRASS-user] question on r.out.gdal in a loop

On 07/11/2018 02:20 AM, Gabriel Cotlier wrote: Dear grass users, I'm using grass 7.4.1 trying to export  all raster layers in the current mapset out of grass to a folder as GeoTiff files

Re: [GRASS-user] r.mask vs g.region

Hello Rich, g.region creates a rectangular computational region while r.mask can create an irregular mask. If you use r.mask, any cells outside the mask will be ignored, so you have to make sure that your mask is big enough to cover your basin. In terms of computational time, r.mask can be faster

Re: [GRASS-user] 7.5svn: create location from existing map

Hello Rich, Create a new directory and copy & rename biota to PERMANENT in this folder. Copy PERMANENT/WIND to PERMANENT/DEFAULT_WIND. Now this folder is a new location with the biota data as a PERMANENT mapset. Best, Huidae On Fri, Jun 29, 2018 at 5:19 PM, Rich Shepard wrote: > The