Hi Markus,
Great, thanks!
Feel free to add those outputs for which I used GRASS:
https://doi.org/10.31223/osf.io/vqgx4
https://doi.org/10.17605/OSF.IO/38W4T
https://dx.doi.org/10.3390/w10020207
Laurent
Le dim. 9 déc. 2018 à 10:24, Markus Neteler a écrit :
> Hi,
>
> for fun, to easily track
a écrit :
> Laurent C. wrote
> > Hello,
> >
> > I am extracting cross-sections from raster maps using r.transect (E,N,
> > azimuth, distance).
> > I would like to create maps where I can show where those xsections are
> > taken using a vector line.
> >
Hello Frank,
It is not very straightforward, but it is possible.
Here is an example on how I do it:
import os
import grass.script as gscript
from grass.pygrass.gis.region import Region
# Set general env
os.environ['GRASS_RENDER_IMMEDIATE'] = "cairo"
os.environ['GRASS_RENDER_FILE_COMPRESSION'] =
Hello,
I am extracting cross-sections from raster maps using r.transect (E,N,
azimuth, distance).
I would like to create maps where I can show where those xsections are
taken using a vector line.
It seems that v.edit only works with absolute coordinates, not bearing and
distance.
Is there a way
Hi Rich,
In the debate kriging vs deterministic rainfall interpolation, the
following article may be of interest:
https://doi.org/10.1016/j.proeng.2016.07.595
Laurent
Le lun. 24 sept. 2018 à 16:32, Rich Shepard a
écrit :
>I want to determine whether GRASS or R is best suited to
>
v.what.vect is doing exactly what I was searching for.
The last module on the list ;)
Laurent
Le sam. 28 juil. 2018 à 18:04, Laurent C. a écrit :
> Hello all,
>
> I have two vectors maps:
> - Map A with only points
> - Map B with polygons representing continents (with a name col
Hello all,
I have two vectors maps:
- Map A with only points
- Map B with polygons representing continents (with a name column)
I want to add a column in map A that contains the continent name on which
the point is located.
I thought of using v.overlay, but it doesn't work with points.
Which
Ken,
You need to allow the access of those mapsets with g.mapsets (plural).
Best regards,
Laurent
2017-12-05 6:56 GMT-06:00 Ken Mankoff :
> Dear GRASS List,
>
> I'm working with some STRDS and STVDS and having trouble accessing an
> STRDS located in a different mapset.
>
> I
Hello all,
I have a similar problem with Fedora 27:
$ grass72 --config path
ERROR: Please install the GRASS GIS development package
Exiting...
$ dnf list installed *grass-devel*
Installed Packages
grass-devel.x86_64 7.2.1-2.fc27 @fedora
I believe it is a
2017-09-21 14:22 GMT-05:00 Markus Metz :
>
> The JPL/NASA website also links to
> http://dds.cr.usgs.gov/srtm/
> where you get the JPL/NASA SRTM products version 1 and 2.1
>
> JPL/NASA SRTM version 3 (completely void-filled version 2.1) is available at
>
2017-09-20 4:07 GMT-05:00 Markus Metz <markus.metz.gisw...@gmail.com>:
>
>
> On Tue, Sep 19, 2017 at 9:10 PM, Markus Metz <markus.metz.gisw...@gmail.com>
> wrote:
>>
>>
>>
>> On Tue, Sep 19, 2017 at 5:53 PM, Laurent C. <lrn...@gmail.com> wr
On Tue, Sep 19, 2017 at 1:41 AM, Laurent C. <lrn...@gmail.com> wrote:
>
> Hello all,
>
> Open-access global DEMs like SRTM are a staple in our GIS community.
> Last year the JAXA released a new product, the ALOS World 3d 30m
> (AW3D30).
> Using GRASS and itzï we'
Hello all,
Open-access global DEMs like SRTM are a staple in our GIS community.
Last year the JAXA released a new product, the ALOS World 3d 30m
(AW3D30).
Using GRASS and itzï we've carried out an assessment of AW3D30 from
the point of view of flood modelling, and compared it to lidar, SRTM
and
Hello Marco,
Itzï (www.itzi.org) is built upon GRASS and employs a partial inertia
simplification of the Saint-Venant equations.
It is certainly not as hydraulically accurate as the finite-volume
approach of ANUGA, but it gives quite good results as a flood model,
and is likely to be much faster.
-08-01 10:44 GMT-05:00 Laurent C. <lrn...@gmail.com>:
> Hello all,
>
> The version 17.7 of the Itzï flood model is out. It fixes some bugs
> and add new features.
> The most important news is that now it includes the SWMM drainage
> model, allowing coupled surface-drainage
Hello all,
The version 17.7 of the Itzï flood model is out. It fixes some bugs
and add new features.
The most important news is that now it includes the SWMM drainage
model, allowing coupled surface-drainage simulations. The results of
the drainage model are exported as GRASS STVDS.
The full
time
> =='2001-01-01'")? Have you tried that?
I did, without success. Actually, I wrote it the exact same way as in
v.db.select, and it worked.
Regards,
Laurent
>
> Best,
> Vero
>
>
> El 14 jun. 2017 5:03 p.m., "Laurent C." <lrn...@gmail.com> escribió:
Hello,
I'm running into a problem with t.rast.what when trying to select a
specific point with the "where" option:
/usr/lib/grass72/scripts/t.rast.what points=stations@chiapas
strds=boqueron_60m_qy@chiapas_results where=name=='Boquerón'
Traceback (most recent call last):
File
Hello,
I just discover this thread. It seems to be the same issue I've
reported two months ago:
https://trac.osgeo.org/grass/ticket/3338
Best regards,
Laurent
2017-03-22 16:44 GMT-06:00 Markus Neteler :
> On Wed, Mar 22, 2017 at 8:21 PM, Andres Solarte
>
Hi Rich,
When I want to modify a raster, I usually draw a vector and use
v.to.rast, then r.mapcalc if necessary.
Regards,
Laurent
2017-06-08 18:57 GMT-05:00 Rich Shepard :
> On Thu, 8 Jun 2017, Mark Seibel wrote:
>
>> Maybe if-then statements with mapcalc are
Hello,
Is there a module that simplifies the reprojection of a complete STDS
from a location to another?
Regards,
Laurent
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user
Hello all,
A new open-access, peer-reviewed journal article has been published
about the use of GRASS and Itzï for the dynamic simulation of floods:
Courty, L. G., Pedrozo-Acuña, A., and Bates, P. D.: Itzï (version
17.1): an open-source, distributed GIS model for dynamic flood
simulation,
Hello Joseph,
There is a solution, but it is not straight forward:
https://grasswiki.osgeo.org/wiki/Working_with_GRASS_without_starting_it_explicitly
Another option is to start your editor from the GRASS command line.
Regards,
Laurent
2017-04-04 9:12 GMT-05:00 Joseph Kariuki
2017-03-30 16:02 GMT-06:00 Rich Shepard <rshep...@appl-ecosys.com>:
> On Thu, 30 Mar 2017, Laurent C. wrote:
>
>> In that case, you will not see the accumulation upstream the culvert that
>> should normally occur if the pipe is actually higher than the ground.
>
>
2017-03-30 12:38 GMT-06:00 Rich Shepard <rshep...@appl-ecosys.com>:
> On Thu, 30 Mar 2017, Laurent C. wrote:
>
>> You can create a new vector map and create a line that link the two
>> points. However, r.carve will interpolate the elevation between the start
&g
2017-03-30 12:09 GMT-06:00 Rich Shepard :
> On Thu, 30 Mar 2017, Rich Shepard wrote:
>
>> Would the process be to use v.edit on a copy (called 'stream') of the
>> culvert point map to generate a straight line between the two points? That
>> would provide the vector input
2017-03-30 12:00 GMT-06:00 Rich Shepard <rshep...@appl-ecosys.com>:
> On Thu, 30 Mar 2017, Laurent C. wrote:
>
>> If you want to 'open' the culvert, you'll have to cut through the road.
>> You can do this with r.carve, for example.
>
>
> Certainly worth doing.
Rich,
If I understand well, you have a 'raw' DEM and you want to open the
culvert in that DEM.
If it is an under-road culvert, your 'raw' DEM will show the road
elevation, and the culvert will not be shown. i.e the culvert is
blocked.
If you want to 'open' the culvert, you'll have to cut through
2017-03-23 12:49 GMT-06:00 Rich Shepard <rshep...@appl-ecosys.com>:
> On Thu, 23 Mar 2017, Laurent C. wrote:
>
>> Is this t.info output taken from the STRDS created by Itzï?
>
>
> Laurent,
>
> I assume so since the registered maps are Itzi output.
Itzï is cr
Rich,
Is this t.info output taken from the STRDS created by Itzï?
The current version of Itzï does not record the initial state.
So if you tell Itzï that the simulation starts at 2013-11-15 00:00:00
and the record_step is 24:00:00, then the first registered map will be
at 2013-11-16 00:00:00.
2017-03-14 16:01 GMT-06:00 Rich Shepard <rshep...@appl-ecosys.com>:
> On Mon, 13 Mar 2017, Laurent C. wrote:
>
>> The water level of 3m seems to be due to accumulation on the edge of the
>> domain. The water flows according to topography. If there is a lower
2017-03-13 15:06 GMT-06:00 Rich Shepard <rshep...@appl-ecosys.com>:
> On Mon, 13 Mar 2017, Laurent C. wrote:
>
>> Additionally, you can check the statistics CSV file to see if there is
>> any numerical instability that could artificially create volume.
>
>
> Lauren
2017-03-13 10:31 GMT-06:00 Laurent C. <lrn...@gmail.com>:
> Copy of my off-list answer.
>
> 2017-03-13 10:30 GMT-06:00 Laurent C. <lrn...@gmail.com>:
>> Rich,
>>
>> Answers below.
>>
>> 2017-03-13 10:16 GMT-06:00 Rich Shepard <rshep...@appl-eco
2017-03-13 10:26 GMT-06:00 Rich Shepard :
> On Mon, 13 Mar 2017, Rich Shepard wrote:
>
>> One model run has completed and I have several questions for which I do
>> not find answers in either the grass or itzi docs.
>
>
> 4.) The water surface elevation (wse) output
Copy of my off-list answer.
2017-03-13 10:30 GMT-06:00 Laurent C. <lrn...@gmail.com>:
> Rich,
>
> Answers below.
>
> 2017-03-13 10:16 GMT-06:00 Rich Shepard <rshep...@appl-ecosys.com>:
>> One model run has completed and I have several questions for which I
Hi Nikos,
What I would like is the possibility to set a percentage of opacity. But if
I understand well, it is what is doing pnmcomp.
Thanks for your appreciation.
Regards,
Laurent
On 11 Mar 2017 03:35, "Nikos Alexandris" <n...@nikosalexandris.net> wrote:
Laurent C.:
Hell
Hello,
I'm using a Python script to generate maps for publication by setting
GRASS_RENDER_FILE and running the different d.* commands. Is it
possible to add opacity to some of the layers?
Regards,
Laurent
___
grass-user mailing list
Hello,
You should make sure that the point you give to r.water.outlet is
located on a main stream. The accumulation map from r.watershed is
useful for this task.
If you're not doing so, you might end up with a tiny watershed.
Regards,
Laurent
2017-03-06 16:50 GMT-06:00 Rich Shepard
> You need to specify the name of the stvds in t.register.
>
> Am 24.02.2017 02:48 schrieb "Laurent C." <lrn...@gmail.com>:
>>
>> Hello,
>>
>> I created an empty STVDS that looks fine:
>>
>> t.info input=water_levels_20150801@kol
Hello,
I created an empty STVDS that looks fine:
t.info input=water_levels_20150801@kolkata type=stvds
+ Space Time Vector Dataset -+
||
+ Basic
2017-02-23 13:13 GMT-06:00 Rich Shepard :
> On Thu, 23 Feb 2017, Sören Gebbert wrote:
>
>> This seems to be a bug.
>
>
> Sören,
>
> Just my luck! :-(
>
>> However, in your last email you performed this task in the topology
>> mapset.
>
>
> That's true.
>
>> Why do you
Rich,
If you want to create a raster map with an uniform value, just use mapcalc:
r.mapcalc exp='manning=0.03'
And that's it.
Laurent
2017-02-22 16:40 GMT-06:00 Rich Shepard :
> For the defined region I want to create a map of the Manning's n friction
> value
2017-02-22 9:28 GMT-06:00 Rich Shepard <rshep...@appl-ecosys.com>:
> On Tue, 21 Feb 2017, Laurent C. wrote:
>
>> In your case, it seems you want to have uniform rainfall. So you can set
>> the region, and then use mapcalc like so: r.mapcalc
>> exp='rainfall_2013115=0.
Rich,
I would use v.in.ascii, v.to.rast and then r.patch.
Laurent
2017-02-21 13:42 GMT-06:00 Rich Shepard :
> What is the appropriate r.mapcalc command to change the value of the
> raster cell at a specified coordinate? Cell resolution is 1m x 1m.
>
> For a
Rich,
I'm not familiar with the 'file' option in r.mapcalc. Maybe removing
the quotes will do?
Laurent
2017-02-21 12:31 GMT-06:00 Rich Shepard <rshep...@appl-ecosys.com>:
> On Tue, 21 Feb 2017, Laurent C. wrote:
>
>> So you can set the region, and then use mapcalc like so
2017-02-15 18:38 GMT-06:00 Rich Shepard <rshep...@appl-ecosys.com>:
> On Wed, 15 Feb 2017, Laurent C. wrote:
>
>> The h_max output map of Itzï will correspond to the maximum water level
>> reached during the simulation.
>
>
> Laurent,
>
> I did this to re
Hello Rich,
The h_max output map of Itzï will correspond to the maximum water
level reached during the simulation.
To simulate the blockage, you might need to change the DEM locally.
But in the case of a under-road culvert, just leaving the measured
elevation of the road might be enough. The
sion not found
>
> Is it intended to install with g.extension?
>
> thanks much!
> Vero
>
> 2017-01-21 1:04 GMT+01:00 Laurent C. <lrn...@gmail.com>:
>>
>> Hello all,
>>
>> I wrote a module that automatically download, import and register GPM
>> IMERG
people working with those data.
Regards,
Laurent
2015-03-19 10:50 GMT-06:00 Laurent C. <lrn...@gmail.com>:
> Hi all,
>
> I'm coming back to that subject because the GPM datas are finally available
> since December 2014. The original data are stored in HDF5, which is
> supported
Hello Sören,
The problem occurs when Itzï is calling:
"tgis.SpaceTimeRasterDataset(name).is_in_db()"
tgis.init() is called just before. "name" is a a formatted map or
strds ID (i.e including mapset)
It is difficult to understand what is happening because I can't
reproduce the problem.
Regards,
2016-10-12 13:52 GMT-05:00 Rich Shepard :
> On Wed, 12 Oct 2016, Rich Shepard wrote:
>
>> Hmm-m-m. I'll read the v.surf.idw manual page.
>
>
> I don't see how v.surf.idw is helpful here. The current map has only
> category values and the label for each station.
Rich,
You can interpolate the data using v.surf.idw to get a raster map.
Literature says that a power parameter superior to 2 (between 3 to 5)
might get better result for rainfall.
Laurent
2016-10-12 12:20 GMT-05:00 Rich Shepard <rshep...@appl-ecosys.com>:
> On Wed, 12 Oct 2016,
Hello Rich,
Sure, raster time series of rainfall intensities is a very good way to
take advantage of GRASS temporal framework in hydrologic modelling.
Regards,
Laurent
2016-10-12 10:45 GMT-05:00 Rich Shepard :
> On Tue, 11 Oct 2016, Rich Shepard wrote:
>
>> I've no
hep...@appl-ecosys.com>:
> On Tue, 11 Oct 2016, Laurent C. wrote:
>
>> - Your rainfall rate is very small.
>
>
> Laurent,
>
> Changed rain to 10 mm/h, time to 2 hrs with 5 min steps. Specified mapset
> for DEM map (be_cropped) as it's not in the analyses mapset.
>
11:18 GMT-05:00 Rich Shepard <rshep...@appl-ecosys.com>:
> On Mon, 10 Oct 2016, Laurent C. wrote:
>
>> Note that Itzï does not need an output to run. If you set a simulation
>> with only the DEM, an uniform friction coefficient and an uniform
>> rainfall, th
Hello Rich,
According to the r.watershed manual:
"Output accumulation map contains the absolute value of each cell in
this output map is the amount of overland flow that traverses the
cell. This value will be the number of upland cells plus one if no
overland flow map is given. If the overland
Hello Rich,
When using r.water.outlet you need to be careful setting the outlet on the
exact flow line. Using the flow accumulation map is useful for that task.
If the outlet point is just one cell on the side, the generated watershed
could be tiny.
Hope this helps.
Laurent
El oct. 2, 2016 6:17
Hi Rich,
The rainfall unit indeed mm/h. The units used for all entry and output maps
are listed in the user manual on the website. In the tutorial, the example
output map is in meters after 100 minutes of simulation (map #20 with a
record step of 5 minutes). Another way to know a map's time stamp
with the units given in the manual.
Regards,
Laurent
El sept. 26, 2016 18:31, "Rich Shepard" <rshep...@appl-ecosys.com> escribió:
> On Thu, 22 Sep 2016, Laurent C. wrote:
>
> If it is the case, you might give Itzï a try:
>> https://www.itzi.org/
>>
>
>
Hello Rich,
For the infiltration, if you can determine the soil texture, then you
could estimate the Green-Ampt parameters using [1]
What do you want to do in that area? Perform hydrological/overland
flow analysis?
If it is the case, you might give Itzï a try:
https://www.itzi.org/
Some of the
2016-09-16 0:14 GMT-05:00 Eddison Araya :
> El 15/9/2016 9:00 a. m., "José Anderson"
> escribió:
>>
>> Dear all,
>>
>> I'm in need to start a study in a urban area in which the DTM
>> (https://en.wikipedia.org/wiki/Digital_elevation_model) is
Hello José,
There is basically three global DEM freely available:
- Aster DEM, which is agreed to be garbage
- SRTM, which is the reference
- the newly released ALOS DEM, which seems to be the best quality of all 3.
They are surface models, so trees and buildings are affecting the altitude
2016-09-08 12:32 GMT-05:00 Veronica Andreo :
> Hello,
>
> You can use r.series for that and pass the list of maps that you get with
> g.list directly into r.series using back ticks, or you can also pass a file
> with all map names. Check the manual for more examples:
>
2016-08-16 13:30 GMT-05:00 Vaclav Petras <wenzesl...@gmail.com>:
>
> On Sat, Aug 13, 2016 at 1:09 PM, Laurent C. <lrn...@gmail.com> wrote:
>>
>> I've tried but never manage to make it work. See:
>> https://lists.osgeo.org/pipermail/grass-dev/2016-March/079481.
roblem. For some reason, pip is not adding the
local bin directory to the PATH when installing with --user, at least
in Ubuntu.
So you might need to add $HOME/.local/bin to your PATH.
Laurent
>
>
>
> On Fri, Aug 12, 2016 at 7:09 PM, Laurent C. <lrn...@gmail.com> wrote:
>>
>&g
2016-08-12 22:21 GMT-05:00 Anna Petrášová <kratocha...@gmail.com>:
> On Fri, Aug 12, 2016 at 7:09 PM, Laurent C. <lrn...@gmail.com> wrote:
>> Hello,
>>
>> This message to let you know that I'll give a talk called "A GRASS GIS
>> module for 2D superficial
Hello,
Then pass quiet=true to run_command()
Laurent
El ago. 13, 2016 2:39, "Micha Silver" <mi...@arava.co.il> escribió:
> -- Original Message -- Subject: Re: [GRASS-user] GRASS_VERBOSE=-1
> Date: Fri, 12 Aug 2016 18:21:55 -0500 To: Micha Silver, Grass User
Hello,
You try to pass a string as an argument to run_command, which will not work.
If you want to change the input, change that string only:
for m in maps:
grass.run_command('r.in.gdal', input=m, output=m)
or, as in your example:
grass.run_command('r.in.gdal', input=sys.argv[5],
Hello,
One solution is to set the environment variable like this:
os.environ['GRASS_VERBOSE'] = '3'
the other is to pass the verbosity as an argument to run_command(verbose=True)
Cheers,
Laurent
2016-08-12 9:56 GMT-05:00 Micha Silver :
> How do I use the VERBOSE=-1 option
Hello,
This message to let you know that I'll give a talk called "A GRASS GIS
module for 2D superficial flow simulations" at the 12th International
Conference on Hydroinformatics that will be held from 21 to 26 of
august in Incheon, Korea. The goal is basically to present Itzï
Eddison,
First align your region to the destination raster with 'g.region
raster=rastername'.
Then resample the other raster.
Laurent
2016-08-02 11:48 GMT-05:00 Eddison Araya :
> Hello , nice to say hello , I have 2 raster : A raster is 20 meters and
> another 30 meters
Hello all,
I'd like to know if there is a way to detect buildings using
superficial lidar data.
I have two types of data:
- 5m DSM raster
- 5m DTM raster
Many described methods use multi-spectral images, but I only have
access to Landsat, which I guess will be too coarse to properly detect
2016-05-24 4:36 GMT-05:00 Gareth Grewcock :
> Again - thanks for all the help so far, it is much appreciated.
>
> Laurent - I ran the read2_command and it returned the following:
>
> Traceback (most recent call last):
> File "", line 254, in run_nodebug
> File
ppets
> but can't get r.shade to work with Python either externally or through the
> UI, I really have no idea what the problem is as all of the other modules
> work fine.
>
> Cheers,
>
> Gareth
>
> On Mon, May 23, 2016 at 3:38 PM, Laurent C. <lrn...@gmail.com> wro
Hello,
Try the .splitlines() method on the string.
Cheers,
Laurent
El abr. 5, 2016 19:52, "Leonardo Hardtke" escribió:
> Dear list members,
> I know this is more a python question than a GRASS question, bot for sure
> some of you knows the ansewer and is able to help
Dimitris,
How did you install GRASS ? From the Suse repositories?
What happen if you try to launch the GUI manually with g.gui (possibly
with --v flag)?
Regards,
Laurent
2016-03-18 7:42 GMT-06:00 ddarlac :
> I am a Qgis user, and wanted to do some 3D videos - I read Nviz
Hello,
I'm working with files that represent the variation of a value over
time for a single point.
Looks like that:
0 1200
0.03125 1380
0.125 1500
0.3125 1620
0.65625 1740
I need to import them as a raster space-time dataset.
What I'm doing now is using a python script that does that
Hi all,
I'm coming back to that subject because the GPM datas are finally available
since December 2014. The original data are stored in HDF5, which is
supported by GDAL.
However, I had a problem with the geo-referencing of the data. It seems
that the lat/long coordinates are flipped.
The problem
Hello all,
I'm trying to write a new area to a new vector map, without success.
If I do the following :
new_boundary = Boundary(points=coordinates)
new_centroid = Point(point_x,point_y)
new_area = Area(boundary = new_boundary, centroid=new_centroid)
new_vect_map.write(new_area)
I get this
again !
Regards,
Laurent
2014-12-03 0:28 GMT-06:00 Pietro peter.z...@gmail.com:
Hi Laurent,
On Tue, Dec 2, 2014 at 4:41 PM, Laurent C. lrn...@gmail.com wrote:
With pygrass I managed to set DB links and import the attributes in
different tables.
But the only way I've found to change the layer
Hello all,
I'm writing a GRASS module with pygrass to import data from a text file
into a new vector map.
Those data have topological relationships: they are nodes and links of a
sewer network. But those objects may have different types of attributes.
For example a link could be a pipe or a pump.
81 matches
Mail list logo