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

2018-07-11 Thread Micha Silver
Thanks for reporting I didn't realize that exporting could fail due to a  file extension. On 07/12/2018 01:31 AM, Gabriel Cotlier wrote: solved with correcting spelling mistake at file extension, now the for loop worked for me as: 

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

2018-07-11 Thread Francois Chartier
Hi Moritz, Thanks for your email. What you described is the concept for developing most geological model by correlating top of layers (rock, soils, formation etc.) which is the correct approach when geological contacts are clear. In my case (as part of a masters/engineering degree), I have so

[GRASS-user] About Spatial PCA in GRASS

2018-07-11 Thread Shiva Khanal
Hi, I am interested to run the spatial-mode PCA on raster time series. In R it is implemented as an ‘EOF’ function in spacetime Package. The function transposes the raster stack and applies ‘prcomp’ function. I would like to know if there is any possible way to implement i.pca similar to this

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

2018-07-11 Thread Gabriel Cotlier
solved with correcting spelling mistake at file extension, now the for loop worked for me as: Gabriel On Wed, Jul 11, 2018 at 6:39 PM, Gabriel Cotlier wrote: > Dear Nikos and Micha, > > I tried both and worked ways for constructing the list and worked well: > > > > and also: > > > > but

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

2018-07-11 Thread Gabriel Cotlier
Dear Nikos and Micha, I tried both and worked ways for constructing the list and worked well: and also: but still have an error problem with the loop why could be this happening? Thanks a lot. Gabriel On Wed, Jul 11, 2018 at 5:00 AM, Nikos Alexandris wrote: > * Nikos Alexandris

Re: [GRASS-user] Working with large attribute table [RESOLVED]

2018-07-11 Thread Rich Shepard
On Wed, 11 Jul 2018, Rich Shepard wrote: I'm trying to obtain attribute data for the two stations within the project area by clicking on the point using the 'query item' button on the map display window. I used gawk to extract the mappable attributes and am using those. Time series and

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

2018-07-11 Thread Rich Shepard
On Wed, 11 Jul 2018, Markus Metz wrote: If 'gdalinfo -proj4' provides a reasonable string, you can just as well use - grass75 -c geofile GISDBASE/LOCATION_NAME/MAPSET - the location wizard using this georeferenced data file - from within GRASS 'g.proj georef= location=' - from with GRASS

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

2018-07-11 Thread Markus Metz
On Wed, Jul 11, 2018 at 2:44 PM, Rich Shepard wrote: > > On Wed, 11 Jul 2018, Nikos Alexandris wrote: > >> The grass --help does not mention support for "Proj" strings. Only >> geo(-referenced)file and EPSG. > > > Nikos, > > Yep. That's why I asked. So, on the infrequent occasions when source

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

2018-07-11 Thread Markus Metz
On Wed, Jul 11, 2018 at 3:55 AM, Huidae Cho wrote: > > 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

[GRASS-user] Working with large attribute table

2018-07-11 Thread Rich Shepard
I have a point map of weather stations and records since about 2005. The attribute table has 113,570 rows; each row has the station name and date as the primary key and there are 20 attribute columns. The reason I imported all 20 attributes is to learn how to turn this map into a

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

2018-07-11 Thread Nikos Alexandris
* Rich Shepard [2018-07-11 09:24:01 -0700]: On Wed, 11 Jul 2018, Nikos Alexandris wrote: By the way, the -gtext option is not relevant. It merely shows the welcome/logo screen when launching a GRASS GIS session, then jumps back to the command line. Nikos, Not with 7.5svn. When I use

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

2018-07-11 Thread Rich Shepard
On Wed, 11 Jul 2018, Nikos Alexandris wrote: By the way, the -gtext option is not relevant. It merely shows the welcome/logo screen when launching a GRASS GIS session, then jumps back to the command line. Nikos, Not with 7.5svn. When I use '-gtext' I get the location wizard. After I run

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

2018-07-11 Thread Rich Shepard
On Wed, 11 Jul 2018, Moritz Lennert wrote: You can have a vector map with features that have a cat value, but no corresponding entry in the attribute table, and you can have entries in the attribute table with no corresponding features in the map. Moritz, Yesterday, initially, there were

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

2018-07-11 Thread Nikos Alexandris
* Rich Shepard [2018-07-11 08:24:27 -0700]: On Wed, 11 Jul 2018, Nikos Alexandris wrote: I see. Maybe we should file an enhancement ticket, to provide this option through the command line too? Nikos, I don't know how common this is. This project is the first time I've had DEM source

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

