code_01.R

library(sf)

## Linking to GEOS 3.8.1, GDAL 3.1.3, PROJ 7.1.1

library(rnaturalearth)

# Get data
pol = ne_countries(country = "israel", scale = 10, returnclass = "sf")
pol = st_geometry(pol)
pol

## Geometry set for 1 feature 
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: 34.24835 ymin: 29.48969 xmax: 35.88807 ymax: 33.40672
## CRS:            +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0

## POLYGON ((34.24835 31.21145, 34.2644 31.22419, ...

# Plot
plot(pol)