library(sf)
## Linking to GEOS 3.8.1, GDAL 3.1.3, PROJ 7.1.1
# Read
dat = st_read("/home/michael/Downloads/גבולות נפות/
(ArcView)/
.shp", options = "ENCODING=UTF-8")
## options: ENCODING=UTF-8
## Reading layer `
' from data source `/home/michael/Downloads/גבולות נפות/
(ArcView)/
.shp' using driver `ESRI Shapefile'
## Simple feature collection with 15 features and 4 fields
## geometry type: POLYGON
## dimension: XY
## bbox: xmin: 130175.7 ymin: 377832.4 xmax: 284068.5 ymax: 804530.3
## projected CRS: Israel 1993 / Israeli TM Grid
# Plot
plot(st_geometry(dat), border = "grey")
text(st_coordinates(st_centroid(dat)), dat$Nafa)
## Warning in st_centroid.sf(dat): st_centroid assumes attributes are constant over
## geometries of x
