Questions tagged [grass.script]

The Python scripting library for GRASS

Filter by
Sorted by
Tagged with
1
vote
1answer
55 views

Error while performing actions in GRASS GIS in Python

While performing actions in python console as described in the module of GRASS GIS,am getting following error: import grass.script as gs r.in.gdal -e input="D:...\ASTGTMV003_N19E073_dem" ...
3
votes
0answers
39 views

CalledModuleError: Module run None g.gisenv -n ended with error

I am trying to use r.horizon GRASS module(version GRASS 7.8) in Python script. But I couldn't able to access the grass modules. Everytime, I am getting the following error. File "C:\OSGEO4~1\apps\...
2
votes
0answers
33 views

GRASS GIS: r.out.gdal raises an exception and hangs if it's called multiple consecutive times

I'm using a small python routine (grass.script, pygrass) to handle several operations on sentinel2 imagery. I'm calling r.out.gdal inside the script and a very annoying problem occurs: when it is ...
0
votes
1answer
23 views

GRASS GIS 7.8 trigometric function of vector attribtue

I need to calculate the sin and cos of a vector attribute column, but these functions do not appear to be available. In Python I've tried using v.db.update: gscript.run_command('v.db.update', map='...
0
votes
0answers
48 views

GRASS GIS 7.8 Delete mapset from command line

Within GRASS GIS 7.8, a new mapset can be created using g.mapset -c mapset=mapset location=location, but is it possible to delete a mapset from the command line? I know that a mapset can be removed ...
0
votes
1answer
12 views

How to display graphics in IOS when connected remotely with Linux server?

