I’m trying to overlay a voronoi tesselation onto a secondary data set. The primary dataset include the gps coordinates of trees. The secondary dataset include soil depth scans with GPS coordinates. I’ve done the voronoi tesselations for the individual trees(using deldir() package). How can I assign soil depths (soil depth scan coordinates), to the respective trees/tiles. .