<div dir="ltr">Hi Micha,<div>Your answer makes perfect sense. Thank you. I was expecting to get the database when only I had imported raster files! </div><div>Anyway, when I tried to import a vector file with the following command, I got an error:</div><div><br></div><div><b>execGRASS("v.import", parameters=list(<br> input = basename(inputWM),<br> output = "WM"<br> ),<br> flags=c("overwrite"))<br></b></div><div><br></div><div>Here is the error: </div><div><br></div><div><b>In system(cmd0, intern = TRUE) :<br> running command 'v.import.bat --interface-description' had status 1<br>Error in parseGRASS(cmd, legacyExec = legacyExec) : v.import not parsed<br></b></div><div><br></div><div><br></div><div>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)</div><div><br></div><div>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.<br></div><div><br></div><div>Any idea what is happening here? </div><div><br></div><div>Thanks for your help,</div><div><br></div><div>Mehrdad</div><div> <br></div><div><br></div></div><div hspace="streak-pt-mark" style="max-height:1px"><img alt="" style="width:0px;max-height:0px;overflow:hidden" src="https://mailfoogae.appspot.com/t?sender=admFyZWRpQGdtYWlsLmNvbQ%3D%3D&type=zerocontent&guid=ea09090a-3525-4d19-8c5e-e0143cdd604c"><font color="#ffffff" size="1">ᐧ</font></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Aug 28, 2019 at 6:24 AM Micha Silver <<a href="mailto:tsvibar@gmail.com">tsvibar@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div style="direction:ltr" bgcolor="#FFFFFF">
<br>
<div class="gmail-m_7062285347713815583moz-cite-prefix">On 28/08/2019 9:56, Mehrdad Varedi
wrote:<br>
</div>
<blockquote type="cite">
<div dir="ltr">
<div dir="ltr">Hi Everyone,
<div>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.</div>
<div><br>
</div>
<div>How can I have the new SQlite database? </div>
<div><br>
</div>
</div>
</div>
</blockquote>
<p><br>
</p>
<p>If you import some vector you'll have the sqlite.db setup for
you:</p>
<p>(Following the example from before)</p>
<p><br>
</p>
<p><tt>> input_vect = "Israel_border.gpkg"</tt><tt><br>
</tt><tt>> grass_vect = "isr_border"</tt><tt><br>
</tt><tt>> execGRASS("v.import", input=input_vect,
output=grass_vect, flags="o")</tt><tt><br>
</tt><tt>Over-riding projection check</tt><tt><br>
</tt><tt>Check if OGR layer <Israel_border> contains
polygons...</tt><tt><br>
</tt><tt> 100%</tt><tt><br>
</tt><tt>Creating attribute table for layer
<Israel_border>...</tt><tt><br>
</tt><tt>Default driver / database set to:</tt><tt><br>
</tt><tt>driver: sqlite</tt><tt><br>
</tt><tt>database:
$GISDBASE/$LOCATION_NAME/$MAPSET/sqlite/sqlite.db</tt><tt><br>
</tt><tt>Importing 1 features (OGR layer <Israel_border>)...</tt><tt><br>
</tt><tt> 100%</tt></p>
<p><tt><br>
</tt></p>
<p>Then in a separate terminal:</p>
<p><tt>micha@TP480:~$ ll /tmp/tmp_location/tmp_mapset/sqlite/</tt><tt><br>
</tt><tt>total 20</tt><tt><br>
</tt><tt>drwxr-xr-x 2 micha micha 4096 Aug 28 13:19 ./</tt><tt><br>
</tt><tt>drwxr-xr-x 6 micha micha 4096 Aug 28 13:19 ../</tt><tt><br>
</tt><tt>-rw-r--r-- 1 micha micha 12288 Aug 28 13:19 sqlite.db</tt><tt><br>
</tt><tt><br>
</tt></p>
<p><br>
</p>
<blockquote type="cite">
<div dir="ltr">
<div dir="ltr">
<div>Thanks for your help.</div>
<div><br>
</div>
</div>
<br>
<div class="gmail_quote">
<div dir="ltr" class="gmail_attr">On Sun, Apr 28, 2019 at
10:55 AM Micha Silver <<a href="mailto:tsvibar@gmail.com" target="_blank">tsvibar@gmail.com</a>> wrote:<br>
</div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div style="direction:ltr" bgcolor="#FFFFFF">
<blockquote type="cite">
<div dir="ltr">
<div>
<div>Maybe this will help:<br>
</div>
</div>
</div>
</blockquote>
<p>(You must have GRASS installed, of course)</p>
<p><br>
</p>
<p><font size="+1"><tt># Parameters for the GRASS
installation and temporary GRASS mapset</tt></font></p>
<p><font size="+1"><tt>GISBase = "/usr/lib/grass76"</tt></font><br>
<font size="+1"><tt><font size="+1"><tt><span style="font-size:14px"># Set these as you like</span></tt></font></tt></font></p>
<p><font size="+1"><tt>GISDbase = "/tmp"</tt><tt><br>
</tt><tt>Location = "tmp_location"</tt><tt><br>
</tt><tt>Mapset = "tmp_mapset"</tt><tt><br>
</tt><tt>georef = "EPSG:2039"</tt><tt><br>
</tt><tt><br>
</tt><tt>mapset_path = file.path(GISDbase, Location,
Mapset)</tt></font></p>
<p><tt><font size="+1"># Now run the 'grass' command
within R and </font></tt><br>
</p>
<p><font size="+1"><tt># use the -c and -e flags to create
the temp mapset, then exit<br>
</tt></font></p>
<p><font size="+1"><tt>setup_grass_cmd = paste("grass",
"-c", georef, "-e", mapset_path)</tt></font></p>
<p><font size="+1"><tt>system(setup_grass_cmd)</tt></font></p>
<p><span style="font-size:14px"><font face="monospace"><font size="+1"><br>
</font></font></span></p>
<p><span style="font-size:14px"><font face="monospace"><font size="+1"># Load the R grass interface and
initialize GRASS within R, </font></font></span></p>
<p><span style="font-size:14px"><font face="monospace"><font size="+1"># using the temporary mapset from above</font></font></span><br>
</p>
<p><font size="+1"><tt>library(rgrass7)</tt><tt><br>
</tt><tt>initGRASS(home=tempdir(), </tt><tt><br>
</tt><tt> gisBase = GISBase,</tt><tt><br>
</tt><tt> gisDbase = GISDbase,</tt><tt><br>
</tt><tt> location = Location,</tt><tt><br>
</tt><tt> mapset = Mapset,</tt><tt><br>
</tt><tt> remove_GISRC = TRUE)</tt><tt><br>
</tt></font></p>
<p><font size="+1"><tt><br>
</tt></font></p>
<p><font size="+1"><tt># Try some GRASS commands</tt></font></p>
<p><font size="+1"><tt>input_tif = "isrlat12.tif" <br>
</tt></font></p>
<p><font size="+1"><tt>grass_rast = "isrlat12"</tt><tt><br>
</tt><tt>execGRASS("r.in.gdal", flags =
c("o","overwrite"),</tt><tt><br>
</tt><tt> input = input_tif,</tt><tt><br>
</tt><tt> output = grass_rast)</tt><tt><br>
</tt><tt>execGRASS("g.region", flags="p", raster =
grass_rast)</tt><tt><br>
</tt></font><br>
</p>
<p>I did not add a command to remove the temporary
mapset. Something like:</p>
<p> <tt>unlink(file.path(GISDbase, Location), recusive =
TRUE)</tt> might be necessary</p>
<p>Micha Silver<br>
</p>
<pre class="gmail-m_7062285347713815583gmail-m_-7760652768775255896moz-signature" cols="72">Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918</pre>
</div>
</blockquote>
</div>
<br clear="all">
<div><br>
</div>
-- <br>
<div dir="ltr" class="gmail-m_7062285347713815583gmail_signature">Mehrdad Varedi</div>
</div>
<div hspace="streak-pt-mark" style="max-height:1px"><img alt="" style="width: 0px; max-height: 0px; overflow: hidden;" src="https://mailfoogae.appspot.com/t?sender=admFyZWRpQGdtYWlsLmNvbQ%3D%3D&type=zerocontent&guid=21afc582-8b6f-47c9-ba1d-ef54adaa0dfb"><font size="1" color="#ffffff">ᐧ</font></div>
</blockquote>
<pre class="gmail-m_7062285347713815583moz-signature" cols="72">--
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918</pre>
</div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature">Mehrdad Varedi</div>