[R-sig-Geo] RQGIS release 0.2.0

Ali Santacruz amsantac at hotmail.com
Sat Jan 21 20:43:22 CET 2017


Dear Patrick,


Thank you very much for your reply.


You're right. It must be an issue with my file. After your email, I tried the example with the dem dataset and it worked well,


Thanks for the other suggestions. I'll keep them in mind,


Best regards,


Ali


________________________________
De: Patrick Schratz <patrick.schratz at gmail.com>
Enviado: sábado, 21 de enero de 2017 12:16:10 p. m.
Para: Ali Santacruz
Cc: Jannes Münchow; r-sig-geo at r-project.org
Asunto: Re: [R-sig-Geo] RQGIS release 0.2.0

Dear Ali,

a general note: Please write a new mail if you have a specific problem with a package instead of replying to the last mail related to the package.

Regarding your issue: I think this error relates to your file and its missing projection and extent information.

I just tried the code using the RQGIS "dem" dataset and everything worked fine. Also, I´m not sure if "GRASS_REGION_PARAMETER" and "GRASS_REGION_CELLSIZE_PARAMETER" need to be set at all or you should better go with the defaults.

If you want to quickly check the algorithm help you can use `open_help("grass7:r.neighbors")`.

library(RQGIS)
library(raster)
r2 <- data("dem")

my_env <- set_env()

params <- get_args_man(alg = "grass7:r.neighbors", qgis_env = my_env)
params$input  <- dem
params$method  <- 2 # "Mode"
params$output  <- file.path(tempdir(), "rst_out.tif")

out <- run_qgis(alg = "grass7:r.neighbors", params = params,
                load_output = params$output, qgis_env = my_env)
out

------

macOS, QGIS 2.18.2, GRASS 7.0

You may also want to post such issues at the Github repository<https://github.com/jannes-m/RQGIS/issues> or at Stackoverflow/GISexchange using the tag 'rqgis' in the future.
Issues like this which feature large code blocks and are related to a specific R package can be much better addressed there than by mail.

Best, Patrick





2017-01-21 15:51 GMT+01:00 Ali Santacruz <amsantac at hotmail.com<mailto:amsantac at hotmail.com>>:

[[elided Hotmail spam]]


I've been able to run QGIS and SAGA algorithms without issues. However, I have some problems running GRASS algorithms,

I'm trying to run r.neighbors from GRASS 7, but the resolution, extent and values of the output raster don't match those of the input,


Here is the code:


# Load raster
library(raster)

r1 <- raster("pred_m1_3000_recl5.tif")
r1


class       : RasterLayer
dimensions  : 837, 849, 710613  (nrow, ncol, ncell)
resolution  : 1, 1  (x, y)
extent      : 0, 849, -837, 0  (xmin, xmax, ymin, ymax)
coord. ref. : NA
data source : E:\amsantac\dstl\working\pred\pred_m1_3000_recl5.tif
names       : pred_m1_3000_recl5
values      : 5, 110  (min, max)


# Load RQGIS
library(RQGIS)
my_env <- set_env("C:/OSGeo4W64/")
dir_tmp <- tempdir()


# Set parameters for grass7:r.neighbors
params <- get_args_man(alg = "grass7:r.neighbors", qgis_env = my_env)
params$input  <- r1
params$method  <- 2 # "Mode"
params$output  <- file.path(dir_tmp, "rst_out.tif")

params$GRASS_REGION_PARAMETER <- "0,849,-837,0"    # as xmin, xmax, ymin, ymax
params$GRASS_REGION_CELLSIZE_PARAMETER <- "1.0"

out <- run_qgis(alg = "grass7:r.neighbors", params = params, load_output = params$output, qgis_env = my_env)
out


class       : RasterLayer
dimensions  : 1, 7, 7  (nrow, ncol, ncell)
resolution  : 0.1519435, 1  (x, y)
extent      : 0, 1.063604, 0, 1  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0
data source : C:\Users\Guest\AppData\Local\Temp\RtmpwJu31Y\rst_out.tif
names       : rst_out



Note the output only has 1 row, 7 cols (vs. 837 rows, 849 cols in input). The CRS has been set as long-lat, while the CRS is NA in the input,


Also note that the resolution has been changed from 1 to 0.1519,


I wonder if there is an issue on how the params$GRASS_REGION_PARAMETER is defined. I have tried different ways eg., "0,849,-837,0" (as xmin, xmax, ymin, ymax), "0,-837,849,0" (as xmin, ymin, xmax, ymax), extent(r1) (which gives error).


