Добавление легенды к диаграмме Вороного со значениями полей

Пытаюсь создать карту со схемой Вороного некоторых климатических станций. Я уже создал диаграмму и раскрасил ее области подходящим цветом. Теперь мне нужна непрерывная легенда, которая расскажет мне, что означают цвета в этих областях.

Вот что уже работает:

library(terra)

f <- system.file("ex/lux.shp", package = "terra")
e <- vect(f)
x <- centroids(e)
v <- voronoi(x, e)
v <- terra::crop(v, e)

FUN <- colorRamp(c("black", "lightgrey"))
vpop <- (v$POP - min(v$POP)) / diff(range(v$POP))
vcols <- rgb(FUN(vpop), maxColorValue = 256)

plot(v, lwd = 2, col = vcols, 
     main = "Voronoi Diagram", 
     background = "beige"
     )
points(x, col = "red")

Мне нужна легенда с правой или нижней стороны графика, которая сопоставляет численность населения территорий с используемыми там цветами. Поскольку исходные данные содержат 50 точек, мне нужна непрерывная легенда. Это то, что я пробовал:

l <- c(5.8, 6.4, 49.35, 49.38)

plot(v, lwd = 2, col = vcols, 
     main = "Voronoi Diagram", 
     background = "beige",
     type = "continuous",
     range = c(min(v$POP, max(v$POP)))
     plg = list(
       ext = l,
       loc = "bottom",
       title = "Population"
       )
     )
points(x, col = "red")

Я также пробовал переместить type и range в список аргумента plg, но это тоже не дало результата. Я также пробовал использовать функцию add_legend из terra.

add_legend("right",
           legend = c(min(v$POP), max(v$POP)),
           col = vcols)

Это также приводит к тому, что на карте плавают только два числа.