Hfrhyu

What is this scratchy sound on the acoustic guitar called?

Why can't we say "I have been having a dog"?

What is the difference between 'contrib' and 'non-free' packages repositories?

Why was Sir Cadogan fired?

Send out email when Apex Queueable fails and test it

Can compressed videos be decoded back to their uncompresed original format?

Does Dispel Magic work on Tiny Hut?

How to show a landlord what we have in savings?

Avoiding the "not like other girls" trope?

Why didn't Boeing produce its own regional jet?

Different meanings of こわい

Can someone clarify Hamming's notion of important problems in relation to modern academia?

Did 'Cinema Songs' exist during Hiranyakshipu's time?

Getting extremely large arrows with tikzcd

What does the same-ish mean?

Why do I get negative height?

Unlock My Phone! February 2018

When handwriting 黄 (huáng; yellow) is it incorrect to have a disconnected 草 (cǎo; grass) radical on top?

Is there a rule of thumb for determining the amount one should accept for a settlement offer?

Machine learning testing data

Finding the reason behind the value of the integral.

Does int main() need a declaration on C++?

How could sorcerers who are able to produce/manipulate almost all forms of energy communicate over large distances?

Do creatures with a speed 0ft., fly 30ft. (hover) ever touch the ground?

Find borders between all unique country pairs with R

How to split polygons by overlapping feature?Find all census tracts N miles from a lat, lon pointExtract all unique values from a RGB rasterFind intersections between spatial lines and points using RUnique SpatialPolygons object with multiple polygons slots in RPostGIS: Extracting borders between contiguous (but not perfectly aligned) polygonsCountry border coordinates with GADM and RMapping Sahel with R and country names?Find overlay between polygons excluding boundaryCreating lines between pairs of pointsRecording all vertices on shortest path with igraph?

0

I would like to obtain a border as a line between all unique country pairs that are contiguous. Here is how to do it:

library(spData)
library(sf)
library(rnaturalearth)
library("rmapshaper")

#load countries dataset:
countries_sp <- ne_countries(scale = 10, continent= "Europe")

# find borders
border_sf <- st_as_sf(countries_sp) %>% 
 ms_innerlines() %>% 
 as_tibble() %>% 
 st_as_sf()

# create ids to visualise lines found by ms_innerlines
border_sf$ids <- rep(1:length(border_sf$geometry),1)

# plot
tm_shape(countries_sp, bbox = c(c(-8.750803, 36.14111), c(40.15954, 70.07531))) +
 tm_fill() +
 tm_shape(border_sf) +
 tm_lines(col = "ids")

Example shows how to find have one line for the border between France and Germany, one for France and Switzerland, France and Italy for each contiguous country pair in the sample.

r

share|improve this question

edited 2 days ago

S Front

asked 2 days ago

S FrontS Front

165

• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  I guess my code already provides the answer
  
  – S Front
  2 days ago
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  So why ask the question? Or at least edit your question and remove the answer then put the answer in an answer below.
  
  – Spacedman
  2 days ago
  

  
  
  
  

add a comment | 

0

I would like to obtain a border as a line between all unique country pairs that are contiguous. Here is how to do it:

library(spData)
library(sf)
library(rnaturalearth)
library("rmapshaper")

#load countries dataset:
countries_sp <- ne_countries(scale = 10, continent= "Europe")

# find borders
border_sf <- st_as_sf(countries_sp) %>% 
 ms_innerlines() %>% 
 as_tibble() %>% 
 st_as_sf()

# create ids to visualise lines found by ms_innerlines
border_sf$ids <- rep(1:length(border_sf$geometry),1)

# plot
tm_shape(countries_sp, bbox = c(c(-8.750803, 36.14111), c(40.15954, 70.07531))) +
 tm_fill() +
 tm_shape(border_sf) +
 tm_lines(col = "ids")

Example shows how to find have one line for the border between France and Germany, one for France and Switzerland, France and Italy for each contiguous country pair in the sample.

r

share|improve this question

edited 2 days ago

S Front

asked 2 days ago

S FrontS Front

165

• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  I guess my code already provides the answer
  
  – S Front
  2 days ago
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  So why ask the question? Or at least edit your question and remove the answer then put the answer in an answer below.
  
  – Spacedman
  2 days ago
  

  
  
  
  

add a comment | 

I would like to obtain a border as a line between all unique country pairs that are contiguous. Here is how to do it:

library(spData)
library(sf)
library(rnaturalearth)
library("rmapshaper")

#load countries dataset:
countries_sp <- ne_countries(scale = 10, continent= "Europe")

# find borders
border_sf <- st_as_sf(countries_sp) %>% 
 ms_innerlines() %>% 
 as_tibble() %>% 
 st_as_sf()

# create ids to visualise lines found by ms_innerlines
border_sf$ids <- rep(1:length(border_sf$geometry),1)

# plot
tm_shape(countries_sp, bbox = c(c(-8.750803, 36.14111), c(40.15954, 70.07531))) +
 tm_fill() +
 tm_shape(border_sf) +
 tm_lines(col = "ids")

Example shows how to find have one line for the border between France and Germany, one for France and Switzerland, France and Italy for each contiguous country pair in the sample.

r

share|improve this question

edited 2 days ago

S Front

asked 2 days ago

S FrontS Front

165

library(spData)
library(sf)
library(rnaturalearth)
library("rmapshaper")

#load countries dataset:
countries_sp <- ne_countries(scale = 10, continent= "Europe")

# find borders
border_sf <- st_as_sf(countries_sp) %>% 
 ms_innerlines() %>% 
 as_tibble() %>% 
 st_as_sf()

# create ids to visualise lines found by ms_innerlines
border_sf$ids <- rep(1:length(border_sf$geometry),1)

# plot
tm_shape(countries_sp, bbox = c(c(-8.750803, 36.14111), c(40.15954, 70.07531))) +
 tm_fill() +
 tm_shape(border_sf) +
 tm_lines(col = "ids")

share|improve this question

edited 2 days ago

S Front

asked 2 days ago

S FrontS Front

165

share|improve this question

edited 2 days ago

S Front

asked 2 days ago

S FrontS Front

165

asked 2 days ago

S FrontS Front

165

asked 2 days ago

S FrontS Front

165