1

I would like to label Voronoi polygons with point attributes (in this case, station name and station number). To create Voronoi polygons, I used the function from this topic.

My station data file (shapefile) is in here.

pts <- readOGR(dsn = ".", layer="rainfallStations")

thiessen.pol <- voronoipolygons(pts)
proj <- CRS('+proj=longlat +ellps=WGS84')
proj4string(thiessen.pol) <- proj
Community
  • 1
  • 1
Nahm
  • 143
  • 1
  • 5

1 Answers1

1

You could text with coordinates, e.g.:

plot(thiessen.pol)
text(coordinates(thiessen.pol), labels=pts$STA_NUMBER, cex=0.9)

enter image description here

jbaums
  • 27,115
  • 5
  • 79
  • 119