Dan Turenne | 1 Jul 23:47 2016

determining values for nmax, nmin and maxdist

Hello listers,

As part of the research I am doing I'm attempting to fit some kriging models, and I was wondering if there are
any "rule of thumb" guidelines for selecting the values for the nmax, nmin and maxdist parameters? 
Intuitively I would set maxdist equal to the range of the variogram but I'm not sure about the other
parameters.  Any advice on this subject would be appreciated.


Daniel Turenne

University of Manitoba

	[[alternative HTML version deleted]]
Fischbach, Anthony | 1 Jul 20:48 2016

projectRaster across the antimeridian

The solution to the previous post addressing projection of rasters across
the antimeridian
no-longer functions
(as of 24 June 2016).

Below I offer permutations of the solution offered in the previous post,
to no avail.

Is this a bug in the raster package, or may my code be adjusted to produce
a projected raster that includes both the western and eastern hemisphere of
the original image?


## Define projections
prj.geo<-CRS('+proj=longlat +datum=WGS84')
prj.StudyArea <- CRS("+proj=aeqd +lat_0=70 +lon_0=-170")
prj.StudyAreaPositive <- CRS("+proj=aeqd +lat_0=70 +lon_0=190")
## Download geoJPEG
Yesterday <- Sys.Date() -1
cat(URL<- format(Yesterday, '
cat(DestFile<-format(Yesterday,'BeringSea.%Y%j.terra.721.2km.jpg' ,
(Continue reading)

Dan Turenne | 30 Jun 01:42 2016

Re: Spatio Temporal kriging in Gstat

Thank you very much for the tip Edzer, that was exactly what I was looking for.

Daniel Turenne

From: Edzer Pebesma <edzer.pebesma <at> uni-muenster.de>
Sent: June 29, 2016 6:17 PM
To: Dan Turenne
Subject: Re: [R-sig-Geo] Spatio Temporal kriging in Gstat

On 29/06/16 11:52, Dan Turenne wrote:
> Hello Ben,
> Thank you for the tip about "tunit", it worked perfectly as you described.  As for my second question, the
object returned by krigeST is indeed of the form STFDF.  After some digging I found the following
information regarding the data from an STFDF object on inside-r.org:
> data<http://inside-r.org/r-doc/utils/data>:
data {utils} | inside-R | A Community Site for R<http://inside-r.org/r-doc/utils/data>
Arguments... literal character strings or names. list a character vector. package a character vector
giving the package(s) to look in for data sets, or NULL.

> Object of class data.frame<http://inside-r.org/r-doc/base/data.frame>, which holds the measured
values; space index cycling first, time order preservedI would interpret this as follows: if I had x
(Continue reading)

Lars-Daniel Weber | 29 Jun 23:02 2016

analysis of point patterns

Hi there!

I've got a road network (polylines) with points on it as geodata. The points have quality measurements of
the paving within a distance of 10 to 20 metres. I'd like to run a analysis on these points to detect
autocorrelation and hot spots with good and bad paving quality: high-high (etc). Local Moran's I might be
a good solution from my point of view. 

I thought of creating 5 classes of quality from 1 (bad) to 6 (perfect) and put it into Local Moran's I. Is this a
good approach to get some analysis out of the data?

Which packages would you suggest? Or would it be better to write it on my own? I can export the point data to
LAT;LON;CLASS to do easy processing in R.

Best regards,
Lars-Daniel Weber
Isaque Daniel | 29 Jun 03:26 2016

How to create an ellipse buffer using SpatialPointsDataFrame?

Hi everyone,

I have been working in a model for spatial analysis and my problem now is generate some ellipses around of SpatialPointsDataFrame.

I looking for some similar gBuffer of rgeos, but until now I can't find a solution.

Any suggestion?

Thanks in advance


Agronomist engineer
Master in Remote Sensing - National  Institute for Space Research (INPE) - Brazil
PHD Student in Transport - Brazilia University (UNB)

	[[alternative HTML version deleted]]
Andy Teucher | 28 Jun 19:45 2016

FYI - bug in GDAL 2.1.0: FID detection in GeoJSON results wrong feature count

There is a small known bug in GDAL 2.1.0 that I was bitten by when using
rgdal::readOGR to read in GeoJSON data.

The bug is with the GeoJSON format, where if you have multiple polygons
with the same "id" property, only one will be recognized (The actual bug is
more complex than that, but that was my circumstance).

The bug report is here: https://trac.osgeo.org/gdal/ticket/6538, and my
demonstration of it with rgdal is here:

It looks like it has been fixed in GDAL 2.1.1. I tested it with GDAL 2.0.1
and the current development version (2.2.x) and the issue is not there, but
I think 2.1.0 is common on a lot of systems right now. So if you are
working with GeoJSON and have that GDAL version, keep your eye out for it.

Andy Teucher

	[[alternative HTML version deleted]]
Maurizio Marchi | 28 Jun 17:16 2016

Absences dataset

Hello everybody,
I have a question about species distribution modelling and the biomod2
I have been working with the biomod2 package many times and often using the
PA selection function to create more than one pseudo-absences dataset. This
was mainly due to the lack of available information about true absences.
Now I'm i a  very good position because I have to work wit both presences
and absences but a new problem arose: the number of absences is much higher
than the presences (e.g 700 presences versus 6000 absences). So I would
like to calculate more than one run to average the results, similarly to
what I'm use to do with the pseudo-absences... So, my question is: is there
the possibility to do it automatically with the biomod2 package or should I
have to create a cycle to repeat the calculation and average the resulting
maps manually?

Many thanks, maurizio


Maurizio Marchi
Calenzano (FI) - Italy
ID Skype: maurizioxyz
Ubuntu 14.04 LTS
linux user 552742

	[[alternative HTML version deleted]]
Csany22 csany | 28 Jun 12:07 2016

spatial anova

Hi R-sig-geo people,

what would be the best way to run spatial anova in R? I have a dependent (continuous), independent
(categorical) variables and lat and long fo reach of the cell(pixel) based observations. I need to run
anova. as would like to compare the means (Tukey) of the independent variable. 
	[[alternative HTML version deleted]]
Edoardo Baldoni | 26 Jun 15:10 2016

The spml function does not work

Dear R users,

I am trying to benchmark my code for estimating spatial panels with the
functions in the splm package. The problem is that I cannot make the spml
functions (from package splm) work. Not even the ones in the standard
examples provided. Here is one of the examples that can be found in the R
help page ( and then the
error that comes out:

data(Produc, package = "plm")
fm <- log(gsp) ~ log(pcap) + log(pc) + log(emp) + unemp
fespaterr <- spml(fm, data = Produc, listw = mat2listw(usaww),
                   model="within", spatial.error="b", Hess = FALSE)

Error in solve(crossprod(xco), crossprod(xco, yco)) :
  error in evaluating the argument 'a' in selecting a method for function
'solve': Error in crossprod(xco) :
  could not find symbol "..." in environment of the generic function

Do you have any idea of what is happening here ? Thanks


	[[alternative HTML version deleted]]
Jefferson Ferreira-Ferreira | 24 Jun 15:45 2016

[off-topic] Raster Distance

Dear listers!

Sorry by this off-topic.
I have an open question on GIS StackExchange: Is there any difference
between gdal_proximity.py and ArcGIS Euclidean Distance?
If someone can help to answer it, I would be grateful.

Thanks a lot.

*Jefferson Ferreira-Ferreira, **PhD Candidate *


*Ecosystem Dynamics Observatory <http://tscanada.wix.com/ecodyn> -
*Department of **Geography *
*Institute of Geosciences and Exact Sciences** (IGCE)*
*São Paulo State University (UNESP)*
*Rio Claro, SP - Brazil*

	[[alternative HTML version deleted]]

R-sig-Geo mailing list
R-sig-Geo <at> r-project.org
Debora Nozza | 24 Jun 16:22 2016

NA data values

Good morning,

I am working on creating a STFDF object using the function meteo2STFDF(). I
am trying to  to reproducing (with my data) the same example provided in

However, after computing meteo2STFDF() function, the data <at> tempc is full of
NA values. It seems that I am doing exactly the same thing of the example
so I am not able to understand where is the problem.

Could someone has a solution?

Best Regards,

Here my variable:

Object of class STFDF
 with Dimensions (s, t, attr): (2, 8761, 1)
Object of class SpatialPointsDataFrame
           min        max
lon 14406.0000 14476.0000
lat    45.6737    45.8771
Is projected: FALSE
proj4string :
[+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0]
Number of points: 2
Data attributes:
(Continue reading)