[GRASS-user] Help; Running Grass from R without installing GRASS

GRASS GRASS user
,
1

Hi Everyone,
I have detail explanations of connecting to an existing dataset in GRASS from R in this address: https://grasswiki.osgeo.org/wiki/R_statistics/rgrass7

What I hope I could do is to create a GRASS dataset in R on the fly and without the need to install GRASS or create a dataset on it. Is there any way to do this?

i.e. I have a vector /raster file and use that to identify the region and details required to set up a dataset and use GRASS functions just by including rgrass library.

Thanks for your help and attention,

Mehrdad

···

Mehrdad Varedi

2
···

On 27/04/2019 22:05, Mehrdad Varedi wrote:

Hi Everyone,
I have detail explanations of connecting to an existing dataset in GRASS from R in this address: https://grasswiki.osgeo.org/wiki/R_statistics/rgrass7

What I hope I could do is to create a GRASS dataset in R on the fly and without the need to install GRASS or create a dataset on it. Is there any way to do this?

Maybe this will help:

(You must have GRASS installed, of course)

# Parameters for the GRASS installation and temporary GRASS mapset

GISBase = "/usr/lib/grass76"
# Set these as you like

GISDbase = "/tmp"
Location = "tmp_location"
Mapset = "tmp_mapset"
georef = "EPSG:2039"

mapset_path = file.path(GISDbase, Location, Mapset)

# Now run the 'grass' command within R and


`# use the -c and -e flags to create the temp mapset, then exit`

`setup_grass_cmd = paste("grass", "-c", georef, "-e", mapset_path)`

`system(setup_grass_cmd)`

# Load the R grass interface and initialize GRASS within R,

# using the temporary mapset from above

`library(rgrass7)`
`initGRASS(home=tempdir(),`
`gisBase = GISBase,`
`gisDbase = GISDbase,`
`location = Location,`
`mapset = Mapset,`
`remove_GISRC = TRUE)`
``

`# Try some GRASS commands`

`input_tif = "isrlat12.tif"`

`grass_rast = "isrlat12"`
`execGRASS("r.in.gdal", flags = c("o","overwrite"),`
`input = input_tif,`
`output = grass_rast)`
`execGRASS("g.region", flags="p", raster = grass_rast)`

I did not add a command to remove the temporary mapset. Something like:

`unlink(file.path(GISDbase, Location), recusive = TRUE)` might be necessary

One other note: you probably will get better suggestions on the R-sig-geo list...

> i.e. I have a vector /raster file and use that to identify the region and details required to set up a dataset and use GRASS functions just by including rgrass library.
> 
> Thanks for your help and attention,
> 
> Mehrdad
> 
> --
> 
> Mehrdad Varedi
> 
> <img alt="" style="width:0px;max-height:0px;overflow:hidden" src="https://mailfoogae.appspot.com/t?sender=admFyZWRpQGdtYWlsLmNvbQ%3D%3D&amp;type=zerocontent&amp;guid=dbe7ab29-5a21-45c9-97bb-082f4f48224a" moz-do-not-send="true">ᐧ
> 
> ```
> _______________________________________________
> grass-user mailing list
> [grass-user@lists.osgeo.org](mailto:grass-user@lists.osgeo.org)
> [https://lists.osgeo.org/mailman/listinfo/grass-user](https://lists.osgeo.org/mailman/listinfo/grass-user)
> ```

```
-- 
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918
```

</details>
3

Thank you for your time and sharing your knowledge Micha for the second time.

I followed your code and found that initGRASS searches for libraries and binaries of GRASS in gisBase folder. If we address a folder that doesn’t have the code it comes back with this error message: “The code execution cannot proceed because libgrass_gis.7.4.1.dll was not found. Reinstalling the program may fix this problem.” it gives the same error even if I copy that file into GISBase folder.

I want to run my code on a shiny server which I can’t install GRASS although I can copy files and libraries. Can I solve this issue by copying a certain folder of GRASS binaries with my app or I am going the wrong path?

Please advise,

Thank you once again for all the time you share to answer my questions,

Mehrdad

ᐧ

···

Mehrdad Varedi

4
···

On 28/04/2019 19:30, Mehrdad Varedi wrote:

Thank you for your time and sharing your knowledge Micha for the second time.

