Re: [R-sig-Geo] gIntersection returns error "TopologyException: no outgoing dirEdge found at"

2015-09-25 Thread Edzer Pebesma
"rgeos_intersection") : >>>>> TopologyException: no outgoing dirEdge found at 0.5 0 >>>> >>>> The error message is coming from GEOS - you are welcome to investigate >>>> further. If you use gUnaryUnion() on the objects first, there is no >>>> error: >>>> >

Re: [R-sig-Geo] Take mean of list of raster stacks

2015-12-02 Thread Edzer Pebesma
0 >>>>> data source : in memory >>>>> names : layer.1, layer.2, layer.3, layer.4, layer.5, layer.6, layer.7, >>>>> layer.8, layer.9 >>>>> min values : 211.6625, 185.8265, 187.7064, 187.3369, 186.3985, >>>>> 149.3203, >>>>> 156.6462, 153.4485,

Re: [R-sig-Geo] Take mean of list of raster stacks

2015-12-02 Thread Edzer Pebesma
+ (3 + 4)) on a julia interpreter. We need something > like this for R, especially for spatial objects. > > Regardless, thanks for your explanation. It's very helpful. > Cordially, > Vijay. > > On Wed, Dec 2, 2015 at 12:50 PM, Edzer Pebesma > <edzer.pebe...@uni-muenst

Re: [R-sig-Geo] Extract points from SpatialLines objects along with the corresponding Lines' ID

2015-12-05 Thread Edzer Pebesma
> If possible I would like to avoid type conversions and get a data.frame > directly though I could indeed recreate one with these matrices, without > information loss I assume. > > Thank you in advance for any tips. > > Best regards, > Bastien > > [[alterna

Re: [R-sig-Geo] SpatialMultiPoints suggestions

2015-12-05 Thread Edzer Pebesma
s"))) > df = do.call(rbind, lapply(dots, function(x) x@data)) > SpatialMultiPointsDataFrame(sp, df) > } > > > HTH, > Kent > > [[alternative HTML version deleted]] > > _______ > R-sig-Geo mailing list &

Re: [R-sig-Geo] Excessive whitespace in saved images from sp::plot()

2015-12-16 Thread Edzer Pebesma
posting here. Since these are simple polygons with no >>>> attributes, I'd prefer to avoid using spplot(). Also, I've tried this >>>> on 2 machines: a Mac and a PC running Ubuntu, in both cases running >>>> the latest version of R and the sp package. >>>> >

Re: [R-sig-Geo] spml package and impacts calcultion

2015-12-23 Thread Edzer Pebesma
; https://stat.ethz.ch/mailman/listinfo/r-sig-geo > -- Edzer Pebesma Institute for Geoinformatics (ifgi), University of Münster Heisenbergstraße 2, 48149 Münster, Germany; +49 251 83 33081 Journal of Statistical Software: http://www.jstatsoft.org/ Computers & Geosciences: http://elsevi

Re: [R-sig-Geo] gstat::krige() - regression kriging vs. kriging with external drift

2015-11-25 Thread Edzer Pebesma
al linear model y = X beta + e, with a spatial covariance function describing the covariance of e, and in that perspective both methods are equivalent. I'm not the person to tell what regression kriging is and what it is not, but it might include the above models as well as those where trend fitting and res

Re: [R-sig-Geo] Match polygon and dataframe IDs after raster::extract

2015-11-19 Thread Edzer Pebesma
ataFrame() doesn't. (I made a mistake in my initial > example, I want to create a SpatialPolygonsDataFrame, not a > SpatialPointsDataFrame). as documented, SpatialPolygonsDataFrame has argument match.ID default to TRUE. > > Cheers, > Loïc > > > On 11/19/2015 03:08 PM, E

Re: [R-sig-Geo] space-time - efficient creation of a TracksCollection

2015-11-18 Thread Edzer Pebesma
his as a way to append a new STI to > appropriate Tracks and to the collection? > > Any hint welcome. > > Thanks, > > Martin > > -- Edzer Pebesma Institute for Geoinformatics (ifgi), University of Münster, Heisenbergstraße 2, 48149 Münster, Germany; +49 251 83

Re: [R-sig-Geo] Match polygon and dataframe IDs after raster::extract

2015-11-19 Thread Edzer Pebesma
ys) > # [1] "1" "2" "3" > getPolyID(polys[-2]) > # [1] "1" "3" > > unique(extract(b, polys[-2], df = TRUE)$ID) > # [1] 1 2 > > > Any suggestions? > > Thanks, > Loïc > > _

Re: [R-sig-Geo] Regarding installation of R package rgdal

2016-01-13 Thread Edzer Pebesma
help > *** installing help indices > ** building package indices > ** installing vignettes > ** testing if installed package can be loaded > * DONE (rgdal) > > The command I used for installing the package is "install.package('rgdal', > repos='http://cran.r-project.org')". I don't know why R package

Re: [R-sig-Geo] Universal Block Kriging covariate definition for krige in gstat

2016-01-13 Thread Edzer Pebesma
, shapefile, mean) > > Thanks in advance, > I hope I could explain it properly, but I will give more details if > necessary. > Kind regards, > Antonio > > [[alternative HTML version deleted]] > > ___ > R-sig-Geo

