水彩ボロノイ図

概要

こちらは水彩画のようなボロノイ図を作成する記事です。
ここではそのようなボロノイ図を便宜上、"水彩ボロノイ図"と呼ぶことにします。
何か不備がありましたら指摘していただければ幸いです。

コード

0.5秒おきに水彩ボロノイ図が自動生成されるプログラムです。

Main.cpp

/*水彩ボロノイ図を生成*/

#include "DxLib.h"

//横画面サイズ
#define MAP_X 512

//縦画面サイズ
#define MAP_Y 512

//原点サイズ
#define PX_SIZE 16

//各位置のID
#define LU VoronoiMap[im][jm]
#define U VoronoiMap[i][jm]
#define RU VoronoiMap[ip][jm]
#define L VoronoiMap[im][j]
#define C VoronoiMap[i][j]
#define R VoronoiMap[ip][j]
#define LD VoronoiMap[im][jp]
#define D VoronoiMap[i][jp]
#define RD VoronoiMap[ip][jp]

//各位置のX座標
#define LU_X WorldX[im][jm]
#define U_X WorldX[i][jm]
#define RU_X WorldX[ip][jm]
#define L_X WorldX[im][j]
#define C_X WorldX[i][j]
#define R_X WorldX[ip][j]
#define LD_X WorldX[im][jp]
#define D_X WorldX[i][jp]
#define RD_X WorldX[ip][jp]

//各位置のY座標
#define LU_Y WorldY[im][jm]
#define U_Y WorldY[i][jm]
#define RU_Y WorldY[ip][jm]
#define L_Y WorldY[im][j]
#define C_Y WorldY[i][j]
#define R_Y WorldY[ip][j]
#define LD_Y WorldY[im][jp]
#define D_Y WorldY[i][jp]
#define RD_Y WorldY[ip][jp]


/*原点を生成する*/
void DrawPX(int x, int y, int l, int id, unsigned char VoronoiMap[][MAP_Y], unsigned short WorldX[][MAP_Y], unsigned short WorldY[][MAP_Y]) {
    for (int i = x; i < x + l; i++) {
        for (int j = y; j < y + l; j++) {
            VoronoiMap[i][j] = id;
            WorldX[i][j] = i;
            WorldY[i][j] = j;
        }
    }
    return;
}


