R График заполненных ячеек сетки долготы-широты на карте

У меня есть кадр данных, содержащий ряд (x, y, z) точек данных, (x, y) правая нижняя координата ячейки долгота-широта размером w(например, сетка с шагом 1 градус). Значение zбыло усреднено по этой ячейке.

Я хотел бы нарисовать эти точки в R так, чтобы вся ячейка сетки была заполнена некоторым цветом, полученным из z.

Результат будет выглядеть примерно так:

Сама проекция (например, конформная коническая Ламберта, равнопрямоугольная) не важна, важно только построение ячеек сетки.

Мои данные скудны: не каждая ячейка долготы-широты будет иметь связанные с ней данные.

Я надеюсь, что решение будет похоже на

library(maps)
map("state")
grid_points(my_data,c("x","y","z"),0.5)

, где 0,5 — это разрешение сетки выше, указывающее на ячейку с углом 0,5 градуса.

Есть мысли?

Спасибо!