Re: [R-sig-Geo] multivariate sequential gaussian cosimulation using gstat in R

2016-06-06 Thread Edzer Pebesma
Which command generates which warning?), and is not reproducible. The way you specify conditioning data is correct, so it is the basis for conditioning your simulations. The warnings should not be ignored, though: have you tried plotting the variograms & cross variograms with the (wrongly?) fitt

Re: [R-sig-Geo] Overlay between polygons and their intersecting lines

2016-05-28 Thread Edzer Pebesma
n (read the help of rgeos::gRelate) needs to be at least one, so touching (a point: 0-dimensional) is ignored, and polygon A would be returned. minDimension is now 0 by default; other values might make more sense, depending on the geometries at hand. A problem is that this increases the computa

Re: [R-sig-Geo] Overlay between polygons and their intersecting lines

2016-05-27 Thread Edzer Pebesma
# The green (1st) segment intersects with the black (1st) polygon => good > > over(Inter[2], PolyD, byid=T) > > 'Attribute > > 11black' > > # The blue(2nd) segment intersects with the red(2nd) polygon => NOT GOOD ! > > over(Inter[1], PolyD, byid=T, returnList = T

Re: [R-sig-Geo] experimental covariogram

2016-06-22 Thread Edzer Pebesma
> R-sig-Geo mailing list > R-sig-Geo@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-sig-geo > -- Edzer Pebesma Institute for Geoinformatics (ifgi), University of Münster Heisenbergstraße 2, 48149 Münster, Germany; +49 251 83 33081 Journal

Re: [R-sig-Geo] Spatio Temporal kriging in Gstat

2016-06-21 Thread Edzer Pebesma
help. > Anything that avoids sending large files to all of the 3500 subscribers of this list: direct email is fine, sharing on dropbox or similar too. > > Daniel Turenne > > University of Manitoba > > > > From: R-sig-Geo <r-sig-geo-

Re: [R-sig-Geo] Spatio Temporal kriging in Gstat

