Edzer Pebesma

Thanks! That should now work.

Thanks, @henningte ! I think this commit addresses all of them. The `to_sftime` was indeed an unfinished bit, and had to be set on depending on the class of `newdata`....

You may want to look into the docs & examples of package `spacetime`, for that.

spurious.

Potentially great to have, but depending on packages not on CRAN...

Cool idea; I've been thinking about the idea of handling tile collections quite a bit, lately. A set of tiles for [sentinel 2](https://sentinel.esa.int/web/sentinel/missions/sentinel-2/data-products) is loaded by: ```r tiles = st_read("https://sentinel.esa.int/documents/247904/1955685/S2A_OPER_GIP_TILPAR_MPC__20151209T095117_V20150622T000000_21000101T000000_B00.kml")...

Something along these lines? ```r > ring = rbind(c(0,0), c(1,1), c(0,1), c(0,0)) > library(sf) Linking to GEOS 3.7.0, GDAL 2.3.2, PROJ 5.2.0 > st_is_valid(st_polygon(list(ring, ring))) [1] FALSE > st_is_valid(st_polygon(list(ring, ring)),...

It works if you add a small value to all direct variogram sills, by ```r meuse.fit = fit.lmc(x, meuse.i, correct.diagonal=1.01) ``` If we don't do this, we get a case...

Yes, it looks like `st_intersects` could do this.

The error message you get clearly points to a bug in gstat, but what you're trying to do is essentially trying to compute residuals from a regression model with two...