[GRASS-user] v.in.csv available?

Hi all, I was wondering if the addon v.in.csv is available or not? I want to import a simple point file (*.csv) which is unfortunately in another projection then my location. So the tool v.in.csv looks handy. This tool is listed here https://grass.osgeo.org/grass78/manuals/addons/ but not

Re: [GRASS-user] Merge spatially connected features

On Thu, Mar 12, 2020 at 10:48 PM Markus Metz wrote: > > > On Thu, Mar 12, 2020 at 12:34 PM Johannes Radinger < > johannesradin...@gmail.com> wrote: > > > > Thank you Markus, > > indeed your approach looks like what I need..The hint with > v.net.com

Re: [GRASS-user] v.edit break lines at coordinates

Hi all, FYI: I just issued a new ticket in github ( https://github.com/OSGeo/grass/issues/415) reporting that issue and accompanied with some more information. /Johannes On Fri, Mar 13, 2020 at 11:53 AM Johannes Radinger < johannesradin...@gmail.com> wrote: > Hi GRASS users, > >

[GRASS-user] v.edit break lines at coordinates

Hi GRASS users, maybe someone of you can explain me following problem (of understanding) I am facing: I want to break a single line into several. More specifically, I am using v.edit with the option break and two coordinate pairs to break a single line. I'd expect that breaking a single line at

Re: [GRASS-user] v.in.ogr layer name

Hi Markus, On Thu, Mar 12, 2020 at 10:19 PM Markus Metz wrote: > > > On Thu, Mar 12, 2020 at 12:43 PM Johannes Radinger < > johannesradin...@gmail.com> wrote: > > > > Hi all, > > I am using v.in.ogr to import a shape file (line vector) that has only &

[GRASS-user] Meaning of flag 'i' in d.vect

Hi all, I am just wondering what exactly the i-flag in d.vector means? It is described as "Use values from 'cats' option as feature id". What is the effect of this in the display? Furthermore, I was wondering how the query tool in the map display works. For example, I display a vector map with 3

[GRASS-user] v.in.ogr layer name

Hi all, I am using v.in.ogr to import a shape file (line vector) that has only one layer. Although I am specifying the output parameter in v.in.ogr the layer name of the imported map still remains the original; it is only the output map and the associated table that gets a new name as defined by

Re: [GRASS-user] Merge spatially connected features

er (v.extract) > identify connected lines (v.net + v.net.components) > update a new attribute of the original lines with the comp attribute of > the output of v.net.components plus some offset to separate different > stream orders > > HTH, > > Markus M > > > On Tue,

Re: [GRASS-user] Merge spatially connected features

So...no also with GRASS-user as recipient... On 05.03.20 16:21, Micha Silver wrote: On 3/5/20 10:47 AM, Johannes Radinger wrote: Hi Micha, hi all, sorry for my late response...however, just today I managed to try your approach of building polylines to connect "touching stream

Re: [GRASS-user] Merge spatially connected features

Hi Micha, hi all, sorry for my late response...however, just today I managed to try your approach of building polylines to connect "touching stream lines"...but... On 24.02.20 16:48, Micha Silver wrote: On 24/02/2020 10:45, Johannes Radinger wrote: Hi all, I have a large river netwo

Re: [GRASS-user] Merge spatially connected features

Hi Micha, hi all, sorry for my late response...however, just today I managed to try your approach of building polylines to connect "touching stream lines"...but... On 24.02.20 16:48, Micha Silver wrote: On 24/02/2020 10:45, Johannes Radinger wrote: Hi all, I have a large riv

[GRASS-user] Splitting lines by other lines that overlap

Hi all, I have two line vectors: LV1) a complete stream network LV2) sections of the same stream network representing impoundments. LV2 is a subset of LV1 and fully spatially overlapping. However, the full stream network consists of polylines with start/end nodes that do not correspond to the

[GRASS-user] Merge spatially connected features

Hi all, I have a large river network dataset (lines). Now I'd to assign unique categories to each group of connected lines that have an attribute in common. For example, my rivers are categorized based on some kind of stream order. I want to group all rivers that belong to stream order 2 and are

[GRASS-user] Error when importing gml file