2016-06-21 Thread Edzer Pebesma
> > Can anyone see what I am doing wrong or give me any pointers? This error is > rather cryptic and I'm not quite sure what I'm doing wrong. Any help would > be appreciated. > > > Many Thanks, > > Daniel Turenne > > University of Manitoba > > > [[alternative

Re: [R-sig-Geo] Using GStat for kriging with Anisotropy

2016-06-15 Thread Edzer Pebesma
> University of Manitoba > > [[alternative HTML version deleted]] > > ___ > R-sig-Geo mailing list > R-sig-Geo@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-sig-geo > -- Edzer Pebesma Institute for Geoinf

Re: [R-sig-Geo] Anisotropy for Cokriging

2016-01-10 Thread Edzer Pebesma
versity > P.O.Box 46417-76489 > Noor > Mazandaran > Iran > Mobile Phone Number 00989111580097 > Tel:0098(122)6253101-3 > http://www.modares.ac.ir/en/Schools/nat/Academic_Staff/~j.alavi > [[alternative HTML version deleted]] > > ________

Re: [R-sig-Geo] Universal Block Kriging covariate definition for krige in gstat

2016-01-13 Thread Edzer Pebesma
2$kr1 [1] 1.00048 > My guess is that the difference can be attributed to how the area is discretized (see ?predict.gstat) > > kind regards > > Antonio Manuel Moreno Rodenas > > /Marie Curie Early Stage Researcher/ > /PhD Candidate/ > > *T**U **Delft / Section Sanitary En

Re: [R-sig-Geo] Universal Block Kriging covariate definition

2016-01-14 Thread Edzer Pebesma
On 14/01/16 12:35, Bruin, Sytze de wrote: > I believe the residual variogram should then be computed using the covariate > data at block support. > > Sytze de Bruin > Wageningen University > Laboratory of Geo-Information Science and Remote Sensing why? -- Edzer

[R-sig-Geo] sp and gstat dev moved to github

2016-01-18 Thread Edzer Pebesma
I moved the development trees of sp and gstat to github. Feel free to send issues or pull requests there. For automatically created windows binaries, I will update r-forge upon request only. https://github.com/edzer/sp https://github.com/edzer/gstat Best regards, -- Edzer Pebesma Institute

Re: [R-sig-Geo] Ordinary kriging variance of the prediction with error structure

2016-02-05 Thread Edzer Pebesma
;> measure >>> at that location? >>> >>> Kind regards, >>> >>> Antonio >>> >>> [[alternative HTML version deleted]] >>> >>> ___ >>> R-sig-Geo mailing list &

Re: [R-sig-Geo] extract mean raster value from polygon

2016-03-11 Thread Edzer Pebesma
> > return(blpi) > } > > my_blpi <- create_perp_buffers(18, 331500, 181000, 332500, 100, 100) > > data(meuse.grid) > coordinates(meuse.grid) <- ~x+y > gridded(meuse.grid) <- TRUE > proj4string(meuse.grid) <- > CRS(paste("+init=epsg:2

Re: [R-sig-Geo] readOGR not reading in prj file?

2016-04-06 Thread Edzer Pebesma
On 06/04/16 19:33, Vinh Nguyen wrote: > > On Apr 6, 2016 10:27 AM, "Edzer Pebesma" <edzer.pebe...@uni-muenster.de > <mailto:edzer.pebe...@uni-muenster.de>> wrote: >> >> What are the contents of the 2013Spring_Flood_Drawn.prj file? >>

Re: [R-sig-Geo] FW: proj4string error on data.frame generated from NCEP data

2016-04-08 Thread Edzer Pebesma
t;> >> library(RNCEP) >> air.sig995.extent <- NCEP.gather(variable='air.sig995', >> level='surface', months.minmax=c(7), >> years.minmax=c(2011), lat.southnorth=c(60,89.9), >> lon.westeast=c(-179,180), >> reanalysis2=FALSE, return.units=TRUE) >> air.sig995.a

Re: [R-sig-Geo] unionSpatialPolygons from the command-line

2016-04-08 Thread Edzer Pebesma
;> Invalid geometry, may only be applied to polygons >>> Calls: example ... eval -> eval -> unionSpatialPolygons -> >>> Execution halted >> >> Is this not possible? >> >> Thank you, >> >> Joseph Stachelek >> >> >&g

Re: [R-sig-Geo] Find a circle center with spatial points

2016-03-19 Thread Edzer Pebesma
take the mean >> coordinate it could fail horribly. See the paper I linked to for an >> example. They use the median to get an initial "robust" estimate of >> x,y,R, and then use some specialised optimisation to improve the >> estimate - you can't just throw it into "optim"! &

Re: [R-sig-Geo] split() function on SpatialPolygonsDataFrame increases file size

2016-03-22 Thread Edzer Pebesma
is that when I try to process this on multiple >> cores, I max out my memory. >> >> Much thanks, >> Luke >> >> ___ >> R-sig-Geo mailing list >> R-sig-Geo@r-project.org >> https://stat.ethz.ch/mailman/

Re: [R-sig-Geo] Tutorials

2016-03-21 Thread Edzer Pebesma
a. I added an "Events" tab to the http://r-spatial.org/ blog, and added these two events. Best regards, -- Edzer Pebesma Institute for Geoinformatics (ifgi), University of Münster Heisenbergstraße 2, 48149 Münster, Germany; +49 251 83 33081 Journal of Statistical Software: http:

Re: [R-sig-Geo] best practice for reading large shapefiles?

2016-04-27 Thread Edzer Pebesma
for 3 levels of granularity. I went with the smallest. I just > realized that the mid-size one would be sufficient for now, which > results in dbf=138mb and shp=501mb. Attempting to read this in now (~ > 30 minutes), which I assume will read in fine after some time. Will > respond to

Re: [R-sig-Geo] spatialpoints: each dot represents 100 individuals?

2016-04-25 Thread Edzer Pebesma
ps://stat.ethz.ch/mailman/listinfo/r-sig-geo > > _______ > R-sig-Geo mailing list > R-sig-Geo@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-sig-geo > -- Edzer Pebesma Institute for Geoinformatics (ifgi), University of Münster H

Re: [R-sig-Geo] spsample bug?

2016-05-11 Thread Edzer Pebesma
dry.grid.u,n=n,type="random") > } > > [[alternative HTML version deleted]] > > ___ > R-sig-Geo mailing list > R-sig-Geo@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-sig-geo > -- Edzer Pebesma Inst

Re: [R-sig-Geo] getting mean values for each raster in a stack

2016-05-16 Thread Edzer Pebesma
t only the overall mean with >> lapply(vals, FUN=mean) >> >> I wonder if there is a way to get the mean (and other parameters) for each >> layer using lapply function. >> >> Thanks in advance. >> >> Antonio Olinto >> >> > >

Re: [R-sig-Geo] Creating a grid

2016-05-18 Thread Edzer Pebesma
n do two things: 1. define a meaningful cell size in degrees, e.g. cellsize = 0.1 2. project the polygon using spTransform to a reference system where the units of coordinates are in metres, and then use cellsize = 1000 > > > > #creating the grid > > gridded(Pointlocationcoun

Re: [R-sig-Geo] Problems with rbind(list(), makeUniqueIDs=T)

2016-07-22 Thread Edzer Pebesma
sh_United States.1252 LC_CTYPE=English_United > States.1252 > [3] LC_MONETARY=English_United States.1252 LC_NUMERIC=C > [5] LC_TIME=English_United States.1252 > > attached base packages: > [1] stats graphics grDevices utils datasets methods base > > other atta

Re: [R-sig-Geo] storing data frame as SpatialPixelsDataFrame with coordinate intervals not constant

2016-07-18 Thread Edzer Pebesma
your case, they are not nearly on a grid. Hth, > > Cheers,Yaara > > > > > [[alternative HTML version deleted]] > > _______ > R-sig-Geo mailing list > R-sig-Geo@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-sig-geo > -- Edzer P

[R-sig-Geo] Simple features for R, part 2

2016-07-18 Thread Edzer Pebesma
There's a second blog post on the ISC project "Simple features for R", describing progress up to date and planned future steps, at: http://r-spatial.org/r/2016/07/18/sf2.html Comments are welcome! -- Edzer Pebesma Institute for Geoinformatics (ifgi), University of Münster Heisenbe

Re: [R-sig-Geo] Problems with installing rgdal on Ubuntu "checking gdal: linking with --libs only... no"

2016-08-04 Thread Edzer Pebesma
configure make sudo make install you can also do a sudo make clean which should remove everything (including e.g. the header files now still in /usr/local/include) -- Edzer Pebesma Institute for Geoinformatics (ifgi), University of Münster Heisenbergstraße 2, 48149 Münster, Germany; +49 251

Re: [R-sig-Geo] Problems with installing rgdal on Ubuntu "checking gdal: linking with --libs only... no"

2016-08-04 Thread Edzer Pebesma
>> `kmlengine::Href::Parse(std::__cxx11::basic_string<char, std:: >>> char_traits, std::allocator > const&)' >>> /usr/local/lib/libgdal.so: undefined reference to >>> `kmlengine::KmzFile::ReadKmlAndGetPath(std::__cxx11::basic_string<char, >>> std:: char_traits, st

Re: [R-sig-Geo] Problems with rbind(list(), makeUniqueIDs=T)

2016-07-23 Thread Edzer Pebesma
length(m[[1]]))) > > spChFIDs(m[[2]]) <- paste0("B", seq(length(m[[2]]))) > > spChFIDs(m[[3]]) <- paste0("C", seq(length(m[[3]]))) > > m <- do.call(rbind, m) > > class(m) > [1] "SpatialPolygonsDataFrame" > attr(,"package&

Re: [R-sig-Geo] analysis of point patterns

2016-06-29 Thread Edzer Pebesma
R. > > Best regards, > Lars-Daniel Weber > > ___ > R-sig-Geo mailing list > R-sig-Geo@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-sig-geo > -- Edzer Pebesma Institute for Geoinformatics (ifgi), University of Münster Heisenbergstraße 2, 481

Re: [R-sig-Geo] "invalid" geometries in shape data

2017-02-08 Thread Edzer Pebesma
ble with code but links to existing emails and online posts are > welcome too. Feel free to construct examples from scratch that reproduce a > particular warning/error. > > If you like, you can use the Github Issues mechanism here, or just email > them in reply. > > http

Re: [R-sig-Geo] Subset a SpLineDF

2017-02-22 Thread Edzer Pebesma
> "TRONCON_TOPO_DSAV_coord_unique_test", : > number of objects mismatch > > > It seems that I have not the same number of objects in the different > shapefile slots. How to update these objets after my subset ? > > > Cheers ! > > " > -- Edze

Re: [R-sig-Geo] How to draw a legend when plotting sf (simple feature) objects

2017-02-26 Thread Edzer Pebesma
(aes(fill = SID79)) Trying this with the gpkg you read in above fails for me; it looks like ggplot now wrongly assumes that geometry columns are always called `geometry' (despite the docs), which is not the case. -- Edzer Pebesma Institute for Geoinformatics (ifgi), University of Münster Heisen

[R-sig-Geo] EGU: R’s deliberate role in Earth sciences

2017-03-01 Thread Edzer Pebesma
This year's EGU meeting has a session with 16 presentation on "R’s deliberate role in Earth sciences"; for the program see: http://meetingorganizer.copernicus.org/EGU2017/pico/24971 -- Edzer Pebesma Institute for Geoinformatics (ifgi), University of Münster Heisenbergstraße 2, 481

Re: [R-sig-Geo] write.tgrass (spacetime) does not work for any RasterStack

2016-09-06 Thread Edzer Pebesma
Indeed: x = setZ(tempmean_new_rs, getZ(test)) write.tgrass(x, "x.tar.gz") seems to work. On 06/09/16 21:31, Veronica Andreo wrote: > Hello Edzer > > Thank you very much for your answer :) > > Yes, sure! Find it attached. > > Best, > Vero > >

Re: [R-sig-Geo] rpostgis announcement

2016-09-02 Thread Edzer Pebesma
ieve I got it, now, and both are supported (but not auto-detected). > > Best, > Mathieu. > > > Le 29/08/2016 16:19, Edzer Pebesma a écrit : >> Nice development! >> >> In the simple features for R project, https://github.com/edzer/sfr , I >> have been wri

Re: [R-sig-Geo] rpostgis announcement

2016-08-29 Thread Edzer Pebesma
gt; ___ > R-sig-Geo mailing list > R-sig-Geo@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-sig-geo > -- Edzer Pebesma Institute for Geoinformatics (ifgi), University of Münster Heisenbergstraße 2, 48149 Münster, Germany; +49 251 83

Re: [R-sig-Geo] A Qgis Map Graphics Device for R

2016-08-18 Thread Edzer Pebesma
ight pop up a map canvas with some scottish regions. Or crash R, or > give a warning about things not being found on your system... > > Successes and failure reports (with error messages, system > configuration, current phase of the moon etc) to the gitlab issue > tracker please! > >

Re: [R-sig-Geo] spatiotemporal co-kriging?

2016-11-10 Thread Edzer Pebesma
, > > Santiago Beguería > CSIC > Spain > > ___ > R-sig-Geo mailing list > R-sig-Geo@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-sig-geo > -- Edzer Pebesma Institute for Geoinformatics (ifgi), University of

Re: [R-sig-Geo] Write sf object to spatialite database

2016-11-10 Thread Edzer Pebesma
e > con <- dbConnect(RSQLite::SQLite(), '~/sandbox/db.sqlite') > > st_write_db(con, sf, "meuse_tbl", dropTable = FALSE) > > _______ > R-sig-Geo mailing list > R-sig-Geo@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-sig

Re: [R-sig-Geo] how to fit 3d variogram in gstat

2016-11-28 Thread Edzer Pebesma
int(v) > > > Thanks very much! > > > > Maogui Hu, Ph.D. > State Key Laboratory of Resources and Environmental Information System > Institute of Geographic Sciences and Natural Resources Research, CAS > > [[alternative HTML version deleted]] > > ___

Re: [R-sig-Geo] Error with gstat::predict [SEC=UNCLASSIFIED]

2016-10-31 Thread Edzer Pebesma
this email, you > acknowledge and accept these risks. > --------- > > > [[alternative HTML version deleted]] > > ___ >

Re: [R-sig-Geo] Error with gstat::predict [SEC=UNCLASSIFIED]

2016-10-31 Thread Edzer Pebesma
On 31/10/16 09:12, Roger Bivand wrote: > On Mon, 31 Oct 2016, Edzer Pebesma wrote: > >> >> >> On 31/10/16 05:09, Li Jin wrote: >>> Hi All, >>> >>> I need to use the predict{gstat} function in one of my functions for >>> a R package. I

[R-sig-Geo] Simple features now on CRAN

2016-11-04 Thread Edzer Pebesma
s, or want to participate in shaping this, please participate; contributions on this mailing list are also welcome. Kudos go to Jeroen Ooms for the Windows port, Roger Bivand for migrating most of the CRAN build farm to GDAL 2.0, and of course the CRAN team itself! -- Edzer Pebesma Institute for Geoinformat

Re: [R-sig-Geo] Add layers to a lattice plot via a for loop

2016-11-02 Thread Edzer Pebesma
loop failed to overlay layers? >> >> Thanks, >> David >> >>[[alternative HTML version deleted]] >> >> _______ >> R-sig-Geo mailing list >> R-sig-Geo@r-project.org >> https://stat.ethz.ch/mailman/listinfo/r

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

2016-12-07 Thread Edzer Pebesma
(or in general: spatial) geometries. -- Edzer Pebesma Institute for Geoinformatics (ifgi), University of Münster Heisenbergstraße 2, 48149 Münster, Germany; +49 251 83 33081 Journal of Statistical Software: http://www.jstatsoft.org/ Computers & Geosciences: http://elsevier.com/locate/cageo/

Re: [R-sig-Geo] Specifying Schema for PostGIS Layer with Simple Features

2016-12-09 Thread Edzer Pebesma
ll step, these days. > > [[alternative HTML version deleted]] > > ___ > R-sig-Geo mailing list > R-sig-Geo@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-sig-geo > -- Edzer Pebesma Institute for Geoinformat

Re: [R-sig-Geo] Simulating or bootstrapping of Spatial Data

2017-01-09 Thread Edzer Pebesma
ted]] > > ___ > R-sig-Geo mailing list > R-sig-Geo@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-sig-geo > -- Edzer Pebesma Institute for Geoinformatics (ifgi), University of Münster Heisenbergstraße 2, 48149 Münster, Germa

Re: [R-sig-Geo] sp::spDists*

2017-01-13 Thread Edzer Pebesma
A new version of gstat, fixing the same bug, is now also on CRAN. On 09/01/17 11:39, Roger Bivand wrote: > New versions of spdep and spgwr fixing the same bug are now on CRAN. > > Roger > > On Thu, 8 Dec 2016, Edzer Pebesma wrote: > >> While building support for geos

Re: [R-sig-Geo] PostGIS layer to R using OGR driver imports unwanted raster column as String

2016-12-01 Thread Edzer Pebesma
e.g. > layer = "SCHEMA.TABLE") > > Some more ideas? Thank you. > > > > > > -- Původní zpráva -- > Od: Roger Bivand > Komu: Roger Bivand , Vilem Ded , Help > R-Sig_Geo > Datum: 30. 11. 2016 12:53:32 > Předmět: Re: [R-sig-Geo

Re: [R-sig-Geo] simple features : st_segmentize question

2016-12-20 Thread Edzer Pebesma
_different_ number of segments for each > line ( and use the value in column 'n') > I tried: > st_segmentize(mls,1/mls$n) > but this throws an error: > "Error in CPL_gdal_segmentize(st_geometry(x), dfMaxLength) : >expecting a single value" > > Is there any way to achi

Re: [R-sig-Geo] Does over() in sp work with "fun=sum"?

2016-12-20 Thread Edzer Pebesma
t; > ___ > R-sig-Geo mailing list > R-sig-Geo@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-sig-geo > -- Edzer Pebesma Institute for Geoinformatics (ifgi), University of Münster Heisenbergstraße 2, 48149 Münster, Germany;

Re: [R-sig-Geo] sp::spDists*

2016-12-28 Thread Edzer Pebesma
The problem described below has been corrected ins sp 1.2-4 which is now on CRAN; new releases of gstat, spgwr and spdep will follow to correct for the same bug. On 08/12/16 15:48, Edzer Pebesma wrote: > While building support for geosphere distance functions in sf and > comparing result

Re: [R-sig-Geo] Scaling coarse raster to finer raster using overlapping area weights

2016-12-21 Thread Edzer Pebesma
ingen Economic Research > > |Schlossplatz 1 - A-2361 Laxenburg, Austria, room s-138 > S +43 2236 807 537 > Skype: michiel.van.dijk > S https://www.researchgate.net/profile/Michiel_Van_Dijk/ > / http://nl.linkedin.com/pub/michiel-van-dijk/1a/805/346 > > www.iiasa.ac.at > www.w

Re: [R-sig-Geo] Building a vector of random sf st_points

2016-12-21 Thread Edzer Pebesma
NA proj4string:NA POINT(1 4) POINT(2 5) POINT(3 6) > ptsdf <- apply(PTS, 1, st_point) > > would work but it doesn't. Why? apply wants to to mold the output into a matrix; it can, and so it does. We want a list of POINT objects, so need to use lapply. -- Edzer Pebesma Inst

Re: [R-sig-Geo] Maintain SRID with st_write to postgis db

2017-03-16 Thread Edzer Pebesma
. > I don't know > whether this will affect swapping in and out of sf <->sp, but include it > for completeness. You can safely ignore it. If `sf' would import `sp' it would go away, but there is no need to do so. Not importing packages makes it easier to understand where t

Re: [R-sig-Geo] Maintain SRID with st_write to postgis db

2017-03-16 Thread Edzer Pebesma
with -lgeos_c... yes > checking for lwgeom_version in -llwgeom... no > configure: error: in > `/tmp/RtmpKLmhyt/devtools9fd11eb9c3/edzer-sfr-d620f3b': > configure: error: lwgeom test failed (--without-readline to disable) > See `config.log' for more details > ERROR:

Re: [R-sig-Geo] xlsx or rJava issue

2017-03-17 Thread Edzer Pebesma
gt; > is it a particular JAVA_HOME setting or a PATH issue? > > I added this in the PATH, re-installed rJava, re-installed xlsx ... nope! > thanks, > > Didier > > on R version 3.3.2 (2016-10-31) > -- Edzer Pebesma Institute for Geoinformatics (ifgi), University of

Re: [R-sig-Geo] Maintain SRID with st_write to postgis db

2017-03-16 Thread Edzer Pebesma
r attached packages: > [1] sf_0.3-4 > > loaded via a namespace (and not attached): > [1] DBI_0.6 tools_3.3.1 units_0.4-2 Rcpp_0.12.9 > udunits2_0.13 grid_3.3.1 > > > (I realize st_write_db is an option, but not finding myself able to > specify a schema (it just

Re: [R-sig-Geo] Ordinary co-kriging with gstat - computing time/syntax error?

2017-04-10 Thread Edzer Pebesma
ing in the predict command… I do > not need to include the training data in it? > I have been reading previous posts from R-sig-Geo, but I still did not find > the solution. I would appreciate any advice you could give me. > > Thanks, > > Mercedes Roman > > INRA >

Re: [R-sig-Geo] Help for STFDF Object Creation

2017-04-23 Thread Edzer Pebesma
GS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0" > ..@ time :An ‘xts’ object on 2017-01-02 01:00:00/2017-03-22 containing: > Data: int [1:1896, 1] 1 2 3 4 5 6 7 8 9 10 ... > - attr(*, "dimnames")=List of 2 > ..$ : NULL > ..$ : chr "timeIndex&qu

Re: [R-sig-Geo] Error in bbox(newdata) : object not a >= 2-column array

2017-03-14 Thread Edzer Pebesma
; https://stat.ethz.ch/mailman/listinfo/r-sig-geo >>> >> >> >> >> -- >> David Warner >> Research Fisheries Biologist >> U.S.G.S. Great Lakes Science Center >> 1451 Green Road >> Ann Arbor, MI 48105 >> 734-214-9392 >> >>

Re: [R-sig-Geo] Create Nodes at lines' segments extremities.

2017-03-03 Thread Edzer Pebesma
ts ? > > > Does somebody has already done /wrote a similar script to process that ?? > > > Thank you in advance for the helpfull advices. > > Cheers. > > Tristan Bourgeois > Hydromorphologist Engineer > Seine Normandie Water Agency > Rouen, Fra

Re: [R-sig-Geo] How many polygons does a single polygon intersect with?

2017-03-08 Thread Edzer Pebesma
quot;C:\\Users\\rameshv\\Terrestrial_dissolved_multipart.shp") > Poly_join <- gIntersects(Poly) > > Regards, > Vijay > -- Edzer Pebesma Institute for Geoinformatics (ifgi), University of Münster Heisenbergstraße 2, 48149 Münster, Germany; +49 251 83 33081 Journal of Statistical

Re: [R-sig-Geo] how to get weights for inverse distance weighting

2017-04-07 Thread Edzer Pebesma
will > in general be 20, but in my real problem it will be less than the number of > times. > > Thank you, > Scott Waichler > Pacific Northwest National Laboratory > Richland, WAUSA > > _______ > R-sig-Geo mailin

Re: [R-sig-Geo] Warning in st_intersection

2017-07-06 Thread Edzer Pebesma
r standing for "attribute geometry relationship"). More here: https://cran.r-project.org/web/packages/sf/vignettes/sf1.html#how-attributes-relate-to-geometries > > Thanks > -- Edzer Pebesma Institute for Geoinformatics (ifgi), University of Münster Heisenbergstr

Re: [R-sig-Geo] How does sf do rgeos::gUnaryUnion or maptools::unionSpatialPolygons?

2017-05-09 Thread Edzer Pebesma
h 2 MULTIPOLYGON and 98 POLYGON. You could st_cast that to MULTIPOLYGON. Alternatively: library(dplyr) bind_cols(nc90, ids=ids) %>% group_by(ids) %>% summarise(ST=head(ST,1), CO=head(CO,1), do_union=FALSE) converts straight into MULTIPOLYGON. I'll update the sp -> sf migr

Re: [R-sig-Geo] units for spatial range in stvariogram gstat

2017-05-17 Thread Edzer Pebesma
nce > Cairo University > > > Tel:(202)35728055-(202)35728116-(202)35736608-(202)35736605 > Fax:(202)35711020 > Follow us on twitter:https://twitter.com/fepsnews > > > ___ > R-sig-Geo mailing list > R-sig-Geo@r-

Re: [R-sig-Geo] Nested Variogram Model for 3D data using gstat

2017-05-18 Thread Edzer Pebesma
ogy and Resources Management > > in the Tropics and Subtropics > > *Technology* *Arts** **Sciences* > *TH Köln* - University of Applied Sciences > > > > T: +49 221 8275-2059 > > F: +49 221 8275-2736 > E : abm.fi...@fh-koeln.de <mailto:abm.fi...@fh-

Re: [R-sig-Geo] Is a simple feature -friendly version of spdep being developed?

2017-05-17 Thread Edzer Pebesma
reference for dataframes, S3 >>>>>>> objects, list columns, and pipeline workflows would be well-suited to >>>> the >>>>>>> field of spatial analysis. Are there some fundamental reasons why the >>>>>>> `spdep`

Re: [R-sig-Geo] Issues installing rgeos on linux server

2017-06-21 Thread Edzer Pebesma
checking for gcc... /usr/bin/gcc -std=gnu99 >>> checking whether the C compiler works... yes >>> checking for C compiler default output file name... a.out >>> checking for suffix of executables... >>> checking whether we are cross compiling... configure: err

Re: [R-sig-Geo] Kriging prediction variances on gstat

2017-06-21 Thread Edzer Pebesma
ure that the sill value (partial sill + nugget) > is used. > > Thank you > Carlo > [[alternative HTML version deleted]] > > ___ > R-sig-Geo mailing list > R-sig-Geo@r-project.org > https://stat.ethz.ch/mailman/listinfo/

Re: [R-sig-Geo] Using gIntersection to split polygons with lines

2017-06-26 Thread Edzer Pebesma
gt; Thanks, > Frederico > ᐧ > > [[alternative HTML version deleted]] > > ___ > R-sig-Geo mailing list > R-sig-Geo@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-sig-geo > -- Edzer Pebesma Institute f

[R-sig-Geo] r-spatial meetup @ dinner on Jul 5, Brussels

2017-06-19 Thread Edzer Pebesma
.com/spreadsheets/d/10YikeTkQzo708yyrNWDxwBPS56kEksPe94obSAeFxIY/edit?usp=sharing Many regards, -- Edzer Pebesma Institute for Geoinformatics (ifgi), University of Münster Heisenbergstraße 2, 48149 Münster, Germany; +49 251 83 33081 Journal of Statistical Software: http://www.jstatsoft

Re: [R-sig-Geo] st_area / st_transform / crs

2017-09-16 Thread Edzer Pebesma
Thank you advance for your help, Thanks for bringing this up - I understand it is confusing, and we might be able to make this simpler from the sf side. > > Jean-Luc Dupouey > -- Edzer Pebesma Institute for Geoinformatics (ifgi), University of Münster Heisenbergstraße 2, 4

Re: [R-sig-Geo] Compile sf in windows with liblwgeom

2017-09-16 Thread Edzer Pebesma
nverting to sp and back, like st_as_sf(as(sp1, "Spatial")) might do the trick, though. > > Thanks > > [[alternative HTML version deleted]] > > ___ > R-sig-Geo mailing list > R-sig-Geo@r-project.org > htt

Re: [R-sig-Geo] Compile sf in windows with liblwgeom

2017-09-16 Thread Edzer Pebesma
On 16/09/17 18:22, Barry Rowlingson wrote: > > > On Sat, Sep 16, 2017 at 3:45 PM, Edzer Pebesma > <edzer.pebe...@uni-muenster.de <mailto:edzer.pebe...@uni-muenster.de>> > wrote: > > > > On 16/09/17 15:20, Eduardo Diez wrote: > > Dear lis

Re: [R-sig-Geo] Compile sf in windows with liblwgeom

2017-09-17 Thread Edzer Pebesma
On 17/09/17 00:07, Edzer Pebesma wrote: > > > On 16/09/17 18:22, Barry Rowlingson wrote: >> >> >> On Sat, Sep 16, 2017 at 3:45 PM, Edzer Pebesma >> <edzer.pebe...@uni-muenster.de <mailto:edzer.pebe...@uni-muenster.de>> >> wrote: &g

Re: [R-sig-Geo] Errors at installing rgdal on Debian Buster

2017-09-12 Thread Edzer Pebesma
org https://stat.ethz.ch/mailman/listinfo/r-sig-geo >> >> >> [[alternative HTML version deleted]] >> >> ___ >> R-sig-Geo mailing list >> R-sig-Geo@r-project.org >> https://stat.ethz.ch/mailman/

Re: [R-sig-Geo] Compatibility between sf and raster

2017-08-29 Thread Edzer Pebesma
es, fdef, mtable) : > unable to find an inherited method for function ‘raster’ for signature > ‘"sf"’ > > > -- Edzer Pebesma Institute for Geoinformatics (ifgi), University of Münster Heisenbergstraße 2, 48149 Münster, Germany; +49 251 83 33081 Journal of Statistical Software:

Re: [R-sig-Geo] Why is the covariance in Universal Kriging modeled this way in lectures by Prof. Edzer Pebesma?

2017-11-22 Thread Edzer Pebesma
On 11/22/2017 12:09 AM, Joelle k. Akram wrote: > Greetings, > > There is code for Universal Kriging from Prof. Edzer Pebesma in > GitHub<https://github.com/edzer/mstp/blob/master/lec7.Rmd>. https://edzer.github.io/mstp/lec7.html gives you the rendered version. > >

Re: [R-sig-Geo] Why is the covariance in Universal Kriging modeled this way in lectures by Prof. Edzer Pebesma?

2017-11-22 Thread Edzer Pebesma
be preferred. > > > thank you, > > Chris Akram > > > > > -------- > *From:* R-sig-Geo <r-sig-geo-boun...@r-project.org> on behalf of Edzer > Pebesma <edzer.pebe...@uni-muenster.de> > *Sent:* November 22, 2017 6:16 AM > *To:* r-sig-geo@r-project.org > *Subjec

Re: [R-sig-Geo] is the nugget parameter included in the "VariogramLine" function in gstat?

2017-12-04 Thread Edzer Pebesma
2),covariance=TRUE) [,1] [,2] [1,] 2.000 0.3678794 [2,] 0.3678794 2.000 > > > thank you. > > Bernard > > [[alternative HTML version deleted]] > > _______ > R-sig-Geo mailing list > R-sig-Geo@r-project.or

