Yes, it's all done by a script I'm currently building.
This one is actually half resolution, so it has around 1500 or so regions. Takes a few minutes to calculate. The problem is quadratic increases in time and memory, especially in the voronoi calculations, preventing me from going much higher. I tried 4096x4096 but it doesn't work out. Crashed with a segfault after about 15 minutes of calculation.