Generación de una red hexagonal usando Spatstat

Estoy analizando el patrón de crecimiento de ciertas partículas y quiero comparar el patrón de puntos con el de una red perfectamente hexagonal con la misma intensidad (mismo número de puntos por unidad de área). He escrito una función que hace esto, pero tiene un error inherente y no estoy seguro de dónde se origina. Esencialmente, después de que la función ha seguido su curso, produce un patrón de puntos perfectamente hexagonal que NO tiene el número correcto de partículas; por lo general, está desviado entre un 1 y un 4 %. Si copia y pega el siguiente código en R, verá esto (para este ejemplo en particular, el error es del 11.25 %): el patrón de puntos original tiene 71 partículas y el patrón de puntos perfectamente hexagonal que se generó tiene 80 partículas. Esto parece muy extraño, ya que el código, como verá, está diseñado para colocar las partículas a una distancia específica entre sí, creando así una ventana del mismo tamaño que la original con el mismo número de partículas.

El siguiente es el código de la función que escribí para generar la red hexagonal.

library(spatstat)

data(swedishpines)

swedishpines.df <- as.data.frame(swedishpines)

MaxXValue <- max(swedishpines.df[1])
MaxYValue <- max(swedishpines.df[2])
#The above two lines dictate the window size
NumberOfParticles <- nrow(swedishpines.df[1])
#Number of entries = number of particles

#calculate delta
intensity <- NumberOfParticles / (MaxXValue*MaxYValue)
#Intensity ---> particles / unit^2
    #Area = ( sqrt(3) / 2 ) * delta^2
    #Now - in substituting intensity in for area, it is key to recognize 
    #that there are 3 particles associated with each hexagonal tile.
    #This is because each particle on the border of the tile is really 1/3 of a 
    #a particle due to it being associated with 3 different hexagonal tiles.  
    #Since intensity = 3 Particles / Area, 
delta <- sqrt(2/(intensity*(sqrt(3))))
#This is derived from the equation for the area of a regular hexagon. 

#xcoords and ycoords represent the x and y coordintes of all of the generated points.  The 'x' and 'y' are temporary holders for the x and y coordinates of a single horizontal line of points (they are appended to xcoords and ycoords at the end of each while loop).  

xcoords <- c()
ycoords <- c()

#The following large while loop calculates the coordinates of the first set of points that are vertically aligned with one another. (alternating horizontal lines of particles)  This is shown in the image below.

y_shift <- 0
while (y_shift < MaxYValue) {

    x <- c(0)
    x_shift <- 0 + delta
    count <- 0

    while (x_shift < MaxXValue) {
        x <- append(x, x_shift)
        x_shift <- x_shift + delta
        count <- count + 1
    }

    y <- c(y_shift)

    for (i in seq(0,(count-1))) {
        y <- append(y, y_shift)
    }

    y_shift <- y_shift + sqrt(3)*delta
    xcoords <- append(xcoords,x)
    ycoords <- append(ycoords,y)

}

#The following large while loop calculates the coordinates of the second set of points that are vertically aligned with one another. This is shown in the image below. 

y_shift <- 0 + (delta*(sqrt(3)))/2
while (y_shift < MaxYValue) {

    x <- c(0 + (1/2)*delta)
    x_shift <- 0 + (1/2)*delta + delta
    count <- 0

    while (x_shift < MaxXValue) {
        x <- append(x, x_shift)
        x_shift <- x_shift + delta
        count <- count + 1
    }

    y <- c(y_shift)

    for (i in seq(0,(count-1))) {
        y <- append(y, y_shift)
    }

    y_shift <- y_shift + sqrt(3)*delta
    xcoords <- append(xcoords,x)
    ycoords <- append(ycoords,y)

}

hexplot <- ppp(xcoords, ycoords, window=owin(c(0,MaxXValue),c(0,MaxYValue)))

Ahora, soy relativamente nuevo en R, por lo que puede ser un error de sintaxis en alguna parte del código que haya causado este error. O puede ser que tenga alguna falla en mi línea de pensamiento con este proceso. Sin embargo, no creo que eso sea probable, ya que mis resultados han estado bastante cerca de lo que he estado intentando (solo un error del 1 al 4 % la mayoría de las veces es bastante bueno).

En resumen, en lo que me gustaría recibir ayuda es en cómo tomar un patrón de puntos y crear otro patrón de puntos del mismo tamaño de ventana con el mismo número de partículas, pero un patrón de puntos perfectamente hexagonal.

Por favor, no dude en pedirme que aclare cualquier cosa si siente que algo no está claro.

¡Gracias!