How to properly plot sf object

Question

I have data like the following:

These data are obtained from a Stata dataset (.dta exported to .csv) and are intended to plot a map of shapefiles. I import it into R, drop rows with missing data, and convert to sf (I am not sure if this is the best way to do it). I then plot.

coords <- read.csv("filename.csv") %>% drop_na()
coords2 <- st_as_sf(x=coords, coords = c("X_X", "X_Y"))
ggplot(coords2) + geom_sf()

This is what I get:

How can I instead get a map that has lines instead of the points? And how can I shade the individual map areas based on another variable, such as X_ID?

Alvaro Morales · Accepted Answer · 2022-01-24T15:53:23.057

1

You can use ggplot2::geom_path(). No need to convert into sf object, just use the coordinates.

edited Jan 24 '22 at 15:53

answered Jan 23 '22 at 15:35

Alvaro Morales

1,845
3
12
21

score 0 · Answer 2 · answered Jan 23 '22 at 01:57

0

You might need to use points_to_lines from the sp package.

library(sp)
map <- points_to_line(data = coords, 
                          long = "X_X", 
                          lat = "Y_Y", 
                          id_field = "X_id", 
                          sort_field = "shape_pt_sequence")

answered Jan 23 '22 at 01:57

Pep

33
6

Thanks. Which package is points_to_line in? – bill999 Jan 23 '22 at 04:14
You can find it in the sp package... – Pep Jan 23 '22 at 04:30

How to properly plot sf object

2 Answers2