如何标记Thiessen（Voronoi）多边形

Question

我想用点属性标记Voronoi多边形（在这种情况下，站名和站号）。为了创建Voronoi多边形，我使用了this topic中的函数。

我的电台数据文件（shapefile）位于here。

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

thiessen.pol <- voronoipolygons(pts)
proj <- CRS('+proj=longlat +ellps=WGS84')
proj4string(thiessen.pol) <- proj

Answer 1

您text coordinates，例如：

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