Создавайте квадратные сетки и экспортируйте их как шейп-файл или таблицу

Это не удается, потому что вы используете зону 32 utm. Вам нужно использовать зону, основанную на долготе страны. Вы можете увидеть их здесь

Вы можете автоматизировать поиск зоны с помощью ceiling((longitude+180)/6)

library(raster)
s <- getData(country = "FRA", level = 0)

Получите центроид. В этом случае вы можете сделать

centr <- coordinates(s)

Если есть несколько полигонов, вы можете сделать что-то вроде этого

centr <- apply(coordinates(s), 2, mean)

Вычислить зону UTM. (обратите внимание, что у вас было 32 для Франции, что не очень хорошо)

zone <- ceiling((centr[1] + 180)/6)
zone
#[1] 31

А затем используйте это так

crs <- paste0("+proj=utm +datum=WGS84 +unit=m +zone=", zone)
st <- spTransform(s, crs)

Для Таиланда вы получите

s <- getData(country = "THA", level = 0)
centr <- apply(coordinates(s), 2, mean)
zone <- ceiling((centr[1] + 180)/6)
zone
#[1] 47

Однако это не подход, который подойдет для всех стран. Зоны UTM имеют ширину 6 градусов, и многие страны охватывают несколько зон (Россия занимает первое место с 28 зонами). Поэтому, в зависимости от ваших целей, вы можете использовать другую систему координат (crs).

После этого альтернативным способом получения квадратных многоугольников является создание RasterLayer с размером s и выбранным разрешением. Но я сомневаюсь, что это лучший способ получить данные из GEE. Вместо этого я бы посоветовал загрузить схему страны.

r <- raster(st, res=10000)
r <- rasterize(st, r, 1)
x <- as(r, "SpatialPolygons")

# write to file
shapefile(x, "test.shp")

# view
plot(x)