I followed your code and found that initGRASS searches for libraries and binaries of GRASS in gisBase folder. If we address a folder that doesn’t have the code it comes back with this error message: “The code execution cannot proceed because libgrass_gis.7.4.1.dll was not found. Reinstalling the program may fix this problem.” it gives the same error even if I copy that file into GISBase folder.

I want to run my code on a shiny server which I can’t install GRASS although I can copy files and libraries. Can I solve this issue by copying a certain folder of GRASS binaries with my app or I am going the wrong path?

I have no idea.

So you’re trying to run an application that’s not installed? I can’t imagine how that would work :-(. After all the libraries (*.dll’s if this is windows) are searched for in a certain OS defined path. And many libraries depend on others, so the whole chain has to be found.

What are the limitations on this server? Can you compile a binary into your home directory? If so, you might be able to install grass with the GISBASE in your home directory, and run from there. But again, GRASS has lot’s of dependencies: proj4, gdal, and many others. Any you need access to the compile toolchain…

Maybe others will have some more insight.

Regards, Micha

Please advise,

Thank you once again for all the time you share to answer my questions,

Mehrdad

ᐧ

On Sun, Apr 28, 2019 at 10:55 AM Micha Silver <tsvibar@gmail.com> wrote:

On 27/04/2019 22:05, Mehrdad Varedi wrote:

Hi Everyone,
I have detail explanations of connecting to an existing dataset in GRASS from R in this address: https://grasswiki.osgeo.org/wiki/R_statistics/rgrass7

What I hope I could do is to create a GRASS dataset in R on the fly and without the need to install GRASS or create a dataset on it. Is there any way to do this?

Maybe this will help:

(You must have GRASS installed, of course)

# Parameters for the GRASS installation and temporary GRASS mapset

GISBase = "/usr/lib/grass76"
# Set these as you like

GISDbase = "/tmp"
Location = "tmp_location"
Mapset = "tmp_mapset"
georef = "EPSG:2039"

mapset_path = file.path(GISDbase, Location, Mapset)

# Now run the 'grass' command within R and

# use the -c and -e flags to create the temp mapset, then exit

setup_grass_cmd = paste("grass", "-c", georef, "-e", mapset_path)

system(setup_grass_cmd)

Load the R grass interface and initialize GRASS within R,

using the temporary mapset from above

library(rgrass7)
initGRASS(home=tempdir(),
gisBase = GISBase,
gisDbase = GISDbase,
location = Location,
mapset = Mapset,
remove_GISRC = TRUE)

# Try some GRASS commands

input_tif = "isrlat12.tif"

grass_rast = "isrlat12"
execGRASS("r.in.gdal", flags = c("o","overwrite"),
input = input_tif,
output = grass_rast)
execGRASS("g.region", flags="p", raster = grass_rast)

I did not add a command to remove the temporary mapset. Something like:

unlink(file.path(GISDbase, Location), recusive = TRUE) might be necessary

One other note: you probably will get better suggestions on the R-sig-geo list…

i.e. I have a vector /raster file and use that to identify the region and details required to set up a dataset and use GRASS functions just by including rgrass library.

Thanks for your help and attention,

Mehrdad

–

Mehrdad Varedi

ᐧ

_______________________________________________
grass-user mailing list
[grass-user@lists.osgeo.org](mailto:grass-user@lists.osgeo.org)
[https://lists.osgeo.org/mailman/listinfo/grass-user](https://lists.osgeo.org/mailman/listinfo/grass-user)

-- 
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918

–

Mehrdad Varedi

-- 
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918
5

Hi Mehrdad,

If this is a permission issue, and the dependencies (and a build environment) are installed, you could try compile GRASS (and install it in your home directory if necessary with the --prefix configure option).

svn checkout https://svn.osgeo.org/grass/grass/ grass_trunk

cd grass_trunk

./configure

make

make install

The make install step is not required for starting GRASS though. The make commando creates a bin…/ directory within the grass_trunk dir, where the grass startupscript is located…

Hope that helps.

Cheers,

Stefan

···

Thank you for your time and sharing your knowledge Micha for the second time.

I followed your code and found that initGRASS searches for libraries and binaries of GRASS in gisBase folder. If we address a folder that doesn’t have the code it comes back with this error message: “The code execution cannot proceed because libgrass_gis.7.4.1.dll was not found. Reinstalling the program may fix this problem.” it gives the same error even if I copy that file into GISBase folder.

I want to run my code on a shiny server which I can’t install GRASS although I can copy files and libraries. Can I solve this issue by copying a certain folder of GRASS binaries with my app or I am going the wrong path?

Please advise,

Thank you once again for all the time you share to answer my questions,

Mehrdad

ᐧ

On Sun, Apr 28, 2019 at 10:55 AM Micha Silver <tsvibar@gmail.com> wrote:

On 27/04/2019 22:05, Mehrdad Varedi wrote:

Hi Everyone,

I have detail explanations of connecting to an existing dataset in GRASS from R in this address: https://grasswiki.osgeo.org/wiki/R_statistics/rgrass7

What I hope I could do is to create a GRASS dataset in R on the fly and without the need to install GRASS or create a dataset on it. Is there any way to do this?

Maybe this will help:

(You must have GRASS installed, of course)

# Parameters for the GRASS installation and temporary GRASS mapset

GISBase = "/usr/lib/grass76"
# Set these as you like

GISDbase = "/tmp"
Location = "tmp_location"
Mapset = "tmp_mapset"
georef = "EPSG:2039"

mapset_path = file.path(GISDbase, Location, Mapset)

# Now run the 'grass' command within R and

# use the -c and -e flags to create the temp mapset, then exit

setup_grass_cmd = paste("grass", "-c", georef, "-e", mapset_path)

system(setup_grass_cmd)

Load the R grass interface and initialize GRASS within R,

using the temporary mapset from above

library(rgrass7)
initGRASS(home=tempdir(),
gisBase = GISBase,
gisDbase = GISDbase,
location = Location,
mapset = Mapset,
remove_GISRC = TRUE)

# Try some GRASS commands

input_tif = "isrlat12.tif"

grass_rast = "isrlat12"
execGRASS("r.in.gdal", flags = c("o","overwrite"),
input = input_tif,
output = grass_rast)
execGRASS("g.region", flags="p", raster = grass_rast)

I did not add a command to remove the temporary mapset. Something like:

unlink(file.path(GISDbase, Location), recusive = TRUE) might be necessary

One other note: you probably will get better suggestions on the R-sig-geo list…

i.e. I have a vector /raster file and use that to identify the region and details required to set up a dataset and use GRASS functions just by including rgrass library.

Thanks for your help and attention,

Mehrdad

–

Mehrdad Varedi

ᐧ

_______________________________________________

grass-user mailing list

[grass-user@lists.osgeo.org](mailto:grass-user@lists.osgeo.org)

[https://lists.osgeo.org/mailman/listinfo/grass-user](https://lists.osgeo.org/mailman/listinfo/grass-user)

-- 

Micha Silver

Ben Gurion Univ.

Sde Boker, Remote Sensing Lab

cell: +972-523-665918

–

Mehrdad Varedi

6

Thank you Micha,

I was under the impression that rgrass7 library has all the functions in GRASS. Apparently, all it does is calling modules from GRASS and GRASS should be completely installed on the same machine that the R code is running. I wish at least we could connect to a GRASS copy on a web-based VM.
I enjoyed the part of your letter that explained how to created a dataset outside of GRASS using R code and thanks for that!

I wish it was possible to read the functions directly from a web host/cloud instead of installing it everywhere (Even if it was a paid service it could be an attractive option). ESRI has some API, I wish we could have it for GRASS or I could help to write the API for GRASS.

Best regards,

Mehrdad

ᐧ

···

Mehrdad Varedi

7

Hi Stefan,

Thanks for your reply. I will try it in case of a permission issue.

Kind regards,

Mehrdad

ᐧ

···

Mehrdad Varedi

8

Hi Everyone,
I used the code from Micha and managed to run grass in R and create a new dataset from R. I can’t see the SQLite database in the created folders like the ones I usually get when creating new datasets from GRASS.

How can I have the new SQlite database?

Thanks for your help.

···

Mehrdad Varedi

9
···

On 28/08/2019 9:56, Mehrdad Varedi wrote:

Hi Everyone,
I used the code from Micha and managed to run grass in R and create a new dataset from R. I can’t see the SQLite database in the created folders like the ones I usually get when creating new datasets from GRASS.

How can I have the new SQlite database?

If you import some vector you’ll have the sqlite.db setup for you:

(Following the example from before)

> input_vect = "Israel_border.gpkg"
> grass_vect = "isr_border"
> execGRASS("v.import", input=input_vect, output=grass_vect, flags="o")
Over-riding projection check
Check if OGR layer <Israel_border> contains polygons...
100%
Creating attribute table for layer <Israel_border>...
Default driver / database set to:
driver: sqlite
database: $GISDBASE/$LOCATION_NAME/$MAPSET/sqlite/sqlite.db
Importing 1 features (OGR layer <Israel_border>)...
100%

Then in a separate terminal:

micha@TP480:~$ ll /tmp/tmp_location/tmp_mapset/sqlite/
total 20
drwxr-xr-x 2 micha micha 4096 Aug 28 13:19 ./
drwxr-xr-x 6 micha micha 4096 Aug 28 13:19 ../
-rw-r--r-- 1 micha micha 12288 Aug 28 13:19 sqlite.db

Thanks for your help.

On Sun, Apr 28, 2019 at 10:55 AM Micha Silver <tsvibar@gmail.com> wrote:

Maybe this will help:

(You must have GRASS installed, of course)

# Parameters for the GRASS installation and temporary GRASS mapset

GISBase = "/usr/lib/grass76"
# Set these as you like

GISDbase = "/tmp"
Location = "tmp_location"
Mapset = "tmp_mapset"
georef = "EPSG:2039"

mapset_path = file.path(GISDbase, Location, Mapset)

# Now run the 'grass' command within R and

# use the -c and -e flags to create the temp mapset, then exit

setup_grass_cmd = paste("grass", "-c", georef, "-e", mapset_path)

system(setup_grass_cmd)

Load the R grass interface and initialize GRASS within R,

using the temporary mapset from above

library(rgrass7)
initGRASS(home=tempdir(),
gisBase = GISBase,
gisDbase = GISDbase,
location = Location,
mapset = Mapset,
remove_GISRC = TRUE)

# Try some GRASS commands

input_tif = "isrlat12.tif"

grass_rast = "isrlat12"
execGRASS("r.in.gdal", flags = c("o","overwrite"),
input = input_tif,
output = grass_rast)
execGRASS("g.region", flags="p", raster = grass_rast)

I did not add a command to remove the temporary mapset. Something like:

unlink(file.path(GISDbase, Location), recusive = TRUE) might be necessary

Micha Silver

Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918

–

Mehrdad Varedi

ᐧ

-- 
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918
10

Hi Micha,
Your answer makes perfect sense. Thank you. I was expecting to get the database when only I had imported raster files!
Anyway, when I tried to import a vector file with the following command, I got an error:

execGRASS(“v.import”, parameters=list(
input = basename(inputWM),
output = “WM”
),
flags=c(“overwrite”))

Here is the error:

In system(cmd0, intern = TRUE) :
running command ‘v.import.bat --interface-description’ had status 1
Error in parseGRASS(cmd, legacyExec = legacyExec) : v.import not parsed

I traced the code and found that the parseGRASS adds the .bat extension for the v.import and the lagacyExec has a TRUE value (default for Windows OS)

FYI, I couldn’t run GRASS from Rstudio unless ran the RStudio from the “OSGeo4W” environment. It helps RStudio to find the functions in GRASS. (I learned it from a previous post). So it works and importing the raster files had no problem (I imported a virtual raster .vrt with “r.in.gdal” command). Also, I processed the raster file with no problem. The problem only shows up when I try to import the vector files with r.import command.

Any idea what is happening here?

Thanks for your help,

Mehrdad

ᐧ

···

Mehrdad Varedi

11
···

On 28/08/2019 21:12, Mehrdad Varedi wrote:

Hi Micha,
Your answer makes perfect sense. Thank you. I was expecting to get the database when only I had imported raster files!
Anyway, when I tried to import a vector file with the following command, I got an error:

execGRASS(“v.import”, parameters=list(
input = basename(inputWM),
output = “WM”
),
flags=c(“overwrite”))

Here is the error:

In system(cmd0, intern = TRUE) :
running command ‘v.import.bat --interface-description’ had status 1
Error in parseGRASS(cmd, legacyExec = legacyExec) : v.import not parsed

I traced the code and found that the parseGRASS adds the .bat extension for the v.import and the lagacyExec has a TRUE value (default for Windows OS)

FYI, I couldn’t run GRASS from Rstudio unless ran the RStudio from the “OSGeo4W” environment. It helps RStudio to find the functions in GRASS. (I learned it from a previous post). So it works and importing the raster files had no problem (I imported a virtual raster .vrt with “r.in.gdal” command). Also, I processed the raster file with no problem. The problem only shows up when I try to import the vector files with r.import command.

Any idea what is happening here?

No, sorry.

I expect others on the list with more experience with Windows installations will be able to help.

Thanks for your help,

Mehrdad

ᐧ

On Wed, Aug 28, 2019 at 6:24 AM Micha Silver <tsvibar@gmail.com> wrote:

On 28/08/2019 9:56, Mehrdad Varedi wrote:

Hi Everyone,
I used the code from Micha and managed to run grass in R and create a new dataset from R. I can’t see the SQLite database in the created folders like the ones I usually get when creating new datasets from GRASS.

How can I have the new SQlite database?

If you import some vector you’ll have the sqlite.db setup for you:

(Following the example from before)

> input_vect = "Israel_border.gpkg"
> grass_vect = "isr_border"
> execGRASS("v.import", input=input_vect, output=grass_vect, flags="o")
Over-riding projection check
Check if OGR layer <Israel_border> contains polygons...
100%
Creating attribute table for layer <Israel_border>...
Default driver / database set to:
driver: sqlite
database: $GISDBASE/$LOCATION_NAME/$MAPSET/sqlite/sqlite.db
Importing 1 features (OGR layer <Israel_border>)...
100%

Then in a separate terminal:

micha@TP480:~$ ll /tmp/tmp_location/tmp_mapset/sqlite/
total 20
drwxr-xr-x 2 micha micha 4096 Aug 28 13:19 ./
drwxr-xr-x 6 micha micha 4096 Aug 28 13:19 ../
-rw-r--r-- 1 micha micha 12288 Aug 28 13:19 sqlite.db

Thanks for your help.

On Sun, Apr 28, 2019 at 10:55 AM Micha Silver <tsvibar@gmail.com> wrote:

Maybe this will help:

(You must have GRASS installed, of course)

# Parameters for the GRASS installation and temporary GRASS mapset

GISBase = "/usr/lib/grass76"
# Set these as you like

GISDbase = "/tmp"
Location = "tmp_location"
Mapset = "tmp_mapset"
georef = "EPSG:2039"

mapset_path = file.path(GISDbase, Location, Mapset)

# Now run the 'grass' command within R and

# use the -c and -e flags to create the temp mapset, then exit

setup_grass_cmd = paste("grass", "-c", georef, "-e", mapset_path)

system(setup_grass_cmd)

Load the R grass interface and initialize GRASS within R,

using the temporary mapset from above

library(rgrass7)
initGRASS(home=tempdir(),
gisBase = GISBase,
gisDbase = GISDbase,
location = Location,
mapset = Mapset,
remove_GISRC = TRUE)

# Try some GRASS commands

input_tif = "isrlat12.tif"

grass_rast = "isrlat12"
execGRASS("r.in.gdal", flags = c("o","overwrite"),
input = input_tif,
output = grass_rast)
execGRASS("g.region", flags="p", raster = grass_rast)

I did not add a command to remove the temporary mapset. Something like:

unlink(file.path(GISDbase, Location), recusive = TRUE) might be necessary

Micha Silver

Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918

–

Mehrdad Varedi

ᐧ

-- 
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918

–

Mehrdad Varedi

-- 
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918
12

Hi,

pls write to grass-stats list where you can get more answers. What version of rgrass7 are you using? The last one includes the possibility to use either sp or sf for vector data. When I run readVECT for example, I first run either use_sp() or use_sf(). Maybe that has something to do with the error? Just guessing…

HTH,
Vero

El mié., 28 ago. 2019 a las 20:33, Micha Silver (<tsvibar@gmail.com>) escribió:

On 28/08/2019 21:12, Mehrdad Varedi wrote:

Hi Micha,
Your answer makes perfect sense. Thank you. I was expecting to get the database when only I had imported raster files!
Anyway, when I tried to import a vector file with the following command, I got an error:

execGRASS(“v.import”, parameters=list(
input = basename(inputWM),
output = “WM”
),
flags=c(“overwrite”))

Here is the error:

In system(cmd0, intern = TRUE) :
running command ‘v.import.bat --interface-description’ had status 1
Error in parseGRASS(cmd, legacyExec = legacyExec) : v.import not parsed

I traced the code and found that the parseGRASS adds the .bat extension for the v.import and the lagacyExec has a TRUE value (default for Windows OS)

FYI, I couldn’t run GRASS from Rstudio unless ran the RStudio from the “OSGeo4W” environment. It helps RStudio to find the functions in GRASS. (I learned it from a previous post). So it works and importing the raster files had no problem (I imported a virtual raster .vrt with “r.in.gdal” command). Also, I processed the raster file with no problem. The problem only shows up when I try to import the vector files with r.import command.

Any idea what is happening here?

No, sorry.

I expect others on the list with more experience with Windows installations will be able to help.

Thanks for your help,

Mehrdad

ᐧ

On Wed, Aug 28, 2019 at 6:24 AM Micha Silver <tsvibar@gmail.com> wrote:

On 28/08/2019 9:56, Mehrdad Varedi wrote:

Hi Everyone,
I used the code from Micha and managed to run grass in R and create a new dataset from R. I can’t see the SQLite database in the created folders like the ones I usually get when creating new datasets from GRASS.

How can I have the new SQlite database?

If you import some vector you’ll have the sqlite.db setup for you:

(Following the example from before)

> input_vect = "Israel_border.gpkg"
> grass_vect = "isr_border"
> execGRASS("v.import", input=input_vect, output=grass_vect, flags="o")
Over-riding projection check
Check if OGR layer <Israel_border> contains polygons...
100%
Creating attribute table for layer <Israel_border>...
Default driver / database set to:
driver: sqlite
database: $GISDBASE/$LOCATION_NAME/$MAPSET/sqlite/sqlite.db
Importing 1 features (OGR layer <Israel_border>)...
100%

Then in a separate terminal:

micha@TP480:~$ ll /tmp/tmp_location/tmp_mapset/sqlite/
total 20
drwxr-xr-x 2 micha micha 4096 Aug 28 13:19 ./
drwxr-xr-x 6 micha micha 4096 Aug 28 13:19 ../
-rw-r--r-- 1 micha micha 12288 Aug 28 13:19 sqlite.db

Thanks for your help.

On Sun, Apr 28, 2019 at 10:55 AM Micha Silver <tsvibar@gmail.com> wrote:

Maybe this will help:

(You must have GRASS installed, of course)

# Parameters for the GRASS installation and temporary GRASS mapset

GISBase = "/usr/lib/grass76"
# Set these as you like

GISDbase = "/tmp"
Location = "tmp_location"
Mapset = "tmp_mapset"
georef = "EPSG:2039"

mapset_path = file.path(GISDbase, Location, Mapset)

# Now run the 'grass' command within R and

# use the -c and -e flags to create the temp mapset, then exit

setup_grass_cmd = paste("grass", "-c", georef, "-e", mapset_path)

system(setup_grass_cmd)

Load the R grass interface and initialize GRASS within R,

using the temporary mapset from above

library(rgrass7)
initGRASS(home=tempdir(),
gisBase = GISBase,
gisDbase = GISDbase,
location = Location,
mapset = Mapset,
remove_GISRC = TRUE)

# Try some GRASS commands

input_tif = "isrlat12.tif"

grass_rast = "isrlat12"
execGRASS("r.in.gdal", flags = c("o","overwrite"),
input = input_tif,
output = grass_rast)
execGRASS("g.region", flags="p", raster = grass_rast)

I did not add a command to remove the temporary mapset. Something like:

unlink(file.path(GISDbase, Location), recusive = TRUE) might be necessary

Micha Silver

Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918

–

Mehrdad Varedi

ᐧ

-- 
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918

–

Mehrdad Varedi

-- 
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918

grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

13

The problem only shows up when I try to import the vector files >with

r.import command.

Vector import works with v.import.

-----
best regards
Helmut
--
Sent from: http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html

14

Hi Veronica,
I think the error should be because of the point you mentioned somehow.
When I use v.import, I can’t import a vector file, although if I use v.in.ogr, it works very well and as Micha had mentioned the sqlite.db is created.

BTW, the version of rgrass7 on my machine is 0.1.12. Also, I tried it on ubuntu and everything worked with no problem following the code provided by Micha.

Thank you all,

Mehrdad

ᐧ

···

Mehrdad Varedi