Hi all, I am trying to import a *.gml file which I received from a governmental agency. However when trying v.import I get an error message (see below). I can successfully import the file in QGIS (and display the attribute table). Other gml files from the same region and from the same provider

[GRASS-user] r.sim.water simulation considering dams

Dear GRASS users, maybe someone of you has already worked with r.sim.water while considering dams as impermeable structures that retain water, i.e. create impoundments/reservoirs. I am wondering how to best specify dams within r.sim.water: There is the parameter flow_control which is a raster

[GRASS-user] Create line perpendicular to other line at specific point

Hi all, I've vector lines representing rivers and a set of points. Now, I'd like to create short lines (e.g. 150 m) that are perpendicular to the river lines and cross at the coordinates of the points vector. This is basically to create transects similar to the tool v.transect. However v.transect

Re: [GRASS-user] Import CORDEX data in GRASS

(if rotated pole is not possible with GRASS/gdal) I could try to transform using CDO and then import into a GRASS location /Johannes On Mon, Oct 21, 2019 at 10:39 PM Markus Neteler wrote: > Hi Johannes, > > On Fri, Oct 18, 2019 at 2:57 PM Johannes Radinger > wrote: > > >

[GRASS-user] Import CORDEX data in GRASS

Hi all, I was wondering if anybody is also working with CORDEX data (Regionally downscaled climate data - Europe) that are provided via the ESGF platforms as netcdf files? I want to import such files in GRASS but first want to create a location with the projection of these data. AFAIK the data

Re: [GRASS-user] v.stream.* questions

Hi Rich, AFAIK, the tool v.stream.order was developed to work with grass vector files only (no raster input needed). So, the input needed for v.stream.order is the vector lines representing the streams (lines must be connected) and points that define the outlet of each river network that

[GRASS-user] Transparent points in ps.map output

Hi all, I'm using ps.map to produce nice GRASS maps that combine some vector lines (vlines) and points (vpoints). For the points I want to distinguish between two types (presence and absence of species at a location) by using different colouring/filling: One set of vpoints should be printed as

[GRASS-user] ps.map with interactive map instructions using python

Hi all, I want to use ps.map in python to create *.eps maps from a GRASS vector line using some interactive mapping instructions, i.e. I don't want to save my mapping instructions into a *.txt file but rather provide the instructions as a python string. For example my mapping instructions:

Re: [GRASS-user] v.extract and layer/table settings with layer=2