2018-07-11 Thread Moritz Lennert
On 11/07/18 16:48, Rich Shepard wrote: On Wed, 11 Jul 2018, Nikos Alexandris wrote: Please post (also) the exact `v.edit` command. Nikos, Very strange. Looking at the map now only one point is visible (cat #3) while both cat #1 and #3 are shown in the attribute table. This means it's

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

2018-07-11 Thread Rich Shepard
On Wed, 11 Jul 2018, Nikos Alexandris wrote: I see. Maybe we should file an enhancement ticket, to provide this option through the command line too? Nikos, I don't know how common this is. This project is the first time I've had DEM source files in ESRI GRID (*.adf) format and this format

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

2018-07-11 Thread Nikos Alexandris
* Rich Shepard [2018-07-11 05:44:47 -0700]: On Wed, 11 Jul 2018, Nikos Alexandris wrote: The grass --help does not mention support for "Proj" strings. Only geo(-referenced)file and EPSG. Nikos, Yep. That's why I asked. So, on the infrequent occasions when source data does not have a

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

2018-07-11 Thread Rich Shepard
On Wed, 11 Jul 2018, Moritz Lennert wrote: See also the general raster info: https://grass.osgeo.org/grass74/manuals/rasterintro.html which explains that: 1.Raster output maps have their bounds and resolution equal to those of the current computational region. 2.Raster input maps are

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

2018-07-11 Thread Rich Shepard
On Wed, 11 Jul 2018, Nikos Alexandris wrote: Please post (also) the exact `v.edit` command. Nikos, Very strange. Looking at the map now only one point is visible (cat #3) while both cat #1 and #3 are shown in the attribute table. This means it's time to g.remove that map and start over.

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

2018-07-11 Thread Moritz Lennert
On 11/07/18 05:46, Huidae Cho wrote: 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

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

2018-07-11 Thread Rich Shepard
On Wed, 11 Jul 2018, Nikos Alexandris wrote: The grass --help does not mention support for "Proj" strings. Only geo(-referenced)file and EPSG. Nikos, Yep. That's why I asked. So, on the infrequent occasions when source data does not have a useful EPSG code but 'gdalinfo -proj4' provies the

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

2018-07-11 Thread Rich Shepard
On Wed, 11 Jul 2018, Huidae Cho wrote: 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. Huidae, Thank you for another complete answer.

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

2018-07-11 Thread Micha Silver
On 07/11/2018 11:46 AM, marion-bru...@sfr.fr wrote: ‌Thank you for this clarification, one last question: the cell size is nsres*ewres right? Yes, exactly -- Micha Silver Ben Gurion Univ. Sde Boker, Remote Sensing Lab

Re: [GRASS-user] creating a new mapset in an existing location

2018-07-11 Thread Ken Mankoff
Right. Sorry for the silly question. I thought that "-c" was only for creating locations, but it clearly says it creates mapsets too. Thank you! -k. On 7/11/18, Moritz Lennert wrote: > On 11/07/18 11:24, Ken Mankoff wrote: >> Dear List, >> >> Is there a way to create a new mapset from the

Re: [GRASS-user] creating a new mapset in an existing location

2018-07-11 Thread Moritz Lennert
On 11/07/18 11:24, Ken Mankoff wrote: Dear List, Is there a way to create a new mapset from the command line? My use case is that I'd like to run parallel jobs (with GNU parallel). But, ...it's recommended to launch each job in its own mapset within the location. according to

[GRASS-user] creating a new mapset in an existing location

2018-07-11 Thread Ken Mankoff
Dear List, Is there a way to create a new mapset from the command line? My use case is that I'd like to run parallel jobs (with GNU parallel). But, > ...it's recommended to launch each job in its own mapset within the > location. according to

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

2018-07-11 Thread marion-brunet
‌Thank you for this clarification, one last question: the cell size is nsres*ewres right? De : "Micha Silver" A : marion-bru...@sfr.fr Envoyé: mardi 10 juillet 2018 21:06 Objet : Re: [GRASS-user] calculate the volume of a DEM body p { margin-bottom: 0cm; margin-top: 0pt; } On 07/10/2018

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

2018-07-11 Thread marion-brunet
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,

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

2018-07-11 Thread Nikos Alexandris
* Nikos Alexandris [2018-07-11 09:14:47 +0200]: Gabriel Cotlier wrote: 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...

Re: [GRASS-user] script option

2018-07-11 Thread Nikos Alexandris
* Frank David [2018-07-10 15:44:51 +0200]: 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

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

2018-07-11 Thread Nikos Alexandris
* Micha Silver [2018-07-10 22:34:53 +0300]:   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

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

2018-07-11 Thread Nikos Alexandris
* Rich Shepard [2018-07-10 15:11:32 -0700]: 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

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

2018-07-11 Thread Nikos Alexandris
Gabriel Cotlier wrote: 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