Search Archives

An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20130727/50bf67bc/attachment.pl>

An embedded and charset-unspecified text was scrubbed... Name: not available Url: https://stat.ethz.ch/pipermail/r-help/attachments/20080126/45df19d2/attachment.pl

For intersecting two polygon layers, is there any advantage of using gIntersection (from the rgeos package) instead of intersect (from the raster package)? After a few tests, it appears that intersect is a little bit faster. Moreover, intersect automatically calculates...

Dear Roger Bivand, I read your message about the rgdal retirement. It was a nice package and an important step in the development of R-spatial. I looked for alternatives for the functions I was using, but I didn't...

Dear forumites, Using R, I want to regularly subsample a raster (for example, one point over two), and get a raster. Usual R subsampling works with raster layers: MyRaster[seq(1,nrow(MyRaster),2),seq(1,ncol(MyRaster),2))] but...

Dear all, is it possible to read Geoconcept files (.gcm and .GCR) in R? ogrDrivers()$name indicates that a Geoconcept driver exists, but I was not able to read any data. I get a "Cannot open data source" error: > ogrDrivers...

Dear all, I am a new user of sf. I have the following error: "lay" is an sf object, composed of POLYGONS and MULTIPOLYGONS in a projected coordinates system. > st_area(lay) # works well and returns the area of polygons...

Dear members, Not a technical question: The number of threads in this mailing list, following a long period of increase, has been regularly and strongly decreasing since 2010, passing from more than 40K threads to less than 11K threads last...

Dear R helpers, I try to use the 'delaunayn' function in the 'geometry' package for Delaunay triangulation in 2 dimensions. For the four following points, I get a warning message : > coord=matrix(ncol=2,byrow=TRUE,c(622633,7073452, + 621228...

Why gIntersection returns the ""TopologyException: no outgoing dirEdge found at" error in the following very simple case: > #construct a spatial layer with two adjacent squares, each of 10 x 10 units: > > Pol1=rbind(c(0,0),c(0,10),c...

I would like to redraw a plot from the .SavedPlots object (list), using R code (not arrows on the keyboard). > windows(record=TRUE) > > x=runif(100) > plot(x,col="blue") > hist(x,col="red") > plot(x,col="green") > #When I...

I faced the same question recently: how to tranform a multipart polygon layer into single polygons? Here is a little bit more general function to do that. It keeps the initial ID of the polygons, adding a sequence number. The...

I get an incorrect match between row names of coordinates and data in SpatialPointsDataFrame: > Mycoord <- data.frame(row.names=letters[1:5], x=rnorm(5), y=rnorm(5)) > > Mycoord x y a -0.80138451 1.9291882 b -1.58229979 -0...

Thanks for your answer. Unfortunately, I have no time to spend debugging the GEOS code. Just making some tests under R, I found an even worse situation, where the program seems to work, but gives a completely false result: library...

Dear forumites, I try to apply a gridded transformation of coordinates in |sf|, but it does not work. Is there something else to do apart from calling |sf_proj_network(TRUE)|? Here is the 7 lines of code I used...