Em softwares como o ArcMap, pode-se criar centróides para polígonosdentro um polígono. Em casos como o mostrado abaixo, isso é necessário.

NoR é possível calcular centróides de polígonos espaciais comrgeos::gCentroid(). No entanto, não há como forçar o cálculo de centróides dentro do polígono.

library(rgdal)
library(rgeos)

x <- readWKT("POLYGON ((1441727.5096940901130438 6550163.0046194596216083, 
             1150685.2609429201111197 6669225.7427449300885201, 
             975398.4520359700545669 6603079.7771196700632572, 
             866257.6087542800232768 6401334.5819626096636057, 
             836491.9242229099618271 6106985.0349301798269153, 
             972091.1537546999752522 5835786.5758665995672345, 
             1547561.0546945100650191 5782869.8033663900569081, 
             1408654.5268814601004124 5600968.3978968998417258, 
             720736.4843787000281736 5663807.0652409195899963, 
             598366.4479719599476084 6001151.4899297598749399, 
             654590.5187534400029108 6341803.2128998702391982, 
             869564.9070355399744585 6784981.1825891500338912, 
             1451649.4045378800947219 6788288.4808704098686576, 
             1441727.5096940901130438 6550163.0046194596216083))")
plot(x)

Este é o polígonox

gCentroid() cria um centróide que, neste caso específico, está localizado fora do polígono. Apesar de serem geometricamente corretos, alguns aplicativos exigem centróides dentro do polígono, pois podem ser calculados pelo ArcMap.

xCent <- gCentroid(x, byid = TRUE)
points(xCent, col = "red", pch = 16)

Uma saída desejada (do ArcMap) é semelhante a esta:

Existe alguma possibilidade de gerar centróides como este em R?

EDITAR:

Após algunsescavação, acontece que o ArcMap escolhe um ponto aleatório dentro do polígono:

"Para um polígono de entrada: o ponto de saída estará dentro do polígono."

Portanto, a pergunta deve ser: existe uma função que cria um ponto em qualquer posição aleatória DENTRO dos polígonos?