The input image (550 kb) can be downloaded from this link: https://www.dropbox.com/s/do95t48foepw62e/pred_m1_3000_recl5.tif?dl=0



Thanks for your help,


Best regards,


Ali



sessionInfo()


R version 3.3.2 (2016-10-31)
Platform: x86_64-w64-mingw32/x64 (64-bit)


locale:
[1] LC_COLLATE=English_United States.1252  LC_CTYPE=English_United States.1252    LC_MONETARY=English_United States.1252
[4] LC_NUMERIC=C                           LC_TIME=English_United States.1252

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base

other attached packages:
[1] RQGIS_0.2.0.9000     raster_2.5-8         sp_1.2-3             RevoUtilsMath_10.0.0

loaded via a namespace (and not attached):
[1] RevoUtils_10.0.2 rgdal_1.1-10     tools_3.3.2      Rcpp_0.12.7      grid_3.3.2       lattice_0.20-34


________________________________
De: R-sig-Geo <r-sig-geo-bounces at r-project.org<mailto:r-sig-geo-bounces at r-project.org>> en nombre de "Jannes Münchow" <malNamalJa at gmx.de<mailto:malNamalJa at gmx.de>>
Enviado: jueves, 19 de enero de 2017 06:31:34 a. m.
Para: r-sig-geo at r-project.org<mailto:r-sig-geo at r-project.org>
Asunto: [R-sig-Geo] RQGIS release 0.2.0

Dear list members,

today we are happy to announce a new version of RQGIS! RQGIS establishes an interface between R and QGIS, i.e. it allows the user to access the more than 1000 QGIS geoalgorithms from within R. The most important news is that you can run again the most recent QGIS releases (>=2.18.2) with RQGIS without having to modify Python source code. This modification became necessary due to a bug introduced in QGIS 2.16 (see https://github.com/jannes-m/RQGIS#rqgis-modifications-v-216-2181).

[https://avatars0.githubusercontent.com/u/9986952?v=3&s=400]<https://github.com/jannes-m/RQGIS#rqgis-modifications-v-216-2181>

GitHub - jannes-m/RQGIS: RQGIS - integrating R with QGIS<https://github.com/jannes-m/RQGIS#rqgis-modifications-v-216-2181>
github.com<http://github.com>
RQGIS - integrating R with QGIS ... Now that we have a spatial object, we can move on to using RQGIS. First of all, we need to specify all the paths necessary to run ...




This release fixes also several bugs which we detected thanks to user feedback. To learn more about these, please refer to https://github.com/jannes-m/RQGIS/releases/tag/v0.2.0.

[https://avatars0.githubusercontent.com/u/9986952?v=3&s=400]<https://github.com/jannes-m/RQGIS/releases/tag/v0.2.0>

jannes-m/RQGIS<https://github.com/jannes-m/RQGIS/releases/tag/v0.2.0>
github.com<http://github.com>
RQGIS - integrating R with QGIS




For more information on RQGIS, how to install it and examples how to use RQGIS, please refer to https://github.com/jannes-m/RQGIS and my https://jannesm.wordpress.com/category/r/.

Posts about R on jannesm<https://jannesm.wordpress.com/category/r/>
jannesm.wordpress.com<http://jannesm.wordpress.com>
Posts about R written by jannesm



[https://avatars0.githubusercontent.com/u/9986952?v=3&s=400]<https://github.com/jannes-m/RQGIS>

GitHub - jannes-m/RQGIS: RQGIS - integrating R with QGIS<https://github.com/jannes-m/RQGIS>
github.com<http://github.com>
RQGIS - integrating R with QGIS ... Now that we have a spatial object, we can move on to using RQGIS. First of all, we need to specify all the paths necessary to run ...




All the best,

Jannes

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo at r-project.org<mailto:R-sig-Geo at r-project.org>
https://stat.ethz.ch/mailman/listinfo/r-sig-geo

R-sig-Geo Info Page - ETH Zurich<https://stat.ethz.ch/mailman/listinfo/r-sig-geo>
stat.ethz.ch<http://stat.ethz.ch>
R-sig-Geo -- R Special Interest Group on using Geographical data and Mapping About R-sig-Geo




        [[alternative HTML version deleted]]


_______________________________________________
R-sig-Geo mailing list
R-sig-Geo at r-project.org<mailto:R-sig-Geo at r-project.org>
https://stat.ethz.ch/mailman/listinfo/r-sig-geo


	[[alternative HTML version deleted]]



More information about the R-sig-Geo mailing list