/*ボロノイ図を生成する*/
void DrawVoronoi(unsigned char VoronoiMap[][MAP_Y]) {

    //マップを初期化
    for (unsigned int i = 0; i < MAP_X; i++) {
        for (unsigned int j = 0; j < MAP_Y; j++) {
            C = 0;
        }
    }

    //マップの座標を格納する
    unsigned short WorldX[MAP_X][MAP_Y] = { 0 };
    unsigned short WorldY[MAP_X][MAP_Y] = { 0 };

    //生成する原点の個数を決める
    int mapCount = 1 + GetRand(12);

    //原点を生成
    for (int i = 1; i <= mapCount; i++) {
        DrawPX(GetRand(MAP_X - PX_SIZE), GetRand(MAP_Y - PX_SIZE), PX_SIZE, i, VoronoiMap, WorldX, WorldY);
    }

    //上下左右の座標
    unsigned int im;
    unsigned int ip;
    unsigned int jm;
    unsigned int jp;

    //原点からの距離の2乗
    unsigned int sizeXY;


    for (unsigned int k = MAP_X; k > 0; k >>= 1) {

        for (unsigned int i = 0; i < MAP_X; i++) {

            im = (MAP_X + i - k) & (MAP_X - 1);
            ip = (i + k) & (MAP_X - 1);

            for (unsigned int j = 0; j < MAP_Y; j++) {

                jm = (MAP_Y + j - k) & (MAP_Y - 1);
                jp = (j + k) & (MAP_Y - 1);

                //原点と中心点の距離の2乗を代入
                sizeXY = (C_X - i)*(C_X - i) + (C_Y - j)*(C_Y - j);

                //上と下を比較
                if (U && U == D && (!C || sizeXY >(D_X - i)*(D_X - i) + (D_Y - j)*(D_Y - j))) {
                    C = D;
                    C_X = D_X;
                    C_Y = D_Y;
                }

                //左と右を比較
                if (L && L == R && (!C || sizeXY > (R_X - i)*(R_X - i) + (R_Y - j)*(R_Y - j))) {
                    C = R;
                    C_X = R_X;
                    C_Y = R_Y;
                }

                //左上と右下を比較
                if (LU && LU == RD && (!C || sizeXY > (RD_X - i)*(RD_X - i) + (RD_Y - j)*(RD_Y - j))) {
                    C = RD;
                    C_X = RD_X;
                    C_Y = RD_Y;
                }

                //右上と左下を比較
                if (RU && RU == LD && (!C || sizeXY > (LD_X - i)*(LD_X - i) + (LD_Y - j)*(LD_Y - j))) {
                    C = LD;
                    C_X = LD_X;
                    C_Y = LD_Y;
                }

                //色を分けて線画
                switch (VoronoiMap[i][j])
                {
                case 0:DrawPixel(i, j, GetColor(255, 255, 255)); break;
                case 1:DrawPixel(i, j, GetColor(255, 0, 0)); break;
                case 2:DrawPixel(i, j, GetColor(0, 255, 0)); break;
                case 3:DrawPixel(i, j, GetColor(0, 0, 255)); break;
                case 4:DrawPixel(i, j, GetColor(255, 255, 0)); break;
                case 5:DrawPixel(i, j, GetColor(255, 0, 255)); break;
                case 6:DrawPixel(i, j, GetColor(0, 255, 255)); break;
                case 7:DrawPixel(i, j, GetColor(127, 127, 127)); break;
                case 8:DrawPixel(i, j, GetColor(127, 0, 0)); break;
                case 9:DrawPixel(i, j, GetColor(0, 127, 0)); break;
                case 10:DrawPixel(i, j, GetColor(0, 0, 127)); break;
                case 11:DrawPixel(i, j, GetColor(127, 127, 0)); break;
                case 12:DrawPixel(i, j, GetColor(127, 0, 127)); break;
                case 13:DrawPixel(i, j, GetColor(0, 127, 127)); break;
                }

            }

        }
    }

}

/*メイン関数*/
int WINAPI WinMain(HINSTANCE, HINSTANCE, LPSTR, int) {

    //log出力停止・ウィンドウモード変更・初期化・裏画面設定
    SetOutApplicationLogValidFlag(FALSE), ChangeWindowMode(TRUE), DxLib_Init(), SetDrawScreen(DX_SCREEN_BACK);

    //画面サイズの決定
    SetGraphMode(MAP_X, MAP_Y, 32);

    //画面サイズが偶数ではない場合
    if (MAP_X & 1 || MAP_Y & 1) {
        //DXライブラリの終了処理
        DxLib_End();
        return -1;
    }

    //マップ
    unsigned char VoronoiMap[MAP_X + 1][MAP_Y] = { 0 };

    //時間
    int time = 0;

    //ボロノイ図を生成
    DrawVoronoi(VoronoiMap);

    while (ScreenFlip() == 0 && ProcessMessage() == 0) {
        //0.5秒おきにボロノイ図を生成
        if (time == 30) {
            DrawVoronoi(VoronoiMap);
            time = 0;
        }

        //ESCキーで終了
        if (CheckHitKey(KEY_INPUT_ESCAPE)) break;

        //時間を増やす
        time++;
    }

    //DXライブラリの終了処理
    DxLib_End();
    return 0;
}

結果

実行してみたらこんな感じ

アニメーションしてみたらこんな感じ

ボロノイ図と違う点は、境界線が凸凹している点ですね。
この性質を利用して例えば、自然地形を線画すると自然な出来になります。

ソースコードのライセンス

These codes are licensed under CC0.

ソースコードは自由に使用してください。

参考

ＤＸライブラリ　関数リファレンスページ
新・C言語 ～ゲームプログラミングの館
ボロノイ図
ボロノイ図を作る
ボロノイ図とその３つの性質
ボロノイ分割とは？
GPUでボロノイ図を描画する
ボロノイ図の描画アルゴリズムについて