Re: [R-sig-Geo] +towgs84 in st_write

2017-12-12 Thread Edzer Pebesma
deral Department of Economic Affairs, Education and Research EAER > Agroscope > > Reckenholzstrasse 191, CH-8046 Z�rich > Ph. +41 58 468 75 98 > Fax +41 58 468 72 01 > manuel.schnei...@agroscope.admin.ch<mailto:manuel.schnei...@art.admin.ch> > https://www.agroscope.a

Re: [R-sig-Geo] +towgs84 in st_write

2017-12-12 Thread Edzer Pebesma
On 12/12/2017 11:27 PM, Roger Bivand wrote: > Does sf use morphToESRI()? No. -- Edzer Pebesma Institute for Geoinformatics Heisenbergstrasse 2, 48151 Muenster, Germany Phone: +49 251 8333081 ___ R-sig-Geo mailing list R-sig-Geo@r-project.org ht

Re: [R-sig-Geo] split lines at vertices with sf

2017-11-01 Thread Edzer Pebesma
estamp.com/email-install?utm_source=promotion_medium=signature_campaign=get_your_own> > -- Edzer Pebesma Institute for Geoinformatics Heisenbergstrasse 2, 48151 Muenster, Germany Phone: +49 251 8333081 ___ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo

<    1   2   3   4   5   6   7   8   >