Questions tagged [sf]

us this tag when asking questions about sf R package (sf: Simple Features for R)

Filter by
Sorted by
Tagged with
4
votes
1answer
366 views

fit two sf polygons seamlessly

Copied from my SO question The problem Suppose we have two shapefiles that should border seamlessly. Only, they don't. Is there a way to force them to stick to one another without gaps? The ...
8
votes
2answers
1k views

Selecting features “above” or “below” a line using R

Given a line and a set of points, I can't figure out how to use sf to identify which side of the line each point falls. A small reproducible example follows, adapted from a different question # Load ...
8
votes
2answers
3k views

Breaking long line segments into shorter ones in R using sf

I have a shapefile of river segments that I've imported into R. The segments range in length from 5-115km. I want to split the segments that are over 10km into 10km chunks. So I use st_segmentize. But ...
0
votes
2answers
806 views

How to apply log scale when plotting sf data?

I want to plot some sf data, similar to this example: library (sf) fishnet <- st_make_grid () n <- length (fishnet) vals <- rlnorm (n) fishnet <- st_sf (fishnet) fishnet <- cbind (...
6
votes
1answer
2k views

Applying filter with the sf package not working

Question: I have a sf-object in R and want to select a few polygons according to their location name. I would like to apply the dpylr approach using the filter function (analog to this example). ...
1
vote
0answers
247 views

sf::st_read_db() considerably slower than postGIStools::get_postgis_query()

I have a large (extent not size) polygon layer of the German federal state borders on my PostgreSQL data base. There is also an index on the geom column. CREATE INDEX germany_bld_geom_idx ON ...
0
votes
2answers
901 views

Get length of components of a multilinestring

I am looking to get not the length of the overall multilinestring but rather the length of the component parts. Here is what I have so far: library(sf) s1 <- rbind(c(0,3),c(0,4),c(1,5),c(2,5)) s2 ...
1
vote
1answer
744 views

lines of a SpatialLinesDataFrame corresponding to points

I have a spatial lines data frame which can be projected on to a map. The projection is defined as WGS84. I have a Spatial Points Data Frame which is also in the same location. How can I isolate the ...
0
votes
2answers
300 views

How to get a sf with line segments from a data table?

Given the following code: library(data.table) x<-runif(10) y<-runif(10) dt<-data.table( someAtt=sample(c("a","b","c"),10,replace = T), long=x, lat=y, prevLong=shift(x), prevLat=shift(y) ...
0
votes
1answer
83 views

Coercing to sp and dealing with duplicate field names

I am dealing with a very similar situation as described at: rgdal::writeOGR collision in column names results in NULL data I have an sf object with duplicate field names (in a case-insensitive way). ...
2
votes
1answer
554 views

Is there a way to create dot density maps with R package sf?

That is, is there a way to place dots at random positions inside a polygon with sf? It's easy enough to do with maptools, using dotsinPolys, but since I like sf's "tidyness", I'd love to be able to ...
2
votes
1answer
1k views

sf lines to polygons with holes

With sp package I generate automatically boundaries for different sets of spatial points. I'm trying to perform the same workflow with sf package but i had problems for detecting holes in the way sp ...
1
vote
0answers
423 views

Using st_read_db in R to read postGIS rasters

I use the st_read_db function in the sf package to read vector data from PostGIS into R. Can i use it to query rasters? ideally extract a raster value at a point: require(sf) # create connection to ...
1
vote
1answer
1k views

How to read a csv dataset with multipoint column in sf

i have a csv with a multipoint WGS84 column: head(df, n=3) CUP the_geom <chr> <chr> 1 ...
15
votes
3answers
6k views

How to combine sfc objects from R package sf

Using the R package sf, how does one combine sfc objects? For example, given the following code, how would one create a single sfc object sfc12 that includes the geometries from both sfc1 and sfc2? (...
8
votes
3answers
2k views

Return end-vertices of an sf linestring

I am trying to return the vertices of a linestring using the sf package. So far I have: ls <- st_sfc(st_linestring(rbind(c(0,0),c(0,1))), st_linestring(rbind(c(0,0),c(10,0)))) ...
8
votes
1answer
3k views

Calculating distances between two geometry columns using R?

I'd like to calculate the pairwise distance between two geometry columns. That is, I want to calculate the distance between several pairs of geometries without calculating the full distance matrix ...
2
votes
1answer
342 views

CRS for US Census 1999 ShapeFile

I'm working with the shapefiles for the US census 1999 msa files which can be downloaded here: https://www2.census.gov/geo/tiger/PREVGENZ/ma/ma99/cm_sa_99_shp.zip Unfortunately, there is no CRS ...
3
votes
1answer
189 views

edit sf point conditionally

I have an sf point dataframe, with a single incorrect point. I know what it's coordinates should be, but I'm not sure how best to directly edit them. I'm aiming for a conditional edit like the ...
7
votes
2answers
7k views

How to calculate distance from POINT to LINESTRING in R using `sf` library and get all POINT features from LINESTRING where distances were calculated?

I'm trying to calculate distances from a set of POINT features (1, 2, 3, 4, 5, 6 & 7) to a LINESTRING feature using sf library in R. I could calculate all the distances using the function ...
18
votes
1answer
10k views

Simplify polygons of sf object

How do I simplify an sf polygon without introducing gaps and slivers? With a shapefile, for example, I would use rmapshaper::ms_simplify(): library("pryr") library("rgdal") library("rmapshaper") ...
1
vote
3answers
2k views

Is there an equivalent of dplyr::bind_rows (or the old plyr::rbind.fill) for spatial objects?

bind_rows can append dataframes even when they do not have identical columns. It fills any non-matching columns with NAs. Is there an equivalent for spatial point data frames in R? I have two ...
2
votes
2answers
2k views

Using simple features (sf) in R, how do I erase polygons overlapping with another layer

With a regular spatialDataFrame I could use erase(x, y) or x-y. How can I perform the equivalent using sf objects? As an example: ##Create layers based on this answer: http://stackoverflow.com/...
34
votes
2answers
26k views

R convert sf object back to SpatialPolygonsDataFrame

How do I convert an sf object back to a SpatialPolygonsDataFrame? I am running R 3.4.0 and sf version 0.4.1. The sf documentation mentions an as_Spatial function, but this isn't available in my ...
1
vote
0answers
618 views

How to execute st_read function including multibyte character?

how to execute st_read function including multibyte character? library(sf) nc <- st_read(system.file("shape/nc.shp", package="sf")) st_write(obj = nc, dsn = enc2utf8("データ"),layer = "nc", driver = ...
1
vote
1answer
260 views

Is it possible to read tables from a specific schema with sf::st_read_db()

How can I specify the schema from where to read a table in the st_read_db() function in the sf package in R? When I try require(RPostgreSQL) require(sf) drv = dbDriver("PostgreSQL") con = dbConnect(...
7
votes
3answers
3k views

Split a line when it crosses a polygon in R

This should be a simple problem but I can't seem to find a solution. I have a line which crosses several polygons and I want to split it into multiple lines at the intersection of each polygon using ...
3
votes
1answer
254 views

Can't find function st_join in R sf package

I am trying to run the st_join example from the sf user manual https://cran.r-project.org/web/packages/sf/sf.pdf (p. 43). a = st_sf(a = 1:3, geom = st_sfc(st_point(c(1,1)), st_point(c(2,2)), st_point(...
22
votes
4answers
13k views

Crop simple features object in R

Is there a function to crop sf map object, similar to maptools::pruneMap(lines, xlim= c(4, 10), ylim= c(10, 15)) used for SpatialPolygon or SpatialLine? I am considering st_intersection() but there ...
5
votes
1answer
1k views

Creating an sf object from the maps package

I am trying to use data from the maps package in an sf workflow. I'm not sure if I have the operation correct (see below). It seems to work until state 20 (massachusetts). What am I doing wrong? Maybe ...
10
votes
1answer
10k views

Create a circle of defined radius around a point and then find the overlapping area on a shapefile

My use case is related to calculating the area of a polygon in a shapefile that is covered by another polygon. Specifically, I want to do the following: Create circles that have a 100m radius, from ...
2
votes
1answer
773 views

left_join breaking sf

Trying to add a data column to a simplefeatures point object using dplyr::left_join, but it appears to drop the sf class: class(sf) [1] "sf" "data.frame" class(df) [1] "tbl_df" "tbl" "...
57
votes
4answers
26k views

Extracting data.frame from simple features object in R

Is there an sf-native (i.e. "correct") way of extracting everything except the geometry column from a simple features object? This works df <- dplyr::select(as.data.frame(sf), -geometry) but the ...
4
votes
1answer
593 views

Iterate through shapefile and output an image for each feature using sf in R

I would like to create some statistical profiles for a range of features; specifically, 2016 Census profiles for the Census Divisions (CD) within British Columbia (BC), Canada. This would require a ...
12
votes
1answer
3k views

Writing multiple layers to GeoPackage using writeOGR() in R?

I'm trying to write multiple layers to the same GeoPackage in R, but I get an error, Creation of output file failed. I have tried to search for documentation of reading and writing to .gpkg files with ...
28
votes
1answer
17k views

Lon-Lat to Simple Features (sfg and sfc) in R

How can I convert lon-lat points to simple features (sfg), and then put them in a simple feature collection (sfc)? Here is a MWE that doesn't work but is the closest I have gotten to. library(data....

1
4 5 6 7
8