Thank you Moritz for clarification... ...I also issued an enhancement ticket for v.extract to include a new output_layer parameter (https://trac.osgeo.org/grass/ticket/3700#ticket) /J On Wed, Nov 28, 2018 at 8:26 AM Moritz Lennert wrote: > > > Am 27. November 2018 20:51:37 MEZ schrieb

[GRASS-user] v.extract and layer/table settings with layer=2

Hi all, this is maybe I quite easy routine for somebody that is using v.extract regularly. Here my challenge: I have a vector map (initially created by v.net tools) that has several layers and associated tables to each layer. Now for further processing I want to extract only some features from

Re: [GRASS-user] Cut a floodplain polygon with 1meter wide polygons perpendicular to its river channel

It seems that would would like to have your floodplain cut based on the columns of the underlying raster-resolution (and not perpendicular to the river line as initially stated). You could use mapcalc to generate a column raster map (r.mapcalc expression="colmap = col()") which you could then

Re: [GRASS-user] Cut a floodplain polygon with 1meter wide polygons perpendicular to its river channel

Shane, have a look at the GRASS v.transect add-on: https://grass.osgeo.org/grass74/manuals/addons/v.transects.html The manual says: "Creates transect lines or quadrilateral areas at regular intervals perpendicular to a polyline." HTH Johannes On Sun, Oct 14, 2018 at 11:28 PM Shane Carey wrote:

[GRASS-user] Identify duplicate points

Hi all, this is maybe a very easy question: But how is it possible to identify duplicate points (that have the same pair of coordinates) of a set of points. I'd like to get the categories of all non-single (i.e.duplicate/triplicate etc.) points. I know that I could use the 'rmdupl' tool from

[GRASS-user] v.net manual describing 'arc_layer'

Hi all, this is probably a simple question, but it is unclear to me when reading the manual of the tool v.net: The parameter 'arc_layer' (same holds true for 'arc_node') is described as "Vector features can have category values in different layers. This number determines which layer to use. When

Re: [GRASS-user] Elevation above a river

rote: > Hey Johannes, > > Thanks for your reply. How does r.grow work if let's say the height above > the river reaches 1m at 3meters away from the river. And in an other area > it reaches the 1meter height at 2meters away from the river. Is it able to > follow that line? >

Re: [GRASS-user] Elevation above a river

To me this looks like a flooding-related question, i.e. to extract the shore lines of a river if it's water level is raised by 1m or 3m? Maybe (1) extract the raster cells of the elevation map that represents the river channel, (2) then apply r.grow and (3) then r.mapcalc to subtract the grown

Re: [GRASS-user] Not matching categories/line IDs

ector_poly_clean_tmp", tool="break,rmdupl,rmline,rmsa,rmdangle", threshold="0,0,0,0,0.001") Maybe this helps to find out what happens here. /Johannes On Thu, Sep 13, 2018 at 10:26 PM Markus Metz wrote: > > > On Thu, Sep 13, 2018 at 10:50 AM Jo

[GRASS-user] Defect v.edit tool=break with multiple coordinates

Hi all, I presumably detected an problem/bug with v.edit and its tool "break". I wanted to break a line at multiple coordinates. But the result is an overlay of multiple lines instead of the original line that is broken into segments. Here a small reproducible example using the NC dataset: from

[GRASS-user] Not matching categories/line IDs

Hi all, I am a little bit puzzled about the actual number of lines/categories of a specific vector map (river network). As discussed yesterday, the output of v.category with the report option provides a column 'count' which is a feature count. Running v.category option=report on my vector

Re: [GRASS-user] v.category report cats

Thanks for your quick answer, Moritz! I filed a ticket for this (https://trac.osgeo.org/grass/ticket/3643) /Johannes On Wed, Sep 12, 2018 at 2:58 PM Moritz Lennert wrote: > On 12/09/18 14:40, Johannes Radinger wrote: > > Hi all, I was wondering about the report output of 'v.category &

[GRASS-user] v.category report cats

Hi all, I was wondering about the report output of 'v.category option=report'. This output contains a column for each layer that is called 'count'. Does this 'count' refer to the number of vector objects or the number of unique category numbers in that layer (I guess it is the first)? This might

Re: [GRASS-user] GRASS-way to delete points from a map and attribute table

ina.no> wrote: > So, removing records from the attribute table would be a job for > db.execute of db.select I guess… > > > > *From:* Johannes Radinger [mailto:johannesradin...@gmail.com] > *Sent:* torsdag 9. november 2017 10.17 > *To:* Stefan Blumentrath <stefan.b

Re: [GRASS-user] GRASS-way to delete points from a map and attribute table

> > https://grass.osgeo.org/grass72/manuals/v.edit.html > > Cheers > > Stefan > > > > *From:* grass-user [mailto:grass-user-boun...@lists.osgeo.org] *On Behalf > Of *Johannes Radinger > *Sent:* torsdag 9. november 2017 10.04 > *To:* GRASS user list <grass-user@l

Re: [GRASS-user] t.register: register 4 consecutive rasters maps per month for long time series

m means, > you can use r.series [1] with every 48 maps for yearly means and, the whole > list of maps for the long term mean. > > hth, > Vero > > [0] https://grass.osgeo.org/grass72/manuals/t.register.html > [1] https://grass.osgeo.org/grass72/manuals/r.series.html > >

[GRASS-user] GRASS-way to delete points from a map and attribute table

Hi, what would be the most GRASS-like way to delete several points from a points vector map. I know the specific cats of the points I would like to delete (10 of 20) and I want to delete the points from the map and, of course, also from the corresponding attribute table. I already tried

[GRASS-user] t.register: register 4 consecutive rasters maps per month for long time series

Hi, I am very new to the temporal functionalities of GRASS. Specifically I have a netCDF file that contains many layers where each layer represents a raster map of a drought index (http://monitordesequia.csic.es/map/) for a timepoint of a timeseries. The time series is from 1/1961 to 12/2015 with

[GRASS-user] Create a point along a line in a specified distance from another point

Hi, is there an easy straight forward command to create a new point on a line that is in a specified distance (e.g. 1000 m) from another point (on the same line)? Just an example (of course from river science): I have a sampling site on a river line. I'd like to use these new point(s) to break

Re: [GRASS-user] Extract subnet from grass vector network based on selected nodes

extract I guess to get > your subnetwork… > > > > Or if you continue with network analysis you could just use: > > http://igraph.org/python/doc/igraph.GraphBase-class.html#subgraph_edges > > > > Hope that is somehow useful… > > > > Cheers > > S

[GRASS-user] Version GRASS for Mac OSX

Hi all, the binaries of GRASS for Mac OSX available via http://grassmac.wikidot.com/downloads are already older than a year (from 14.09.2016). Are there any plans to upload and provide a newer version of GRASS as binary download in near future? Best regards, Johannes

Re: [GRASS-user] Extract subnet from grass vector network based on selected nodes

;> On Wed, Sep 27, 2017 at 4:07 PM, Moritz Lennert < > mlenn...@club.worldonline.be <mailto:mlenn...@club.worldonline.be> > <mailto:mlenn...@club.worldonline.be <mailto:mlenn...@club.worldonline.be>>> > wrote: > >> >> > > >> >>

[GRASS-user] Extract subnet from grass vector network based on selected nodes

Hi, I have a GRASS vector network that represents a river network (with many first order tributaries) and that has additional connected nodes that represent sampling sites. I'd like to extract a minimum subnetwork of the full network that still connects a set of selected nodes (e.g. identified

Re: [GRASS-user] Merge lines into common categories based an attribute information

calculating sinuosity. /J On Sun, Sep 17, 2017 at 3:53 PM, Moritz Lennert < mlenn...@club.worldonline.be> wrote: > > > Le 17 septembre 2017 15:16:02 GMT+02:00, Johannes Radinger < > johannesradin...@gmail.com> a écrit : > >Hi, > > > >how can I merge mult

[GRASS-user] Merge lines into common categories based an attribute information

Hi, how can I merge multiple lines into one category based on some attribute information. For example I have a line vector consisting of 10 lines and several lines belong to three different groups (lines belonging to the same group share the same value in a specific attribute column). Now I'd

[GRASS-user] Running GRASS within R (Rstudio) on MacOS

Hi, I try to run GRASS73 within R (actually R-studio) using the rgrass7 library on MacOS 10.12.6. I followed the instructions stated in the wiki ( https://grasswiki.osgeo.org/wiki/R_statistics/rgrass7), however some of the commands were not working for me. For instance, "grass70 --config path"

[GRASS-user] Image segmentation to separate object from background

Hi GRASS users, has anyone of you tried to use GRASS image tools (e.g. segmentation etc) to identify an object in a picture. For example I have multiple photos of fish and would like to separate the fish from its background in an automatized way. The images look like:

[GRASS-user] Single centroid/point for multiple polygons sharing the same category

Hi GRASS users, I am working on a polygon map of France that shows the single areas belonging to a specific postal code. I want to extract for each area the centroid (using v.extract), transform it to a point map (using v.type) and finally add the X and Y coordinates to the attribute table for

Re: [GRASS-user] Obtain land use within 50 km buffers around vector points

AS CLC_forest, SUM(CLC_others) AS CLC_others, SUM(CLC_intertidal_coastal) AS CLC_intertidal_coastal, SUM(CLC_water_bodies) AS CLC_water_bodies, SUM(CLC_sea) AS CLC_sea FROM {} AS A LEFT JOIN {} AS B ON A.buffer_cat_2=B.buffer_cat_2 GROUP BY buffer_cat_1".format(clc_summary_table,correspondance

[GRASS-user] Obtain land use within 50 km buffers around vector points

Hi, I have a vector map of 50 km buffers surrounding 4000 points in entire Europe. The single buffers where calculated using v.buffer with the -t flag. Thus, many of my buffers overlap (i.e. geometries of buffers are not merged but split up). Now, I want to obtain for each 50 km buffer the

[GRASS-user] Indentify cell with maximum value per category of a cover map

Hi, I'd like to identify the cell of maximum flow accumulation for each subbasin. For example, I have a cover raster map representing all my subbasins (each subbasin has its own cat value). Additionally, I have a flow accumulation map. Now, I'd like to get a map with cells that are actually

Re: [GRASS-user] Column names of sqlite attribute tables while using v.import

Thank you Moritz, that's what I thought. So I'll try the 'columns' option in v.in.ogr to rename the columns during import. /johannes On Tue, May 16, 2017 at 11:50 PM, Moritz Lennert < mlenn...@club.worldonline.be> wrote: > > > Le 16 mai 2017 22:16:03 GMT+02:00, Jo

[GRASS-user] Column names of sqlite attribute tables while using v.import

Hi all, I tried to import a polygon shape, specifically the HydroBASINS from hydrosheds.org however I got following error: v.import input=/.../Hydrosheds/hybas_na_lev00_v1c/hybas_na_lev00_v1c.shp layer=hybas_na_lev00_v1c output=hybas_na_lev00_v1c WARNING: All available OGR layers will be

Re: [GRASS-user] [snap pour point for r.water.outlet]

Sorry I ment r.stream.snap... Original message From: Ang Sherpa <angsherpa...@gmail.com> Date:30/03/2017 07:32 (GMT+01:00) To: Johannes Radinger <johannesradin...@gmail.com> Cc: Helmut Kudrnovsky <hel...@web.de>, GRASS user list <grass-user@lists.os

Re: [GRASS-user] [snap pour point for r.water.outlet]

Dear Ang Sherpa, r.stream.snap is an add-on. So you first need to install it using e.g g.extension. Have you installed r.stream.order before calling it? /j Original message From: Ang Sherpa Date:30/03/2017 05:34 (GMT+01:00) To: Helmut Kudrnovsky

Re: [GRASS-user] Change elevation value for entire DEM raster cell [UPDATE 2]

Rich, Maybe also this conversation on setting a specific raster value using r.mapcalc might help:  http://osgeo-org.1560.x6.nabble.com/Set-raster-value-for-specific-cell-td5052853.html /j Original message From: Rich Shepard Date:23/03/2017 22:07

[GRASS-user] Auto-detection of GCP for rectifying image

Hi GRASS users! I have an image of a small object (e.g. photograph of blue square 15x5cm) which I want to automatically georeference. In addition, the image contains 6 points where I know the exact coordinates, i.e. these points can be use as ground control points (GCP) for rectifying the image.

Re: [GRASS-user] r.flow: define contributing area

Rich, Maybe have a look at r.lake:  https://grass.osgeo.org/grass70/manuals/r.lake.html. This module fills the area upstream a dam or a blocked stream. Cheers, Johannes Original message From: Thomas Adams Date:04/10/2016 04:32 (GMT+01:00) To: Rich

[GRASS-user] Update categories/attribute table after v.clean with option break

Hi, I have a vector line map (river network) which I want to clean and break the lines at intersections. Therefore, I use the tool v.clean with the option break. This of course increases the number of vector lines as several lines got broken into two or more separate lines. Now there is a

Re: [GRASS-user] Select specific layer when copying vector map

Indeed, v.extract using the layer option did the job! Thank you Moritz! Probably, I can then use v.category option="chlayer" to change the layer Nr. 3 to Nr. 1!? On Tue, Aug 9, 2016 at 2:23 PM, Moritz Lennert <mlenn...@club.worldonline.be > wrote: > On 09/08/16 13:48, Joha

[GRASS-user] Select specific layer when copying vector map

Hi, I'd like to create from an existing GRASS vector map that has 3 layers (with three attribute tables) a new vector map that has only one layer and its associated attribute table (here the information related to layer 3). What I am doing so far (in python): # First copy the map to a new one

[GRASS-user] Updating v.net-report results directly to attribute table

Hi all, is there an easy way to transfer the information reported by v.net (operation=report) directly to the attribute table of the vector network (arc-table)? The approach I am using so far is to create an temporary database table (using python) where I store the output from v.net

[GRASS-user] New GRASS GIS addon: v.stream.order

Dear GRASS Users, As some of you might already have noticed, there is a new GRASS GIS add-on available, called v.stream.order: https://grass.osgeo.org/grass70/manuals/addons/v.stream.order.html This module computes various types of stream order (Strahler, Shreve, Scheidegger Drwal) of stream

Re: [GRASS-user] Average raster value of the upslope contributing area,

Hi Pierluigi, I think here you could use the tool r.watershed ( https://grass.osgeo.org/grass70/manuals/r.watershed.html) with your underlying DEM as input. This tool will provide you an accumulation map with increasing values in downstream direction. I am not exactly sure if this is your task,

Re: [GRASS-user] Unique IDs for network segments

Moritz, thank you! Your proposed approach solved my problem! On Thu, Dec 3, 2015 at 3:04 PM, Moritz Lennert <mlenn...@club.worldonline.be > wrote: > On 03/12/15 14:19, Johannes Radinger wrote: > >> Hi Moritz, >> >> The two-step v.category approach (del +

Re: [GRASS-user] Unique IDs for network segments

nn...@club.worldonline.be > wrote: > On 02/12/15 15:40, Johannes Radinger wrote: > >> Hi all, >> >> I built a GRASS7 network using v.net <http://v.net>: >> >> First, I used the 'connect' option to snap some points to my network >> (using the -s fl

[GRASS-user] Unique IDs for network segments

Hi all, I built a GRASS7 network using v.net: First, I used the 'connect' option to snap some points to my network (using the -s flag) Second, I used that first network and also added all nodes with the 'nodes' option to get a final network. With the 'report' option of v.net I can report all

Re: [GRASS-user] Error starting GRASS71 trunk: module catalog missing?

Thanks Martin, It seems that deleting the installation directory helped! /Johannes On Tue, Nov 24, 2015 at 2:58 PM, Martin Landa <landa.mar...@gmail.com> wrote: > Hi, > > 2015-11-24 14:39 GMT+01:00 Johannes Radinger <johannesradin...@gmail.com>: > > from d

[GRASS-user] Error starting GRASS71 trunk: module catalog missing?

Hi, I tried to start the most recent (SVN up) self-compiled version of GRASS71 on a Ubuntu-machine. GRASS starts and I can selected a location/mapset from the GUI but after the selection I get following error and the GUI does not start: Launching GUI in the background, please wait... GRASS

Re: [GRASS-user] river and dams in netwok analysis

Hi Etienne, Maybe you could use v.net.allpairs with negative/positve costs for up- and downstream direction and then select the smallest positive and negative distance for each dam to get the two neighbours. Here what I found in the manual of v.net.distance which might also be useful for your

Re: [GRASS-user] Problems running g.gui.gcp

On Fri, Oct 23, 2015 at 10:17 AM, Moritz Lennert < mlenn...@club.worldonline.be> wrote: > On 22/10/15 16:25, Johannes Radinger wrote: > >> Dear GRASS users, >> >> Today I wanted to try the georectification tool g.gui.gcp in GRASS 7 for >> the first time. I was

Re: [GRASS-user] Problems running g.gui.gcp

On Fri, Oct 23, 2015 at 11:19 AM, Moritz Lennert < mlenn...@club.worldonline.be> wrote: > On 23/10/15 10:56, Johannes Radinger wrote: > >> >> >> On Fri, Oct 23, 2015 at 10:17 AM, Moritz Lennert >> <mlenn...@club.worldonline.be <mailto:mlenn...@club.worldo

Re: [GRASS-user] Problems running g.gui.gcp

On Fri, Oct 23, 2015 at 1:37 PM, Anna Petrášová <kratocha...@gmail.com> wrote: > > On Oct 23, 2015 6:04 AM, "Johannes Radinger" <johannesradin...@gmail.com> > wrote: > > > > > > > > On Fri, Oct 23, 2015 at 11:19 AM, Moritz Lennert < >

[GRASS-user] Problems running g.gui.gcp

Dear GRASS users, Today I wanted to try the georectification tool g.gui.gcp in GRASS 7 for the first time. I was basically following the description on the GRASS wiki (https://grasswiki.osgeo.org/wiki/Georeferencing). In particular, I have a photograph from my camera with known reference points

Re: [GRASS-user] Need help in reprojecting raster images

Hi Uttam, if I understand your problem correctly, you can use the -p flag of r.proj to print input map's bounds in the current projection. With this information you can set the region boundaries appropriately using g.region. Afterwards you can 'import' your raster into your laea-projection

Re: [GRASS-user] Stream extraction from pseudo-elevation map using r.watershed

(e.g. 1000) before using it with r.watershed. Cheers, Johannes On Thu, Jul 2, 2015 at 2:18 PM, Johannes Radinger johannesradin...@gmail.com wrote: Hi, I tried something that was working at least several months ago, but somehow r.watershed has changed. What I am trying to do is to extract

[GRASS-user] Stream extraction from pseudo-elevation map using r.watershed

Hi, I tried something that was working at least several months ago, but somehow r.watershed has changed. What I am trying to do is to extract a stream network from a pseudo elevation map. This pseudo elevation map is created by getting the distance of each stream cell from to outlet in upstream

[GRASS-user] Dataset file for add-on

Hi, I'd like to use a HDF5 datafile (multidimensional matrix with parameter values) to be used and called from an python add-on. Thus, I thought about providing the data file together with the python add-on an load the data via a relative path (relative to the path of the add-on). Is it generally

[GRASS-user] Copy raster map from one unprojected location to another

Hi, I'd like to import a raster map from an unprojected location (XY-location) to my actual location (also unprojected XY). Because these are two different locations I can't use g.copy. So I tried r.proj. Here I get the error that the projection is unknown, which makes sense. Of course, I could

Re: [GRASS-user] Stream order for vector networks

provide you with some sample code in R if that would be of interest. Cheers Stefan *From:* grass-user-boun...@lists.osgeo.org [mailto: grass-user-boun...@lists.osgeo.org] *On Behalf Of *Johannes Radinger *Sent:* 18. mars 2015 10:30 *To:* GRASS user list *Subject:* [GRASS-user] Stream

Re: [GRASS-user] Copy raster map from one unprojected location to another

wrote: On Thu, Mar 19, 2015, at 11:25, Roy wrote: Hi, Il 19/03/2015 11:04, Johannes Radinger ha scritto: Of course, I could export the map to geotiff and import them in the other location; IMHO this is the simpler way, Bad idea. You will lose all associated metadata: colour

Re: [GRASS-user] Copy raster map from one unprojected location to another

(filename, r).readlines() IOError: [Errno 2] No such file or directory: 'PROJ_INFO' On Thu, Mar 19, 2015 at 12:26 PM, Nikos Alexandris n...@nikosalexandris.net wrote: On 19.03.2015 13:15, Johannes Radinger wrote: Thanks for your ideas. I tried the pack-unpack approach. However, it seems that I

[GRASS-user] Stream order for vector networks

Hi, I am interested in calculating stream order (Strahler, Shreve) for stream networks based on GRASS vector networks (e.g. created by v.net) without using the raster approach of the GRASS add-on r.stream.order. Therefore, I came across following paper: Gleyzer, A., Denisyuk, M., Rimmer, A. and

Re: [GRASS-user] sampling points in a grid

Maybe you can use a database selection approach for example in sqlite. First update your points with your grid cell ID and then randomly select 7 rows from the selection where the grid cell ID == 42. In sqlite there exists the function random(). Here you can find an approach to select random rows

[GRASS-user] Choosing mapset from command line

Hi, I tried to access my machine with GRASS GIS via SSH and wanted to get the text based welcome screen of GRASS 7 (RC in Ubuntu) to select a mapset. Usually I am starting the GUI selection menu so I am not familiar with the text based startup. I thought that the command grass70 -gtext should

Re: [GRASS-user] Choosing mapset from command line

On Fri, Feb 13, 2015 at 4:58 PM, Johannes Radinger johannesradin...@gmail.com wrote: Maybe I expressed myself unclear or mixed things up: I want to get the text based menu to select e.g. a location and a mapset, which I think is called text-based location wizzard [1]. I mixed that up

Re: [GRASS-user] Choosing mapset from command line

/grass70/manuals/grass7.html On Fri, Feb 13, 2015 at 3:43 PM, Moritz Lennert mlenn...@club.worldonline.be wrote: On 13/02/15 15:33, Johannes Radinger wrote: Hi, I tried to access my machine with GRASS GIS via SSH and wanted to get the text based welcome screen of GRASS 7 (RC in Ubuntu

[GRASS-user] Transform single raster cells to 3d bars

Hi, I have a 2D raster map with scattered single cells (e.g. created by r.random.cells) with lots of empty space (NaNs) between the cells. Each cell has a value that represents actually the height of this cell. Now I'd like to transform the 2D raster cells into a 3D raster where each cell should

Re: [GRASS-user] temporary files of grass.script.array

Hi, I just tested in trunk, and it works: The temporary file that is create when assigning a raster to an python array is deleted once e.g. the python session is closed again. So this would be a nice thing also to have implemented in the next RC of GRASS7. /johannes On Tue, Feb 3, 2015 at

[GRASS-user] r.in.onearth in GRASS 7

Hi, I just try to follow some examples r3.out.vtk to get some nice 3D plots in paraview. http://grass.osgeo.org/grass70/manuals/r3.out.vtk.html In the manual there is an example that uses the r.in.onearth add-on (http://trac.osgeo.org/grass/browser/grass-addons/grass6/raster/r.in.onearth ) that

[GRASS-user] temporary files of grass.script.array

Hi, In GRASS 7 (RC1) python it is possible to export a raster map to an array which can be used e.g. by Numpy. For example: import grass.script.array as garray x1 = garray.array() x1.read(my_raster) It seems that this creates a temporary file in the .tmp folder of the respective location. I am

Re: [GRASS-user] how to change all cats in a vector line

Hi Michael, I have not read the entire thread, but maybe follwowing might work if you want to extract this upstream-downstream watercourse in raster format: 1) Use r.cost with your outlet as starting point and your cell resolution as cost surface. So the output is a pseudo-elevation map with

Re: [GRASS-user] A question before I embark on a programming exercise

Hi Tom, what comes into my mind is following approach (not tested): 1) calculate flow direction using r.watershed for each pixel of your map 2) transforming raster to vector points (r.to.vect), adding some columns for your results (X, Y, ID). This step should provide you already with unique IDs

[GRASS-user] Overlapping points in ps.map

Hi, I am trying to use ps.map to create a map of many vector points which should be colored based on the values (0-1) in a specific column in the attribute table. If I understood the manual correctly, I need to create an rgbcolumn indicating the color for each point. Is there an example

[GRASS-user] Upper case letters for defining columns via v.db.addtable

Hi, when I want to specify a column in v.db.addtable, all upper case letters in the column name get converted to lower case. I'm working on GRASS 71 (trunk). Is this an intended behaviour? v.db.addtable map=myMap table=test layer=3 columns=My_COL INT cheers, /Johannes

[GRASS-user] Categories in new layers of vector maps

Hi, I am slightly confused about creating new layers for vector maps and adding a category value. I want to add a layer and subsequently I want to use v.what.rast to query a raster map and update that layer. Hence, I need categories in the new layer. However, it seems that adding a new table and

Re: [GRASS-user] my first script

Hi,  It seems to raise an error related to the flag(s) you are specifying in oneof your grass modules you are using in your script. Flags in grass python scripts should be used like: flags=a or for two flags: flags=ab. Note this does not need a - like in the bash command. Best, /johannes

[GRASS-user] Null-value in grass python array

Hi, I am using the GRASS-numpy functionality [1] to read a python numpy.array from a GRASS raster map (GRASS7): import grass.script.array as garray a = garray.array() a.read(map) Here [2] it says that also the null value can be specified. So what is the actual null value in GRASS so that

Re: [GRASS-user] raster exchange between GRASS and R with nodata

On 4 September 2014 14:27, Johannes Radinger johannesradin...@gmail.com wrote: Hi all, of course it is possible to load the raster maps directly via spgrass6. However, we use this work flow also to exchange some of the maps between different users (e.g. via email

  1   2   3   4   >