data-raw/glottolog/scripts/hexagons.R

In cysouw/qlcData: Processing Data for Quantitative Language Comparison

# check this on hexagon plotting!
library(rmapshaper)
# https://rpubs.com/RobWHickman/testhex


hex_number <- length(cart)
size = 6.3
guessed_hex_centre <- spsample(cart, type = "hexagonal", cellsize = size)
guessed_hex <- HexPoints2SpatialPolygons(guessed_hex_centre, dx = size)
print(length(guessed_hex));print(size)

while( hex_number - length(guessed_hex) > 0 | hex_number - length(guessed_hex) < -20) {
  ifelse(hex_number > length(guessed_hex), size <- size * 0.999, size <- size * 1.001)
  guessed_hex_centre <- spsample(cart, type = "hexagonal", cellsize = size)
  guessed_hex <- HexPoints2SpatialPolygons(guessed_hex_centre, dx = size)
  print(length(guessed_hex));print(size)
}

whole_shape <- ms_dissolve(guessed_hex)
plot(whole_shape)


proj4string(guessed_hex) <- CRS(mollweide)
proj4string(cart) <- CRS(mollweide)
hex <- spTransform(guessed_hex, CRS(mollweide))
ids <- data.frame("ID" = 1:length(hex@polygons))
hex <- SpatialPolygonsDataFrame(hex, ids, match.ID = FALSE)

library(rgeos)
link <- Matrix(gIntersects(hex,cart, byid = T), sparse = T)
table(rowSums(link))
table(colSums(link))