[R-sig-Geo] Store sf as 'line-by-line' GeoJSON

)) This writes every feature on a single line. However, it doesn't store them as a FeatureCollection. Suggestions are welcome. Best regards, ir. Thierry Onkelinx Statisticus / Statistician Vlaamse Overheid / Government of Flanders INSTITUUT VOOR NATUUR- EN BOSONDERZOEK / RESEARCH INSTITUTE FOR NA

Re: [R-sig-Geo] GRTS sampling - 2-level design

is a reimplementation of my GRTS package (https://github.com/ThierryO/grts). That package has a vignette describing GRTS via the Reversed Randomized Quadrant-Recursive Raster strategy ( https://doi.org/10.1007/s00267-005-0199-x). Best regards, ir. Thierry Onkelinx Statisticus / Statistician

Re: [R-sig-Geo] GRTS sampling - 2-level design

his allows for repeated measures for stable locations while taking into account the changes in population. Best regards, Thierry ir. Thierry Onkelinx Statisticus / Statistician Vlaamse Overheid / Government of Flanders INSTITUUT VOOR NATUUR- EN BOSONDERZOEK / RESEARCH INSTITUTE FOR NATURE AND FOREST Team

Re: [R-sig-Geo] Impute missing values along a spatial network

that problem. Best regards, ir. Thierry Onkelinx Statisticus / Statistician Vlaamse Overheid / Government of Flanders INSTITUUT VOOR NATUUR- EN BOSONDERZOEK / RESEARCH INSTITUTE FOR NATURE AND FOREST Team Biometrie & Kwaliteitszorg / Team Biometrics & Quality Assurance thierry.onkel...

[R-sig-Geo] Importing layers from Google My Maps

Dear all, Can I import the layers from a Google My Maps based on it's link? I know I can go to the link, export the layers as KML and import that in R using rgdal. But that is tedious when you have a lot of links. Any suggestions? Best regards, ir. Thierry Onkelinx Statisticus / Statistician

Re: [R-sig-Geo] Worked out example of spatiotemporal interpolation of daily precipitation?

Dear Tomislav, Blangiardo and Cameletti (2015) Spatial and Spatio-temporal Bayesian Models with R - INLA models the Parana state rainfall data (chapter 8). See https://sites.google.com/a/r-inla.org/stbook/ Not kriging but maybe useful for you. Best regards, ir. Thierry Onkelinx Statisticus

Re: [R-sig-Geo] How to fit Spatial logistic regression models to DHS data

Have a look at Zuur et al (2017) Beginner's Guide to Spatial, Temporal and Spatial-Temporal Ecological Data Analysis with R-INLA Best regards, ir. Thierry Onkelinx Statisticus / Statistician Vlaamse Overheid / Government of Flanders INSTITUUT VOOR NATUUR- EN BOSONDERZOEK / RESEARCH INSTITUTE

Re: [R-sig-Geo] testing nlme models for spatial correlation

for the SPDE approach in INLA. I would recommend that you consult a (local) statistician. Best regards, ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest team Biometrie & Kwaliteitszorg / team Biometrics & Quality Assurance Kliniekstraat

Re: [R-sig-Geo] testing nlme models for spatial correlation

Dear Javier, Your problem is hard to understand without a reproducible example. You only gives the code, not the data nor the error message. Best regards, ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest team Biometrie & Kwaliteits

Re: [R-sig-Geo] testing nlme models for spatial correlation

with R-INLA. Best regards, ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest team Biometrie & Kwaliteitszorg / team Biometrics & Quality Assurance Kliniekstraat 25 1070 Anderlecht Belgium To call in the statistician after the ex

Re: [R-sig-Geo] Problem with Proj when loading 'rgdal' on Super Computer Facility

Dear Ryan, You **first** need to install the dependencies of rgdal. sudo apt-get update && sudo apt-get install libgdal-dev libproj-dev Then you should be able to install rgdal. For more details google "rgdal install linux" Best regards, ir. Thierry Onkelinx Institu

Re: [R-sig-Geo] Problems With Cartography - Help

Dear Andres, You'll need spTransform() to reproject the layer into WGS84 carto2012x <- unionSpatialPolygons(carto2012, IDs=carto2012$DPA_CANTON) spTransform(carto2012x, CRS("+proj=longlat")) Best regards, ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Resea

Re: [R-sig-Geo] How to fit a spatial hurdle model with the INLA package ?

I wonder whether it would be possible to do something with zero-inflated models of type 0: see http://www.math.ntnu.no/inla/r-inla.org/doc/likelihood/zeroinflated.pdf, first equation on page 1. All the best Renaud Le 17/02/2017 à 15:29, Thierry Onkelinx a écrit : > Dear Renaud, > > IMHO you c

Re: [R-sig-Geo] How to fit a spatial hurdle model with the INLA package ?

Dear Renaud, IMHO you can't. INLA currently only fits zero-inflated distributions (with a single zero-inflation parameter). A hurdle model would require to fit a zero-truncated distribution. Best regards, ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research Institute

Re: [R-sig-Geo] Plotting simple features sf with ggplot2

regards, ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest team Biometrie & Kwaliteitszorg / team Biometrics & Quality Assurance Kliniekstraat 25 1070 Anderlecht Belgium To call in the statistician after the experiment is done may be no m

Re: [R-sig-Geo] Converting csv files to SpatialPolygonsDataFrame

Dear Ahmed, A SpatialPolygons object is more complex than just a set of coordinates. I suggest that you read the helpfile of SpatialPolygonsDataFrame. The example in the helpfile for the SpatialPolygonsDataFrame class should be clear enough. Best regards, ir. Thierry Onkelinx Instituut voor

Re: [R-sig-Geo] Installing packages of R 3.2.3 in Linux cluster computers

Dear John, You'll need to install raster from source. Binaries are only available for the previous, current and development versions of R. Another option is that you are missing dependencies. We need to full error message to see what's wrong. Best regards, ir. Thierry Onkelinx Instituut voor

Re: [R-sig-Geo] Questions about spatio-temporal modelling

Dear Luca, Have a look at Blangiarde & Cameletti (2015) Spatial and Spatio-temporal Bayesian Models with R - INLA ISBN: 978-1-118-32655-8 They describe how you can tackle this problem with mixed models with correlated random effects. Best regards, ir. Thierry Onkelinx Instituut voor natuur

Re: [R-sig-Geo] geocoding street addresses within R, not using any web-based services

Dear Christopher, Have a look at the rgeos and spgrass6 packages. rgeos provides several GIS operations. spgrass6 allows to move data from R to GRASS GIS (and vice versa). You can run GRASS operations from within R. Best regards, Thierry ir. Thierry Onkelinx Instituut voor natuur- en

Re: [R-sig-Geo] help with 2d- Density plots (ggplot2)

Dear Amao, It's not clear how the "hundreds of species" come into the play. Do you want a map for each species? Or rather a single map that combines the species info. ggplot2() has AFAIK no options to restrict the output of stat_bin2d() to a polygon. Best regards, ir. Thierr

Re: [R-sig-Geo] is this a clash between wrong version of rgdal and sp?

Dear Frede, Why not installing INLA on a recent R version? I'm running INLA on R 3.2.2. Best regards, ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest team Biometrie & Kwaliteitszorg / team Biometrics & Quality Assurance Klinieks

Re: [R-sig-Geo] Rasterize SpatialPolygons in R is ignoring small values

Please don't post in HTML and don't cross post (http://stackoverflow.com/questions/32652814/rasterize-spatialpolygons-in-r-is-ignoring-small-values). ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest team Biometrie & Kwaliteitszorg /

Re: [R-sig-Geo] polygon neighbor matrix

Maybe gTouches() from the rgeos package. ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest team Biometrie & Kwaliteitszorg / team Biometrics & Quality Assurance Kliniekstraat 25 1070 Anderlecht Belgium To call in the statistici

Re: [R-sig-Geo] Modeling areal data with lots of holes, islands

- C - D - E. The correlation between two neighbours in rho. If C has no data, then INLA will estimate it's effect so that the correlation between B-C and C-D is rho. As a result, the correlation between B-D will be rho ^ 2. Best regards, ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek

Re: [R-sig-Geo] Problem in geoR

It would be easier for us to help you if the example is reproducible. ​ ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest team Biometrie Kwaliteitszorg / team Biometrics Quality Assurance Kliniekstraat 25 1070 Anderlecht Belgium To call

Re: [R-sig-Geo] variogram problem

You'll need to provide a reproducible example. Your code is not reproducible without the data. ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest team Biometrie Kwaliteitszorg / team Biometrics Quality Assurance Kliniekstraat 25 1070

Re: [R-sig-Geo] Calculate shortest distance between points belonging to different polygons

(x)[apply(x, 1, which.min)] ) }) Best regards, Thierry ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest team Biometrie Kwaliteitszorg / team Biometrics Quality Assurance Kliniekstraat 25 1070 Anderlecht Belgium To call in the statistician

Re: [R-sig-Geo] Calculate shortest distance between points belonging to different polygons

Can you provide a reproducible example? The output of as.data.frame(nndist(randp.ppp, by = randp.ppp$Polygon_ID)) should be a data.frame with one row per point and one column per polygon. Best regards, Thierry ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research Institute