I am trying to work using GRASS GIS for my data remotely from Mac. Here is what I have done: XQuartz was installed and run; Terminal was run. I typed: ssh -Y root@IPaddress Then I typed: grass64 (to ...
0
votes
0answers
2k views

eval and cats is not recognized as an internal or external command, operable program or batch file in GRASS GIS

I'm trying to do the r.lfp and r.accumulate examples for multiple longest flow paths, except I keep getting the "eval" or "cats" "is not recognized as an internal or external ...
0
votes
1answer
29 views

r.green error installation on GRASS

I am using GRASS GIS 7.0.3 on linux, however I've been trying to install the r.green extension on GRASS and always get the error message: "ImportError: cannot import name set_path" "...
0
votes
0answers
26 views

Create a loop to run and store the r.stats output of several viewsheds?

I have 15 viewsheds (view 0-14) from which I would like to export the cell count. I can do this for one viewshed through: count=[] i=0 viewshedSize=grass.read_command("r.stats", overwrite=True,flags=...
1
vote
1answer
26 views

run r.random through script in python shell

Trying to run r.random using the script but receiving an error message. It works when the module is run outside of the python shell so it must be te code, I'm new to this so I can't figure out what is ...
1
vote
0answers
39 views

v.overlay with empty layer [closed]

I'm writing a Grass Python Script (using GRASS GIS 7.6.1) for a project for university. Sadly, I never learned Python, so it's a bit tricky but so far I got everything done I wanted to. The Script ...
0
votes
1answer
30 views

Creating reclass-rules-textfile within Simple Python Editor of GRASS GIS

We are writing a script with the simple Python editor (GRASS GIS). We try to set the rules for a reclass command within the script (without creating a txt.file manually). Is it possible to set the ...
1
vote
1answer
81 views

GRASS GIS r.in.gdal, import multiple band ncdf

I am trying to import a ncdf file with 20000 bands (time steps of precipitation) to grass gis using r.in.gdal. My problem is that i dont need all the bands, only a portion of it, e.g : band 1 to band ...
0
votes
1answer
60 views

GRASS GIS: creating raster with a certain gradient and extent

Constraints - Use only GRASS GIS - Automated using pyGRASS Input DEM River raster OR river as vector (both are computed) gradient of the valley (as a float) vector line with direction and magnitude ...
1
vote
1answer
68 views

Creating cumulative friction map

I am looking at bottom friction (Manning's n) values in relation to storm surge as large volumes of water flow over a surface. In order to calculate total head loss occurring as water passes from the ...
0
votes
0answers
107 views

GRASS GIS: From multiple lines to single line

This seems like a simple process but somehow I am not able to solve this. The situation is as follows: I have three separate lines which form visually 1 line. As I want to further process these lines ...
1
vote
0answers
126 views

GRASS r.sun returns “ERROR: insol_time and incidout are incompatible options”

I'm trying to measure irradiance on a DEM using the r.sun tool in QGIS 3.6. I have grids for the elevation (EPSG:32735 - WGS 84 / UTM zone 35S), the slope and aspect (EPSG:4326 - WGS 84). I've tried ...
0
votes
0answers
46 views

GRASS v.hull for many point sets

I have map "my_pointmap" with 3D points with ids=(1,2,3...N). I have also list of subsets of these point ids , S = [(1,2,3,4),(2,4,10,12), ...]. I would love to get GRASS GIS vector map with |S| 3D ...
1
vote
1answer
142 views

Adding GRASS module to QGIS 2.18.21

I did a module in grass wich I called r.out.haralickFeatures.py. I did a .qgm. I can't get the parameters rigth. #!/usr/bin/env python #%module #% description: Compute asm, contrast, var, idm and ...
0
votes
1answer
308 views

Loading ADF file from National Elevation Dataset into GRASS with Python?

I'm trying to load data from the National Elevation Dataset into GRASS/Python, using this Python script (outside of GRASS): # Imports import os import grass_session from grass.script import setup as ...
0
votes
1answer
68 views

Entry point not found Error initiating GRASS GIS from Python

I cannot run any GRASS tools from within Python 2.7.15. While the GRASS environment, grass.script etc. are loading correctly and I can run basic checks like gscript.read_command('g.gisenv') gscript....
2
votes
1answer
634 views

Watershed Delineation Loop - r.water.outlet resulting in NAN area while snapped

I am attempting to create multiple watershed delineations for stations within a CSV file. I have files for flow direction (raster with a cell size of 30m), flow accumulation(raster with a cell size of ...
1
vote
0answers
53 views

GRASS GIS Error Flags

I am attempting to run a GRASS command from the terminal and eventually a bash script, in which it links rasters in a folder to the Mapset previously defined. I was able to successfully create the ...
2
votes
1answer
197 views

grass GIS + spyder python script error

I'm getting crazy with grass gis and python. I want to call a grass function from an external python script (using spyder for instance). I've installed QGIS + GRASS with the OSGeo4W64 installation ...
1
vote
1answer
203 views

Setting variable from input map in GRASS python script?

I'm currently working on a GRASS script that takes maps from user input and computes them into a mapcalc expression. The issue i'm stucked with that when i ask for the map i can't set a variable name ...
0
votes
1answer
135 views

How to execute a GRASS module developed in Python?

I am following the FOSS4G-E 2015 tutorial on how to develop GRASS modules with Python. The minimum template given in this tutorial is this script: %%file r.viewshed.points.py #!/usr/bin/env python #%...
6
votes
2answers
4k views

How to clip raster in GRASS GIS?

I am using GRASS GIS 7.2.0. The intent is to clip my DEM raster. I have searched available raster tools but didn't found anything allow me to clip my raster using polygon vector or interactively. How ...
1
vote
0answers
32 views

(d.*) of GRASS GIS does not execute from standalone script

I am in windows 10. All my module import works well and grass env for standalone script is fine. I tried to execute following commands from standalone script. Neither wx1 display gets executed nor it ...
1
vote
2answers
424 views

How to calculate area of the output of r.water.outlet?

I wrote the following code to use r.water.outlet module in GRASS GIS: import os import sys #set up GRASS environment variables sys.path.append(os.path.join(os.environ['GISBASE'], 'etc', 'python')) ...
0
votes
1answer
210 views

Extracting raster information within/inside polygon (for each contour)

I have raster data of NDVI and DEM I have already created contour from DEM. Now I'm trying to find following: for each contour find the raster information within the contour. Can anyone suggest ...
0
votes
1answer
137 views

Possible to convert v.net.spanningtree to QGIS custom script so don't need GRASS installed?

Is it possible (and if yes then how or has this been done) to convert v.net.spanningtree in to a QGIS/python (custom) script that can be accessed in the QGIS graphical modeler in a way that GRASS is ...
2
votes
1answer
1k views

Connecting Python script external to Grass GIS 7 program in Windows 10?

I've tried to find resources online to link and run my Python script externally from Grass GIS 7, however, I have not been successful. I found resources here that provide instructions (e.g., https://...
1
vote
0answers
454 views

Use Function V.clean of the grass tool in a Python script

I need to integrate V.clean functionality of Grass in a script python. this functionnality is accessible in Qgis (Processing->tools). I found that I should use Grass API but I don't know the code ...
5
votes
3answers
3k views

GRASS 7 environment setup for Python

I would like to use some GRASS functionality from a Python programme. With GRASS 6 this required the set up of various environment variables, as detailed in the GRASS Wiki. I am trying a similar ...
1
vote
1answer
47 views

Problem with numeric variable in i.topo.corr

I have a script to do topographic correction on landsat images. Following the manual, I first have to create the illumination model. I do it with the following line: i.topo.corr -i --overwrite ...
2
votes
0answers
194 views

GRASS GIS: how to run g.run_command(“v.extract”) in Python with sqlite, smooth-spline in GRASS GIS?

I would like to solve several problems by using a python script in GRASS GIS: I already uploaded collected GPS-waypoints and tracks ("v.in.org"), patched several GPS-vector data ("v.patch") and ...
2
votes
1answer
114 views

Pygrass reading raster issue

I'm tring to read a raster map from my grassData. I'm using this code import os import grass.script.setup as gsetup from grass.pygrass import raster gisbase = os.environ['GISBASE'] gisdbase = os....
4
votes
1answer
108 views

GRASS GIS temporal raster individual layer statistical summary

I have created a temporal raster in GRASS GIS which has 90 layers. I would like to extract the mean of each raster layer and export it as text file. Right now I have to export the rasters as GeoTiffs ...
1
vote
1answer
248 views

Solar radiation flux at 5 minute interval using GRASS GIS r.sun

I have been trying to estimate instantaneous total solar irradiance(w/m2) at five minute interval for a day (eg. 21 sept). I can do this in SAGA GIS by using "moment, mode 1" and running the solar ...
2
votes
1answer
786 views

Capturing Output Messages from GRASS when Running Python Scripts

Using documentation here: https://grasswiki.osgeo.org/wiki/Working_with_GRASS_without_starting_it_explicitly#Python:_GRASS_GIS_7 I've figured out how to access GRASS modules in Python without ...
1
vote
1answer
160 views

Making GRASS GIS 7 Ubuntu bash file work on Windows 7

I have created this bash file (.sh) while working on watershed analysis stochastic simulation in GRASS GIS 7.0.2RC1, which is running on Ubuntu. It worked perfectly and I received the correct results. ...
3
votes
2answers
1k views

Missing v.in.ogr module of GRASS in QGIS Processing toolbox?

I am a new user to QGIS/GRASS and I am having a problem utilizing GRASS. I have several files that I generated in QGIS that I want to import into GRASS so I can use some of the Grass modules. ...
1
vote
1answer
65 views

Set region with rounded values/coordinates

How do I set a region update with rounded coordinates. I usually set my region from a puffer around my points of interest and round up the north and east coordinates and round down the south and west ...
2
votes
1answer
99 views

Return variable with GRASS GIS on a Windows machine

How do I return a variable in GRASS GIS on a Windows machine? For example I set a list of files to define my region. I already figured out that I can't use the bash command export, I need to use set ...
4
votes
1answer
2k views

Using Spyder IDE to develop python for Grass GIS, get 'ImportError: No module named grass.script'

I am trying to use Spyder IDE to develop a python script using python library from grass-6.4.3. I run the script below, and get ImportError: No module named grass.script I am using python 2.7, ...
2
votes
2answers
4k views

Extract Elevation of Points from DEM using Python QGIS or GRASS

I am utterly new to GIS, so please bear with me if I am providing incomplete information here. I have a DEM file (.dt2). What I am looking to do is as follows: take a series of points specified as ...
3
votes
0answers
173 views

Extract representative slope profiles from hillslope polygons

With the help of the Whitebox Geospatial Analysis Tools I was able to easily extract hill slopes from a river network in combination with a flow direction map. Now I would like to extract ...
1
vote
1answer
355 views

How to fix this “Sorry <1> is not a valid option” in GRASS-Ubuntu using r.water.outlet?

I am working with this .sh code (a while loop, extracted from: http://www.surfaces.co.il/) to create watersheds from drainage raster and outlet point of every watershed, but I obtain this error: "...
1
vote
0answers
75 views

Python script split up

I have a python script which works with Grass 6.4. At some point during the project I decided to split this big python script into several scripts. I have made all necessary variables global, created ...
1
vote
1answer
575 views

GRASS scripting for producing multi-band raster?

I am trying to loop the r.sun function within a Python script. Anyway, I seem not to understand how I can stack the raster maps each loop would produce into a single multi-band raster (preferably in ...