【发布时间】:2014-12-07 07:20:23
【问题描述】:
我想用点属性(在本例中为站名和站号)标记 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
【问题讨论】: