#include "igs_os_type.h"
#include <stdio.h>
#include <math.h>
#include <stdint.h>
#include <string.h>
#include <stdarg.h>
#include <time.h>
#include <assert.h>
#include <stdlib.h>
#include <iostream>
#include <stdexcept>
namespace
{

#ifndef _pri_h_
#define _pri_h_

/* Windowsではstdint.hが見付からない */
#if defined _WIN32
typedef int int32_t;
#else
#include <stdint.h> /* for int32_t */
#endif

#ifdef __cplusplus
extern "C" {
#endif

extern void pri_funct_cv_start(int32_t i32_ys);
extern void pri_funct_cv_run(int32_t i32_y);
extern void pri_funct_cv_end(void);

extern void pri_funct_set_cp_title(char *cp_title);
extern void pri_funct_msg_ttvr(const char *fmt, ...);
extern void pri_funct_msg_vr(const char *fmt, ...);
extern void pri_funct_err_bttvr(const char *fmt, ...);

#ifdef __cplusplus
}
#endif

#endif /* !_pri_h_ */
#include <stdio.h>

#include "igs_line_blur.h" // "pri.h"

static int32_t pri_param_i32_ysize,
	pri_param_i32_pos_before;

/* カウントダウン開始 */
void pri_funct_cv_start(int32_t i32_ys)
{
	/* 目盛表示 */
	(void)fprintf(stdout, "0%%  10   20   30   40   50   60   70   80   90   100%%\n");
	(void)fprintf(stdout, "....v....v....v....v....v....v....v....v....v....v\n");

	/* パラメータ初期設定 */
	pri_param_i32_ysize = i32_ys;
	pri_param_i32_pos_before = 0;
}
/* カウントダウン実行中 */
void pri_funct_cv_run(int32_t i32_y)
{
	int32_t i32_pos_current, ii;

	/* 目盛上の現在位置(1...50) */
	i32_pos_current = 50 * (i32_y + 1) / pri_param_i32_ysize;

	/* 一回で45度くるくるカウンター */
	/******switch (i32_y%4) {
	case 0: (void)fprintf(stdout, "-\b"  ); break;
	case 1: (void)fprintf(stdout, "\\\b" ); break;
	case 2: (void)fprintf(stdout, "|\b"  ); break;
	case 3: (void)fprintf(stdout, "/\b"  ); break;
	}
	(void)fflush(stdout);******/

	/* 一回前の目盛位置と同じなら抜ける */
	if (i32_pos_current == pri_param_i32_pos_before)
		return;

	/* 目盛が進んだとき */
	for (ii = pri_param_i32_pos_before; ii < i32_pos_current; ++ii) {
		/* 一目盛以上進んだときのすき間 */
		if ((ii + 1) < i32_pos_current) {
			(void)fprintf(stdout, ".");
		}
		/* 一目盛分 */
		else {
			(void)fprintf(stdout, "^");
		}
	}
	(void)fflush(stdout);

	/* 現在位置を記憶する */
	pri_param_i32_pos_before = i32_pos_current;
}
/* カウントダウン終了 */
void pri_funct_cv_end(void)
{
	/* 改行 */
	(void)fprintf(stdout, "\n");
}

#include <stdio.h>
#include <string.h>
#include <stdarg.h>
#include <time.h>

#include "igs_line_blur.h" // "pri.h"

/* Windowsではvsnprintf()の頭にアンダーバーが付く!!! */
#if defined _WIN32
#define vsnprintf(buf, len, fmt, ap) _vsnprintf(buf, len, fmt, ap)
#endif

static char *pri_param_cp_com_name = "#";

void pri_funct_set_cp_title(char *cp_title)
{
	pri_param_cp_com_name = cp_title;
}

/* print  Valiable_argument, Return_code */
void pri_funct_msg_vr(const char *fmt, ...)
{
	int i_ret;
	va_list ap;
	char buf[FILENAME_MAX]; /* redhat9 : stdio.h --> bits/stdio_lim.h */

	va_start(ap, fmt);
	i_ret = vsnprintf(buf, FILENAME_MAX, fmt, ap);
	va_end(ap);
	if (i_ret < 0) {
		(void)strcpy(buf, "bad argument for fprintf stdout");
	}

	/* 可変引数,改行 */
	(void)fprintf(stdout, "%s\n", buf);
	(void)fflush(stdout);
}

/* print  Title, Time_stamp, Valiable_argument, Return_code */
void pri_funct_msg_ttvr(const char *fmt, ...)
{
	time_t tt;
	char *cc;

	int i_ret;
	va_list ap;
	char buf[FILENAME_MAX]; /* redhat9 : stdio.h --> bits/stdio_lim.h */

	tt = time(NULL);
	cc = asctime(localtime(&tt)); /* 注意：ccはstatic変数 */
	cc[24] = '\0';				  /* 26-character Example : "Fri Sep 13 00:00:00 1986\n\0" */

	va_start(ap, fmt);
	i_ret = vsnprintf(buf, FILENAME_MAX, fmt, ap);
	va_end(ap);
	if (i_ret < 0) {
		(void)strcpy(buf, "bad argument for fprintf stdout");
	}

	/* ベル,タイトル,日時,可変引数,改行 */
	(void)fprintf(stdout, "%s  %s  %s\n", pri_param_cp_com_name, cc, buf);
	(void)fflush(stdout);
}
/* print  Bell, Title, Time_stamp, Valiable_argument, Return_code */
void pri_funct_err_bttvr(const char *fmt, ...)
{
	time_t tt;
	char *cc;

	int i_ret;
	va_list ap;
	char buf[FILENAME_MAX];

	tt = time(NULL);
	cc = asctime(localtime(&tt)); /* 注意：ccはstatic変数 */
	cc[24] = '\0';				  /* 26-character Example : "Fri Sep 13 00:00:00 1986\n\0" */

	va_start(ap, fmt);
	i_ret = vsnprintf(buf, FILENAME_MAX, fmt, ap);
	va_end(ap);
	if (i_ret < 0) {
		(void)strcpy(buf, "bad argument for fprintf stderr");
	}

	/* ベル,タイトル,日時,可変引数,改行 */
	(void)fprintf(stderr, "\007%s  %s  %s\n", pri_param_cp_com_name, cc, buf);
	(void)fflush(stderr);
}

#ifndef _list_node_h_
#define _list_node_h_

#include <stdio.h>

class list_node
{
public:
	list_node()
	{
		this->_clp_previous = NULL;
		this->_clp_next = NULL;
	}

	list_node *get_clp_next(void) { return this->_clp_next; }
	list_node *get_clp_previous(void) { return this->_clp_previous; }

	void set_clp_next(list_node *clp) { this->_clp_next = clp; }
	void set_clp_previous(list_node *clp) { this->_clp_previous = clp; }

private:
	list_node *_clp_previous,
		*_clp_next;
};

#endif /* !_list_node_h_ */
#ifndef _list_root_h_
#define _list_root_h_

#include <stdio.h>

/* Windowsではstdint.hが見つからない */
#if defined _WIN32
typedef int int32_t;
typedef unsigned short uint16_t;
#else
#include <stdint.h> /* for int32_t, uint16_t */
#endif

#include "igs_line_blur.h" // "list_node.h"

class list_root
{
public:
	list_root()
	{
		this->_clp_first = NULL;
		this->_clp_last = NULL;
		this->_i32_count = 0;
	}

	list_node *push(list_node *clp_previous, list_node *clp_new);
	list_node *pop(list_node *clp_old);

	list_node *get_clp_first(void) { return this->_clp_first; }
	list_node *get_clp_last(void) { return this->_clp_last; }

	int32_t get_i32_count(void) { return this->_i32_count; }

private:
	list_node *_clp_first,
		*_clp_last;
	int32_t _i32_count;

	void _set_clp_first(list_node *clp_) { this->_clp_first = clp_; }
	void _set_clp_last(list_node *clp_) { this->_clp_last = clp_; }

	void _connect(list_node *clp_previous, list_node *clp_next);
};

#endif				/* !_list_root_h_ */
#include <assert.h> /* assert() */

#include "igs_line_blur.h" // "list_root.h"

/* ピクセルリストを前後でつなげるメソッド */
void list_root::_connect(list_node *clp_previous, list_node *clp_next)
{
	/* リスト初期化 */
	if ((NULL == clp_previous) && (NULL == clp_next)) {
		this->_clp_first = NULL;
		this->_clp_last = NULL;
	}
	/* (clp_previous)をリストの最後部に追加 */
	else if ((NULL != clp_previous) && (NULL == clp_next)) {
		this->_clp_last = clp_previous;
		clp_previous->set_clp_next(NULL);
	}
	/* (clp_next)をリストの最前部に挿入 */
	else if ((NULL == clp_previous) && (NULL != clp_next)) {
		this->_clp_first = clp_next;
		clp_next->set_clp_previous(NULL);
	}
	/* (clp_previous)と、(clp_next)をつなげる */
	else if ((NULL != clp_previous) && (NULL != clp_next)) {
		clp_next->set_clp_previous(clp_previous);
		clp_previous->set_clp_next(clp_next);
	}
}

list_node *list_root::push(list_node *clp_previous, list_node *clp_new)
{
	/* あってはならないプログラムバグのチェック */
	assert(NULL != clp_new);

	/* リスト上で、clp_newの後ろを接続 */
	if (NULL != clp_previous) { /* clp_previousの後ろに接続 */
		this->_connect(clp_new, clp_previous->get_clp_next());
	} else { /* clp_previousがNULLなら、リストの頭に挿入 */
		this->_connect(clp_new, this->get_clp_first());
	}

	/* リスト上で、clp_newの前を接続 */
	this->_connect(clp_previous, clp_new);

	/* 子リストの数のカウント */
	++(this->_i32_count);

	return clp_new;
}
list_node *list_root::pop(list_node *clp_old)
{
	/* あってはならないプログラムバグのチェック */
	assert(NULL != clp_old);

	/* はずすノードの前後をつなげる */
	this->_connect(clp_old->get_clp_previous(), clp_old->get_clp_next());

	/* 子リストの数のカウント */
	--(this->_i32_count);

	/* あってはならないプログラムバグのチェック */
	assert(0 <= (this->_i32_count));

	return clp_old;
}

#ifndef _brush_curve_blur_h_
#define _brush_curve_blur_h_

#include <stdio.h>

/* Windowsではstdint.hが見つからない */
#if defined _WIN32
typedef int int32_t;
#else
#include <stdint.h> /* for int32_t */
#endif

#ifndef ON
#define ON (1)
#endif
#ifndef OFF
#define OFF (0)
#endif

#ifndef OK
#define OK (0)
#endif
#ifndef NG
#define NG (-1)
#endif

#ifndef CHANNEL_COUNT
#define CHANNEL_COUNT (4)
#endif

class brush_curve_blur
{
public:
	/* constructer */
	brush_curve_blur(void)
	{
		int32_t ii;

		this->_i_mv_sw = OFF;
		this->_i_pv_sw = OFF;
		this->_i_cv_sw = OFF;

		this->_i32_count = 51;
		this->_i32_subpixel_divide = 2;
		this->_d_effect_area_radius = 25.0;
		this->_d_power = 1.0;

		this->_dp_ratio = NULL;
		this->_dp_linepixels = NULL;
		this->_dp_xp = NULL;
		this->_dp_yp = NULL;
		this->_dp_subpixel = NULL;
		for (ii = 0; ii < CHANNEL_COUNT; ++ii) {
			this->_da_pixel[ii] = 0.0;
		}
	}
	/* constructer */
	~brush_curve_blur(void)
	{
		this->mem_free();
	}

	/* パラメータ設定 */
	void set_i_mv_sw(int ii) { this->_i_mv_sw = ii; }
	void set_i_pv_sw(int ii) { this->_i_pv_sw = ii; }
	void set_i_cv_sw(int ii) { this->_i_cv_sw = ii; }

	/* ぼかし線のポイント数 */
	void set_i32_count(int32_t ii) { this->_i32_count = ii; }
	int32_t get_i32_count(void) { return this->_i32_count; }

	void set_i32_subpixel_divide(int32_t ii)
	{
		this->_i32_subpixel_divide = ii;
	}
	int32_t get_i32_subpixel_divide(void)
	{
		return this->_i32_subpixel_divide;
	}

	void set_d_effect_area_radius(double dd)
	{
		this->_d_effect_area_radius = dd;
	}
	double get_d_effect_area_radius(void)
	{
		return this->_d_effect_area_radius;
	}

	void set_d_power(double dd)
	{
		this->_d_power = dd;
	}
	double get_d_power(void)
	{
		return this->_d_power;
	}

	/* get for using */
	double *get_dp_linepixels(void) { return this->_dp_linepixels; }
	double *get_dp_xp(void) { return this->_dp_xp; }
	double *get_dp_yp(void) { return this->_dp_yp; }
	double *get_dp_pixel(void) { return this->_da_pixel; }

	void mem_free(void);
	int mem_alloc(void);
	void init_ratio_array(void);
	void set_subpixel_value(int32_t i32_x_sub, int32_t i32_y_sub);
	void set_pixel_value(void);

	int save(double d_xp, double d_yp, char *cp_fname);
	void debug_print(void);

private:
	int _i_mv_sw, /* Method    Verbose */
		_i_pv_sw, /* Parameter Verbose */
		_i_cv_sw; /* Counter   Verbose */

	int32_t _i32_count;
	int32_t _i32_subpixel_divide;
	double _d_effect_area_radius;
	double _d_power;
	double *_dp_ratio,
		*_dp_linepixels,
		*_dp_xp,
		*_dp_yp,
		*_dp_subpixel,
		_da_pixel[CHANNEL_COUNT];
};

#endif				/* !_brush_curve_blur_h_ */
#include <math.h>   /* pow() */
#include <stdlib.h> /* free(), calloc() */
#include <assert.h> /* assert() */

#include "igs_line_blur.h" // "pri.h" "brush_curve_blur.h"

/* メモリ開放 */
void brush_curve_blur::mem_free(void)
{
	if (NULL != this->_dp_ratio) {
		if (ON == this->_i_mv_sw) {
			pri_funct_msg_ttvr("brush_curve_blur::mem_free()");
		}
		free(this->_dp_ratio);
		this->_dp_ratio = NULL;
		this->_dp_linepixels = NULL;
		this->_dp_xp = NULL;
		this->_dp_yp = NULL;
		this->_dp_subpixel = NULL;
	}
}

/* メモリ確保 */
int brush_curve_blur::mem_alloc(void)
{
	/* 以前のメモリが残っていたら開放する */
	this->mem_free();

	/* ユーザー指定がゼロなら動作キャンセル */
	if (this->_i32_count <= 0) {
		return OK;
	}

	/* 処理ごとのメッセージ */
	if (ON == this->_i_mv_sw) {
		pri_funct_msg_ttvr("brush_curve_blur::mem_alloc()");
	}
	if (ON == this->_i_pv_sw) {
		pri_funct_msg_ttvr(
			"alloc brush_curve_blur memory ((%d * %d) + %d) x %d bytes",
			this->_i32_count, (1 + CHANNEL_COUNT + 2), /* ratio,accum,xp,yp */
			this->_i32_subpixel_divide *
				this->_i32_subpixel_divide * CHANNEL_COUNT, /* pixel(rgba) */
			sizeof(double));
	}

	this->_dp_ratio = (double *)calloc(
		/* ratio,linepixels,xp,yp */
		this->_i32_count * (1 + CHANNEL_COUNT + 2) +
			this->_i32_subpixel_divide *
				this->_i32_subpixel_divide * CHANNEL_COUNT, /* pixel(rgba) */
		sizeof(double));

	if (NULL == this->_dp_ratio) {
		pri_funct_err_bttvr("Error : calloc(-) returns NULL.");
		return NG;
	}

	this->_dp_linepixels = this->_dp_ratio + this->_i32_count;
	this->_dp_xp = this->_dp_linepixels + this->_i32_count * CHANNEL_COUNT;
	this->_dp_yp = this->_dp_xp + this->_i32_count;
	this->_dp_subpixel = this->_dp_yp + this->_i32_count;

	return OK;
}

/********************************************************************/

void brush_curve_blur::init_ratio_array(void)
{
	int32_t ii;
	double d_total;

	/* ユーザー指定がゼロなら動作キャンセル */
	if (this->_i32_count <= 0) {
		return;
	}

	/* 前半の変化、ゼロからリニア増 */
	for (ii = 1; ii < this->_i32_count / 2; ++ii) {
		this->_dp_ratio[ii] = ii;
	}
	/* 後半の変化、ゼロへリニア減 */
	for (ii = this->_i32_count / 2; ii < this->_i32_count; ++ii) {
		this->_dp_ratio[ii] = this->_i32_count - ii;
	}

	/* 正規化(0...1)する */
	/****for (ii = 0; ii < this->_i32_count; ++ii) {
		this->_dp_ratio[ii] /=
			(double)(this->_i32_count - this->_i32_count/2);
	}****/
	/* 影響の強さの設定 */
	for (ii = 0; ii < this->_i32_count; ++ii) {
		this->_dp_ratio[ii] = pow(this->_dp_ratio[ii], 1.0 / this->_d_power);
	}

	/* 総計 */
	d_total = 0.0;
	for (ii = 0; ii < this->_i32_count; ++ii) {
		d_total += this->_dp_ratio[ii];
	}
	/* 総計で割って、総計を1にする */
	for (ii = 0; ii < this->_i32_count; ++ii) {
		this->_dp_ratio[ii] /= d_total;
	}
}

/********************************************************************/

/* 線分から取り込んだピクセル値をサブピクセル値として合成 */
void brush_curve_blur::set_subpixel_value(int32_t i32_x_sub, int32_t i32_y_sub)
{
	int32_t ii, zz;
	double d_accum;

	/* ユーザー指定がゼロなら動作キャンセル */
	if (this->_i32_count <= 0) {
		return;
	}

	/* チャネルごと */
	for (zz = 0; zz < CHANNEL_COUNT; ++zz) {
		/* 加算の前にゼロ初期化 */
		this->_dp_subpixel[(this->_i32_subpixel_divide * i32_y_sub + i32_x_sub) *
							   CHANNEL_COUNT +
						   zz] = 0.0;
		/* ピクセル値を加算 */
		d_accum = 0.0;
		for (ii = 0; ii < this->_i32_count; ++ii) {
			/* 画像の範囲外 */
			if (this->_dp_linepixels[ii * CHANNEL_COUNT + zz] < 0.0) {
				continue;
			}
			/* ピクセル値に比を掛けて加算 */
			this->_dp_subpixel[(this->_i32_subpixel_divide * i32_y_sub + i32_x_sub) *
								   CHANNEL_COUNT +
							   zz] += this->_dp_linepixels[ii * CHANNEL_COUNT + zz] *
									  this->_dp_ratio[ii];
			d_accum += this->_dp_ratio[ii];
		}
		/* 画像をスキャンしているのだから、
		画像の範囲となる値はかならず一つ以上ある */
		assert(0.0 < d_accum);

		/* 値の計算 */
		this->_dp_subpixel[(this->_i32_subpixel_divide * i32_y_sub + i32_x_sub) *
							   CHANNEL_COUNT +
						   zz] /= d_accum;
	}
}

/* サブピクセル値を平均してピクセル値を得る */
void brush_curve_blur::set_pixel_value(void)
{
	int32_t xx, yy, zz;

	/* ユーザー指定がゼロなら動作キャンセル */
	if (this->_i32_count <= 0) {
		return;
	}

	for (zz = 0; zz < CHANNEL_COUNT; ++zz) {
		this->_da_pixel[zz] = 0.0;
		for (yy = 0; yy < this->_i32_subpixel_divide; ++yy) {
			for (xx = 0; xx < this->_i32_subpixel_divide; ++xx) {
				this->_da_pixel[zz] +=
					this->_dp_subpixel[(this->_i32_subpixel_divide * yy + xx) *
										   CHANNEL_COUNT +
									   zz];
			}
		}
		this->_da_pixel[zz] /= this->_i32_subpixel_divide *
							   this->_i32_subpixel_divide;
	}
}

#include <assert.h> /* assert() */

#include "igs_line_blur.h" // "pri.h" "brush_curve_blur.h"

int brush_curve_blur::save(double d_xp, double d_yp, char *cp_fname)
{
	FILE *fp;
	int32_t ii;

	/* ファイル開く */
	fp = fopen(cp_fname, "w");
	if (NULL == fp) {
		pri_funct_err_bttvr("Error : fopen(%s,w) returns NULL", cp_fname);
		return NG;
	}

	/* 選択数保存 */
	if (fprintf(fp, "# curve blur count %d\n", this->get_i32_count()) < 0) {
		pri_funct_err_bttvr(
			"Error : fprintf(# curve blur count) returns minus");
		fclose(fp);
		return NG;
	}

	for (ii = 0; ii < this->get_i32_count(); ++ii) {

		/* ピクセル位置から近点位置保存 */
		if (fprintf(fp, "%g %g\n",
					d_xp + this->_dp_xp[ii],
					d_yp + this->_dp_yp[ii]) < 0) {
			pri_funct_err_bttvr(
				"Error : fprintf(avarage x y) returns minus");
			fclose(fp);
			return NG;
		}
	}

	/* ファイル閉じる */
	fclose(fp);

	return OK;
}

#ifndef __brush_smudge_circle_h__
#define __brush_smudge_circle_h__

#include <stdio.h>

/* Windowsではstdint.hが見つからない */
#if defined _WIN32
typedef int int32_t;
typedef unsigned short uint16_t;
#else
#include <stdint.h> /* for int32_t, uint16_t */
#endif

#ifndef ON
#define ON (1)
#endif
#ifndef OFF
#define OFF (0)
#endif

#ifndef OK
#define OK (0)
#endif
#ifndef NG
#define NG (-1)
#endif

class brush_smudge_circle
{
public:
	brush_smudge_circle()
	{
		this->_i_mv_sw = OFF;
		this->_i_pv_sw = OFF;
		this->_i_cv_sw = OFF;

		this->_i32_size_by_pixel = 7; /* 画像上の線の最大幅 */
		this->_i32_subpixel_divide = 4;
		this->_d_ratio = 0.85;

		this->_dp_brush = NULL;
		this->_dp_subpixel_image = NULL;
		this->_dp_pixel_image = NULL;
	}
	~brush_smudge_circle()
	{
		this->mem_free();
	}

	/* パラメータ設定 */
	void set_i_mv_sw(int ii) { this->_i_mv_sw = ii; }
	void set_i_pv_sw(int ii) { this->_i_pv_sw = ii; }
	void set_i_cv_sw(int ii) { this->_i_cv_sw = ii; }

	void set_i32_size_by_pixel(int32_t ii)
	{
		this->_i32_size_by_pixel = ii;
	}
	void set_i32_subpixel_divide(int32_t ii)
	{
		this->_i32_subpixel_divide = ii;
	}
	void set_d_ratio(double dd) { this->_d_ratio = dd; }

	/* パラメータを得る */
	int32_t get_i32_size_by_pixel(void)
	{
		return this->_i32_size_by_pixel;
	}
	int32_t get_i32_subpixel_divide(void)
	{
		return this->_i32_subpixel_divide;
	}
	double get_d_ratio(void) { return this->_d_ratio; }

	void get_dp_area(
		double d_xp, double d_yp,
		double *dp_x1, double *dp_y1, double *dp_x2, double *dp_y2);

	/* メモリ確保 */
	int mem_alloc(void);

	/* メモリへのポインターを得る */
	double *get_dp_brush(void) { return this->_dp_brush; }
	double *get_dp_subpixel_image(void)
	{
		return this->_dp_subpixel_image;
	}
	double *get_dp_pixel_image(void)
	{
		return this->_dp_pixel_image;
	}

	/* 実行 */
	void set_brush_circle(void);
	void copy_to_brush_from_image(void);
	void exec(void);
	void to_subpixel_from_pixel(
		double d_x1, double d_y1, double d_x2, double d_y2);
	void to_pixel_from_subpixel(
		double d_x1, double d_y1, double d_x2, double d_y2);

	/* メモリ開放 */
	void mem_free(void);

private:
	int _i_mv_sw, /* Method    Verbose */
		_i_pv_sw, /* Parameter Verbose */
		_i_cv_sw; /* Counter   Verbose */

	int32_t _i32_size_by_pixel,
		_i32_subpixel_divide;
	double _d_ratio;

	double *_dp_brush,
		*_dp_subpixel_image,
		*_dp_pixel_image;
};

#endif				/* !__brush_smudge_circle_h__ */
#include <stdlib.h> /* free(), calloc() */
#include <math.h>   /* sqrt() */
#include <string.h> /* memset() */
#include <assert.h> /* assert() */

#include "igs_line_blur.h" // "pri.h" "brush_smudge_circle.h"

/* メモリ開放 */
void brush_smudge_circle::mem_free(void)
{
	if (NULL != this->_dp_brush) {
		if (ON == this->_i_mv_sw) {
			pri_funct_msg_ttvr("brush_smudge_circle::mem_free()");
		}
		free(this->_dp_brush);
		this->_dp_brush = NULL;
		this->_dp_pixel_image = NULL;
		this->_dp_subpixel_image = NULL;
	}
}

/* データ設定とメモリ確保 */
int brush_smudge_circle::mem_alloc(void)
{
	int32_t i32_sz;

	/* もし以前のメモリがあるなら開放する */
	this->mem_free();

	/* メモリ確保のためのサイズの1単位 */
	i32_sz = this->_i32_size_by_pixel * this->_i32_subpixel_divide;

	/* 処理ごとのメッセージ */
	if (ON == this->_i_mv_sw) {
		pri_funct_msg_ttvr("brush_smudge_circle::mem_alloc()");
	}

	if (ON == this->_i_pv_sw) {
		pri_funct_msg_ttvr(
			"calloc((%d x %d + %d x %d + %d x %d) x (%d x %d))",
			i32_sz, i32_sz,
			i32_sz, i32_sz,
			(this->_i32_size_by_pixel + 1), (this->_i32_size_by_pixel + 1),
			sizeof(double), 5);
	}

	this->_dp_brush = (double *)calloc(
		i32_sz * i32_sz +
			i32_sz * i32_sz +
			(this->_i32_size_by_pixel + 1) * (this->_i32_size_by_pixel + 1),
		sizeof(double) * 5);
	if (NULL == this->_dp_brush) {
		pri_funct_err_bttvr("Error : calloc(-) returns NULL.");
		return NG;
	}

	this->_dp_subpixel_image = this->_dp_brush +
							   i32_sz * i32_sz * 5;
	this->_dp_pixel_image = this->_dp_subpixel_image +
							i32_sz * i32_sz * 5;

	return OK;
}

/********************************************************************/

/* 画像上に置いたブラシの範囲
	+-------+-------+ --> *dp_y2
	|       |       |
	|   +   |   +   |
	|       |       |
	+-------+-------+ <-- d_yp
	|       |       |
	|   +   |   +   |
	|       |       |
	+-------+-------+ --> *dp_y2
	|       ^       |
	v       |       v
	*dp_x1  d_xp    *dp_x2
	    0       1     <---- d_xp,d_yp is pixel position
	                     |
	                     | +0.5
	                     v
	0       1       2 <---- *dp_x1,*dp_y1,*dp_x2,*dp_y2 is image position
*/
void brush_smudge_circle::get_dp_area(double d_xp, double d_yp, double *dp_x1, double *dp_y1, double *dp_x2, double *dp_y2)
{
	*dp_x1 = d_xp + 0.5 - this->_i32_size_by_pixel / 2.0;
	*dp_x2 = *dp_x1 + this->_i32_size_by_pixel;
	*dp_y1 = d_yp + 0.5 - this->_i32_size_by_pixel / 2.0;
	*dp_y2 = *dp_y1 + this->_i32_size_by_pixel;
}

/********************************************************************/

/* ブラシの形を、有効フラグの設定で決める
   例
	pixel_size      3 x 3
	subpixel_divide 2 x 2
	+-------+-------+-------+   ^
	| +   + | +   + | +   + | 5 |
	|       |       |       |   |
	| +   + | +   + | +   + | 4 |
	+-------+-------+-------+   |
	| +   + | +   + | +   + | 3 |
center->|       |       |       |  i32_size(3*2)
	| +   + | +   + | +   + | 2 |
	+-------+-------+-------+   |
	| +   + | +   + | +   + | 1 |
	|       |       |       |   |
	| +   + | +   + | +   + | 0 |
	+-------+-------+-------+ ^ v
	  0   1   2   3   4   5 <-+---- subpixel count
	0   1   2   3   4   5   6  <----- position
	<---- i32_size(3*2) ---->
	            ^
	            |
	           center
*/
void brush_smudge_circle::set_brush_circle(void)
{
	double *dp_brush, d_radius, d_tmp;
	int32_t xx, yy,
		i32_size;

	dp_brush = this->_dp_brush;
	i32_size = this->_i32_size_by_pixel * this->_i32_subpixel_divide;

	d_radius = i32_size / 2.0;

	for (yy = 0; yy < i32_size; ++yy) {
		for (xx = 0; xx < i32_size; ++xx, dp_brush += 5) {
			d_tmp = sqrt(
				(xx + 0.5 - d_radius) * (xx + 0.5 - d_radius) +
				(yy + 0.5 - d_radius) * (yy + 0.5 - d_radius));
			if (d_tmp < d_radius) {
				dp_brush[4] = 1.0;
			} else {
				dp_brush[4] = 0.0;
			}
		}
	}
}

/* 切り取った画像をブラシにセット */
void brush_smudge_circle::copy_to_brush_from_image(void)
{
	double *dp_brush, *dp_image;
	int32_t xx, yy, zz, i32_size;

	dp_brush = this->_dp_brush;
	dp_image = this->_dp_subpixel_image;
	i32_size = this->_i32_size_by_pixel * this->_i32_subpixel_divide;

	for (yy = 0; yy < i32_size; ++yy) {
		for (xx = 0; xx < i32_size; ++xx, dp_brush += 5, dp_image += 5) {
			if (0.0 < dp_brush[4]) {
				/* rgbaの4チャンネルをコピーする */
				for (zz = 0; zz < 4; ++zz) {
					dp_brush[zz] = dp_image[zz];
				}
			}
		}
	}
}

/* こする
	以前の場所での画像と指定の比率で差分を画像に足し、その画像をブラシに保存
 */
void brush_smudge_circle::exec(void)
{
	double *dp_brush, *dp_image;
	int32_t xx, yy, zz, i32_size;

	dp_brush = this->_dp_brush;
	dp_image = this->_dp_subpixel_image;
	i32_size = this->_i32_size_by_pixel * this->_i32_subpixel_divide;

	for (yy = 0; yy < i32_size; ++yy) {
		for (xx = 0; xx < i32_size; ++xx, dp_brush += 5, dp_image += 5) {
			if (0.0 < dp_brush[4]) {
				/* rgbaの4チャンネルで実行 */
				for (zz = 0; zz < 4; ++zz) {
					dp_image[zz] += (dp_brush[zz] - dp_image[zz]) *
									this->_d_ratio;
					dp_brush[zz] = dp_image[zz];
				}
			}
		}
	}
}

/********************************************************************/

void brush_smudge_circle::to_subpixel_from_pixel(double d_x1, double d_y1, double d_x2, double d_y2)
{
	double d_subsize;
	int32_t i32_xsize;
	double *dp_image,
		*dp_save;
	double d_x, d_y,
		d_x1floor, d_y1floor,
		d_xsave, d_ysave;
	int32_t zz;

	d_subsize = 1.0 / this->_i32_subpixel_divide;

	/* 保存(復元)位置 */
	d_x1floor = floor(d_x1 + d_subsize / 2.0);
	d_y1floor = floor(d_y1 + d_subsize / 2.0);
	i32_xsize = (int32_t)floor(d_x2 - d_subsize / 2.0) - (int32_t)d_x1floor + 1;

	dp_image = this->_dp_subpixel_image;
	dp_save = this->_dp_pixel_image;

	/* d_x,d_yでループ、d_xsave,d_ysaveで位置 */
	for (d_y = d_y1 + d_subsize / 2.0; d_y < d_y2; d_y += d_subsize) {
		for (d_x = d_x1 + d_subsize / 2.0; d_x < d_x2; d_x += d_subsize,
			dp_image += 5) {

			d_xsave = d_x - d_x1floor;
			d_ysave = d_y - d_y1floor;

			assert(0 <= (int32_t)d_xsave);
			assert(0 <= (int32_t)d_ysave);
			assert((int32_t)d_xsave < (this->_i32_size_by_pixel + 1));
			assert((int32_t)d_ysave < (this->_i32_size_by_pixel + 1));

			for (zz = 0; zz < 5; ++zz) {
				dp_image[zz] = dp_save[(int32_t)d_ysave * i32_xsize * 5 +
									   (int32_t)d_xsave * 5 +
									   zz];
			}
		}
	}
}

/* 切取ったサブピクセル画像を、(元画像へ)保存のためピクセルサイズに縮小
	例えば、元画像の2 x 2ピクセルの任意の位置から切り取ってあるなら、
	復元に3 x 3ピクセルサイズ必要となる。
	+-------+-------+-------+
	|     +-|-----+-|-----+ |
	|   + | |   + | |   + | |
	|     | |     | |     | |
	+-------+-------+-------+
	|     +-|-----+ |-----+ |
	|   + | |   + | |   + | |
	|     | |     | |     | |
	+-------+-------+-------+
	|     +-|-----+-|-----+ |
	|   +   |   +   |   +   |
	|       |       |       |
	+-------+-------+-------+
*/
void brush_smudge_circle::to_pixel_from_subpixel(double d_x1, double d_y1, double d_x2, double d_y2)
{
	double d_subsize;
	int32_t i32_xsize;
	double *dp_image,
		*dp_save;
	double d_x, d_y,
		d_x1floor, d_y1floor,
		d_xsave, d_ysave;
	int32_t zz;
	int32_t ii, jj;

	d_subsize = 1.0 / this->_i32_subpixel_divide;

	/* 初期化 */
	dp_save = this->_dp_pixel_image;
	for (ii = 0; ii < this->_i32_size_by_pixel + 1; ++ii) {
		for (jj = 0; jj < this->_i32_size_by_pixel + 1; ++jj, dp_save += 5) {
			for (zz = 0; zz < 5; ++zz) {
				dp_save[zz] = 0.0;
			}
		}
	}

	/* 保存(復元)位置 */
	d_x1floor = floor(d_x1 + d_subsize / 2.0);
	d_y1floor = floor(d_y1 + d_subsize / 2.0);
	i32_xsize = (int32_t)floor(d_x2 - d_subsize / 2.0) - (int32_t)d_x1floor + 1;

	dp_image = this->_dp_subpixel_image;
	dp_save = this->_dp_pixel_image;

	/* d_x,d_yでループ、d_xsave,d_ysaveで位置 */
	for (d_y = d_y1 + d_subsize / 2.0; d_y < d_y2; d_y += d_subsize) {
		for (d_x = d_x1 + d_subsize / 2.0; d_x < d_x2; d_x += d_subsize,
			dp_image += 5) {

			d_xsave = d_x - d_x1floor;
			d_ysave = d_y - d_y1floor;

			assert(0 <= (int32_t)d_xsave);
			assert(0 <= (int32_t)d_ysave);
			assert((int32_t)d_xsave < (this->_i32_size_by_pixel + 1));
			assert((int32_t)d_ysave < (this->_i32_size_by_pixel + 1));

			for (zz = 0; zz < 5; ++zz) {
				dp_save[(int32_t)d_ysave * i32_xsize * 5 +
						(int32_t)d_xsave * 5 +
						zz] += dp_image[zz];
			}
		}
	}

	/* ピクセル値に */
	dp_save = this->_dp_pixel_image;
	for (ii = 0; ii < this->_i32_size_by_pixel + 1; ++ii) {
		for (jj = 0; jj < this->_i32_size_by_pixel + 1; ++jj, dp_save += 5) {
			for (zz = 0; zz < 5; ++zz) {
				dp_save[zz] /= this->_i32_subpixel_divide *
							   this->_i32_subpixel_divide;
			}
		}
	}
}

#ifndef __calculator_geometry_h__
#define __calculator_geometry_h__

class calculator_geometry
{
public:
	double get_d_radian(double d_xv, double d_yv);
	double get_d_radian_by_2_vector(double d_xv1, double d_yv1, double d_xv2, double d_yv2);
	void get_dd_rotate(double d_xp1, double d_yp1, double d_radian, double *dp_xp2, double *dp_yp2);
	void get_dd_mirror(double d_xp1, double d_yp1, double d_mirror_xc, double d_mirror_yc, double d_mirror_radian, double *dp_xp2, double *dp_yp2);
	void get_dd_rotate_by_pos(double d_xp1, double d_yp1, double d_xpos, double d_ypos, double d_radian, double *dp_xp2, double *dp_yp2);

private:
};

#endif			  /* !__calculator_geometry_h__ */
#include <math.h> /* sin(),cos(),atan(), M_PI */

#include "igs_line_blur.h" // "pri.h" "calculator_geometry.h"

/* WindowsではM_PIが見つからない */
#if defined _WIN32
#define M_PI 3.14159265358979323846
#endif

/* ベクトルの角度を"0"から"2*PI"で返す */
double calculator_geometry::get_d_radian(double d_xv, double d_yv)
{
	double d_radian;

	/* ゼロエラー */
	if ((0.0 == d_xv) && (0.0 == d_yv)) {
		pri_funct_err_bttvr(
			"Warning : calculator_geometry::get_d_radian(d_xv,d_yv is zero).");
		return 0.0;
	}
	/* 第1象限 (0 <= angle < 90)*/
	else if ((0.0 < d_xv) && (0.0 <= d_yv)) {
		d_radian = atan(d_yv / d_xv);
	}
	/* 第2象限 (第1象限に置き換えて... 0 <= angle < 90) */
	else if ((d_xv <= 0.0) && (0.0 < d_yv)) {
		d_radian = atan(-d_xv / d_yv) + M_PI / 2.0;
	}
	/* 第3象限 (第1象限に置き換えて... 0 <= angle < 90) */
	else if ((d_xv < 0.0) && (d_yv <= 0.0)) {
		d_radian = atan(d_yv / d_xv) + M_PI;
	}
	/* 第4象限 (第1象限に置き換えて... 0 <= angle < 90) */
	else if ((0.0 <= d_xv) && (d_yv < 0.0)) {
		d_radian = atan(d_xv / -d_yv) + M_PI + M_PI / 2.0;
	}
	return d_radian;
}

/* 2つのベクトルのなす角度を反時計回りに調べる */
double calculator_geometry::get_d_radian_by_2_vector(double d_xv1, double d_yv1, double d_xv2, double d_yv2)
{
	double d_radian_start, d_radian_end;

	if ((0.0 == d_xv1) && (0.0 == d_yv1)) {
		pri_funct_err_bttvr(
			"Warning : calculator_geometry::get_d_radian_by_2_vector(d_xv1,d_yv1 is zero).");
	}
	if ((0.0 == d_xv2) && (0.0 == d_yv2)) {
		pri_funct_err_bttvr(
			"Warning : calculator_geometry::get_d_radian_by_2_vector(d_xv2,d_yv2 is zero).");
	}

	d_radian_start = this->get_d_radian(d_xv1, d_yv1);
	d_radian_end = this->get_d_radian(d_xv2, d_yv2);
	if (d_radian_end < d_radian_start) {
		d_radian_end += M_PI * 2.0;
	}
	return d_radian_end - d_radian_start;
}

/* x,y座標の回転 */
void calculator_geometry::get_dd_rotate(double d_xp1, double d_yp1, double d_radian, double *dp_xp2, double *dp_yp2)
{
	*dp_xp2 = d_xp1 * cos(d_radian) - d_yp1 * sin(d_radian);
	*dp_yp2 = d_xp1 * sin(d_radian) + d_yp1 * cos(d_radian);
}

/* x,y座標を線(鏡)対象に移動する */
void calculator_geometry::get_dd_mirror(double d_xp1, double d_yp1, double d_mirror_xpos, double d_mirror_ypos, double d_mirror_radian, double *dp_xp2, double *dp_yp2)
{
	d_xp1 -= d_mirror_xpos;
	d_yp1 -= d_mirror_ypos;
	this->get_dd_rotate(d_xp1, d_yp1, -d_mirror_radian, &d_xp1, &d_yp1);
	this->get_dd_rotate(d_xp1, -d_yp1, d_mirror_radian, &d_xp1, &d_yp1);
	*dp_xp2 = d_xp1 + d_mirror_xpos;
	*dp_yp2 = d_yp1 + d_mirror_ypos;
}

/* 指定座標を中心に回転する */
void calculator_geometry::get_dd_rotate_by_pos(double d_xp1, double d_yp1, double d_xpos, double d_ypos, double d_radian, double *dp_xp2, double *dp_yp2)
{
	d_xp1 -= d_xpos;
	d_yp1 -= d_ypos;
	this->get_dd_rotate(d_xp1, d_yp1, d_radian, &d_xp1, &d_yp1);
	*dp_xp2 = d_xp1 + d_xpos;
	*dp_yp2 = d_yp1 + d_ypos;
}

#ifndef _pixel_point_node_h_
#define _pixel_point_node_h_

#include <stdio.h>

/* Windowsではstdint.hが見つからない */
#if defined _WIN32
typedef int int32_t;
typedef unsigned short uint16_t;
#else
#include <stdint.h> /* for int32_t, uint16_t */
#endif

#include "igs_line_blur.h" // "list_node.h"

#ifndef OK
#define OK (0)
#endif
#ifndef NG
#define NG (-1)
#endif

#ifndef LINK_NEAR_COUNT
#define LINK_NEAR_COUNT 4
#endif

/* x,yポイント座標のリストノード、画素連結、及び、線分連結、機能付き */
class pixel_point_node : public list_node
{
public:
	pixel_point_node()
	{
		int32_t ii;

		this->_i32_xp = 0;
		this->_i32_yp = 0;
		this->_d_xp_tgt = 0.0;
		this->_d_yp_tgt = 0.0;

		for (ii = 0; ii < LINK_NEAR_COUNT; ++ii) {
			this->_clpa_link_near[ii] = NULL;
		}
		this->_clp_previous_point = NULL;
		this->_clp_next_point = NULL;
	}

	/* 位置のセットと、ゲット */
	void set_i32_xp(int32_t i32) { this->_i32_xp = i32; }
	void set_i32_yp(int32_t i32) { this->_i32_yp = i32; }
	void set_d_xp_tgt(double dd) { this->_d_xp_tgt = dd; }
	void set_d_yp_tgt(double dd) { this->_d_yp_tgt = dd; }

	int32_t get_i32_xp(void) { return this->_i32_xp; }
	int32_t get_i32_yp(void) { return this->_i32_yp; }
	double get_d_xp_tgt(void) { return this->_d_xp_tgt; }
	double get_d_yp_tgt(void) { return this->_d_yp_tgt; }

	/* 第１のリンク : list_nodeの継承による全データ検索用リスト */
	/* void set_clp_next( pixel_point_node *clp ); */
	/* void set_clp_previous( pixel_point_node *clp ); */
	/* pixel_point_node *get_clp_next( void ); */
	/* pixel_point_node *get_clp_previous( void ); */

	/* 第２のリンク : 画像上の(上下左右斜め)ピクセル連結 */
	int link_near(pixel_point_node *clp_);
	pixel_point_node *get_clp_link_near(int32_t i32);

	/* 第３のリンク : 線分を表すポイントリスト */
	void set_clp_next_point(pixel_point_node *clp)
	{
		this->_clp_next_point = clp;
	}
	void set_clp_previous_point(pixel_point_node *clp)
	{
		this->_clp_previous_point = clp;
	}
	pixel_point_node *get_clp_next_point(void)
	{
		return this->_clp_next_point;
	}
	pixel_point_node *get_clp_previous_point(void)
	{
		return this->_clp_previous_point;
	}

	/* 連結情報表示 */
	void print_xy_around(void);

private:
	int32_t _i32_xp, _i32_yp;	/* source(元)画像上の位置 */
	double _d_xp_tgt, _d_yp_tgt; /* target(最終)画像上の位置 */

	pixel_point_node *_clpa_link_near[LINK_NEAR_COUNT];
	pixel_point_node *_clp_previous_point,
		*_clp_next_point;
};

#endif				/* !_pixel_point_node_h_ */
#include <assert.h> /* assert() */

#include "igs_line_blur.h" // "pri.h" "pixel_point_node.h"

int pixel_point_node::link_near(pixel_point_node *clp_)
{
	int32_t ii;

	for (ii = 0; ii < LINK_NEAR_COUNT; ++ii) {
		if (NULL == this->_clpa_link_near[ii]) {
			this->_clpa_link_near[ii] = clp_;
			return ii;
		}
	}

	pri_funct_err_bttvr("Error : no link_near point, over %d.",
						LINK_NEAR_COUNT - 1);
	pri_funct_err_bttvr("this   x %d y %d",
						this->get_i32_xp(),
						this->get_i32_yp());
	for (ii = 0; ii < LINK_NEAR_COUNT; ++ii) {
		assert(NULL != this->_clpa_link_near[ii]);
		pri_funct_err_bttvr(
			"link_near %d x %d y %d",
			ii,
			this->_clpa_link_near[ii]->get_i32_xp(),
			this->_clpa_link_near[ii]->get_i32_yp());
	}

	return NG;
}

pixel_point_node *pixel_point_node::get_clp_link_near(int32_t i32)
{
	if (LINK_NEAR_COUNT <= i32)
		return NULL;
	return this->_clpa_link_near[i32];
}

void pixel_point_node::print_xy_around(void)
{
	int32_t ii;

	pri_funct_msg_ttvr(
		"pixel_point_node::debug_print_xy_around() : self address <0x%lx>",
		this);
	pri_funct_msg_ttvr(
		" self    x %d y %d",
		this->get_i32_xp(),
		this->get_i32_yp());

	for (ii = 0; ii < LINK_NEAR_COUNT; ++ii) {
		if (NULL != this->get_clp_link_near(ii)) {
			pri_funct_msg_ttvr(
				" link_near %d  x %d y %d",
				ii,
				this->get_clp_link_near(ii)->get_i32_xp(),
				this->get_clp_link_near(ii)->get_i32_yp());
		} else {
			pri_funct_msg_ttvr(
				" link_near %ld is not exist(NULL).", ii);
		}
	}
}

#ifndef _pixel_point_root_h_
#define _pixel_point_root_h_

#include "igs_line_blur.h" // "list_root.h" "pixel_point_node.h"

#ifndef ON
#define ON (1)
#endif
#ifndef OFF
#define OFF (0)
#endif

class pixel_point_root : public list_root
{
public:
	pixel_point_root(void)
	{
		this->_i_mv_sw = OFF;
		this->_i_cv_sw = OFF;
		this->_i_pv_sw = OFF;
	}
	~pixel_point_root(void)
	{
		this->mem_free();
	}

	void set_i_mv_sw(int sw) { this->_i_mv_sw = sw; }
	void set_i_cv_sw(int sw) { this->_i_cv_sw = sw; }
	void set_i_pv_sw(int sw) { this->_i_pv_sw = sw; }

	/* 16Bits 1channel画像からゼロ以上のピクセルをプロットする */
	int alloc_mem_and_list_node(int32_t i32_xs, int32_t i32_ys, uint16_t *ui16p_src);

	/* リストの追加 */
	pixel_point_node *append(pixel_point_node *clp_previous);

	/* for debug */
	int save(char *cp_fname);

	void mem_free(void);

private:
	int _i_mv_sw,
		_i_cv_sw,
		_i_pv_sw;

	/* リストの削除 */
	void _remove(pixel_point_node *clp_target);
};

#endif				/* !_pixel_point_root_h_ */
#include <assert.h> /* assert() */

#include "igs_line_blur.h" // "pri.h" "pixel_point_root.h"

pixel_point_node *pixel_point_root::append(pixel_point_node *clp_previous)
{
	pixel_point_node *clp_new;

	clp_new = new pixel_point_node;
	if (NULL == clp_new) {
		pri_funct_err_bttvr("Error : 'new pixel_point_node' returns NULL.");
		return NULL;
	}
	clp_new = (pixel_point_node *)this->push(clp_previous, clp_new);

	return clp_new;
}

void pixel_point_root::_remove(pixel_point_node *clp_old)
{
	assert(NULL != clp_old); /* あってはならないプログラムバグのチェック */
	this->pop(clp_old);
	delete clp_old;
}

void pixel_point_root::mem_free(void)
{
	pixel_point_node *clp_;
	int32_t ii;

	if (NULL != this->get_clp_last()) {
		if (ON == this->_i_mv_sw) {
			pri_funct_msg_ttvr("pixel_point_root::mem_free()");
		}
		ii = 0;
		while (NULL != (clp_ = (pixel_point_node *)this->get_clp_last())) {
			this->_remove(clp_);
			++ii;
		}
		if (ON == this->_i_pv_sw) {
			pri_funct_msg_ttvr("free point node %d", ii);
		}
	}
}

int pixel_point_root::alloc_mem_and_list_node(int32_t i32_xs, int32_t i32_ys, uint16_t *ui16p_src)
{
	int32_t xx, yy;
	pixel_point_node *clp_pp;

	/* 処理実行表示 */
	if (ON == this->_i_mv_sw) {
		pri_funct_msg_ttvr("pixel_point_root::alloc_mem_and_list_node()");
	}

	/* カウントダウン表示始め */
	if (ON == this->_i_cv_sw) {
		pri_funct_cv_start(i32_ys);
	}

	clp_pp = NULL;
	for (yy = 0; yy < i32_ys; ++yy) {
		/* カウントダウン表示中 */
		if (ON == this->_i_cv_sw) {
			pri_funct_cv_run(yy);
		}

		for (xx = 0; xx < i32_xs; ++xx, ++ui16p_src) {
			if (0 < (*ui16p_src)) {
				clp_pp = this->append(clp_pp);
				if (NULL == clp_pp) {
					pri_funct_err_bttvr(
						"Error : this->append(clp_pp) returns NULL.");
					return NG;
				}
				clp_pp->set_i32_xp(xx);
				clp_pp->set_i32_yp(yy);
			}
		}
	}
	/* カウントダウン表示終了 */
	if (ON == this->_i_cv_sw) {
		pri_funct_cv_end();
	}

	/* 処理結果表示 */
	if (ON == this->_i_pv_sw) {
		pri_funct_msg_ttvr("alloc and list %d points", this->get_i32_count());
	}

	return OK;
}

int pixel_point_root::save(char *cp_fname)
{
	FILE *fp;
	int32_t ii;
	pixel_point_node *clp_pp_node;

	/* あってはならないプログラムバグのチェック */
	assert(NULL != cp_fname);

	/* ファイル開く */
	fp = fopen(cp_fname, "w");
	if (NULL == fp) {
		pri_funct_err_bttvr("Error : fopen(%s) returns NULL.", cp_fname);
		return NG;
	}

	/* データ数保存 */
	if (fprintf(fp, "# count %d\n", this->get_i32_count()) < 0) {
		pri_funct_err_bttvr("Error : fprintf(count) returns minus.");
		return NG;
	}

	/* データ保存 */
	for (clp_pp_node = (pixel_point_node *)this->get_clp_first(),
		ii = 0L;
		 (NULL != clp_pp_node) && (ii < this->get_i32_count());
		 clp_pp_node = (pixel_point_node *)clp_pp_node->get_clp_next(),
		++ii) {
		if (fprintf(fp, "%d %d\n",
					clp_pp_node->get_i32_xp(),
					clp_pp_node->get_i32_yp()) < 0) {
			pri_funct_err_bttvr(
				"Error : fprintf(count xp yp) returns minus.");
			return NG;
		}
	}
	/* ファイル閉じる */
	fclose(fp);

	return OK;
}

#ifndef _pixel_line_node_h_
#define _pixel_line_node_h_

#include "igs_line_blur.h" // "calculator_geometry.h" "list_node.h" "pixel_point_root.h"

#ifndef OK
#define OK (0)
#endif
#ifndef NG
#define NG (-1)
#endif

class pixel_line_node : public list_node
{
public:
	pixel_line_node()
	{
		this->_i32_point_count = 0;
		this->_d_same_way_radian_one = 0.0;
		this->_d_same_way_radian_another = 0.0;
		this->_d_bbox_x_min = 0.0;
		this->_d_bbox_x_max = 0.0;
		this->_d_bbox_y_min = 0.0;
		this->_d_bbox_y_max = 0.0;

		this->_clpa_link[0] = NULL;
		this->_clpa_link[1] = NULL;
		this->_clpa_link[2] = NULL;
		this->_clpa_link[3] = NULL;
		this->_clpa_link[4] = NULL;
	}

	/* 値の参照 */
	int32_t get_i32_point_count(void) { return this->_i32_point_count; }
	void set_d_same_way_radian_one(double dd)
	{
		this->_d_same_way_radian_one = dd;
	}
	void set_d_same_way_radian_another(double dd)
	{
		this->_d_same_way_radian_another = dd;
	}
	double get_d_same_way_radian_one(void)
	{
		return this->_d_same_way_radian_one;
	}
	double get_d_same_way_radian_another(void)
	{
		return this->_d_same_way_radian_another;
	}
	double get_d_bbox_x_min(void) { return this->_d_bbox_x_min; }
	double get_d_bbox_x_max(void) { return this->_d_bbox_x_max; }
	double get_d_bbox_y_min(void) { return this->_d_bbox_y_min; }
	double get_d_bbox_y_max(void) { return this->_d_bbox_y_max; }

	/* 対象(pixel_point_node)へのリンク機能 */
	void link_one(pixel_point_node *clp) { this->_clpa_link[0] = clp; }
	void link_another(pixel_point_node *clp) { this->_clpa_link[1] = clp; }

	void link_middle(pixel_point_node *clp) { this->_clpa_link[2] = clp; }
	void link_one_expand(pixel_point_node *clp) { this->_clpa_link[3] = clp; }
	void link_another_expand(pixel_point_node *clp) { this->_clpa_link[4] = clp; }

	pixel_point_node *get_clp_link_one(void) { return this->_clpa_link[0]; }
	pixel_point_node *get_clp_link_another(void) { return this->_clpa_link[1]; }

	pixel_point_node *get_clp_link_middle(void) { return this->_clpa_link[2]; }
	pixel_point_node *get_clp_link_one_expand(void) { return this->_clpa_link[3]; }
	pixel_point_node *get_clp_link_another_expand(void) { return this->_clpa_link[4]; }

	void set_middle(void);
	void int2double_body(void);
	void smooth_body(int32_t i32_smooth_retry);
	void smooth_expand(int32_t i32_smooth_retry);
	void link_line(pixel_point_node *clp_crnt, pixel_point_node *clp_next, int32_t i32_count);
	void get_near_point(double d_xp, double d_yp, int32_t *i32p_pos, pixel_point_node **clpp_point, double *dp_length);

	pixel_point_node *get_next_point_by_count(pixel_point_node *clp_point, int32_t i32_count);
	pixel_point_node *get_prev_point_by_count(pixel_point_node *clp_point, int32_t i32_count);
	void set_bbox(void);

	int save_line(FILE *fp);
	int save_one_point(FILE *fp);
	int save_middle_point(FILE *fp);
	int save_another_point(FILE *fp);

	int expand_line(pixel_point_root *clp_pixel_point_root);

	int save_expand_line(FILE *fp);
	int save_expand_vector(FILE *fp);
	int save_one_expand_point(FILE *fp);
	int save_another_expand_point(FILE *fp);

private:
	int32_t _i32_point_count;
	double _d_same_way_radian_one,
		_d_same_way_radian_another;
	double _d_bbox_x_min, _d_bbox_x_max,
		_d_bbox_y_min, _d_bbox_y_max;

	pixel_point_node *_clpa_link[5];
	calculator_geometry _cl_cal_geom;

	int _expand_line_from_one(pixel_point_root *clp_pp_root, int32_t i32_body_point_count, pixel_point_node *clp_one, pixel_point_node *clp_another, double d_radian);
	int _expand_line_from_another(pixel_point_root *clp_pp_root, int32_t i32_body_point_count, pixel_point_node *clp_one, pixel_point_node *clp_another, double d_radian);
	void _smooth_point(double d_x1, double d_y1, double d_x2, double d_y2, double d_x3, double d_y3, double *dp_x, double *dp_y);

	void _get_link_line_selecter_vector(pixel_point_node *clp_crnt, pixel_point_node *clp_next, double *dp_xv, double *dp_yv, int32_t i32_count);
	pixel_point_node *_get_link_line_selecter(double d_xv, double d_yv, pixel_point_node *clp_crnt, int32_t i32_count);
};

#endif /* !_pixel_line_node_h_ */
#include <assert.h>

#include "igs_line_blur.h" // "pri.h" "pixel_line_node.h"

int pixel_line_node::_expand_line_from_one(pixel_point_root *clp_pp_root, int32_t i32_body_point_count, pixel_point_node *clp_one, pixel_point_node *clp_another, double d_radian)
{
	pixel_point_node *clp_last, *clp_before, *clp_exist;
	double d_xp, d_yp;
	int32_t ii;

	clp_before = clp_one;

	/* ライン開始点(clp_one)から後ろへたどる、
	   始まりは開始点でなくその一つ後ろから */
	for (ii = 1, clp_exist = clp_one->get_clp_next_point();
		 NULL != clp_exist;
		 ++ii, clp_exist = clp_exist->get_clp_next_point()) {
		/* 偽の場合、たぶん無限ループ */
		assert(ii < i32_body_point_count);

		/* 点node生成とデータリンク(free用) */
		clp_last = clp_pp_root->append(
			(pixel_point_node *)(clp_pp_root->get_clp_last()));
		if (NULL == clp_last) {
			pri_funct_err_bttvr(
				"Error : ii %d : this->append(clp_last) returns NULL.",
				ii);
			return NG;
		}

		/* ラインとしてのリンク */
		clp_before->set_clp_previous_point(clp_last);
		clp_last->set_clp_next_point(clp_before);
		this->link_one_expand(clp_last);
		++this->_i32_point_count;

		/* 位置を鏡面反転 */
		this->_cl_cal_geom.get_dd_mirror(
			clp_exist->get_d_xp_tgt(),
			clp_exist->get_d_yp_tgt(),
			clp_one->get_d_xp_tgt(),
			clp_one->get_d_yp_tgt(),
			d_radian, &d_xp, &d_yp);

		/* 位置をセット */
		clp_last->set_d_xp_tgt(d_xp);
		clp_last->set_d_yp_tgt(d_yp);

		/* 接続のため前の点を記憶 */
		clp_before = clp_last;

		/* 終了 */
		if (clp_another == clp_exist) {
			break;
		}
	}
	return OK;
}
int pixel_line_node::_expand_line_from_another(pixel_point_root *clp_pp_root, int32_t i32_body_point_count, pixel_point_node *clp_one, pixel_point_node *clp_another, double d_radian)
{
	pixel_point_node *clp_last, *clp_before, *clp_exist;
	double d_xp, d_yp;
	int32_t ii;

	clp_before = clp_another;

	/* ライン最終点(clp_another)から前へたどる、
	   始まりは最終点でなくその一つ前から */
	for (ii = 1, clp_exist = clp_another->get_clp_previous_point();
		 NULL != clp_exist;
		 ++ii, clp_exist = clp_exist->get_clp_previous_point()) {
		/* 偽の場合、たぶん無限ループ */
		assert(ii < i32_body_point_count);

		/* 点node生成とデータリンク(free用) */
		clp_last = clp_pp_root->append(
			(pixel_point_node *)(clp_pp_root->get_clp_last()));
		if (NULL == clp_last) {
			pri_funct_err_bttvr(
				"Error : ii %d : this->append(clp_last) returns NULL.",
				ii);
			return NG;
		}

		/* ラインとしてのリンク */
		clp_before->set_clp_next_point(clp_last);
		clp_last->set_clp_previous_point(clp_before);
		this->link_another_expand(clp_last);
		++this->_i32_point_count;

		/* 位置を鏡面反転 */
		this->_cl_cal_geom.get_dd_mirror(
			clp_exist->get_d_xp_tgt(),
			clp_exist->get_d_yp_tgt(),
			clp_another->get_d_xp_tgt(),
			clp_another->get_d_yp_tgt(),
			d_radian, &d_xp, &d_yp);

		/* 位置をセット */
		clp_last->set_d_xp_tgt(d_xp);
		clp_last->set_d_yp_tgt(d_yp);

		/* 接続のため前の点を記憶 */
		clp_before = clp_last;

		/* 終了 */
		if (clp_one == clp_exist) {
			break;
		}
	}
	return OK;
}

/********************************************************************/

int pixel_line_node::expand_line(pixel_point_root *clp_pixel_point_root)
{
	pixel_point_node *clp_one, *clp_middle, *clp_another;
	double d_radian, d_radian_one, d_radian_another;
	int32_t i32_body_point_count;

	/* 3点より少ない、短すぎる線分は伸ばすことができない */
	if (this->_i32_point_count < 3) {
		return OK;
	}

	/* 両端点と中点を得る */
	clp_one = this->get_clp_link_one();
	clp_another = this->get_clp_link_another();
	clp_middle = this->get_clp_link_middle();

	assert(clp_one != clp_another);
	assert(clp_another != clp_middle);
	assert(clp_middle != clp_one);

	/* 2つのベクトルの角度を反時計回りに求める */
	assert((clp_one->get_i32_xp() != clp_middle->get_i32_xp()) ||
		   (clp_one->get_i32_yp() != clp_middle->get_i32_yp()));
	assert((clp_another->get_i32_xp() != clp_middle->get_i32_xp()) ||
		   (clp_another->get_i32_yp() != clp_middle->get_i32_yp()));
	d_radian = this->_cl_cal_geom.get_d_radian_by_2_vector(
		clp_one->get_i32_xp() - clp_middle->get_i32_xp(),
		clp_one->get_i32_yp() - clp_middle->get_i32_yp(),
		clp_another->get_i32_xp() - clp_middle->get_i32_xp(),
		clp_another->get_i32_yp() - clp_middle->get_i32_yp());

	/* 半分にして */
	d_radian /= 2.0;

	/* 始点から中間点へのベクトルの角度 */
	d_radian_one = this->_cl_cal_geom.get_d_radian(
		clp_middle->get_i32_xp() - clp_one->get_i32_xp(),
		clp_middle->get_i32_yp() - clp_one->get_i32_yp());

	/* 始点における、線分の線対象な線の角度 */
	d_radian_one -= d_radian;

	/* 終点から中間点へのベクトルの角度 */
	d_radian_another = this->_cl_cal_geom.get_d_radian(
		clp_middle->get_i32_xp() - clp_another->get_i32_xp(),
		clp_middle->get_i32_yp() - clp_another->get_i32_yp());

	/* 終点における、線分の線対象な線の角度 */
	d_radian_another += d_radian;

	/* 始点が端点ならば先へ伸ばす */
	i32_body_point_count = this->_i32_point_count;
	if (NULL == clp_one->get_clp_link_near(1)) { /* 2点目がないなら端点 */
		if (OK != this->_expand_line_from_one(clp_pixel_point_root, i32_body_point_count, this->get_clp_link_one(), this->get_clp_link_another(), d_radian_one)) {
			pri_funct_err_bttvr(
				"Error : this->_expand_line_from_one(-) returns NULL.");
			return NG;
		}
	}

	/* 終点が端点ならば先へ伸ばす */
	if (NULL == clp_another->get_clp_link_near(1)) { /* 2点目がないなら端点 */
		if (OK != this->_expand_line_from_another(clp_pixel_point_root, i32_body_point_count, this->get_clp_link_one(), this->get_clp_link_another(), d_radian_another)) {
			pri_funct_err_bttvr(
				"Error : this->_expand_line_from_another(-) returns NULL.");
			return NG;
		}
	}

	return OK;
}

#include <math.h>   /* sqrt() */
#include <assert.h> /* assert() */

#include "igs_line_blur.h" // "pixel_line_node.h"

void pixel_line_node::get_near_point(double d_xp, double d_yp, int32_t *i32p_pos, pixel_point_node **clpp_point, double *dp_length)
{
	pixel_point_node *clp_loop;
	int32_t ii;
	double d_length;

	assert(NULL != this->get_clp_link_one_expand());

	/* 各ポイントを探索 */
	*dp_length = 1000.0;
	clp_loop = this->get_clp_link_one_expand();
	for (ii = 0; NULL != clp_loop; ++ii,
		clp_loop = clp_loop->get_clp_next_point()) {
		/* 偽の場合、たぶん無限ループ */
		assert(ii < this->_i32_point_count);

		/* 距離 */
		d_length = sqrt(
			(clp_loop->get_d_xp_tgt() - d_xp) *
				(clp_loop->get_d_xp_tgt() - d_xp) +
			(clp_loop->get_d_yp_tgt() - d_yp) *
				(clp_loop->get_d_yp_tgt() - d_yp));
		/* 近いものならセットする */
		if (d_length < (*dp_length)) {
			*i32p_pos = ii;
			*clpp_point = clp_loop;
			*dp_length = d_length;
		}
	}
}

#include <assert.h> /* assert() */

#include "igs_line_blur.h" // "pixel_line_node.h"

void pixel_line_node::int2double_body(void)
{
	pixel_point_node *clp_1;
	int32_t ii;

	/* 浮動小数化 */
	clp_1 = this->get_clp_link_one();
	for (ii = 0; NULL != clp_1;
		 clp_1 = clp_1->get_clp_next_point(), ++ii) {
		/* 偽の場合、たぶん無限ループ */
		assert(ii < this->_i32_point_count);

		/* x,y座標をint(整数値)からdouble(浮動小数点)に変換し格納 */
		clp_1->set_d_xp_tgt((double)clp_1->get_i32_xp());
		clp_1->set_d_yp_tgt((double)clp_1->get_i32_yp());
	}
}

#include <math.h>   /* M_PI */
#include <assert.h> /* assert() */

#include "igs_line_blur.h" // "pri.h" "pixel_line_node.h"

/* WindowsではM_PIが見つからない */
#if defined _WIN32
#define M_PI 3.14159265358979323846
#endif

void pixel_line_node::_get_link_line_selecter_vector(pixel_point_node *clp_crnt, pixel_point_node *clp_next, double *dp_xv, double *dp_yv, int32_t i32_count)
{
	pixel_point_node *clp_start;
	int32_t ii;

	clp_start = clp_crnt;

	/* 次への端点を求める */
	for (ii = 0; (NULL != clp_next) && (ii < i32_count); ++ii) {
		/* 次の点が線分としてリンク済みのとき */
		if ((NULL != clp_next->get_clp_next_point()) || (NULL != clp_next->get_clp_previous_point())) {
			break;
		}

		/* もう一方の端点に到着し終端(リンク2点目がないので) */
		if ((NULL == clp_next->get_clp_link_near(1))) {
			clp_crnt = clp_next;
			break;
		}

		/* 分岐点(リンク3点目があるので) */
		if (NULL != clp_next->get_clp_link_near(2)) {
			clp_crnt = clp_next;
			break;
		}
		/* 次へのリンクをたどる */
		if (clp_crnt == clp_next->get_clp_link_near(0)) {
			clp_crnt = clp_next;
			clp_next = clp_next->get_clp_link_near(1);
		} else if (clp_crnt == clp_next->get_clp_link_near(1)) {
			clp_crnt = clp_next;
			clp_next = clp_next->get_clp_link_near(0);
		}
		/* リンクが壊れている */
		else {
			pri_funct_err_bttvr("Error : bad link");
			assert(0);
		}
	}
	*dp_xv = clp_crnt->get_i32_xp() - clp_start->get_i32_xp();
	*dp_yv = clp_crnt->get_i32_yp() - clp_start->get_i32_yp();
}

pixel_point_node *pixel_line_node::_get_link_line_selecter(double d_xv, double d_yv, pixel_point_node *clp_crnt, int32_t i32_count)
{
	int32_t ii, i32_pos;
	double da_xv[LINK_NEAR_COUNT],
		da_yv[LINK_NEAR_COUNT],
		da_radian[LINK_NEAR_COUNT],
		d_radian;

	/* あってはならないプログラムバグのチェック */
	assert((0.0 != d_xv) || (0.0 != d_yv));
	assert(NULL != clp_crnt);

	/* 各分岐の方向ベクトルを得る */
	for (ii = 0; ii < LINK_NEAR_COUNT; ++ii) {
		if (NULL == clp_crnt->get_clp_link_near(ii)) {
			break;
		}
		this->_get_link_line_selecter_vector(
			clp_crnt, clp_crnt->get_clp_link_near(ii),
			&da_xv[ii], &da_yv[ii], i32_count);
	}

	/* 各分岐線の方向ベクトルから、手前の線との角度を求める */
	for (ii = 0; ii < LINK_NEAR_COUNT; ++ii) {
		if (NULL == clp_crnt->get_clp_link_near(ii)) {
			break;
		}
		if ((0.0 != da_xv[ii]) || (0.0 != da_yv[ii])) {
			da_radian[ii] = this->_cl_cal_geom.get_d_radian_by_2_vector(
				d_xv, d_yv, da_xv[ii], da_yv[ii]);
			if (M_PI < da_radian[ii]) {
				da_radian[ii] = M_PI - (da_radian[ii] - M_PI);
			}
		} else {
			da_radian[ii] = M_PI;
		}
	}

	/* もっとも直線と見なされる(角度の小さい)分岐を返す */
	d_radian = M_PI * 2.0;
	for (ii = 0; ii < LINK_NEAR_COUNT; ++ii) {
		if (NULL == clp_crnt->get_clp_link_near(ii)) {
			break;
		}
		if (da_radian[ii] < d_radian) {
			i32_pos = ii;
			d_radian = da_radian[ii];
		}
	}
	/* もっとも直線と見なされる(角度の小さい)といっても
	90度以上も曲がってたらだめ */
	if ((M_PI / 2.0 <= d_radian) && (d_radian <= M_PI * 3.0 / 2.0))
		return NULL;

	return clp_crnt->get_clp_link_near(i32_pos);
}

void pixel_line_node::link_line(pixel_point_node *clp_crnt, pixel_point_node *clp_next, int32_t i32_count)
{
	int32_t ii;
	pixel_point_node *clp_tmp1, *clp_tmp2;
	double d_xv, d_yv;

	/* あってはならないプログラムバグのチェック */
	assert(NULL != clp_crnt);
	assert(NULL != clp_next);

	/* スタート点をセット */
	this->link_one(clp_crnt);

	++(this->_i32_point_count);

	clp_tmp1 = clp_tmp2 = this->get_clp_link_one();
	d_xv = d_yv = 0.0;

	for (ii = 0; (NULL != clp_next) && (ii < i32_count); ++ii) {
		/* 次の点が線分としてリンク済みのとき */
		if ((NULL != clp_next->get_clp_next_point()) || (NULL != clp_next->get_clp_previous_point())) {
			/* カレント点をゴール点としてセット */
			this->link_another(clp_crnt);
			return;
		}

		/* 線分としてのリンクをする */
		clp_crnt->set_clp_next_point(clp_next);
		clp_next->set_clp_previous_point(clp_crnt);
		++(this->_i32_point_count);

		/* もう一方の端点に到着し終端(リンク2点目がないので) */
		if ((NULL == clp_next->get_clp_link_near(1))) {
			/* 次の点をゴール点としてセット */
			this->link_another(clp_next);
			return;
		}

		/* 分岐点(リンク3点目があるので)から次の点へ */
		if (NULL != clp_next->get_clp_link_near(2)) {

			/* それまでの線分の方向を保存しておく */
			clp_tmp2 = clp_tmp1;
			clp_tmp1 = clp_next;
			assert(clp_tmp1 != clp_tmp2);

			d_xv = clp_tmp1->get_i32_xp() - clp_tmp2->get_i32_xp();
			d_yv = clp_tmp1->get_i32_yp() - clp_tmp2->get_i32_yp();
			assert((0.0 != d_xv) || (0.0 != d_yv));

			/* 各分岐線をたどり、方向の最も合う線を探す */
			clp_crnt = clp_next;
			clp_next = this->_get_link_line_selecter(
				d_xv, d_yv, clp_next, i32_count);

			/* たどる線がすでにたどっているものの場合、
			そこで、終端 */
			if (NULL == clp_next) {
				/* ゴール点をセット */
				this->link_another(clp_crnt);
				return;
			}
		}
		/* 線分構成点(リンク3点目がない)から次の点へ */
		else {
			/* 次へのリンクをたどる */
			if (clp_crnt == clp_next->get_clp_link_near(0)) {
				clp_crnt = clp_next;
				clp_next = clp_next->get_clp_link_near(1);
			} else if (clp_crnt == clp_next->get_clp_link_near(1)) {
				clp_crnt = clp_next;
				clp_next = clp_next->get_clp_link_near(0);
			}
			/* リンクが壊れている */
			else {
				pri_funct_err_bttvr("Error : bad link");
				assert(0);
			}
		}
	}
	/* リンクが不自然に長すぎる */
	pri_funct_err_bttvr("Error : too long link");
	assert(0);
	return;
}

#include "igs_line_blur.h" // "pixel_line_node.h"

/*
pixel_select_curve_blur.cxx
の、
void pixel_select_curve_blur_root::exec( double d_xp, double d_yp, pixel_line_node *clp_line_first, int32_t i32_count, int32_t i32_blur_count )
で使用。
*/

pixel_point_node *pixel_line_node::get_next_point_by_count(pixel_point_node *clp_point, int32_t i32_count)
{
	for (; (0 < i32_count) && (NULL != clp_point); --i32_count) {
		clp_point = (pixel_point_node *)clp_point->get_clp_next_point();
	};
	return clp_point;
}

pixel_point_node *pixel_line_node::get_prev_point_by_count(pixel_point_node *clp_point, int32_t i32_count)
{
	for (; (0 < i32_count) && (NULL != clp_point); --i32_count) {
		clp_point = (pixel_point_node *)clp_point->get_clp_previous_point();
	};
	return clp_point;
}

#include <assert.h>		   /* assert() */
#include "igs_line_blur.h" // "pri.h" "pixel_line_node.h"

int pixel_line_node::save_line(FILE *fp)
{
	pixel_point_node *clp_crnt;
	int32_t ii;

	clp_crnt = this->get_clp_link_one();
	for (ii = 0; NULL != clp_crnt;
		 clp_crnt = clp_crnt->get_clp_next_point(), ++ii) {
		/* 偽の場合、たぶん無限ループ */
		assert(ii < this->_i32_point_count);

		/* 次の点(シード2含む)の座標値保存 */
		if (fprintf(fp, "%d %d\n",
					clp_crnt->get_i32_xp(), clp_crnt->get_i32_yp()) < 0) {
			pri_funct_err_bttvr(
				"Error : point %d : fprintf(%d %d) returns minus",
				ii, clp_crnt->get_i32_xp(), clp_crnt->get_i32_yp());
			return NG;
		}
	}
	return OK;
}

/********************************************************************/

int pixel_line_node::save_one_point(FILE *fp)
{
	pixel_point_node *clp_crnt;

	/* シード1の座標値保存 */
	clp_crnt = this->get_clp_link_one();

	if (NULL != clp_crnt) {
		if (fprintf(fp, "%d %d\n",
					clp_crnt->get_i32_xp(), clp_crnt->get_i32_yp()) < 0) {
			pri_funct_err_bttvr(
				"Error : one point : fprintf(%d %d) returns minus",
				clp_crnt->get_i32_xp(), clp_crnt->get_i32_yp());
			return NG;
		}
	}
	return OK;
}

int pixel_line_node::save_middle_point(FILE *fp)
{
	pixel_point_node *clp_crnt;

	/* シード1の座標値保存 */
	clp_crnt = this->get_clp_link_middle();

	if (NULL != clp_crnt) {
		if (fprintf(fp, "%d %d\n",
					clp_crnt->get_i32_xp(), clp_crnt->get_i32_yp()) < 0) {
			pri_funct_err_bttvr(
				"Error : middle point : fprintf(%d %d) returns minus",
				clp_crnt->get_i32_xp(), clp_crnt->get_i32_yp());
			return NG;
		}
	}
	return OK;
}

int pixel_line_node::save_another_point(FILE *fp)
{
	pixel_point_node *clp_crnt;

	/* シード2の座標値保存 */
	clp_crnt = this->get_clp_link_another();

	if (NULL != clp_crnt) {
		if (fprintf(fp, "%d %d\n",
					clp_crnt->get_i32_xp(), clp_crnt->get_i32_yp()) < 0) {
			pri_funct_err_bttvr(
				"Error : another point : fprintf(%d %d) returns minus",
				clp_crnt->get_i32_xp(), clp_crnt->get_i32_yp());
			return NG;
		}
	}

	return OK;
}

#include <assert.h> /* assert() */

#include "igs_line_blur.h" // "pri.h" "pixel_line_node.h"

int pixel_line_node::save_expand_line(FILE *fp)
{
	pixel_point_node *clp_crnt;
	int32_t ii;

	clp_crnt = this->get_clp_link_one_expand();
	for (ii = 0; NULL != clp_crnt;
		 clp_crnt = clp_crnt->get_clp_next_point(), ++ii) {
		/* 伸長した後なので3倍-2の長さ。偽の場合、たぶん無限ループ */
		/***assert( ii < this->_i32_point_count * 3 - 2 );***/
		assert(ii < this->_i32_point_count);

		/* 次の点(シード2含む)の座標値保存 */
		if (fprintf(fp, "%g %g\n",
					clp_crnt->get_d_xp_tgt(), clp_crnt->get_d_yp_tgt()) < 0) {
			pri_funct_err_bttvr(
				"Error : point %d : fprintf(%g %g) returns minus",
				ii, clp_crnt->get_d_xp_tgt(), clp_crnt->get_d_yp_tgt());
			return NG;
		}
	}
	return OK;
}

/********************************************************************/

int pixel_line_node::save_one_expand_point(FILE *fp)
{
	pixel_point_node *clp_crnt;

	clp_crnt = this->get_clp_link_one_expand();

	if (NULL != clp_crnt) {
		if (fprintf(fp, "%g %g\n",
					clp_crnt->get_d_xp_tgt(), clp_crnt->get_d_yp_tgt()) < 0) {
			pri_funct_err_bttvr(
				"Error : one expand point : fprintf(%g %g) returns minus",
				clp_crnt->get_d_xp_tgt(), clp_crnt->get_d_yp_tgt());
			return NG;
		}
	}
	return OK;
}

int pixel_line_node::save_another_expand_point(FILE *fp)
{
	pixel_point_node *clp_crnt;

	clp_crnt = this->get_clp_link_another_expand();

	if (NULL != clp_crnt) {
		if (fprintf(fp, "%g %g\n",
					clp_crnt->get_d_xp_tgt(), clp_crnt->get_d_yp_tgt()) < 0) {
			pri_funct_err_bttvr(
				"Error : another expand point : fprintf(%g %g) returns minus",
				clp_crnt->get_d_xp_tgt(), clp_crnt->get_d_yp_tgt());
			return NG;
		}
	}

	return OK;
}

#include <assert.h> /* assert() */

#include "igs_line_blur.h" // "pri.h" "pixel_line_node.h"

int pixel_line_node::save_expand_vector(FILE *fp)
{
	pixel_point_node *clp_term, *clp_expa;

	clp_term = this->get_clp_link_one();
	clp_expa = this->get_clp_link_one_expand();

	if ((NULL != clp_term) && (NULL != clp_expa)) {
		(void)fprintf(fp, "%g %g %g %g\n",
					  clp_term->get_d_xp_tgt(),
					  clp_term->get_d_yp_tgt(),
					  clp_expa->get_d_xp_tgt() - clp_term->get_d_xp_tgt(),
					  clp_expa->get_d_yp_tgt() - clp_term->get_d_yp_tgt());
		if (ferror(fp)) {
			pri_funct_err_bttvr(
				"Error : fprintf(one and one_expand xp and yp)");
			return NG;
		}
	}

	clp_term = this->get_clp_link_another();
	clp_expa = this->get_clp_link_another_expand();

	if ((NULL != clp_term) && (NULL != clp_expa)) {
		(void)fprintf(fp, "%g %g %g %g\n",
					  clp_term->get_d_xp_tgt(),
					  clp_term->get_d_yp_tgt(),
					  clp_expa->get_d_xp_tgt() - clp_term->get_d_xp_tgt(),
					  clp_expa->get_d_yp_tgt() - clp_term->get_d_yp_tgt());
		if (ferror(fp)) {
			pri_funct_err_bttvr(
				"Error : fprintf(another and another_expand xp and yp)");
			return NG;
		}
	}

	return OK;
}

#include <assert.h> /* assert() */

#include "igs_line_blur.h" // "pixel_line_node.h"

void pixel_line_node::set_bbox(void)
{
	pixel_point_node *clp_point;
	int32_t ii;

	/* 浮動小数化 */
	clp_point = this->get_clp_link_one_expand();
	if (NULL == clp_point) {
		clp_point = this->get_clp_link_one();
	}
	for (ii = 0; NULL != clp_point;
		 clp_point = clp_point->get_clp_next_point(), ++ii) {
		/* 偽の場合、たぶん無限ループ */
		assert(ii < this->_i32_point_count);

		if (0 == ii) {
			this->_d_bbox_x_min = clp_point->get_d_xp_tgt();
			this->_d_bbox_x_max = clp_point->get_d_xp_tgt();
			this->_d_bbox_y_min = clp_point->get_d_yp_tgt();
			this->_d_bbox_y_max = clp_point->get_d_yp_tgt();
		} else {
			if (clp_point->get_d_xp_tgt() < this->_d_bbox_x_min) {
				this->_d_bbox_x_min = clp_point->get_d_xp_tgt();
			} else if (this->_d_bbox_x_max < clp_point->get_d_xp_tgt()) {
				this->_d_bbox_x_max = clp_point->get_d_xp_tgt();
			}
			if (clp_point->get_d_yp_tgt() < this->_d_bbox_y_min) {
				this->_d_bbox_y_min = clp_point->get_d_yp_tgt();
			} else if (this->_d_bbox_y_max < clp_point->get_d_yp_tgt()) {
				this->_d_bbox_y_max = clp_point->get_d_yp_tgt();
			}
		}
	}
}

#include <assert.h> /* assert() */

#include "igs_line_blur.h" // "pixel_line_node.h"

void pixel_line_node::set_middle(void)
{
	int32_t ii;
	pixel_point_node *clp_middle;

	/* すでにラインデータとしてのリンクがあること */
	assert(NULL != this->get_clp_link_one());

	/* 中間点を得る(nodeが偶数個のときはendよりのnode) */
	for (ii = 0, clp_middle = this->get_clp_link_one();
		 (ii < (this->_i32_point_count) / 2) && (NULL != clp_middle);
		 ++ii) {
		if (NULL != clp_middle->get_clp_next_point()) {
			clp_middle = clp_middle->get_clp_next_point();
		}
	}
	/* 短すぎる線は選択していないはず */
	assert(clp_middle != this->get_clp_link_one());
	assert(clp_middle != this->get_clp_link_another());

	/* 中間点を登録 */
	this->link_middle(clp_middle);
}

#include <assert.h> /* assert() */

#include "igs_line_blur.h" // "pixel_line_node.h"

/* 3点から線分がスムースになるような中点を計算する
	* p1
	|
	|
	|
	|       +
	|
	|   * <-- new p2
	|
	*---------------*
	p2              |p3
		    *   |
			|
		+	|       +
			|
			|   *
			|
			*----------------*
*/
void pixel_line_node::_smooth_point(double d_x1, double d_y1, double d_x2, double d_y2, double d_x3, double d_y3, double *dp_x, double *dp_y)
{
	*dp_x = (d_x2 + ((d_x1 + d_x3) / 2.0)) / 2.0;
	*dp_y = (d_y2 + ((d_y1 + d_y3) / 2.0)) / 2.0;
}

void pixel_line_node::smooth_body(int32_t i32_smooth_retry)
{
	pixel_point_node *clp_1, *clp_2, *clp_3;
	int32_t kk, ii;
	double d_x, d_y, d_x1, d_y1, d_x2 = 0, d_y2 = 0, d_x3, d_y3;

	/* スムース化 */
	for (kk = 0; kk < i32_smooth_retry; ++kk) {
		clp_1 = this->get_clp_link_one();
		clp_2 = clp_3 = NULL;
		for (ii = 0; NULL != clp_1; ++ii) {
			/* 偽の場合、たぶん無限ループ */
			assert(ii < this->_i32_point_count);

			d_x1 = clp_1->get_d_xp_tgt();
			d_y1 = clp_1->get_d_yp_tgt();

			/* スムース化 */
			if (NULL != clp_3) {
				this->_smooth_point(
					d_x1, d_y1, d_x2, d_y2, d_x3, d_y3, &d_x, &d_y);
				clp_2->set_d_xp_tgt(d_x);
				clp_2->set_d_yp_tgt(d_y);
			}

			/* 細線化部分の終点で終り */
			if (this->get_clp_link_another() == clp_1)
				break;

			/* 次の3ポイントへずらす */
			clp_3 = clp_2;
			d_x3 = d_x2;
			d_y3 = d_y2;
			clp_2 = clp_1;
			d_x2 = d_x1;
			d_y2 = d_y1;
			clp_1 = clp_1->get_clp_next_point();
		}
	}
}

void pixel_line_node::smooth_expand(int32_t i32_smooth_retry)
{
	pixel_point_node *clp_1, *clp_2, *clp_3;
	int32_t kk, ii;
	double d_x, d_y, d_x1, d_y1, d_x2, d_y2, d_x3, d_y3;

	/* スムース化 */
	for (kk = 0; kk < i32_smooth_retry; ++kk) {
		clp_1 = this->get_clp_link_one_expand();
		clp_2 = clp_3 = NULL;
		for (ii = 0; NULL != clp_1; ++ii) {
			/* 偽の場合、たぶん無限ループ */
			/***assert( ii < (this->_i32_point_count * 3) );***/
			assert(ii < (this->_i32_point_count));

			d_x1 = clp_1->get_d_xp_tgt();
			d_y1 = clp_1->get_d_yp_tgt();

			/* スムース化 */
			if (NULL != clp_3) {
				this->_smooth_point(
					d_x1, d_y1, d_x2, d_y2, d_x3, d_y3, &d_x, &d_y);
				clp_2->set_d_xp_tgt(d_x);
				clp_2->set_d_yp_tgt(d_y);
			}

			/* 伸ばした終点で終り */
			if (this->get_clp_link_another_expand() == clp_1)
				break;

			/* 次の3ポイントへずらす */
			clp_3 = clp_2;
			d_x3 = d_x2;
			d_y3 = d_y2;
			clp_2 = clp_1;
			d_x2 = d_x1;
			d_y2 = d_y1;
			clp_1 = clp_1->get_clp_next_point();
		}
	}
}

#ifndef _pixel_select_same_way_h_
#define _pixel_select_same_way_h_

#include <stdio.h>

#include "igs_line_blur.h" // "list_root.h" "calculator_geometry.h" "pixel_point_root.h" "pixel_line_root.h"

class pixel_select_same_way_node : public list_node
{
public:
	pixel_select_same_way_node(void)
	{
		this->clp_point_middle = NULL;
		this->clp_point_term = NULL;
		this->clp_point_expand = NULL;
		this->d_length = -1.0;
	}
	void copy(pixel_select_same_way_node *clp)
	{
		this->clp_point_middle = clp->clp_point_middle;
		this->clp_point_term = clp->clp_point_term;
		this->clp_point_expand = clp->clp_point_expand;
		this->d_length = clp->d_length;
	}
	pixel_point_node *clp_point_middle;
	pixel_point_node *clp_point_term;
	pixel_point_node *clp_point_expand;
	double d_length;

private:
};

class pixel_select_same_way_root : public list_root
{
public:
	pixel_select_same_way_root(void)
	{
		this->_i_mv_sw = OFF;
		this->_i_cv_sw = OFF;
		this->_i_pv_sw = OFF;

		this->_i32_count_max = 4;
		this->_d_length_max = 160;
	}
	~pixel_select_same_way_root(void)
	{
		this->mem_free();
	}
	void set_i_mv_sw(int sw) { this->_i_mv_sw = sw; }
	void set_i_cv_sw(int sw) { this->_i_cv_sw = sw; }
	void set_i_pv_sw(int sw) { this->_i_pv_sw = sw; }

	int32_t get_i32_count_max(void) { return this->_i32_count_max; }
	void set_i32_count_max(int32_t ii) { this->_i32_count_max = ii; }
	double get_d_length_max(void) { return this->_d_length_max; }
	void set_d_length_max(double dd) { this->_d_length_max = dd; }

	void exec(pixel_line_node *clp_line_first, int32_t i32_count, pixel_point_node *clp_point_middle, pixel_point_node *clp_point_term, pixel_point_node *clp_point_expand);
	void get_vector(double *dp_xv, double *dp_yv);

	void mem_free(void);

private:
	int _i_mv_sw,
		_i_cv_sw,
		_i_pv_sw;

	int32_t _i32_count_max;
	double _d_length_max;

	calculator_geometry _cl_cal_geom;

	pixel_select_same_way_node *_append(pixel_select_same_way_node *clp_previous);
	void _remove(pixel_select_same_way_node *clp_old);

	int _sort_append(pixel_select_same_way_node *clp_src);
	double _term_length(pixel_point_node *clp_middle1, pixel_point_node *clp_term1, pixel_point_node *clp_middle2, pixel_point_node *clp_term2);
};

#endif				/* !_pixel_select_same_way_h_ */
#include <math.h>   /* sqrt() */
#include <stdlib.h> /* free(), calloc() */
#include <assert.h> /* assert() */

#include "igs_line_blur.h" // "pri.h" "pixel_line_root.h" "pixel_select_same_way.h"

/* WindowsではM_PIが見つからない */
#if defined _WIN32
#define M_PI 3.14159265358979323846
#endif

pixel_select_same_way_node *pixel_select_same_way_root::_append(pixel_select_same_way_node *clp_previous)
{
	pixel_select_same_way_node *clp_new;

	clp_new = new pixel_select_same_way_node;
	if (NULL == clp_new) {
		pri_funct_err_bttvr(
			"Error : 'new pixel_select_same_way_node' returns NULL.");
		return NULL;
	}
	clp_new = (pixel_select_same_way_node *)this->push(clp_previous, clp_new);

	return clp_new;
}

void pixel_select_same_way_root::_remove(pixel_select_same_way_node *clp_old)
{
	assert(NULL != clp_old); /* あってはならないプログラムバグのチェック */
	this->pop(clp_old);
	delete clp_old;
}

void pixel_select_same_way_root::mem_free(void)
{
	pixel_select_same_way_node *clp_;
	int32_t ii;

	if (NULL != this->get_clp_last()) {
		if (ON == this->_i_mv_sw) {
			pri_funct_msg_ttvr(
				"pixel_select_same_way_root::mem_free()");
		}

		ii = 0;
		while (NULL !=
			   (clp_ = (pixel_select_same_way_node *)this->get_clp_last())) {
			this->_remove(clp_);
			++ii;
		}
		if (ON == this->_i_pv_sw) {
			pri_funct_msg_ttvr("free select same way node %d", ii);
		}
	}
}

/********************************************************************/

/* lengthの値を見て、小さい値順になるように追加する
	他の選択より小さいなら先頭
	他の選択より大きいなら最後に追加
	他の選択の値の間入るならそこに追加
*/
int pixel_select_same_way_root::_sort_append(pixel_select_same_way_node *clp_src)
{
	pixel_select_same_way_node *clp_loop, *clp_prev;
	int32_t ii;

	assert(0.0 < clp_src->d_length);

	clp_loop = (pixel_select_same_way_node *)this->get_clp_first();

	/* まだなにも選択されていないとき */
	if (NULL == clp_loop) {
		clp_prev = NULL;
	}
	/* 1つだけ選択されているとき */
	else if (NULL == clp_loop->get_clp_next()) {
		/* 既選択より値が小さいので先頭位置を示す */
		if (clp_src->d_length < clp_loop->d_length) {
			clp_prev = NULL;
		}
		/* 既選択より値が大きいので既選択を示す */
		else {
			clp_prev = clp_loop;
		}
	}
	/* 2つ以上の選択があるとき */
	else {
		/* 既選択より値が大きいときのため最後尾を示しておく */
		clp_prev = (pixel_select_same_way_node *)this->get_clp_last();
		/* 既選択より値が小さい場所があればその手前を示す */
		for (ii = 0; NULL != clp_loop; ++ii,
			clp_loop = (pixel_select_same_way_node *)clp_loop->get_clp_next()) {
			assert(ii < this->get_i32_count());

			if (clp_src->d_length < clp_loop->d_length) {
				clp_prev = (pixel_select_same_way_node *)clp_loop->get_clp_previous();
				break;
			}
		}
	}

	/* 示した場所の後ろに追加 */
	clp_prev = this->_append(clp_prev);
	if (NULL == clp_prev) {
		pri_funct_err_bttvr(
			"Error : this->_append(clp_prev) returns NULL.");
		return NG;
	}

	/* 情報をコピー */
	clp_prev->copy(clp_src);

	return OK;
}

/********************************************************************/

double pixel_select_same_way_root::_term_length(pixel_point_node *clp_middle1, pixel_point_node *clp_term1, pixel_point_node *clp_middle2, pixel_point_node *clp_term2)
{
	double d_length, d_radian;

	/* 2つの端点の距離 */
	assert(clp_term2 != clp_term1);
	d_length = sqrt(
		(clp_term2->get_d_xp_tgt() - clp_term1->get_d_xp_tgt()) *
			(clp_term2->get_d_xp_tgt() - clp_term1->get_d_xp_tgt()) +
		(clp_term2->get_d_yp_tgt() - clp_term1->get_d_yp_tgt()) *
			(clp_term2->get_d_yp_tgt() - clp_term1->get_d_yp_tgt()));

	/* 同方向のものでなければならない
		2つの線分の角度をしらべ...
	*/
	assert(clp_term1 != clp_middle1);
	assert(clp_term2 != clp_middle2);
	d_radian = this->_cl_cal_geom.get_d_radian_by_2_vector(
		clp_term1->get_d_xp_tgt() - clp_middle1->get_d_xp_tgt(),
		clp_term1->get_d_yp_tgt() - clp_middle1->get_d_yp_tgt(),
		clp_term2->get_d_xp_tgt() - clp_middle2->get_d_xp_tgt(),
		clp_term2->get_d_yp_tgt() - clp_middle2->get_d_yp_tgt());
	/* 線分が逆方向ならだめ */
	if ((M_PI / 2.0 <= d_radian) && (d_radian <= M_PI * 3.0 / 2.0))
		return -1.0;

	/* 端点が離れすぎならだめ */
	if (this->_d_length_max <= d_length)
		return -1.0;

	/* 距離を返す */
	return d_length;
}

/********************************************************************/

/* 同方向で、一定の距離の範囲にある端点を選択する */
void pixel_select_same_way_root::exec(pixel_line_node *clp_line_first, int32_t i32_count, pixel_point_node *clp_point_middle, pixel_point_node *clp_point_term, pixel_point_node *clp_point_expand)
{
	pixel_line_node *clp_line;
	int32_t ii;
	double d_length;
	pixel_select_same_way_node cl_select;

	assert(clp_point_middle != clp_point_term);
	assert(clp_point_term != clp_point_expand);
	assert(clp_point_expand != clp_point_middle);

	/* 選択リストをクリア */
	this->mem_free();

	/* 自身以外の全端点と比較 */
	for (
		clp_line = clp_line_first, ii = 0;
		NULL != clp_line;
		clp_line = (pixel_line_node *)clp_line->get_clp_next(), ++ii) {
		assert(ii < i32_count);

		/* 端点1が自分以外のときで、
		   線分が短すぎて(1点以下)真中指定と端点1が同じではない */
		if ((clp_line->get_clp_link_one() != clp_point_term) && (clp_line->get_clp_link_middle() !=
																 clp_line->get_clp_link_one())) {
			assert(NULL != clp_line->get_clp_link_middle());
			assert(NULL != clp_line->get_clp_link_one());
			d_length = this->_term_length(
				clp_point_middle,
				clp_point_term,
				clp_line->get_clp_link_middle(),
				clp_line->get_clp_link_one());
			/* 同方向で一定の距離の範囲内にあるなら */
			if ((0.0 < d_length) && (NULL != clp_line->get_clp_link_one_expand())) {
				cl_select.clp_point_middle =
					clp_line->get_clp_link_middle();
				cl_select.clp_point_term =
					clp_line->get_clp_link_one();
				cl_select.clp_point_expand =
					clp_line->get_clp_link_one_expand();
				cl_select.d_length = d_length;
				/* 登録 */
				this->_sort_append(&cl_select);
			}
		}

		/* 端点2が自分以外のときで、
		   線分が短すぎて(2点以下)真中指定と端点2が同じではない */
		if ((clp_line->get_clp_link_another() != clp_point_term) && (clp_line->get_clp_link_middle() !=
																	 clp_line->get_clp_link_another())) {
			assert(NULL != clp_line->get_clp_link_middle());
			assert(NULL != clp_line->get_clp_link_another());
			/* 端点(another)を調べて */
			d_length = this->_term_length(
				clp_point_middle,
				clp_point_term,
				clp_line->get_clp_link_middle(),
				clp_line->get_clp_link_another());
			/* 同方向で一定の距離の範囲内にあるなら */
			if ((0.0 < d_length) && (NULL != clp_line->get_clp_link_another_expand())) {
				cl_select.clp_point_middle =
					clp_line->get_clp_link_middle();
				cl_select.clp_point_term =
					clp_line->get_clp_link_another();
				cl_select.clp_point_expand =
					clp_line->get_clp_link_another_expand();
				cl_select.d_length = d_length;
				/* 登録 */
				this->_sort_append(&cl_select);
			}
		}
	}
}

/********************************************************************/

/* 距離の近い順(lengthの小さい値)に指定個数分のベクトルを合成して返す
	先頭から
	this->_i32_count_max数個分
	のベクトルを合成する
*/
void pixel_select_same_way_root::get_vector(double *dp_xv, double *dp_yv)
{
	pixel_select_same_way_node *clp_loop;
	int32_t ii;

	clp_loop = (pixel_select_same_way_node *)this->get_clp_first();
	for (ii = 0; (NULL != clp_loop) && (ii < this->_i32_count_max); ++ii,
		clp_loop = (pixel_select_same_way_node *)clp_loop->get_clp_next()) {
		assert(ii < this->get_i32_count());
		assert(NULL != clp_loop->clp_point_expand);
		assert(NULL != clp_loop->clp_point_term);

		*dp_xv +=
			clp_loop->clp_point_expand->get_d_xp_tgt() -
			clp_loop->clp_point_term->get_d_xp_tgt();
		*dp_yv +=
			clp_loop->clp_point_expand->get_d_yp_tgt() -
			clp_loop->clp_point_term->get_d_yp_tgt();
	}
}

#ifndef _pixel_line_root_h_
#define _pixel_line_root_h_

#include <stdio.h>

#include "igs_line_blur.h" // "list_root.h" "calculator_geometry.h" "pixel_point_root.h" "pixel_line_node.h" "pixel_select_same_way.h"

class pixel_line_root : public list_root
{
public:
	pixel_line_root(void)
	{
		this->_i_mv_sw = OFF;
		this->_i_cv_sw = OFF;
		this->_i_pv_sw = OFF;

		this->_d_bbox_x_min = 0.0;
		this->_d_bbox_x_max = 0.0;
		this->_d_bbox_y_min = 0.0;
		this->_d_bbox_y_max = 0.0;

		this->_i32_smooth_retry = 100;
		this->_i_same_way_exec_sw = ON;
	}
	~pixel_line_root(void)
	{
		this->mem_free();
	}
	void set_i_mv_sw(int sw) { this->_i_mv_sw = sw; }
	void set_i_cv_sw(int sw) { this->_i_cv_sw = sw; }
	void set_i_pv_sw(int sw) { this->_i_pv_sw = sw; }

	double get_d_bbox_x_min(void) { return this->_d_bbox_x_min; }
	double get_d_bbox_x_max(void) { return this->_d_bbox_x_max; }
	double get_d_bbox_y_min(void) { return this->_d_bbox_y_min; }
	double get_d_bbox_y_max(void) { return this->_d_bbox_y_max; }

	void set_i32_smooth_retry(int32_t i32) { this->_i32_smooth_retry = i32; }
	int32_t get_i32_smooth_retry(void) { return this->_i32_smooth_retry; }
	void set_i_same_way_exec_sw(int sw) { this->_i_same_way_exec_sw = sw; }

	int exec01020304(pixel_point_root *clp_pixel_point_root);
	void exec05_set_middle(void);
	void exec06_int2double_body(void);
	void exec07_smooth_body(void);
	int exec08_expand_lines(pixel_point_root *clp_pixel_point_root);
	void exec09_same_way_expand(pixel_select_same_way_root *clp_select);
	void exec10_smooth_expand(void);
	void exec11_set_bbox(void);

	int save_not_include(pixel_point_root *clp_pixel_point_root, char *cp_fname);
	int save_lines(char *cp_fname);
	int save_one_point(char *cp_fname);
	int save_middle_point(char *cp_fname);
	int save_another_point(char *cp_fname);

	int save_expand_lines(char *cp_fname);
	int save_one_expand_point(char *cp_fname);
	int save_another_expand_point(char *cp_fname);
	int save_expand_vector(char *cp_fname);

	void mem_free(void);

private:
	int _i_mv_sw,
		_i_cv_sw,
		_i_pv_sw;

	double _d_bbox_x_min,
		_d_bbox_x_max,
		_d_bbox_y_min,
		_d_bbox_y_max;

	int32_t _i32_smooth_retry;
	int _i_same_way_exec_sw;

	calculator_geometry _cl_cal_geom;

	pixel_line_node *_append(pixel_line_node *clp_previous);
	void _remove(pixel_line_node *clp_old);

	int _exec01_link_left_right(pixel_point_root *clp_pixel_point_root);
	int _exec02_link_up_down(pixel_point_root *clp_pixel_point_root);
	int _exec03_link_slant(pixel_point_root *clp_pixel_point_root);
	int _exec04_grouping(pixel_point_root *clp_pixel_point_root);

	double _same_way_expand_radian_diff(pixel_point_node *clp_point_middle, pixel_point_node *clp_point_term, pixel_point_node *clp_point_term_expand, pixel_select_same_way_root *clp_select);
};

#endif					   /* !_pixel_line_root_h_ */
#include "igs_line_blur.h" // "pri.h" "pixel_line_root.h"

int pixel_line_root::exec01020304(pixel_point_root *clp_pixel_point_root)
{
	/* ピクセルノードの左右リンク */
	if (OK != this->_exec01_link_left_right(clp_pixel_point_root)) {
		pri_funct_err_bttvr(
			"Error : this->_exec01_link_left_right() returns NG.");
		return NG;
	}
	/* ピクセルノードの上下リンク */
	if (OK != this->_exec02_link_up_down(clp_pixel_point_root)) {
		pri_funct_err_bttvr(
			"Error : this->_exec02_link_up_down() returns NG.");
		return NG;
	}
	/* ピクセルノードの斜めリンク */
	if (OK != this->_exec03_link_slant(clp_pixel_point_root)) {
		pri_funct_err_bttvr(
			"Error : this->_exec03_link_slant() returns NG.");
		return NG;
	}

	/* リンクしたものを、グループ化 */
	if (OK != this->_exec04_grouping(clp_pixel_point_root)) {
		pri_funct_err_bttvr(
			"Error : this->_exec04_grouping() returns NG.");
		return NG;
	}

	return OK;
}

/********************************************************************/

int pixel_line_root::_exec01_link_left_right(pixel_point_root *clp_pixel_point_root)
{
	pixel_point_node *clp_point,
		*clp_point2;
	int32_t ii;

	if (ON == this->_i_mv_sw) {
		pri_funct_msg_ttvr("pixel_line_root::_exec01_link_left_right()");
	}

	ii = 0;
	clp_point2 = NULL;
	for (
		clp_point = (pixel_point_node *)(clp_pixel_point_root->get_clp_first());
		NULL != clp_point;
		clp_point = (pixel_point_node *)clp_point->get_clp_next()) {
		/* 隣同士の２ポイントで */
		if ((NULL != clp_point2)
			/* 同じスキャンラインで */
			&& (clp_point->get_i32_yp() == clp_point2->get_i32_yp())
			/* 左右隣同士なら */
			&& (((clp_point2->get_i32_xp()) + 1) == (clp_point->get_i32_xp()))) {
			/* 双方向リンクする */
			if (NG == clp_point->link_near(clp_point2)) {
				pri_funct_err_bttvr(
					"Error : count %d : clp_point->link_near() returns NG.", ii);
				return NG;
			}
			if (NG == clp_point2->link_near(clp_point)) {
				pri_funct_err_bttvr(
					"Error : count %d : clp_point2->link_near() returns NG.", ii);
				return NG;
			}
			++ii;
		}
		clp_point2 = clp_point;
	}

	if (ON == this->_i_pv_sw) {
		pri_funct_msg_ttvr(" link left right %d", ii);
	}

	return OK;
}
int pixel_line_root::_exec02_link_up_down(pixel_point_root *clp_pixel_point_root)
{
	pixel_point_node *clp_point,
		*clp_point2;
	int32_t ii;

	if (ON == this->_i_mv_sw) {
		pri_funct_msg_ttvr("pixel_line_root::_exec02_link_up_down()");
	}

	ii = 0;

	for (
		clp_point = (pixel_point_node *)(clp_pixel_point_root->get_clp_first());
		NULL != clp_point;
		clp_point = (pixel_point_node *)clp_point->get_clp_next()) {
		for (
			clp_point2 = (pixel_point_node *)clp_point->get_clp_next();
			(NULL != clp_point2)
			/* １つ上のスキャンラインまでを見る */
			&& (clp_point2->get_i32_yp() <= (1 + clp_point->get_i32_yp()));
			clp_point2 = (pixel_point_node *)clp_point2->get_clp_next()) {
			if (
				/* １つ上のスキャンラインで */
				(clp_point2->get_i32_yp() ==
				 (1 + (clp_point->get_i32_yp())))
				/* 上下隣接同士なら */
				&& (clp_point->get_i32_xp() ==
					clp_point2->get_i32_xp())) {
				/* 双方向リンクする */
				if (NG == clp_point->link_near(clp_point2)) {
					pri_funct_err_bttvr(
						"Error : count %d : clp_point->link_near() returns NG.", ii);
					return NG;
				}
				if (NG == clp_point2->link_near(clp_point)) {
					pri_funct_err_bttvr(
						"Error : count %d : clp_point2->link_near() returns NG.", ii);
					return NG;
				}
				++ii;

				/* 一個つなげたら抜ける(上下一致するのは一か所だけ)*/
				break;
			}
		}
	}

	if (ON == this->_i_pv_sw) {
		pri_funct_msg_ttvr(" link up down %d", ii);
	}

	return OK;
}

int pixel_line_root::_exec03_link_slant(pixel_point_root *clp_pixel_point_root)
{
	pixel_point_node *clp_point, *clp_point2,
		*clp_point_link;
	int32_t ii, jj;

	int i_left_link_sw,
		i_right_link_sw,
		i_up_link_sw;

	if (ON == this->_i_mv_sw) {
		pri_funct_msg_ttvr("pixel_line_root::_exec03_link_slant()");
	}

	ii = 0;

	for (
		clp_point = (pixel_point_node *)(clp_pixel_point_root->get_clp_first());
		NULL != clp_point;
		clp_point = (pixel_point_node *)clp_point->get_clp_next()) {
		/* LINK_NEAR_COUNT個所リンクしているところは、
		   これ以上リンクできないので、次へ */
		if (NULL != clp_point->get_clp_link_near(LINK_NEAR_COUNT - 1)) {
			continue;
		}

		i_left_link_sw = OFF;
		i_right_link_sw = OFF;
		i_up_link_sw = OFF;

		/* 全リンクの調査 */
		for (jj = 0; jj < LINK_NEAR_COUNT; ++jj) {
			clp_point_link = clp_point->get_clp_link_near(jj);
			/* リンクの終了なので抜ける */
			if (NULL == clp_point_link)
				break;
			/* 上につながっている */
			if (((clp_point_link->get_i32_xp()) == (clp_point->get_i32_xp())) && ((clp_point_link->get_i32_yp()) == (clp_point->get_i32_yp() + 1))) {
				i_up_link_sw = ON;
			}
			/* 左につながっている */
			if (((clp_point_link->get_i32_xp()) == (clp_point->get_i32_xp() - 1)) && ((clp_point_link->get_i32_yp()) == (clp_point->get_i32_yp()))) {
				i_left_link_sw = ON;
			}
			/* 右につながっている */
			if (((clp_point_link->get_i32_xp()) == (clp_point->get_i32_xp() + 1)) && ((clp_point_link->get_i32_yp()) == (clp_point->get_i32_yp()))) {
				i_right_link_sw = ON;
			}
		}

		/* 左右両方につながりあるなら、斜めにはつなげず、次へ */
		if ((ON == i_left_link_sw) && (ON == i_right_link_sw))
			continue;
		/* 上につながっているなら、斜めにはつなげず、次へ */
		if (ON == i_up_link_sw)
			continue;

		for (
			clp_point2 = (pixel_point_node *)clp_point->get_clp_next();
			(NULL != clp_point2)
			/* １つ上のスキャンラインまでを見る */
			&& (clp_point2->get_i32_yp() <= (1 + clp_point->get_i32_yp()));
			clp_point2 = (pixel_point_node *)clp_point2->get_clp_next()) {
			/* １つ上のスキャンラインでないなら次へ */
			if ((clp_point2->get_i32_yp()) !=
				(1 + (clp_point->get_i32_yp()))) {
				continue;
			}

			/* 左斜め上につながるものあり */
			if ((OFF == i_left_link_sw) && ((clp_point2->get_i32_xp()) ==
											(clp_point->get_i32_xp() - 1)) &&
				((clp_point2->get_i32_yp()) ==
				 (clp_point->get_i32_yp() + 1))) {
				/* 双方向リンクする */
				if (NG == clp_point->link_near(clp_point2)) {
					pri_funct_err_bttvr(
						"Error : count %d : clp_point->link_near() returns NG.",
						ii);
					return NG;
				}
				if (NG == clp_point2->link_near(clp_point)) {
					pri_funct_err_bttvr(

						"Error : count %d : clp_point2->link_near() returns NG.",
						ii);
					return NG;
				}
				++ii;
			}
			/* 右斜め上につながるものあり */
			if ((OFF == i_right_link_sw) && ((clp_point2->get_i32_xp()) ==
											 (clp_point->get_i32_xp() + 1)) &&
				((clp_point2->get_i32_yp()) ==
				 (clp_point->get_i32_yp() + 1))) {
				/* 双方向リンクする */
				if (NG == clp_point->link_near(clp_point2)) {
					pri_funct_err_bttvr(
						"Error : count %d : clp_point->link_near() returns NG.",
						ii);
					return NG;
				}
				if (NG == clp_point2->link_near(clp_point)) {
					pri_funct_err_bttvr(
						"Error : count %d : clp_point2->link_near() returns NG.",
						ii);
					return NG;
				}
				++ii;
			}
		}
	}

	if (ON == this->_i_pv_sw) {
		pri_funct_msg_ttvr(" link slant %d", ii);
	}

	return OK;
}

int pixel_line_root::_exec04_grouping(pixel_point_root *clp_pixel_point_root)
{
	pixel_point_node *clp_point;
	pixel_line_node *clp_line_before, *clp_line_crnt;
	int32_t ii;

	if (ON == this->_i_mv_sw) {
		pri_funct_msg_ttvr("pixel_line_root::_exec04_grouping()");
	}

	ii = 0;
	clp_line_before = NULL;

	for (
		clp_point = (pixel_point_node *)(clp_pixel_point_root->get_clp_first());
		NULL != clp_point;
		clp_point = (pixel_point_node *)clp_point->get_clp_next()) {
		/* 線分としてリンク済みのは飛ばして次へ */
		if ((NULL != clp_point->get_clp_next_point()) || (NULL != clp_point->get_clp_previous_point())) {
			continue;
		}

		/* 端点をみつけた */
		if ((NULL != clp_point->get_clp_link_near(0)) && (NULL == clp_point->get_clp_link_near(1))) {

			clp_line_crnt = this->_append(clp_line_before);
			if (NULL == clp_line_crnt) {
				pri_funct_err_bttvr(
					"Error : this->_append() returns NULL.");
				return NG;
			}

			/* 線分としてリンクし、両端点をセットする */
			clp_line_crnt->link_line(
				clp_point,
				clp_point->get_clp_link_near(0),
				clp_pixel_point_root->get_i32_count());
			/* リンクした線分は3点以上でないとだめ */
			if (clp_line_crnt->get_i32_point_count() < 3) {
				this->_remove(clp_line_crnt);
			}
			/* 選択が3点以上ならOKなので、次のラインへ */
			else {
				clp_line_before = clp_line_crnt;
			}

			++ii;
		}
	}

	if (ON == this->_i_pv_sw) {
		pri_funct_msg_ttvr(" make %d lines", this->get_i32_count());
	}

	return OK;
}

#include <assert.h> /* assert() */

#include "igs_line_blur.h" // "pri.h" "pixel_line_root.h"

void pixel_line_root::exec05_set_middle(void)
{
	pixel_line_node *clp_line;
	int32_t ii;

	if (ON == this->_i_mv_sw) {
		pri_funct_msg_ttvr("pixel_line_root::exec05_set_middle()");
	}

	for (
		clp_line = (pixel_line_node *)this->get_clp_first(), ii = 0;
		NULL != clp_line;
		clp_line = (pixel_line_node *)clp_line->get_clp_next(), ++ii) {
		assert(ii < this->get_i32_count());

		clp_line->set_middle();
	}
	if (ON == this->_i_pv_sw) {
		pri_funct_msg_ttvr(" set middle point about %d lines", ii);
	}
}

#include <assert.h> /* assert() */

#include "igs_line_blur.h" // "pri.h" "pixel_line_root.h"

void pixel_line_root::exec06_int2double_body(void)
{
	pixel_line_node *clp_line;
	int32_t ii;

	if (ON == this->_i_mv_sw) {
		pri_funct_msg_ttvr("pixel_line_root::exec06_int2double_body()");
	}

	for (
		clp_line = (pixel_line_node *)this->get_clp_first(), ii = 0;
		NULL != clp_line;
		clp_line = (pixel_line_node *)clp_line->get_clp_next(), ++ii) {
		assert(ii < this->get_i32_count());

		clp_line->int2double_body();
	}

	if (ON == this->_i_pv_sw) {
		pri_funct_msg_ttvr(" int to double %d lines", ii);
	}
}
#include <assert.h> /* assert() */

#include "igs_line_blur.h" // "pri.h" "pixel_line_root.h"

void pixel_line_root::exec07_smooth_body(void)
{
	pixel_line_node *clp_line;
	int32_t ii;

	if (ON == this->_i_mv_sw) {
		pri_funct_msg_ttvr("pixel_line_root::exec07_smooth_body()");
	}

	if (ON == this->_i_pv_sw) {
		pri_funct_msg_ttvr(" smooth retry %u", this->_i32_smooth_retry);
	}

	/* カウントダウン表示始め */
	if (ON == this->_i_cv_sw) {
		pri_funct_cv_start(this->get_i32_count());
	}

	for (
		clp_line = (pixel_line_node *)this->get_clp_first(), ii = 0;
		NULL != clp_line;
		clp_line = (pixel_line_node *)clp_line->get_clp_next(), ++ii) {
		assert(ii < this->get_i32_count());

		/* カウントダウン表示中 */
		if (ON == this->_i_cv_sw) {
			pri_funct_cv_run(ii);
		}

		clp_line->smooth_body(this->_i32_smooth_retry);
	}
	/* カウントダウン表示終了 */
	if (ON == this->_i_cv_sw) {
		pri_funct_cv_end();
	}

	if (ON == this->_i_pv_sw) {
		pri_funct_msg_ttvr(" smooth %d lines", ii);
	}
}

void pixel_line_root::exec10_smooth_expand(void)
{
	pixel_line_node *clp_line;
	int32_t ii;

	if (ON == this->_i_mv_sw) {
		pri_funct_msg_ttvr("pixel_line_root::exec10_smooth_expand()");
	}
	if (ON == this->_i_pv_sw) {
		pri_funct_msg_ttvr(" smooth retry %u", this->_i32_smooth_retry);
	}

	/* カウントダウン表示始め */
	if (ON == this->_i_cv_sw) {
		pri_funct_cv_start(this->get_i32_count());
	}

	for (
		clp_line = (pixel_line_node *)this->get_clp_first(), ii = 0;
		NULL != clp_line;
		clp_line = (pixel_line_node *)clp_line->get_clp_next(), ++ii) {
		assert(ii < this->get_i32_count());

		/* カウントダウン表示中 */
		if (ON == this->_i_cv_sw) {
			pri_funct_cv_run(ii);
		}

		clp_line->smooth_expand(this->_i32_smooth_retry);
	}
	/* カウントダウン表示終了 */
	if (ON == this->_i_cv_sw) {
		pri_funct_cv_end();
	}

	if (ON == this->_i_pv_sw) {
		pri_funct_msg_ttvr(" smooth %d lines", ii);
	}
}

#include <assert.h> /* assert() */

#include "igs_line_blur.h" // "pri.h" "pixel_line_root.h"

int pixel_line_root::exec08_expand_lines(pixel_point_root *clp_pixel_point_root)
{
	pixel_line_node *clp_line;
	int32_t ii;

	if (ON == this->_i_mv_sw) {
		pri_funct_msg_ttvr("pixel_line_root::exec08_expand_lines()");
	}

	for (
		clp_line = (pixel_line_node *)this->get_clp_first(), ii = 0;
		NULL != clp_line;
		clp_line = (pixel_line_node *)clp_line->get_clp_next(), ++ii) {
		assert(ii < this->get_i32_count());

		if (OK != clp_line->expand_line(clp_pixel_point_root)) {
			pri_funct_err_bttvr(
				"Error : line number %d : clp_line->expand_line() returns NG",
				ii);
			return NG;
		}
	}

	if (ON == this->_i_pv_sw) {
		pri_funct_msg_ttvr(" expand %d lines", ii);
	}

	return OK;
}

#include <assert.h> /* assert() */

#include "igs_line_blur.h" // "pri.h" "pixel_line_root.h"

double pixel_line_root::_same_way_expand_radian_diff(pixel_point_node *clp_point_middle, pixel_point_node *clp_point_term, pixel_point_node *clp_point_term_expand, pixel_select_same_way_root *clp_select)
{
	double d_xv, d_yv;

	assert(clp_point_middle != clp_point_term);
	assert(clp_point_term != clp_point_term_expand);
	assert(clp_point_term_expand != clp_point_middle);

	/* 自身以外の全端点と比較 */
	clp_select->exec(
		(pixel_line_node *)this->get_clp_first(),
		this->get_i32_count(),
		clp_point_middle,
		clp_point_term,
		clp_point_term_expand);

	if (clp_select->get_i32_count() <= 0) {
		return 0.0;
	}

	/* ベクトルを付加合成 */
	d_xv = 0.0;
	d_yv = 0.0;
	clp_select->get_vector(&d_xv, &d_yv);

	/* ベクトルがゼロの場合
	   --> 何も選択していない
	       --> 角度ゼロを返す
	*/
	if ((0.0 == d_xv) && (0.0 == d_yv)) {
		return 0.0;
	}

	/* 自分自身のベクトル、他より優位生を持たせるため2倍する */
	d_xv += (clp_point_term_expand->get_d_xp_tgt() -
			 clp_point_term->get_d_xp_tgt()) *
			2.0;
	d_yv += (clp_point_term_expand->get_d_yp_tgt() -
			 clp_point_term->get_d_yp_tgt()) *
			2.0;

	/* 元の延長線との角度の差を返す */
	return this->_cl_cal_geom.get_d_radian_by_2_vector(
		clp_point_term_expand->get_d_xp_tgt() -
			clp_point_term->get_d_xp_tgt(),
		clp_point_term_expand->get_d_yp_tgt() -
			clp_point_term->get_d_yp_tgt(),
		d_xv, d_yv);
}

/********************************************************************/

void pixel_line_root::exec09_same_way_expand(pixel_select_same_way_root *clp_select)
{
	pixel_point_node *clp_point_start,
		*clp_point_tmp;
	pixel_line_node *clp_line;
	int32_t ii;
	double d_radian, d_x, d_y;
	int32_t i32_count_one, i32_count_another;

	/* 同方向曲げをしないときはここで抜けてしまう */
	if (ON != this->_i_same_way_exec_sw) {
		return;
	}

	if (ON == this->_i_mv_sw) {
		pri_funct_msg_ttvr("pixel_line_root::exec09_same_way_expand()");
	}
	if (ON == this->_i_pv_sw) {
		pri_funct_msg_ttvr(" select max length %g count %d",
						   clp_select->get_d_length_max(),
						   clp_select->get_i32_count_max());
	}

	/* 各端点(one,anohter)ごと
	   同方向線分となる端点の近点を探し、
	   それらの合成ベクトルによる角度との差を一時保管する */
	for (
		clp_line = (pixel_line_node *)this->get_clp_first(), ii = 0;
		NULL != clp_line;
		clp_line = (pixel_line_node *)clp_line->get_clp_next(), ++ii) {
		assert(ii < this->get_i32_count());

		if (NULL != clp_line->get_clp_link_one_expand()) {
			d_radian = this->_same_way_expand_radian_diff(
				clp_line->get_clp_link_middle(),
				clp_line->get_clp_link_one(),
				clp_line->get_clp_link_one_expand(),
				clp_select);
			clp_line->set_d_same_way_radian_one(d_radian);
		}

		if (NULL != clp_line->get_clp_link_another_expand()) {
			d_radian = this->_same_way_expand_radian_diff(
				clp_line->get_clp_link_middle(),
				clp_line->get_clp_link_another(),
				clp_line->get_clp_link_another_expand(),
				clp_select);
			clp_line->set_d_same_way_radian_another(d_radian);
		}
	}

	/* 各端点(one,anohter)ごと
	   伸ばした線分を、一時保管した角度に曲げる */
	i32_count_one = 0;
	i32_count_another = 0;
	for (
		clp_line = (pixel_line_node *)this->get_clp_first(), ii = 0;
		NULL != clp_line;
		clp_line = (pixel_line_node *)clp_line->get_clp_next(), ++ii) {
		assert(ii < this->get_i32_count());

		if (0.0 < clp_line->get_d_same_way_radian_one()) {
			clp_point_start = clp_line->get_clp_link_one();
			clp_point_tmp = clp_point_start->get_clp_previous_point();
			for (; NULL != clp_point_tmp;
				 clp_point_tmp = clp_point_tmp->get_clp_previous_point()) {
				this->_cl_cal_geom.get_dd_rotate_by_pos(
					clp_point_tmp->get_d_xp_tgt(),
					clp_point_tmp->get_d_yp_tgt(),
					clp_point_start->get_d_xp_tgt(),
					clp_point_start->get_d_yp_tgt(),
					clp_line->get_d_same_way_radian_one(),
					&d_x, &d_y);
				clp_point_tmp->set_d_xp_tgt(d_x);
				clp_point_tmp->set_d_yp_tgt(d_y);
			}
			++i32_count_one;
		}
		if (0.0 < clp_line->get_d_same_way_radian_another()) {
			clp_point_start = clp_line->get_clp_link_another();
			clp_point_tmp = clp_point_start->get_clp_next_point();
			for (; NULL != clp_point_tmp;
				 clp_point_tmp = clp_point_tmp->get_clp_next_point()) {
				this->_cl_cal_geom.get_dd_rotate_by_pos(
					clp_point_tmp->get_d_xp_tgt(),
					clp_point_tmp->get_d_yp_tgt(),
					clp_point_start->get_d_xp_tgt(),
					clp_point_start->get_d_yp_tgt(),
					clp_line->get_d_same_way_radian_another(),
					&d_x, &d_y);
				clp_point_tmp->set_d_xp_tgt(d_x);
				clp_point_tmp->set_d_yp_tgt(d_y);
			}
			++i32_count_another;
		}
	}

	clp_select->mem_free();

	if (ON == this->_i_pv_sw) {
		pri_funct_msg_ttvr(" same way expand  one %d another %d",
						   i32_count_one, i32_count_another);
	}
}

#include <assert.h> /* assert() */

#include "igs_line_blur.h" // "pri.h" "pixel_line_root.h"

void pixel_line_root::exec11_set_bbox(void)
{
	pixel_line_node *clp_line;
	int32_t ii;

	if (ON == this->_i_mv_sw) {
		pri_funct_msg_ttvr("pixel_line_root::exec11_set_bbox()");
	}

	for (
		clp_line = (pixel_line_node *)this->get_clp_first(), ii = 0;
		NULL != clp_line;
		clp_line = (pixel_line_node *)clp_line->get_clp_next(), ++ii) {
		assert(ii < this->get_i32_count());

		clp_line->set_bbox();

		if (0 == ii) {
			this->_d_bbox_x_min = clp_line->get_d_bbox_x_min();
			this->_d_bbox_x_max = clp_line->get_d_bbox_x_max();
			this->_d_bbox_y_min = clp_line->get_d_bbox_y_min();
			this->_d_bbox_y_max = clp_line->get_d_bbox_y_max();
		} else {
			if (clp_line->get_d_bbox_x_min() < this->_d_bbox_x_min) {
				this->_d_bbox_x_min = clp_line->get_d_bbox_x_min();
			} else if (this->_d_bbox_x_max < clp_line->get_d_bbox_x_max()) {
				this->_d_bbox_x_max = clp_line->get_d_bbox_x_max();
			}
			if (clp_line->get_d_bbox_y_min() < this->_d_bbox_y_min) {
				this->_d_bbox_y_min = clp_line->get_d_bbox_y_min();
			} else if (this->_d_bbox_y_max < clp_line->get_d_bbox_y_max()) {
				this->_d_bbox_y_max = clp_line->get_d_bbox_y_max();
			}
		}
	}

	if (ON == this->_i_pv_sw) {
		pri_funct_msg_ttvr(
			" set bbox %d lines : min x %g y %g : max x %g y %g",
			ii,
			this->_d_bbox_x_min, this->_d_bbox_y_min,
			this->_d_bbox_x_max, this->_d_bbox_y_max);
	}
}
#include <assert.h> /* assert() */

#include "igs_line_blur.h" // "pri.h" "pixel_line_root.h"

pixel_line_node *pixel_line_root::_append(pixel_line_node *clp_previous)
{
	pixel_line_node *clp_new;

	clp_new = new pixel_line_node;
	if (NULL == clp_new) {
		pri_funct_err_bttvr("Error : 'new pixel_line_node' returns NULL.");
		return NULL;
	}
	clp_new = (pixel_line_node *)this->push(clp_previous, clp_new);

	return clp_new;
}

void pixel_line_root::_remove(pixel_line_node *clp_old)
{
	assert(NULL != clp_old); /* あってはならないプログラムバグのチェック */
	this->pop(clp_old);
	delete clp_old;
}

void pixel_line_root::mem_free(void)
{
	pixel_line_node *clp_;
	int32_t ii;

	if (NULL != this->get_clp_last()) {
		if (ON == this->_i_mv_sw) {
			pri_funct_msg_ttvr("pixel_line_root::mem_free()");
		}

		ii = 0;
		while (NULL != (clp_ = (pixel_line_node *)this->get_clp_last())) {
			this->_remove(clp_);
			++ii;
		}
		if (ON == this->_i_pv_sw) {
			pri_funct_msg_ttvr("free line node %d", ii);
		}
	}
}

#include "igs_line_blur.h" // "pri.h" "pixel_line_root.h"

int pixel_line_root::save_not_include(pixel_point_root *clp_pixel_point_root, char *cp_fname)
{
	FILE *fp;
	int32_t ii;
	pixel_point_node *clp_point;
	int32_t i32_not_include_line;

	/* ファイル開く */
	fp = fopen(cp_fname, "w");
	if (NULL == fp) {
		pri_funct_err_bttvr("Error : fopen(%s,w) returns NULL", cp_fname);
		return NG;
	}

	/* 全ポイント数保存 */
	if (fprintf(fp, "# all point count %d\n",
				clp_pixel_point_root->get_i32_count()) < 0) {
		pri_funct_err_bttvr("Error : fprintf(# all point count) returns minus");
		fclose(fp);
		return NG;
	}

	i32_not_include_line = 0;
	for (
		clp_point = (pixel_point_node *)(clp_pixel_point_root->get_clp_first()),
	   ii = 0L;
		NULL != clp_point;
		clp_point = (pixel_point_node *)clp_point->get_clp_next(),
	   ++ii) {
		if ((NULL == clp_point->get_clp_next_point()) && (NULL == clp_point->get_clp_previous_point())) {
			if (fprintf(fp, "%d %d\n",
						clp_point->get_i32_xp(),
						clp_point->get_i32_yp()) < 0) {
				pri_funct_err_bttvr("Error : fprintf(%d %d) returns minus",
									clp_point->get_i32_xp(),
									clp_point->get_i32_yp());
				fclose(fp);
				return NG;
			}
			++i32_not_include_line;
		}
	}
	/* 独立ポイント数保存 */
	if (fprintf(fp, "# not include is %d\n",
				i32_not_include_line) < 0) {
		pri_funct_err_bttvr("Error : fprintf(# not include is %d) returns minus",
							i32_not_include_line);
		fclose(fp);
		return NG;
	}

	/* ファイル閉じる */
	fclose(fp);

	return OK;
}

/********************************************************************/

int pixel_line_root::save_lines(char *cp_fname)
{
	FILE *fp;
	int32_t ii;
	pixel_line_node *clp_line;

	/* ファイル開く */
	fp = fopen(cp_fname, "w");
	if (NULL == fp) {
		pri_funct_err_bttvr("Error : fopen(%s,w) returns NULL", cp_fname);
		return NG;
	}

	/* グループ(ライン)数保存 */
	if (fprintf(fp, "# lines count %d\n", this->get_i32_count()) < 0) {
		pri_funct_err_bttvr("Error : fprintf(# group count) returns minus");
		fclose(fp);
		return NG;
	}

	/* データ保存 */
	for (
		clp_line = (pixel_line_node *)this->get_clp_first(),
	   ii = 0;
		NULL != clp_line;
		clp_line = (pixel_line_node *)clp_line->get_clp_next(),
	   ++ii) {
		/* グループ(ライン)番号保存 */
		if (fprintf(fp, "# line number %d  and points count %d\n",
					ii, clp_line->get_i32_point_count()) < 0) {
			pri_funct_err_bttvr(
				"Error : fprintf(# group number %d) returns minus",
				ii);
			fclose(fp);
			return NG;
		}

		/* ポイント保存 */
		if (OK != clp_line->save_line(fp)) {
			pri_funct_err_bttvr(
				"Error : clp_line->save_line() returns NG",
				ii);
			fclose(fp);
			return NG;
		}

		if (NULL != clp_line->get_clp_next()) {
			/* ２行の改行(for gnuplot)保存 */
			if (fprintf(fp, "\n\n") < 0) {
				pri_funct_err_bttvr(
					"Error : fprintf(LF,LF) returns minus");
				fclose(fp);
				return NG;
			}
		}
	}

	/* ファイル閉じる */
	fclose(fp);

	return OK;
}

/********************************************************************/

int pixel_line_root::save_one_point(char *cp_fname)
{
	FILE *fp;
	int32_t ii;
	pixel_line_node *clp_line;

	/* ファイル開く */
	fp = fopen(cp_fname, "w");
	if (NULL == fp) {
		pri_funct_err_bttvr("Error : fopen(%s,w) returns NULL", cp_fname);
		return NG;
	}

	/* グループ(ライン)数保存 */
	if (fprintf(fp, "# line count %d\n", this->get_i32_count()) < 0) {
		pri_funct_err_bttvr("Error : fprintf(# group count) returns minus");
		fclose(fp);
		return NG;
	}

	/* データ保存 */
	for (
		clp_line = (pixel_line_node *)this->get_clp_first(),
	   ii = 0;
		NULL != clp_line;
		clp_line = (pixel_line_node *)clp_line->get_clp_next(),
	   ++ii) {
		/* グループ(ライン)番号保存 */
		if (fprintf(fp, "# line number %d ... start point\n", ii) < 0) {
			pri_funct_err_bttvr(
				"Error : fprintf(# group number %d) returns minus",
				ii);
			fclose(fp);
			return NG;
		}

		/* ポイント保存 */
		if (OK != clp_line->save_one_point(fp)) {
			pri_funct_err_bttvr(
				"Error : clp_line->save_one_point() returns NG",
				ii);
			fclose(fp);
			return NG;
		}
	}

	/* ファイル閉じる */
	fclose(fp);

	return OK;
}

int pixel_line_root::save_middle_point(char *cp_fname)
{
	FILE *fp;
	int32_t ii;
	pixel_line_node *clp_line;

	/* ファイル開く */
	fp = fopen(cp_fname, "w");
	if (NULL == fp) {
		pri_funct_err_bttvr("Error : fopen(%s,w) returns NULL", cp_fname);
		return NG;
	}

	/* グループ(ライン)数保存 */
	if (fprintf(fp, "# line count %d\n", this->get_i32_count()) < 0) {
		pri_funct_err_bttvr("Error : fprintf(# group count) returns minus");
		fclose(fp);
		return NG;
	}

	/* データ保存 */
	for (
		clp_line = (pixel_line_node *)this->get_clp_first(),
	   ii = 0;
		NULL != clp_line;
		clp_line = (pixel_line_node *)clp_line->get_clp_next(),
	   ++ii) {
		/* グループ(ライン)番号保存 */
		if (fprintf(fp, "# line number %d ... start point\n", ii) < 0) {
			pri_funct_err_bttvr(
				"Error : fprintf(# group number %d) returns minus",
				ii);
			fclose(fp);
			return NG;
		}

		/* ポイント保存 */
		if (OK != clp_line->save_middle_point(fp)) {
			pri_funct_err_bttvr(
				"Error : clp_line->save_middle_point() returns NG",
				ii);
			fclose(fp);
			return NG;
		}
	}

	/* ファイル閉じる */
	fclose(fp);

	return OK;
}

int pixel_line_root::save_another_point(char *cp_fname)
{
	FILE *fp;
	int32_t ii;
	pixel_line_node *clp_line;

	/* ファイル開く */
	fp = fopen(cp_fname, "w");
	if (NULL == fp) {
		pri_funct_err_bttvr("Error : fopen(%s,w) returns NULL", cp_fname);
		return NG;
	}

	/* グループ(ライン)数保存 */
	if (fprintf(fp, "# line count %d\n", this->get_i32_count()) < 0) {
		pri_funct_err_bttvr("Error : fprintf(# group count) returns minus");
		fclose(fp);
		return NG;
	}

	/* データ保存 */
	for (
		clp_line = (pixel_line_node *)this->get_clp_first(),
	   ii = 0;
		NULL != clp_line;
		clp_line = (pixel_line_node *)clp_line->get_clp_next(),
	   ++ii) {
		/* グループ(ライン)番号保存 */
		if (fprintf(fp, "# line number %d ... end point\n", ii) < 0) {
			pri_funct_err_bttvr(
				"Error : fprintf(# group number %d) returns minus",
				ii);
			fclose(fp);
			return NG;
		}

		/* ポイント保存 */
		if (OK != clp_line->save_another_point(fp)) {
			pri_funct_err_bttvr(
				"Error : clp_line->save_another_point() returns NG",
				ii);
			fclose(fp);
			return NG;
		}
	}

	/* ファイル閉じる */
	fclose(fp);

	return OK;
}

#include "igs_line_blur.h" // "pri.h" "pixel_line_root.h"

int pixel_line_root::save_expand_lines(char *cp_fname)
{
	FILE *fp;
	int32_t ii;
	pixel_line_node *clp_line;

	/* ファイル開く */
	fp = fopen(cp_fname, "w");
	if (NULL == fp) {
		pri_funct_err_bttvr("Error : fopen(%s,w) returns NULL", cp_fname);
		return NG;
	}

	/* グループ(ライン)数保存 */
	if (fprintf(fp, "# lines count %d\n", this->get_i32_count()) < 0) {
		pri_funct_err_bttvr("Error : fprintf(# group count) returns minus");
		fclose(fp);
		return NG;
	}

	/* データ保存 */
	for (
		clp_line = (pixel_line_node *)this->get_clp_first(),
	   ii = 0;
		NULL != clp_line;
		clp_line = (pixel_line_node *)clp_line->get_clp_next(),
	   ++ii) {
		/* グループ(ライン)番号保存 */
		if (fprintf(fp, "# line number %d  and points count %d\n",
					ii, clp_line->get_i32_point_count()) < 0) {
			pri_funct_err_bttvr(
				"Error : fprintf(# group number %d) returns minus",
				ii);
			fclose(fp);
			return NG;
		}

		/* ポイント保存 */
		if (OK != clp_line->save_expand_line(fp)) {
			pri_funct_err_bttvr(
				"Error : clp_line->save_line() returns NG",
				ii);
			fclose(fp);
			return NG;
		}

		if (NULL != clp_line->get_clp_next()) {
			/* ２行の改行(for gnuplot)保存 */
			if (fprintf(fp, "\n\n") < 0) {
				pri_funct_err_bttvr(
					"Error : fprintf(LF,LF) returns minus");
				fclose(fp);
				return NG;
			}
		}
	}

	/* ファイル閉じる */
	fclose(fp);

	return OK;
}

/********************************************************************/

int pixel_line_root::save_one_expand_point(char *cp_fname)
{
	FILE *fp;
	int32_t ii;
	pixel_line_node *clp_line;

	/* ファイル開く */
	fp = fopen(cp_fname, "w");
	if (NULL == fp) {
		pri_funct_err_bttvr("Error : fopen(%s,w) returns NULL", cp_fname);
		return NG;
	}

	/* グループ(ライン)数保存 */
	if (fprintf(fp, "# line count %d\n", this->get_i32_count()) < 0) {
		pri_funct_err_bttvr("Error : fprintf(# group count) returns minus");
		fclose(fp);
		return NG;
	}

	/* データ保存 */
	for (
		clp_line = (pixel_line_node *)this->get_clp_first(),
	   ii = 0;
		NULL != clp_line;
		clp_line = (pixel_line_node *)clp_line->get_clp_next(),
	   ++ii) {
		/* グループ(ライン)番号保存 */
		if (fprintf(fp, "# line number %d ... one expand point\n", ii) < 0) {
			pri_funct_err_bttvr(
				"Error : fprintf(# group number %d) returns minus",
				ii);
			fclose(fp);
			return NG;
		}

		/* ポイント保存 */
		if (OK != clp_line->save_one_expand_point(fp)) {
			pri_funct_err_bttvr(
				"Error : clp_line->save_one_expand_point() returns NG",
				ii);
			fclose(fp);
			return NG;
		}
	}

	/* ファイル閉じる */
	fclose(fp);

	return OK;
}

int pixel_line_root::save_another_expand_point(char *cp_fname)
{
	FILE *fp;
	int32_t ii;
	pixel_line_node *clp_line;

	/* ファイル開く */
	fp = fopen(cp_fname, "w");
	if (NULL == fp) {
		pri_funct_err_bttvr("Error : fopen(%s,w) returns NULL", cp_fname);
		return NG;
	}

	/* グループ(ライン)数保存 */
	if (fprintf(fp, "# line count %d\n", this->get_i32_count()) < 0) {
		pri_funct_err_bttvr("Error : fprintf(# group count) returns minus");
		fclose(fp);
		return NG;
	}

	/* データ保存 */
	for (
		clp_line = (pixel_line_node *)this->get_clp_first(),
	   ii = 0;
		NULL != clp_line;
		clp_line = (pixel_line_node *)clp_line->get_clp_next(),
	   ++ii) {
		/* グループ(ライン)番号保存 */
		if (fprintf(fp, "# line number %d ... another expand point\n", ii) < 0) {
			pri_funct_err_bttvr(
				"Error : fprintf(# group number %d) returns minus",
				ii);
			fclose(fp);
			return NG;
		}

		/* ポイント保存 */
		if (OK != clp_line->save_another_expand_point(fp)) {
			pri_funct_err_bttvr(
				"Error : clp_line->save_another_expand_point() returns NG",
				ii);
			fclose(fp);
			return NG;
		}
	}

	/* ファイル閉じる */
	fclose(fp);

	return OK;
}

#include "igs_line_blur.h" // "pri.h" "pixel_line_root.h"

int pixel_line_root::save_expand_vector(char *cp_fname)
{
	FILE *fp;
	int32_t ii;
	pixel_line_node *clp_line;

	/* ファイル開く */
	fp = fopen(cp_fname, "w");
	if (NULL == fp) {
		pri_funct_err_bttvr("Error : fopen(%s,w) returns NULL", cp_fname);
		return NG;
	}

	/* グループ(ライン)数保存 */
	if (fprintf(fp, "# lines count %d\n", this->get_i32_count()) < 0) {
		pri_funct_err_bttvr(
			"Error : fprintf(# group count) returns minus");
		fclose(fp);
		return NG;
	}

	/* データ保存 */
	for (
		clp_line = (pixel_line_node *)this->get_clp_first(),
	   ii = 0;
		NULL != clp_line;
		clp_line = (pixel_line_node *)clp_line->get_clp_next(),
	   ++ii) {
		/* グループ(ライン)番号保存 */
		if (fprintf(fp, "# line number %d  and points count %d\n",
					ii, clp_line->get_i32_point_count()) < 0) {
			pri_funct_err_bttvr(
				"Error : fprintf(# line number %d) returns minus", ii);
			fclose(fp);
			return NG;
		}

		/* ポイント保存 */
		if (OK != clp_line->save_expand_vector(fp)) {
			pri_funct_err_bttvr(
				"Error : clp_line->save_expand_vector() returns NG");
			fclose(fp);
			return NG;
		}
	}

	/* ファイル閉じる */
	fclose(fp);

	return OK;
}

#ifndef _pixel_select_curve_blur_h_
#define _pixel_select_curve_blur_h_

#include <stdio.h>

#include "igs_line_blur.h" // "list_root.h" "calculator_geometry.h" "pixel_point_root.h" "pixel_line_root.h"

class pixel_select_curve_blur_node : public list_node
{
public:
	pixel_select_curve_blur_node(void)
	{
		this->clp_line = NULL;
		this->clp_start_point = NULL;
		this->clp_near_point = NULL;
		this->i32_near_point_pos = 0;
		this->d_length = -1.0;
		this->i_reverse_sw = OFF;
	}
	void copy(pixel_select_curve_blur_node *clp)
	{
		this->clp_line = clp->clp_line;
		this->clp_start_point = clp->clp_start_point;
		this->clp_near_point = clp->clp_near_point;
		this->i32_near_point_pos = clp->i32_near_point_pos;
		this->d_length = clp->d_length;
		this->i_reverse_sw = clp->i_reverse_sw;
	}
	pixel_line_node *clp_line;
	pixel_point_node *clp_start_point;
	pixel_point_node *clp_near_point;
	int32_t i32_near_point_pos;
	double d_length;
	int i_reverse_sw;

private:
};

class pixel_select_curve_blur_root : public list_root
{
public:
	pixel_select_curve_blur_root(void)
	{
		this->_i_mv_sw = OFF;
		this->_i_cv_sw = OFF;
		this->_i_pv_sw = OFF;

		this->_i32_count_max = 4;
		this->_d_length_max = 160;
	}
	~pixel_select_curve_blur_root(void)
	{
		this->mem_free();
	}
	void set_i_mv_sw(int sw) { this->_i_mv_sw = sw; }
	void set_i_cv_sw(int sw) { this->_i_cv_sw = sw; }
	void set_i_pv_sw(int sw) { this->_i_pv_sw = sw; }

	int32_t get_i32_count_max(void) { return this->_i32_count_max; }
	void set_i32_count_max(int32_t ii) { this->_i32_count_max = ii; }
	double get_d_length_max(void) { return this->_d_length_max; }
	void set_d_length_max(double dd) { this->_d_length_max = dd; }

	void exec(double d_xp, double d_yp, pixel_line_root *clp_pixel_line_root, int32_t i32_blur_count, double d_effect_length_radius);
	/******void exec( double d_xp, double d_yp, pixel_line_node *clp_line_first, int32_t i32_count, int32_t i32_blur_count, double d_effect_length_radius );******/
	int get_line(int32_t i32_blur_count, double *dp_xv, double *dp_yv);

	int save(double d_xp, double d_yp, int32_t i32_blur_count, char *cp_fname);

	void mem_free(void);

private:
	int _i_mv_sw,
		_i_cv_sw,
		_i_pv_sw;

	int32_t _i32_count_max;
	double _d_length_max;

	calculator_geometry _cl_cal_geom;

	pixel_select_curve_blur_node *_append(pixel_select_curve_blur_node *clp_previous);
	void _remove(pixel_select_curve_blur_node *clp_old);
	int _sort_append(pixel_select_curve_blur_node *clp_src);

	pixel_point_node *_get_next_points(pixel_point_node *clp_point, int32_t i32_count);
	pixel_point_node *_get_prev_points(pixel_point_node *clp_point, int32_t i32_count);

	double _get_line_accum_count(pixel_select_curve_blur_node *clp_select, int32_t i32_blur_count);
};

#endif				/* !_pixel_select_curve_blur_h_ */
#include <stdlib.h> /* exit() */
#include <math.h>   /* M_PI */

/* WindowsではM_PIが見つからない */
#if defined _WIN32
#define M_PI 3.14159265358979323846
#endif

#include <assert.h>		   /* assert() */
#include "igs_line_blur.h" // "pri.h" "pixel_line_root.h" "pixel_select_curve_blur.h"

pixel_select_curve_blur_node *pixel_select_curve_blur_root::_append(pixel_select_curve_blur_node *clp_previous)
{
	pixel_select_curve_blur_node *clp_new;

	clp_new = new pixel_select_curve_blur_node;
	if (NULL == clp_new) {
		pri_funct_err_bttvr(
			"Error : 'new pixel_select_curve_blur_node' returns NULL.");
		return NULL;
	}
	clp_new = (pixel_select_curve_blur_node *)this->push(clp_previous, clp_new);

	return clp_new;
}

void pixel_select_curve_blur_root::_remove(pixel_select_curve_blur_node *clp_old)
{
	assert(NULL != clp_old); /* あってはならないプログラムバグのチェック */
	this->pop(clp_old);
	delete clp_old;
}

void pixel_select_curve_blur_root::mem_free(void)
{
	int32_t ii;

	if (NULL != this->get_clp_last()) {
		if (ON == this->_i_mv_sw) {
			pri_funct_msg_ttvr(
				"pixel_select_curve_blur_root::mem_free()");
		}

		ii = 0;
		while (NULL != this->get_clp_last()) {
			this->_remove(
				(pixel_select_curve_blur_node *)this->get_clp_last());
			++ii;
		}
		if (ON == this->_i_pv_sw) {
			pri_funct_msg_ttvr("free select curve blur node %d", ii);
		}
	}
}

/********************************************************************/

/* lengthの値を見て、距離の近い(小さい値)順になるように追加する
	すべての選択リストより小さいなら先頭
	すべての選択リストより大きいなら最後に追加
	値から選択リストの間に入るならそこに追加
*/
int pixel_select_curve_blur_root::_sort_append(pixel_select_curve_blur_node *clp_src)
{
	pixel_select_curve_blur_node *clp_prev;
	int32_t ii;

	assert(0.0 < clp_src->d_length);

	/* リストの最後から距離を比較し、より小さい場所 */
	clp_prev = (pixel_select_curve_blur_node *)this->get_clp_last();
	for (ii = 0; NULL != clp_prev; ++ii,
		clp_prev = (pixel_select_curve_blur_node *)clp_prev->get_clp_previous()) {
		assert(ii < this->get_i32_count());
		if (clp_prev->d_length < clp_src->d_length) {
			break;
		}
	}

	/* 示した場所の後ろに追加 */
	clp_prev = this->_append(clp_prev);
	if (NULL == clp_prev) {
		pri_funct_err_bttvr(
			"Error : this->_append(clp_prev) returns NULL.");
		return NG;
	}

	/* 情報をコピー */
	clp_prev->copy(clp_src);

	return OK;
}

/********************************************************************/
#define NOT_USE_PARAMETER_VAL (-10000.0)

/* 同方向で、一定の距離の範囲にある線とその位置を選択する */
void pixel_select_curve_blur_root::exec(double d_xp, double d_yp, pixel_line_root *clp_pixel_line_root, int32_t i32_blur_count, double d_effect_area_radius)
{
	pixel_line_node *clp_line;
	int32_t ii,
		i32_pos;
	int i_reverse_sw;
	pixel_point_node *clp_near_point, *clp_start_point;
	double d_length,
		d_radius, d_radius_1st;
	pixel_select_curve_blur_node cl_select;

	/* 選択リストをクリア */
	this->mem_free();

	/* 全体の bbox を見る */
	if ((d_xp < (clp_pixel_line_root->get_d_bbox_x_min() - d_effect_area_radius)) || ((clp_pixel_line_root->get_d_bbox_x_max() + d_effect_area_radius) < d_xp) || (d_yp < (clp_pixel_line_root->get_d_bbox_y_min() - d_effect_area_radius)) || ((clp_pixel_line_root->get_d_bbox_y_max() + d_effect_area_radius) < d_yp)) {
		return;
	}

	d_radius_1st = NOT_USE_PARAMETER_VAL;

	/* 全ライン */
	ii = 0;
	for (
		clp_line = (pixel_line_node *)clp_pixel_line_root->get_clp_first();
		NULL != clp_line;
		clp_line = (pixel_line_node *)clp_line->get_clp_next(), ++ii) {
		/* 無限ループ? */
		assert(ii < clp_pixel_line_root->get_i32_count());
		/* 選択してない? */
		assert(NULL != clp_line->get_clp_link_middle());
		assert(NULL != clp_line->get_clp_link_one());
		assert(NULL != clp_line->get_clp_link_another());
		/* 線が短すぎ? */
		assert(clp_line->get_clp_link_middle() !=
			   clp_line->get_clp_link_one());
		assert(clp_line->get_clp_link_middle() !=
			   clp_line->get_clp_link_another());

		/* 各ラインの bbox を見る */
		if ((d_xp < (clp_line->get_d_bbox_x_min() - d_effect_area_radius)) || ((clp_line->get_d_bbox_x_max() + d_effect_area_radius) < d_xp) || (d_yp < (clp_line->get_d_bbox_y_min() - d_effect_area_radius)) || ((clp_line->get_d_bbox_y_max() + d_effect_area_radius) < d_yp)) {
			continue;
		}

		/* ラインとの近点を探す */
		clp_line->get_near_point(d_xp, d_yp,
								 &i32_pos, &clp_near_point, &d_length);

		/* 距離が遠すぎる */
		if (this->_d_length_max < d_length) {
			continue;
		}

		/* 線分の端で必要な長さがとれない
			画面の端の場合、問題あり!!!!!!
		 */
		if ((i32_pos < (i32_blur_count / 2)) || ((clp_line->get_i32_point_count() - 1 - i32_pos) < (i32_blur_count / 2))) {
			continue;
		}

		/* 逆方向スイッチoffしておく */
		i_reverse_sw = OFF;

		/* 対応する部分のスタート点 */
		clp_start_point = clp_line->get_prev_point_by_count(
			clp_near_point, i32_blur_count / 2);

		/* 近点からスタート点へのベクトル角度 */
		d_radius = this->_cl_cal_geom.get_d_radian(
			clp_start_point->get_d_xp_tgt() - clp_near_point->get_d_xp_tgt(),
			clp_start_point->get_d_yp_tgt() - clp_near_point->get_d_yp_tgt());

		/* 1番始めは線分の方向を記憶しておく */
		if (NOT_USE_PARAMETER_VAL == d_radius_1st) {
			d_radius_1st = d_radius;
		}
		/* 2番目以後は記憶した方向との比較し
		   線の対応部分の、方向と開始点を最設定 */
		else {
			/* 1番目の方向との角度差を出す */
			if (d_radius < d_radius_1st) {
				d_radius += M_PI * 2.0;
			}
			d_radius -= d_radius_1st;
			/* 反対方向なら、逆方向スイッチをonし、スタート点を最設定 */
			if ((M_PI / 2.0 < d_radius) && (d_radius < M_PI * 3.0 / 2.0)) {
				clp_start_point = clp_line->get_next_point_by_count(
					clp_near_point, i32_blur_count / 2);
				i_reverse_sw = ON;
			}
		}

		/* 選択リストに追加(近い順に) */
		cl_select.clp_line = clp_line;
		cl_select.clp_start_point = clp_start_point;
		cl_select.clp_near_point = clp_near_point;
		cl_select.i32_near_point_pos = i32_pos;
		cl_select.d_length = d_length;
		cl_select.i_reverse_sw = i_reverse_sw;
		/* 登録 */
		this->_sort_append(&cl_select);
	}
	/*****if ((NULL != this->get_clp_first()) && (510.0 < d_xp) && (100.0 < d_yp)) {
//if (NULL != this->get_clp_first()) {
 this->save( d_xp, d_yp, i32_blur_count,"tmp19_select.txt" ); exit(1); }
******/
}

/********************************************************************/

double pixel_select_curve_blur_root::_get_line_accum_count(pixel_select_curve_blur_node *clp_select, int32_t i32_blur_count)
{
	int32_t i32_pos, i32_count;
	double d_ratio;

	i32_pos = clp_select->i32_near_point_pos - (i32_blur_count / 2);
	i32_count = clp_select->clp_line->get_i32_point_count() - (i32_blur_count / 2) - (i32_blur_count / 2);

	assert(0 <= i32_pos);
	assert(0 < i32_count);
	assert(i32_pos < i32_count);

	/**if (i32_pos < (i32_count/2)) {
		d_ratio = (double)(i32_pos+1) / (i32_count/2);
	} else {
		d_ratio = (double)(i32_count-i32_pos) / (i32_count/2);
	}**/

	/* ゼロから１の間の値(ゼロと１はだめ) */
	d_ratio = ((double)i32_pos + 0.5) / (double)i32_count;
	/* sinカーブ
		--> 0より大きく1より小さい(d_ratio)
		--> 0-(PI/2)より大きく、(2*PI-(PI/2))より小さい
		--> sin()
		--> -1より大きく1より小さい
		-->  0より大きく2より小さい */
	d_ratio = 1.0 + sin(2.0 * M_PI * d_ratio - M_PI / 2);

	assert(0.0 < d_ratio);

	return d_ratio;
}

/* 距離の近い順(lengthの小さい値)に指定個数分のベクトルを合成して返す
	先頭から
	this->_i32_count_max数個分
	のベクトルを合成する
*/
int pixel_select_curve_blur_root::get_line(int32_t i32_blur_count, double *dp_xv, double *dp_yv)
{
	int32_t ii, jj, i32_count;
	double d_xp, d_yp;
	double d_ratio, d_accum;
	pixel_select_curve_blur_node *clp_select;
	pixel_point_node *clp_point;

	for (ii = 0; ii < i32_blur_count; ++ii) {
		dp_xv[ii] = 0.0;
		dp_yv[ii] = 0.0;
	}

	i32_count = 0;
	d_accum = 0.0;

	clp_select = (pixel_select_curve_blur_node *)this->get_clp_first();
	for (ii = 0; (NULL != clp_select) && (ii < this->_i32_count_max); ++ii,
		++i32_count,
		clp_select = (pixel_select_curve_blur_node *)clp_select->get_clp_next()) {
		assert(ii < this->get_i32_count());
		assert(NULL != clp_select->clp_line);
		assert(NULL != clp_select->clp_start_point);
		assert(NULL != clp_select->clp_near_point);

		d_ratio = this->_get_line_accum_count(clp_select, i32_blur_count);
		d_accum += d_ratio;

		d_xp = clp_select->clp_near_point->get_d_xp_tgt();
		d_yp = clp_select->clp_near_point->get_d_yp_tgt();
		clp_point = clp_select->clp_start_point;

		/* 線分合成(いきなりの変化をへらす) */
		if (OFF == clp_select->i_reverse_sw) {
			for (jj = 0; jj < i32_blur_count; ++jj) {
				dp_xv[jj] += (clp_point->get_d_xp_tgt() - d_xp) * d_ratio;
				dp_yv[jj] += (clp_point->get_d_yp_tgt() - d_yp) * d_ratio;
				clp_point = clp_point->get_clp_next_point();
			}
		} else {
			for (jj = 0; jj < i32_blur_count; ++jj) {
				dp_xv[jj] += (clp_point->get_d_xp_tgt() - d_xp) * d_ratio;
				dp_yv[jj] += (clp_point->get_d_yp_tgt() - d_yp) * d_ratio;
				clp_point = clp_point->get_clp_previous_point();
			}
		}
	}

	if (i32_count <= 0) {
		return NG;
	}

	assert(0.0 < d_accum);
	for (ii = 0; ii < i32_blur_count; ++ii) {
		dp_xv[ii] /= d_accum;
		dp_yv[ii] /= d_accum;
	}

	return OK;
}

#include <assert.h> /* assert() */

#include "igs_line_blur.h" // "pri.h" "pixel_select_curve_blur.h"

int pixel_select_curve_blur_root::save(double d_xp, double d_yp, int32_t i32_blur_count, char *cp_fname)
{
	FILE *fp;
	int32_t ii, jj;
	pixel_select_curve_blur_node *clp_loop;
	pixel_point_node *clp_point;

	/* ファイル開く */
	fp = fopen(cp_fname, "w");
	if (NULL == fp) {
		pri_funct_err_bttvr("Error : fopen(%s,w) returns NULL", cp_fname);
		return NG;
	}

	/* ライン長さ(ポイント数)保存 */
	if (fprintf(fp, "# blur point count %d\n",
				i32_blur_count) < 0) {
		pri_funct_err_bttvr(
			"Error : fprintf(# blur point count) returns minus");
		fclose(fp);
		return NG;
	}
	/* 選択数保存 */
	if (fprintf(fp, "# select count %d in %d\n",
				this->_i32_count_max, this->get_i32_count()) < 0) {
		pri_funct_err_bttvr(
			"Error : fprintf(# select count) returns minus");
		fclose(fp);
		return NG;
	}

	clp_loop = (pixel_select_curve_blur_node *)this->get_clp_first();
	for (ii = 0; (NULL != clp_loop) && (ii < this->_i32_count_max); ++ii,
		clp_loop = (pixel_select_curve_blur_node *)clp_loop->get_clp_next()) {
		assert(ii < this->get_i32_count());
		assert(NULL != clp_loop);
		assert(NULL != clp_loop->clp_line);
		assert(NULL != clp_loop->clp_start_point);
		assert(NULL != clp_loop->clp_near_point);

		/* グループ(ライン)番号保存 */
		if (fprintf(fp, "# selct number %d : reverse sw is %s\n",
					ii, (OFF == clp_loop->i_reverse_sw) ? "off" : "on") < 0) {
			pri_funct_err_bttvr(
				"Error : fprintf(# line number %d) returns minus", ii);
			fclose(fp);
			return NG;
		}

		clp_point = clp_loop->clp_near_point;
		/* ピクセル位置から近点位置保存 */
		if (fprintf(fp, "%g %g\n",
					clp_point->get_d_xp_tgt(),
					clp_point->get_d_yp_tgt()) < 0) {
			pri_funct_err_bttvr(
				"Error : fprintf(near point) returns minus");
			fclose(fp);
			return NG;
		}
		if (fprintf(fp, "%g %g\n\n\n", d_xp, d_yp) < 0) {
			pri_funct_err_bttvr(
				"Error : fprintf(pixel point) returns minus");
			fclose(fp);
			return NG;
		}

		clp_point = clp_loop->clp_start_point;
		if (OFF == clp_loop->i_reverse_sw) {
			for (jj = 0; jj < i32_blur_count; ++jj) {
				/* グループ(ライン)番号保存 */
				if (fprintf(fp, "%g %g\n",
							clp_point->get_d_xp_tgt(),
							clp_point->get_d_yp_tgt()) < 0) {
					pri_funct_err_bttvr(
						"Error : fprintf(x,y) returns minus");
					fclose(fp);
					return NG;
				}
				clp_point = clp_point->get_clp_next_point();
			}
		} else {
			for (jj = 0; jj < i32_blur_count; ++jj) {
				/* グループ(ライン)番号保存 */
				if (fprintf(fp, "%g %g\n",
							clp_point->get_d_xp_tgt(),
							clp_point->get_d_yp_tgt()) < 0) {
					pri_funct_err_bttvr(
						"Error : fprintf(x,y) returns minus");
					fclose(fp);
					return NG;
				}
				clp_point = clp_point->get_clp_previous_point();
			}
		}
		if (fprintf(fp, "\n\n") < 0) {
			pri_funct_err_bttvr(
				"Error : fprintf(EnterEnter) returns minus");
			fclose(fp);
			return NG;
		}
	}

	/* ファイル閉じる */
	fclose(fp);

	return OK;
}

#ifndef __thinnest_ui16_image_h__
#define __thinnest_ui16_image_h__

#include <stdio.h>

/* Windowsではstdint.hが見つからない */
#if defined _WIN32
typedef int int32_t;
typedef unsigned short uint16_t;
#else
#include <stdint.h> /* for int32_t, uint16_t */
#endif

#ifndef UINT16_MAX
#define UINT16_MAX (65535)
#endif

#ifndef OK
#define OK (0)
#endif
#ifndef NG
#define NG (-1)
#endif

#ifndef ON
#define ON (1)
#endif
#ifndef OFF
#define OFF (0)
#endif

class thinnest_ui16_image
{
public:
	thinnest_ui16_image()
	{
		this->_i_mv_sw = OFF;
		this->_i_pv_sw = OFF;
		this->_i_cv_sw = OFF;

		this->_i32_xs = 0;
		this->_i32_ys = 0;
		this->_i32_xd = 1;
		this->_i32_yd = 1;
		this->_ui16_threshold = (uint16_t)0x4000; /* =64*256 */
		this->_i32_exec_loop_count = 100;

		this->_ui16p_channel[0] = NULL;
		this->_ui16p_channel[1] = NULL;
		this->memory_free_this_ = NULL;
	}
	~thinnest_ui16_image()
	{
		this->mem_free();
	}

	/* パラメータ設定 */
	void set_i_mv_sw(int ii) { this->_i_mv_sw = ii; }
	void set_i_pv_sw(int ii) { this->_i_pv_sw = ii; }
	void set_i_cv_sw(int ii) { this->_i_cv_sw = ii; }

	void set_i32_xs(int32_t ii) { this->_i32_xs = ii; }
	void set_i32_ys(int32_t ii) { this->_i32_ys = ii; }
	void set_i32_xd(int32_t ii) { this->_i32_xd = ii; }
	void set_i32_yd(int32_t ii) { this->_i32_yd = ii; }
	void set_ui16_threshold(uint16_t ii) { this->_ui16_threshold = ii; }
	void set_i32_exec_loop_count(int32_t ii) { this->_i32_exec_loop_count = ii; }

	/* パラメータを得る */
	int32_t get_i32_xs(void) { return this->_i32_xs; }
	int32_t get_i32_ys(void) { return this->_i32_ys; }
	int32_t get_i32_xd(void) { return this->_i32_xd; }
	int32_t get_i32_yd(void) { return this->_i32_yd; }
	uint16_t get_ui16_threshold(void) { return this->_ui16_threshold; }
	int32_t get_i32_exec_loop_count(void) { return this->_i32_exec_loop_count; }

	/* メモリ確保 */
	int mem_alloc(void);

	/* 実行 */
	int32_t exec01_fill_noise_pixel(void);
	void exec02_scale_add_edge_pixel(void);
	void exec03_scale_liner(void);
	void exec04_bw(void);
	int exec05_thin(void);

	/* メモリへのポインターを得る */
	uint16_t *get_ui16p_src_channel(void) { return this->_ui16p_channel[0]; }
	uint16_t *get_ui16p_tgt_channel(void) { return this->_ui16p_channel[1]; }

	/* メモリ開放 */
	void mem_free(void);

private:
	int _i_mv_sw, /* Method    Verbose */
		_i_pv_sw, /* Parameter Verbose */
		_i_cv_sw; /* Counter   Verbose */

	int32_t _i32_xs, _i32_ys,	 /* Size */
		_i32_xd, _i32_yd;		  /* sub Divide */
	uint16_t _ui16_threshold;	 /* bw_by_threshold */
	int32_t _i32_exec_loop_count; /* thinnest root */

	uint16_t *memory_free_this_;
	uint16_t *_ui16p_channel[2];

	void _swap_channel(void)
	{
		uint16_t *_ui16p_tmp;

		_ui16p_tmp = this->_ui16p_channel[0];
		this->_ui16p_channel[0] = this->_ui16p_channel[1];
		this->_ui16p_channel[1] = _ui16p_tmp;
	}

	int32_t _exec01_fill_noise_pixel_scanline(uint16_t *ui16p_y1, uint16_t *ui16p_y2, uint16_t *ui16p_y3);
	int32_t _exec01_fill_noise_pixel_pixel(uint16_t *ui16p_x1a, uint16_t *ui16p_x1b, uint16_t *ui16p_x1c, uint16_t *ui16p_x2a, uint16_t *ui16p_x2b, uint16_t *ui16p_x2c, uint16_t *ui16p_x3a, uint16_t *ui16p_x3b, uint16_t *ui16p_x3c);

	int32_t _one_side_thinner(void);
	int32_t _one_side_thinner_scanline(uint16_t *ui16p_src_y1, uint16_t *ui16p_src_y2, uint16_t *ui16p_src_y3, uint16_t *ui16p_tgt);
	int32_t _one_side_thinner_pixel(uint16_t *ui16p_src_x1a, uint16_t *ui16p_src_x1b, uint16_t *ui16p_src_x1c, uint16_t *ui16p_src_x2a, uint16_t *ui16p_src_x2b, uint16_t *ui16p_src_x2c, uint16_t *ui16p_src_x3a, uint16_t *ui16p_src_x3b, uint16_t *ui16p_src_x3c, uint16_t *ui16p_tgt);

	void _rot90_by_clockwork(void);
};

#endif /* !__thinnest_ui16_image_h__ */

#include <assert.h> /* assert() */

#include "igs_line_blur.h" // "pri.h" "thinnest_ui16_image.h"

/* 閾値より大きい白線上の独立点を潰す */
int32_t thinnest_ui16_image::exec01_fill_noise_pixel(void)
{
	uint16_t *ui16p_y1, *ui16p_y2, *ui16p_y3;
	int32_t yy,
		i32_fill_count;

	/* 処理ごとのメッセージ */
	if (ON == this->_i_mv_sw) {
		pri_funct_msg_ttvr("thinnest_ui16_image::exec01_fill_noise_pixel()");
	}
	/* カウントダウン表示始め */
	if (ON == this->_i_cv_sw) {
		pri_funct_cv_start(this->_i32_ys);
	}

	/* 始めのスキャンライン位置(画像から外れているものはNULLをいれる) */
	ui16p_y1 = this->get_ui16p_src_channel() + this->_i32_xs;
	ui16p_y2 = this->get_ui16p_src_channel();
	ui16p_y3 = NULL;

	/* 画像を縦にループ */
	i32_fill_count = 0;
	for (yy = 0; yy < this->_i32_ys; ++yy) {

		/* カウントダウン表示中 */
		if (ON == this->_i_cv_sw) {
			pri_funct_cv_run(yy);
		}

		/* スキャンライン(とその前後のスキャンライン)毎の処理 */
		i32_fill_count += this->_exec01_fill_noise_pixel_scanline(
			ui16p_y1, ui16p_y2, ui16p_y3);

		/* 次のスキャンラインへ進める(画像から外れたらNULLをいれる) */
		ui16p_y3 = ui16p_y2;
		ui16p_y2 = ui16p_y1;
		if ((yy + 1) < (this->_i32_ys)) {
			ui16p_y1 += this->_i32_xs;
		} else
			ui16p_y1 = NULL;
	}

	/* カウントダウン表示終了 */
	if (ON == this->_i_cv_sw) {
		pri_funct_cv_end();
	}

	/* 終了 */
	return i32_fill_count;
}

int32_t thinnest_ui16_image::_exec01_fill_noise_pixel_scanline(uint16_t *ui16p_y1, uint16_t *ui16p_y2, uint16_t *ui16p_y3)
{
	uint16_t *ui16p_x1a, *ui16p_x1b, *ui16p_x1c,
		*ui16p_x2a, *ui16p_x2b, *ui16p_x2c,
		*ui16p_x3a, *ui16p_x3b, *ui16p_x3c;
	int32_t xx,
		i32_fill_count;

	/* 始めのピクセル位置(画像から外れているものはNULLをいれる) */
	ui16p_x1a = NULL;
	ui16p_x1b = NULL;
	ui16p_x1c = NULL;
	ui16p_x2a = NULL;
	ui16p_x2b = NULL;
	ui16p_x2c = NULL;
	ui16p_x3a = NULL;
	ui16p_x3b = NULL;
	ui16p_x3c = NULL;
	if (NULL != ui16p_y1) {
		ui16p_x1a = ui16p_y1 + 1;
		ui16p_x1b = ui16p_y1;
	}
	if (NULL != ui16p_y2) {
		ui16p_x2a = ui16p_y2 + 1;
		ui16p_x2b = ui16p_y2;
	}
	if (NULL != ui16p_y3) {
		ui16p_x3a = ui16p_y3 + 1;
		ui16p_x3b = ui16p_y3;
	}

	/* 画像を横にループ */
	i32_fill_count = 0;
	for (xx = 0; xx < this->_i32_xs; ++xx) {
		i32_fill_count += this->_exec01_fill_noise_pixel_pixel(
			ui16p_x1a, ui16p_x1b, ui16p_x1c,
			ui16p_x2a, ui16p_x2b, ui16p_x2c,
			ui16p_x3a, ui16p_x3b, ui16p_x3c);

		/* 次のピクセルへ進める(画像から外れたらNULLをいれる) */
		ui16p_x1c = ui16p_x1b;
		ui16p_x1b = ui16p_x1a;
		ui16p_x2c = ui16p_x2b;
		ui16p_x2b = ui16p_x2a;
		ui16p_x3c = ui16p_x3b;
		ui16p_x3b = ui16p_x3a;
		if ((xx + 1) < (this->_i32_xs)) {
			if (NULL != (ui16p_x1a))
				++ui16p_x1a;
			if (NULL != (ui16p_x2a))
				++ui16p_x2a;
			if (NULL != (ui16p_x3a))
				++ui16p_x3a;
		} else {
			ui16p_x1a = NULL;
			ui16p_x2a = NULL;
			ui16p_x3a = NULL;
		}
	}

	/* 終了 */
	return i32_fill_count;
}

int32_t thinnest_ui16_image::_exec01_fill_noise_pixel_pixel(uint16_t *ui16p_x1a, uint16_t *ui16p_x1b, uint16_t *ui16p_x1c, uint16_t *ui16p_x2a, uint16_t *ui16p_x2b, uint16_t *ui16p_x2c, uint16_t *ui16p_x3a, uint16_t *ui16p_x3b, uint16_t *ui16p_x3c)
{
	int32_t i32_fill_count;
	int32_t i32_exist_count, i32_eqgt_count;

	/* ui16p_x2bがNULLであってはならない */
	assert(NULL != ui16p_x2b);

	/* パラメータの初期値 */
	i32_fill_count = 0;
	i32_eqgt_count = 0L;
	i32_exist_count = 0L;

	/* 時計回り(右回り)に見ていく */
	/* 始め */
	if (NULL != (ui16p_x1a)) {
		if (this->_ui16_threshold <= (*ui16p_x1a)) {
			++i32_eqgt_count;
		}
		++i32_exist_count;
	}
	/* 次 ... */
	if (NULL != (ui16p_x1b)) {
		if (this->_ui16_threshold <= (*ui16p_x1b)) {
			++i32_eqgt_count;
		}
		++i32_exist_count;
	}
	if (NULL != (ui16p_x1c)) {
		if (this->_ui16_threshold <= (*ui16p_x1c)) {
			++i32_eqgt_count;
		}
		++i32_exist_count;
	}
	if (NULL != (ui16p_x2c)) {
		if (this->_ui16_threshold <= (*ui16p_x2c)) {
			++i32_eqgt_count;
		}
		++i32_exist_count;
	}
	if (NULL != (ui16p_x3c)) {
		if (this->_ui16_threshold <= (*ui16p_x3c)) {
			++i32_eqgt_count;
		}
		++i32_exist_count;
	}
	if (NULL != (ui16p_x3b)) {
		if (this->_ui16_threshold <= (*ui16p_x3b)) {
			++i32_eqgt_count;
		}
		++i32_exist_count;
	}
	if (NULL != (ui16p_x3a)) {
		if (this->_ui16_threshold <= (*ui16p_x3a)) {
			++i32_eqgt_count;
		}
		++i32_exist_count;
	}
	if (NULL != (ui16p_x2a)) {
		if (this->_ui16_threshold <= (*ui16p_x2a)) {
			++i32_eqgt_count;
		}
		++i32_exist_count;
	}

	/* 回りが全部閾値より大きかったら白線上の独立点なので潰す */
	if (i32_exist_count == i32_eqgt_count) {
		*ui16p_x2b = UINT16_MAX;
		++i32_fill_count;
	}

	/* 終了 */
	return i32_fill_count;
}
#include "igs_line_blur.h" // "pri.h" "thinnest_ui16_image.h"

void thinnest_ui16_image::exec02_scale_add_edge_pixel(void)
{
	int32_t xx, yy, i32_tmp1, i32_tmp2;
	uint16_t *ui16p_src, *ui16p_src1, *ui16p_src2,
		*ui16p_tgt, *ui16p_tgt1, *ui16p_tgt2;

	/* 画像サイズが足りないか、分割数が指定されていない、ときは実行キャンセル */
	if ((this->_i32_xs < 2) || (this->_i32_ys < 2) || (this->_i32_xd < 2) || (this->_i32_yd < 2)) {
		return;
	}

	/* 処理ごとのメッセージ */
	if (ON == this->_i_mv_sw) {
		pri_funct_msg_ttvr("thinnest_ui16_image::exec02_scale_add_edge_pixel()");
	}
	/* カウントダウン表示始め */
	if (ON == this->_i_cv_sw) {
		pri_funct_cv_start(this->_i32_ys);
	}

	ui16p_src = this->get_ui16p_src_channel();
	ui16p_tgt1 = this->get_ui16p_tgt_channel() + (this->_i32_xs + 2) + 1;

	/* 縁以外の部分のコピー */
	for (yy = 0; yy < this->_i32_ys; ++yy) {
		/* カウントダウン表示中 */
		if (ON == this->_i_cv_sw) {
			pri_funct_cv_run(yy);
		}

		ui16p_tgt2 = ui16p_tgt1;
		for (xx = 0; xx < this->_i32_xs; ++xx) {
			*ui16p_tgt2 = *ui16p_src;
			++ui16p_src;
			++ui16p_tgt2;
		}
		ui16p_tgt1 += (this->_i32_xs + 2);
	}
	/* カウントダウン表示終了 */
	if (ON == this->_i_cv_sw) {
		pri_funct_cv_end();
	}

	/* 縁の上部分 */
	ui16p_src1 = this->get_ui16p_src_channel();
	ui16p_src2 = this->get_ui16p_src_channel() + this->_i32_xs;
	ui16p_tgt = this->get_ui16p_tgt_channel() + 1;
	for (xx = 0; xx < this->_i32_xs; ++xx) {
		i32_tmp1 = (int32_t)(*ui16p_src1); /* 隣接値 */
		i32_tmp2 = (int32_t)(*ui16p_src2); /* さらに一つ隣の値 */
		i32_tmp1 *= 2;					   /* 隣接値２倍 */
		i32_tmp1 -= i32_tmp2;			   /* カレントの値を得る */
		if (i32_tmp1 < 0) {
			*ui16p_tgt = 0;
		} else if ((int32_t)UINT16_MAX < i32_tmp1) {
			*ui16p_tgt = UINT16_MAX;
		} else {
			*ui16p_tgt = (uint16_t)i32_tmp1;
		}

		++ui16p_src1;
		++ui16p_src2;
		++ui16p_tgt;
	}
	/* 縁の下部分 */
	ui16p_src1 = this->get_ui16p_src_channel() +
				 ((this->_i32_ys - 1) * this->_i32_xs);
	ui16p_src2 = this->get_ui16p_src_channel() +
				 ((this->_i32_ys - 2) * this->_i32_xs);
	ui16p_tgt = this->get_ui16p_tgt_channel() +
				(this->_i32_ys + 2 - 1) * (this->_i32_xs + 2) + 1;
	for (xx = 0; xx < this->_i32_xs; ++xx) {
		i32_tmp1 = (int32_t)(*ui16p_src1); /* 隣接値 */
		i32_tmp2 = (int32_t)(*ui16p_src2); /* さらに一つ隣の値 */
		i32_tmp1 *= 2;					   /* 隣接値２倍 */
		i32_tmp1 -= i32_tmp2;			   /* カレントの値を得る */
		if (i32_tmp1 < 0) {
			*ui16p_tgt = 0;
		} else if ((int32_t)UINT16_MAX < i32_tmp1) {
			*ui16p_tgt = UINT16_MAX;
		} else {
			*ui16p_tgt = (uint16_t)i32_tmp1;
		}

		++ui16p_src1;
		++ui16p_src2;
		++ui16p_tgt;
	}
	/* 縁の左部分 */
	ui16p_src1 = this->get_ui16p_src_channel();
	ui16p_src2 = this->get_ui16p_src_channel() + 1;
	ui16p_tgt = this->get_ui16p_tgt_channel() + (this->_i32_xs + 2);
	for (yy = 0; yy < this->_i32_ys; ++yy) {
		i32_tmp1 = (int32_t)(*ui16p_src1); /* 隣接値 */
		i32_tmp2 = (int32_t)(*ui16p_src2); /* さらに一つ隣の値 */
		i32_tmp1 *= 2;					   /* 隣接値２倍 */
		i32_tmp1 -= i32_tmp2;			   /* カレントの値を得る */
		if (i32_tmp1 < 0) {
			*ui16p_tgt = 0;
		} else if ((int32_t)UINT16_MAX < i32_tmp1) {
			*ui16p_tgt = UINT16_MAX;
		} else {
			*ui16p_tgt = (uint16_t)i32_tmp1;
		}

		ui16p_src1 += this->_i32_xs;
		ui16p_src2 += this->_i32_xs;
		ui16p_tgt += (this->_i32_xs + 2);
	}
	/* 縁の右部分 */
	ui16p_src1 = this->get_ui16p_src_channel() + this->_i32_xs - 1;
	ui16p_src2 = this->get_ui16p_src_channel() + this->_i32_xs - 1 - 1;
	ui16p_tgt = this->get_ui16p_tgt_channel() + (this->_i32_xs + 2) * 2 - 1;
	for (yy = 0; yy < this->_i32_ys; ++yy) {
		i32_tmp1 = (int32_t)(*ui16p_src1); /* 隣接値 */
		i32_tmp2 = (int32_t)(*ui16p_src2); /* さらに一つ隣の値 */
		i32_tmp1 *= 2;					   /* 隣接値２倍 */
		i32_tmp1 -= i32_tmp2;			   /* カレントの値を得る */
		if (i32_tmp1 < 0) {
			*ui16p_tgt = 0;
		} else if ((int32_t)UINT16_MAX < i32_tmp1) {
			*ui16p_tgt = UINT16_MAX;
		} else {
			*ui16p_tgt = (uint16_t)i32_tmp1;
		}

		ui16p_src1 += this->_i32_xs;
		ui16p_src2 += this->_i32_xs;
		ui16p_tgt += (this->_i32_xs + 2);
	}

	/* 画像サイズを縁１ピクセル増やす */
	this->_i32_xs += 2;
	this->_i32_ys += 2;

	/* 処理終了したらsrc,tgt画像交換 */
	this->_swap_channel();
}

#include <math.h> /* floor(), ceil() */

#include "igs_line_blur.h" // "pri.h" "thinnest_ui16_image.h"

void thinnest_ui16_image::exec03_scale_liner(void)
{
	int32_t xx, yy,
		i32_tgt_xs,
		i32_tgt_ys;
	uint16_t *ui16p_src,
		*ui16p_tgt;
	double d_src_xp, d_src_yp,
		d_tgt_xp, d_tgt_yp;
	int32_t i32_src_xpos1, i32_src_ypos1,
		i32_src_xpos2, i32_src_ypos2;
	double d_src_xratio1, d_src_yratio1,
		d_src_xratio2, d_src_yratio2;
	double d_tmp;

	/* 画像サイズが足りないか、分割数が指定されていない、ときは実行キャンセル */
	if ((this->_i32_xs < 2) || (this->_i32_ys < 2) || (this->_i32_xd < 2) || (this->_i32_yd < 2)) {
		return;
	}

	/* 処理ごとのメッセージ */
	if (ON == this->_i_mv_sw) {
		pri_funct_msg_ttvr("thinnest_ui16_image::exec03_scale_liner()");
	}
	if (ON == this->_i_mv_sw) {
		pri_funct_msg_ttvr("thi : Scale %d x %d",
						   this->_i32_xd,
						   this->_i32_yd);
	}
	/* カウントダウン表示始め */
	if (ON == this->_i_cv_sw) {
		pri_funct_cv_start((this->_i32_ys - 2) * this->_i32_yd);
	}

	ui16p_src = this->get_ui16p_src_channel();
	ui16p_tgt = this->get_ui16p_tgt_channel();

	i32_tgt_ys = (this->_i32_ys - 2) * this->_i32_yd;
	i32_tgt_xs = (this->_i32_xs - 2) * this->_i32_xd;

	/* ターゲットの大きさでループ */
	for (yy = 0; yy < i32_tgt_ys; ++yy) {
		for (xx = 0; xx < i32_tgt_xs; ++xx, ++ui16p_tgt) {
			/* カウントダウン表示中 */
			if (ON == this->_i_cv_sw) {
				pri_funct_cv_run(yy);
			}

			/*
	画像scaleのピクセルの対応図
	  (3倍拡大の場合)
	+-------+
	|       |
	|       |
	|       |
	|       |
	|       |
	|   +   |外周付加
	|       |ピクセル
	|       |
	|       |
	|       |
	|       |
	+-------+	+-------+
	|       |	|       |
	|       |  -->  |   +   |
	|       |	|       |
	|       |	+-------+
	|       |	|       |
	|   +   |  -->  |   +   |
	|       |	|       |
	|       |	+-------+
	|       |	|       |
	|       |  -->  |   +   |
	|       |	|       |
	+-------+	+-------+
	|       |	|       |
	|       |  -->  |   +   |
	|       |	|       |
	|       |	+-------+
	|       |	|       |
	|   +   |  -->  |   +   |
	|       |	|       |
	|       |	+-------+
	|       |	|       |
	|       |  -->  |   +   |
	|       |	|       |
	+-------+	+-------+
	|       |
	|       |
	|       |
	|       |
	|       |
	|   +   |外周付加
	|       |ピクセル
	|       |
	|       |
	|       |
	|       |
	+-------+
*/
			/* ターゲット上の正規化された位置 */
			d_tgt_xp = ((double)xx + 0.5) / (double)i32_tgt_xs;
			d_tgt_yp = ((double)yy + 0.5) / (double)i32_tgt_ys;

			/* ソース画像上の正規化された位置 */
			d_src_xp = ((double)(this->_i32_xs) - 2.0) * d_tgt_xp + 0.5;
			d_src_yp = ((double)(this->_i32_ys) - 2.0) * d_tgt_yp + 0.5;

			/* ソース画像上の正規化された上下、左右の位置 */
			i32_src_xpos1 = (int32_t)floor(d_src_xp);
			i32_src_ypos1 = (int32_t)floor(d_src_yp);
			i32_src_xpos2 = (int32_t)ceil(d_src_xp);
			i32_src_ypos2 = (int32_t)ceil(d_src_yp);
			/* ソース画像上の正規化された上下、左右の比率 */
			if (i32_src_xpos1 == i32_src_xpos2) {
				d_src_xratio1 = 0.0;
				d_src_xratio2 = 1.0;
			} else {
				d_src_xratio1 = d_src_xp - floor(d_src_xp);
				d_src_xratio2 = ceil(d_src_xp) - d_src_xp;
			}
			if (i32_src_ypos1 == i32_src_ypos2) {
				d_src_yratio1 = 0.0;
				d_src_yratio2 = 1.0;
			} else {
				d_src_yratio1 = d_src_yp - floor(d_src_yp);
				d_src_yratio2 = ceil(d_src_yp) - d_src_yp;
			}

			/* 比率でその場所の値を計算する */
			d_tmp = ((ui16p_src[this->_i32_xs * i32_src_ypos1 + i32_src_xpos1] *
						  d_src_xratio2 +
					  ui16p_src[this->_i32_xs * i32_src_ypos1 + i32_src_xpos2] *
						  d_src_xratio1) *
						 d_src_yratio2 +
					 (ui16p_src[this->_i32_xs * i32_src_ypos2 + i32_src_xpos1] *
						  d_src_xratio2 +
					  ui16p_src[this->_i32_xs * i32_src_ypos2 + i32_src_xpos2] *
						  d_src_xratio1) *
						 d_src_yratio1);
			if (UINT16_MAX <= d_tmp) {
				ui16p_tgt[0] = (uint16_t)UINT16_MAX;
			} else {
				ui16p_tgt[0] = (uint16_t)d_tmp;
			}
		}
	}
	/* カウントダウン表示終了 */
	if (ON == this->_i_cv_sw) {
		pri_funct_cv_end();
	}

	/* 画像サイズの設定 */
	this->_i32_xs = (this->_i32_xs - 2) * this->_i32_xd;
	this->_i32_ys = (this->_i32_ys - 2) * this->_i32_yd;

	/* 処理終了したらsrc,tgt画像交換 */
	this->_swap_channel();
}

#include "igs_line_blur.h" // "pri.h" "thinnest_ui16_image.h"

/* データB/W化 */
void thinnest_ui16_image::exec04_bw(void)
{
	int32_t xx, yy;
	uint16_t *ui16p_src, *ui16p_tgt;

	/* 処理実行表示 */
	if (ON == this->_i_mv_sw) {
		pri_funct_msg_ttvr("thinnest_ui16_image::exec04_bw()");
	}
	/* パラメータ表示 */
	if (ON == this->_i_pv_sw) {
		pri_funct_msg_ttvr("thi : threshold %u", this->_ui16_threshold);
	}
	/* カウントダウン表示始め */
	if (ON == this->_i_cv_sw) {
		pri_funct_cv_start(this->_i32_ys);
	}

	ui16p_src = this->get_ui16p_src_channel();
	ui16p_tgt = this->get_ui16p_tgt_channel();

	for (yy = 0L; yy < this->_i32_ys; ++yy) {
		/* カウントダウン表示中 */
		if (ON == this->_i_cv_sw) {
			pri_funct_cv_run(yy);
		}

		for (xx = 0L; xx < this->_i32_xs; ++xx, ++ui16p_src, ++ui16p_tgt) {
			if (this->_ui16_threshold <= *ui16p_src) {
				*ui16p_tgt = UINT16_MAX;
			} else {
				*ui16p_tgt = 0;
			}
		}
	}
	/* カウントダウン表示終了 */
	if (ON == this->_i_cv_sw) {
		pri_funct_cv_end();
	}

	/* 処理終了したらsrc,tgt画像交換 */
	this->_swap_channel();
}

#include "igs_line_blur.h" // "pri.h" "thinnest_ui16_image.h"

int thinnest_ui16_image::exec05_thin(void)
{
	int32_t ii, jj,
		i32_pixel_count_total,
		i32_pixel_count_one_round,
		i32_pixel_count_one_side_tmp,
		i32_pixel_count_one_side[4];

	/* 処理ごとのメッセージ */
	if (ON == this->_i_mv_sw) {
		pri_funct_msg_ttvr("thinnest_ui16_image::exec05_thin()");
	}

	/* 削るところがなくなるまでループする */
	i32_pixel_count_total = 0;
	i32_pixel_count_one_side[0] = 0;
	i32_pixel_count_one_side[1] = 0;
	i32_pixel_count_one_side[2] = 0;
	i32_pixel_count_one_side[3] = 0;
	for (ii = 0;; ++ii) {

		if (this->get_i32_exec_loop_count() <= ii) {
			pri_funct_err_bttvr("Error : loop counter over %ld.",
								this->get_i32_exec_loop_count());
			return NG;
		}

		/* 4方向に1pixelずつ細線化処理を行なう */
		i32_pixel_count_one_round = 0;
		for (jj = 0; jj < 4; ++jj) {

			/* 右方向への細線化 */
			i32_pixel_count_one_side_tmp = this->_one_side_thinner();

			/* 排除したpixel数のカウント */
			i32_pixel_count_one_round += i32_pixel_count_one_side_tmp;
			i32_pixel_count_total += i32_pixel_count_one_side_tmp;
			i32_pixel_count_one_side[jj] += i32_pixel_count_one_side_tmp;

			/* 時計方向(右)へ画像を90度回転する */
			this->_rot90_by_clockwork();

			/* カウントダウン表示(上下左右ひとつずつピリオド表示) */
			if (ON == this->_i_cv_sw) {
				(void)fprintf(stdout, ".");
				(void)fflush(stdout);
			}
		}
		/* カウントダウン表示(上下左右やって"....|"となる) */
		if (ON == this->_i_cv_sw) {
			(void)fprintf(stdout, "|");
			(void)fflush(stdout);
			/* 5 x 10カラムで改行 */
			if (9 == (ii % 10)) {
				(void)fprintf(stdout, "\n");
				(void)fflush(stdout);
			}
		}

		if (i32_pixel_count_one_round <= 0)
			break;
	}
	/* カウントダウン表示終了 */
	if (ON == this->_i_cv_sw) {
		(void)fprintf(stdout, "\nthin line ... end.\n");
		(void)fflush(stdout);
	}

	if (ON == this->_i_pv_sw) {
		pri_funct_msg_ttvr(
			"thi : total %d loop, and %ld pixel deleted",
			ii,
			i32_pixel_count_total);
		pri_funct_msg_ttvr(
			"thi : and each r<%ld>+ t<%ld>+ l<%ld>+ b<%ld> pixel deleted",
			i32_pixel_count_one_side[0],
			i32_pixel_count_one_side[1],
			i32_pixel_count_one_side[2],
			i32_pixel_count_one_side[3]);
	}

	/* 正常終了 */
	return OK;
}
#include <stdlib.h> /* free(), calloc() */

#include "igs_line_blur.h" // "pri.h" "thinnest_ui16_image.h"

/* メモリ開放 */
void thinnest_ui16_image::mem_free(void)
{
#if 0
	if (NULL !=	this->_ui16p_channel[0]) {
		if (ON == this->_i_mv_sw) {
			pri_funct_msg_ttvr( "thinnest_ui16_image::mem_free()" );
		}

		free(	this->_ui16p_channel[0]);/* ここで落ちる2014-5-16 */
			this->_ui16p_channel[0] = NULL;
			this->_ui16p_channel[1] = NULL;
	}
#endif
	if (NULL != this->memory_free_this_) {
		if (ON == this->_i_mv_sw) {
			pri_funct_msg_ttvr(
				"thinnest_ui16_image::mem_free() <%x>", this->memory_free_this_);
		}

		free(this->memory_free_this_); /* これだと落ちない2014-5-16 */
		this->memory_free_this_ = NULL;
		this->_ui16p_channel[0] = NULL;
		this->_ui16p_channel[1] = NULL;
	}
}

/* データ設定とメモリ確保 */
int thinnest_ui16_image::mem_alloc(void)
{
	/* 以前のメモリが残っていたら開放する */
	this->mem_free();

	/* 処理ごとのメッセージ */
	if (ON == this->_i_mv_sw) {
		pri_funct_msg_ttvr("thinnest_ui16_image::mem_alloc()");
	}
	if (ON == this->_i_pv_sw) {
		pri_funct_msg_ttvr(
			"alloc ui16_image memory (%d+%d) x (%d+%d) x %d x %d bytes",
			this->_i32_xs * this->_i32_xd, 2,
			this->_i32_ys * this->_i32_yd, 2, 2,
			sizeof(uint16_t));
	}

	this->memory_free_this_ = this->_ui16p_channel[0] = (uint16_t *)calloc(
		(this->_i32_xs * this->_i32_xd + 2) *
			(this->_i32_ys * this->_i32_yd + 2) * 2,
		sizeof(uint16_t));
	if (ON == this->_i_pv_sw) {
		pri_funct_msg_ttvr(
			"thinnest_ui16_image::mem_alloc() memory <%x>", this->memory_free_this_);
	}
	if (NULL == this->_ui16p_channel[0]) {
		pri_funct_err_bttvr("Error : calloc(-) returns NULL.");
		return NG;
	}

	this->_ui16p_channel[1] = this->_ui16p_channel[0] +
							  (this->_i32_xs * this->_i32_xd + 2) *
								  (this->_i32_ys * this->_i32_yd + 2);

	return OK;
}

#include <assert.h> /* assert() */

#include "igs_line_blur.h" // "pri.h" "thinnest_ui16_image.h"

int32_t thinnest_ui16_image::_one_side_thinner(void)
{
	uint16_t *ui16p_src_y1, *ui16p_src_y2, *ui16p_src_y3,
		*ui16p_tgt;
	int32_t yy,
		i32_delete_count;

	/* 処理ごとのメッセージ */
	/******if (ON == this->_i_mv_sw) {
		pri_funct_msg_ttvr( "thinnest_ui16_image::_one_side_thinner()" );
	}******/
	/* カウントダウン表示始め */
	//if (ON == this->_i_cv_sw) { pri_funct_cv_start( this->_i32_ys ); }

	/* 始めのスキャンライン位置(画像から外れているものはNULLをいれる) */
	ui16p_src_y1 = this->get_ui16p_src_channel() + this->_i32_xs;
	ui16p_src_y2 = this->get_ui16p_src_channel();
	ui16p_src_y3 = NULL;
	ui16p_tgt = this->get_ui16p_tgt_channel();

	/* 画像を縦にループ */
	i32_delete_count = 0;
	for (yy = 0; yy < this->_i32_ys; ++yy) {

		/* カウントダウン表示中 */
		//if (ON == this->_i_cv_sw) { pri_funct_cv_run(yy); }

		/* スキャンライン(とその前後のスキャンライン)毎の処理 */
		i32_delete_count += this->_one_side_thinner_scanline(
			ui16p_src_y1, ui16p_src_y2, ui16p_src_y3, ui16p_tgt);

		/* 次のスキャンラインへ進める(画像から外れたらNULLをいれる) */
		ui16p_src_y3 = ui16p_src_y2;
		ui16p_src_y2 = ui16p_src_y1;
		if ((yy + 1) < (this->_i32_ys)) {
			ui16p_src_y1 += this->_i32_xs;
		} else
			ui16p_src_y1 = NULL;
		ui16p_tgt += this->_i32_xs;
	}

	/* カウントダウン表示終了 */
	//if (ON == this->_i_cv_sw) { pri_funct_cv_end(); }

	/* 処理終了したらsrc,tgt画像交換 */
	this->_swap_channel();

	/* 終了 */
	return i32_delete_count;
}

int32_t thinnest_ui16_image::_one_side_thinner_scanline(uint16_t *ui16p_src_y1, uint16_t *ui16p_src_y2, uint16_t *ui16p_src_y3, uint16_t *ui16p_tgt)
{
	uint16_t *ui16p_src_x1a, *ui16p_src_x1b, *ui16p_src_x1c,
		*ui16p_src_x2a, *ui16p_src_x2b, *ui16p_src_x2c,
		*ui16p_src_x3a, *ui16p_src_x3b, *ui16p_src_x3c,
		ui16_src_x2c_before, ui16_src_x2b_before;
	int32_t xx,
		i32_delete_count;

	/* 始めのピクセル位置(画像から外れているものはNULLをいれる) */
	ui16p_src_x1a = NULL;
	ui16p_src_x1b = NULL;
	ui16p_src_x1c = NULL;
	ui16p_src_x2a = NULL;
	ui16p_src_x2b = NULL;
	ui16p_src_x2c = NULL;
	ui16p_src_x3a = NULL;
	ui16p_src_x3b = NULL;
	ui16p_src_x3c = NULL;
	if (NULL != ui16p_src_y1) {
		ui16p_src_x1a = ui16p_src_y1 + 1;
		ui16p_src_x1b = ui16p_src_y1;
	}
	if (NULL != ui16p_src_y2) {
		ui16p_src_x2a = ui16p_src_y2 + 1;
		ui16p_src_x2b = ui16p_src_y2;
	}
	if (NULL != ui16p_src_y3) {
		ui16p_src_x3a = ui16p_src_y3 + 1;
		ui16p_src_x3b = ui16p_src_y3;
	}
	ui16_src_x2c_before = 0;
	ui16_src_x2b_before = 0;

	/* 画像を横にループ */
	i32_delete_count = 0;
	for (xx = 0; xx < this->_i32_xs; ++xx) {
		/* ui16p_src_x2bがNULLであってはならない */
		assert(NULL != ui16p_src_x2b);

		/* 処理前に、次のループで参照のため、値をとっておく */
		ui16_src_x2b_before = (*ui16p_src_x2b);

		/*	ピクセル位置が一番左はじ、あるいは、
			前のピクセルがゼロ値(黒)で、
		カレントピクセルがゼロ以外(白)のとき
		*/
		if (((NULL == ui16p_src_x2c) ||
			 (0 == ui16_src_x2c_before)) &&
			(0 < (*ui16p_src_x2b))) {
			i32_delete_count += this->_one_side_thinner_pixel(
				ui16p_src_x1a, ui16p_src_x1b, ui16p_src_x1c,
				ui16p_src_x2a, ui16p_src_x2b, ui16p_src_x2c,
				ui16p_src_x3a, ui16p_src_x3b, ui16p_src_x3c,
				ui16p_tgt);
		}
		/* 細線化に関係ない部分は単にコピーする */
		else {
			*ui16p_tgt = *ui16p_src_x2b;
		}

		/* 次のピクセルへ進める(画像から外れたらNULLをいれる) */
		ui16_src_x2c_before = ui16_src_x2b_before;
		ui16p_src_x1c = ui16p_src_x1b;
		ui16p_src_x1b = ui16p_src_x1a;
		ui16p_src_x2c = ui16p_src_x2b;
		ui16p_src_x2b = ui16p_src_x2a;
		ui16p_src_x3c = ui16p_src_x3b;
		ui16p_src_x3b = ui16p_src_x3a;
		if ((xx + 1) < (this->_i32_xs)) {
			if (NULL != (ui16p_src_x1a))
				++ui16p_src_x1a;
			if (NULL != (ui16p_src_x2a))
				++ui16p_src_x2a;
			if (NULL != (ui16p_src_x3a))
				++ui16p_src_x3a;
		} else {
			ui16p_src_x1a = NULL;
			ui16p_src_x2a = NULL;
			ui16p_src_x3a = NULL;
		}
		++ui16p_tgt;
	}

	/* 終了 */
	return i32_delete_count;
}

int32_t thinnest_ui16_image::_one_side_thinner_pixel(uint16_t *ui16p_src_x1a, uint16_t *ui16p_src_x1b, uint16_t *ui16p_src_x1c, uint16_t *ui16p_src_x2a, uint16_t *ui16p_src_x2b, uint16_t *ui16p_src_x2c, uint16_t *ui16p_src_x3a, uint16_t *ui16p_src_x3b, uint16_t *ui16p_src_x3c, uint16_t *ui16p_tgt)
{
	int32_t i32_delete_count;
	long l_off_count,
		l_white_count;
	int i_sw, i_sw2;

	/* ui16p_src_x2bがNULLであってはならない */
	assert(NULL != ui16p_src_x2b);

	/* パラメータの初期値 */
	i32_delete_count = 0;
	l_off_count = 0L;
	l_white_count = 0L;
	i_sw = i_sw2 = OFF;

	/* 時計回り(右回り)に見ていく */
	/* 始め */
	if (NULL != (ui16p_src_x1a)) {
		i_sw = (0 < (*ui16p_src_x1a)) ? ON : OFF;
		if (ON == i_sw)
			++l_white_count;
		if ((i_sw != i_sw2) && (OFF == i_sw))
			++l_off_count;
		i_sw2 = i_sw;
	}
	/* 次 ... */
	if (NULL != (ui16p_src_x1b)) {
		i_sw = (0 < (*ui16p_src_x1b)) ? ON : OFF;
		if (ON == i_sw)
			++l_white_count;
		if ((i_sw != i_sw2) && (OFF == i_sw))
			++l_off_count;
		i_sw2 = i_sw;
	}
	if (NULL != (ui16p_src_x1c)) {
		i_sw = (0 < (*ui16p_src_x1c)) ? ON : OFF;
		if (ON == i_sw)
			++l_white_count;
		if ((i_sw != i_sw2) && (OFF == i_sw))
			++l_off_count;
		i_sw2 = i_sw;
	}
	if (NULL != (ui16p_src_x2c)) {
		i_sw = (0 < (*ui16p_src_x2c)) ? ON : OFF;
		if (ON == i_sw)
			++l_white_count;
		if ((i_sw != i_sw2) && (OFF == i_sw))
			++l_off_count;
		i_sw2 = i_sw;
	}
	if (NULL != (ui16p_src_x3c)) {
		i_sw = (0 < (*ui16p_src_x3c)) ? ON : OFF;
		if (ON == i_sw)
			++l_white_count;
		if ((i_sw != i_sw2) && (OFF == i_sw))
			++l_off_count;
		i_sw2 = i_sw;
	}
	if (NULL != (ui16p_src_x3b)) {
		i_sw = (0 < (*ui16p_src_x3b)) ? ON : OFF;
		if (ON == i_sw)
			++l_white_count;
		if ((i_sw != i_sw2) && (OFF == i_sw))
			++l_off_count;
		i_sw2 = i_sw;
	}
	if (NULL != (ui16p_src_x3a)) {
		i_sw = (0 < (*ui16p_src_x3a)) ? ON : OFF;
		if (ON == i_sw)
			++l_white_count;
		if ((i_sw != i_sw2) && (OFF == i_sw))
			++l_off_count;
		i_sw2 = i_sw;
	}
	if (NULL != (ui16p_src_x2a)) {
		i_sw = (0 < (*ui16p_src_x2a)) ? ON : OFF;
		if (ON == i_sw)
			++l_white_count;
		if ((i_sw != i_sw2) && (OFF == i_sw))
			++l_off_count;
		i_sw2 = i_sw;
	}
	/* 調査の開始点のみがゼロのときのため(l_off_count)開始点を再度調査する */
	if (NULL != (ui16p_src_x1a)) {
		i_sw = (0 < (*ui16p_src_x1a)) ? ON : OFF;
		if ((i_sw != i_sw2) && (OFF == i_sw))
			++l_off_count;
		i_sw2 = i_sw;
	}

	/* 回りが全部黒だったら独立点なので消す */
	if (l_white_count <= 0L) {
		*ui16p_tgt = 0;
		++i32_delete_count;
	}
	/* ３近傍以上のときのエッジを消す(１近傍、２近傍は端点として残す) */
	else if ((1L == l_off_count) && (3 <= l_white_count)) {
		*ui16p_tgt = 0;
		++i32_delete_count;
	}
	/* １ピクセル幅の線分上(両端含む)のとき */
	else {
		*ui16p_tgt = *ui16p_src_x2b;
	}

	/* 終了 */
	return i32_delete_count;
}

#include "igs_line_blur.h" // "pri.h" "thinnest_ui16_image.h"

/* 時計回りに90度回転 */
void thinnest_ui16_image::_rot90_by_clockwork(void)
{
	int32_t i32_tmp;
	int32_t xx, yy;
	uint16_t *ui16p_src, *ui16p_tgt_y, *ui16p_tgt_x;

	/* 処理ごとのメッセージ */
	/******if (ON == this->_i_mv_sw) {
		pri_funct_msg_ttvr( "thinnest_ui16_image::_rot90_by_clockwork()" );
	}******/
	/* カウントダウン表示始め */
	//if (ON == this->_i_cv_sw) { pri_funct_cv_start( this->_i32_ys ); }

	ui16p_src = this->get_ui16p_src_channel();
	ui16p_tgt_y = this->get_ui16p_tgt_channel() + (this->_i32_ys - 1);

	for (yy = 0L; yy < this->_i32_ys; ++yy) {
		/* カウントダウン表示中 */
		//if (ON == this->_i_cv_sw) { pri_funct_cv_run(yy); }

		ui16p_tgt_x = ui16p_tgt_y;
		for (xx = 0L; xx < this->_i32_xs; ++xx) {
			*ui16p_tgt_x = *ui16p_src;
			++ui16p_src;
			ui16p_tgt_x += this->_i32_ys;
		}
		--ui16p_tgt_y;
	}
	/* カウントダウン表示終了 */
	//if (ON == this->_i_cv_sw) { pri_funct_cv_end(); }

	/* 横と縦のサイズを交換 */
	i32_tmp = this->_i32_xs;
	this->_i32_xs = this->_i32_ys;
	this->_i32_ys = i32_tmp;

	i32_tmp = this->_i32_xd;
	this->_i32_xd = this->_i32_yd;
	this->_i32_yd = i32_tmp;

	/* 処理終了したらsrc,tgt画像交換 */
	this->_swap_channel();
}
#include <stdexcept>

#include "igs_line_blur.h" // "brush_curve_blur.h" "igs_line_blur.h"

template <class T>
void igs_line_blur_brush_curve_point_put_image_template_(
	double *dp_pixel, int xp, int yp, const int height // no_margin
	,
	const int width // no_margin
	,
	const int channels, T *image_top // no_margin
	)
{
	for (int zz = 0; zz < channels; ++zz) {
		image_top[yp * channels * width +
				  xp * channels +
				  zz] = static_cast<T>(dp_pixel[zz]);
	}
}

/*
int xp, int yp
 --> ピクセル位置
*/
void igs_line_blur_brush_curve_point_put_image_(
	brush_curve_blur &cl_brush_curve_blur, int xp, int yp, const int height // no_margin
	,
	const int width // no_margin
	,
	const int channels, const int bits, void *out // no_margin
	)
{
	if ((xp < 0) && (width <= xp) && (yp < 0) && (height <= yp)) {
		throw std::domain_error(
			"Error : igs::line_blur::_brush_curve_point_put_image(-)");
	}

	if (16 == bits) {
		igs_line_blur_brush_curve_point_put_image_template_(
			cl_brush_curve_blur.get_dp_pixel(), xp, yp, height, width, channels, static_cast<unsigned short *>(out));
	} else if (8 == bits) {
		igs_line_blur_brush_curve_point_put_image_template_(
			cl_brush_curve_blur.get_dp_pixel(), xp, yp, height, width, channels, static_cast<unsigned char *>(out));
	}
}
#include "igs_line_blur.h" // "brush_curve_blur.h"

template <class T>
void igs_line_blur_brush_curve_line_get_image_template_(
	const T *image_top, int height, int width, int channels, int i32_blur_count, double *xp_array, double *yp_array, double *linepixels_array, double xp, double yp)
{
	for (int ii = 0; ii < i32_blur_count; ++ii) {
		const int xx = (int)(xp_array[ii] + xp + 0.5);
		const int yy = (int)(yp_array[ii] + yp + 0.5);
		if ((0 <= xx) && (xx < width) && (0 <= yy) && (yy < height)) {
			for (int zz = 0; zz < channels; ++zz) {
				linepixels_array[ii * CHANNEL_COUNT + zz] =
					(double)(image_top[yy * channels * width +
									   xx * channels +
									   zz]) +
					0.999999;
			}
		} else {
			for (int zz = 0; zz < channels; ++zz) {
				linepixels_array[ii * CHANNEL_COUNT + zz] = -1.0;
			}
		}
	}
}

/*
this->cl_brush_curve_blur.get_dp_xp()
this->cl_brush_curve_blur.get_dp_yp()
double xp, double yp
 --> ピクセルの真中が原点
*/
void igs_line_blur_brush_curve_line_get_image_(
	brush_curve_blur &cl_brush_curve_blur, const void *in // no_margin
	,
	const int height // no_margin
	,
	const int width // no_margin
	,
	const int channels, const int bits, double xp, double yp)
{
	if (16 == bits) {
		igs_line_blur_brush_curve_line_get_image_template_(
			static_cast<const unsigned short *>(in), height, width, channels, cl_brush_curve_blur.get_i32_count(), cl_brush_curve_blur.get_dp_xp(), cl_brush_curve_blur.get_dp_yp(), cl_brush_curve_blur.get_dp_linepixels(), xp, yp);
	} else if (8 == bits) {
		igs_line_blur_brush_curve_line_get_image_template_(
			static_cast<const unsigned char *>(in), height, width, channels, cl_brush_curve_blur.get_i32_count(), cl_brush_curve_blur.get_dp_xp(), cl_brush_curve_blur.get_dp_yp(), cl_brush_curve_blur.get_dp_linepixels(), xp, yp);
	}
}
#include <iostream>
#include <string.h> /* memcpy() */

#include "igs_line_blur.h" // "pri.h" "brush_curve_blur.h" "pixel_select_curve_blur.h" "pixel_line_root.h" "igs_line_blur.h"

int igs_line_blur_brush_curve_blur_subpixel_(
	brush_curve_blur &cl_brush_curve_blur, pixel_select_curve_blur_root &cl_pixel_select_curve_blur_root, pixel_line_root &cl_pixel_line_root

	,
	const void *in // no_margin
	,
	const int height // no_margin
	,
	const int width // no_margin
	,
	const int channels, const int bits, int32_t xx, int32_t yy)
{
	int32_t i32_subpixel,
		xsub, ysub;

	i32_subpixel = cl_brush_curve_blur.get_i32_subpixel_divide();

	for (ysub = 0; ysub < i32_subpixel; ++ysub) {
		for (xsub = 0; xsub < i32_subpixel; ++xsub) {

			/* 近い線分の、近い部分を選択 */
			cl_pixel_select_curve_blur_root.exec(
				(double)xx + (double)xsub / i32_subpixel - 0.5,
				(double)yy + (double)ysub / i32_subpixel - 0.5,
				&(cl_pixel_line_root),
				cl_brush_curve_blur.get_i32_count(),
				cl_brush_curve_blur.get_d_effect_area_radius());

			/* 選択できなかったらサブピクセルループを抜けて次のピクセルへ */
			if (cl_pixel_select_curve_blur_root.get_i32_count() <= 0) {
				return NG;
			}

			/* 合成線を生成 */
			cl_pixel_select_curve_blur_root.get_line(
				cl_brush_curve_blur.get_i32_count(),
				cl_brush_curve_blur.get_dp_xp(),
				cl_brush_curve_blur.get_dp_yp());

			/* 線分の各ピクセル値を取得 */
			igs_line_blur_brush_curve_line_get_image_(
				cl_brush_curve_blur, in, height, width, channels, bits, (double)xx + (double)xsub / i32_subpixel - 0.5, (double)yy + (double)ysub / i32_subpixel - 0.5);
			/* サブピクセル値を計算 */
			cl_brush_curve_blur.set_subpixel_value(xsub, ysub);
		}
	}
	return OK;
}

int igs_line_blur_brush_curve_blur_all_(
	int mv_sw, int pv_sw, int cv_sw, brush_curve_blur &cl_brush_curve_blur, pixel_select_curve_blur_root &cl_pixel_select_curve_blur_root, pixel_line_root &cl_pixel_line_root

	,
	const void *in // no_margin
	,
	const int height // no_margin
	,
	const int width // no_margin
	,
	const int channels, const int bits, void *out // no_margin
	)
{
	/* 処理ごとのメッセージ */
	if (ON == mv_sw) {
		std::cout << "igs::line_blur::_brush_curve_blur_all()" << std::endl;
	}
	if (ON == pv_sw) {
		std::cout
			<< " curve blur points count is " << std::endl
			<< cl_brush_curve_blur.get_i32_count() << std::endl
			<< " power is " << std::endl
			<< cl_brush_curve_blur.get_d_power() << std::endl
			<< " subpixel divide is " << std::endl
			<< cl_brush_curve_blur.get_i32_subpixel_divide() << std::endl
			<< " clip area for speedup is " << std::endl
			<< cl_brush_curve_blur.get_d_effect_area_radius() << std::endl;
	}

	/* ブラシメモリの確保 */
	if (OK != cl_brush_curve_blur.mem_alloc()) {
		throw std::domain_error(
			"Error : cl_brush_curve_blur.mem_alloc() returns NG");
	}

	/* ブラシの線ぼかし変化比率の設定 */
	cl_brush_curve_blur.init_ratio_array();

	/* 画像をinからoutへコピーしておく */
	(void)memcpy(out, in, height * width * channels * ((16 == bits) ? 2 : 1));
	/* カウントダウン表示始め */
	if (ON == cv_sw) {
		pri_funct_cv_start(height);
	}

	for (int yy = 0; yy < height; ++yy) {
		/* カウントダウン表示中 */
		if (ON == cv_sw) {
			pri_funct_cv_run(yy);
		}

		for (int xx = 0; xx < width; ++xx) {
			if (OK == igs_line_blur_brush_curve_blur_subpixel_(
						  cl_brush_curve_blur, cl_pixel_select_curve_blur_root, cl_pixel_line_root, in, height, width, channels, bits, xx, yy)) {
				/* ピクセル値を計算 */
				cl_brush_curve_blur.set_pixel_value();

				/* 結果をピクセルへ置く
		(取ったものと別の画像におくこと) */
				igs_line_blur_brush_curve_point_put_image_(
					cl_brush_curve_blur, xx, yy, height, width, channels, bits, out);
			}
		}
	}

	/* カウントダウン表示終了 */
	if (ON == cv_sw) {
		pri_funct_cv_end();
	}

	/* ブラシメモリの開放 */
	cl_brush_curve_blur.mem_free();

	return OK;
}

#include <math.h> /* floor() */

#include "igs_line_blur.h" // "brush_smudge_circle.h"

template <class T>
void igs_line_blur_brush_smudge_get_image_template_(
	T *in, int height, int width, int channels, double x1, double y1, double x2, double y2, double d_subsize, double *dp_image)
{
	/* 保存(復元)位置 */
	int x1p = (int)floor(x1 + d_subsize / 2.0);
	int x2p = (int)floor(x2 - d_subsize / 2.0);
	int y1p = (int)floor(y1 + d_subsize / 2.0);
	int y2p = (int)floor(y2 - d_subsize / 2.0);

	for (int yy = y1p; yy <= y2p; ++yy) {
		for (int xx = x1p; xx <= x2p; ++xx, dp_image += 5) {
			if ((0 <= xx) && (xx < width) && (0 <= yy) && (yy < height)) {
				for (int zz = 0; zz < 4; ++zz) {
					if (zz < channels) {
						dp_image[zz] = (double)(in[yy * channels * width +
												   xx * channels +
												   zz]) +
									   0.999999;
					} else {
						/* 本来来てはいけない条件なので意味のない値にしとく */
						dp_image[zz] = 0.0;
					}
				}
				dp_image[4] = 1.0;
			} else {
				dp_image[4] = 0.0;
			}
		}
	}
}

void igs_line_blur_brush_smudge_get_image_(
	brush_smudge_circle &cl_brush_smudge_circle, const void *in // no_margin
	,
	const int height // no_margin
	,
	const int width // no_margin
	,
	const int channels, const int bits, double d_xp, double d_yp)
{
	/* 画像上に置いたブラシの範囲 */
	double x1, y1, x2, y2;
	cl_brush_smudge_circle.get_dp_area(
		d_xp, d_yp, &x1, &y1, &x2, &y2);

	if (16 == bits) {
		igs_line_blur_brush_smudge_get_image_template_(
			static_cast<const unsigned short *>(in), height, width, channels, x1, y1, x2, y2, 1.0 / cl_brush_smudge_circle.get_i32_subpixel_divide(), cl_brush_smudge_circle.get_dp_pixel_image());
	} else if (8 == bits) {
		igs_line_blur_brush_smudge_get_image_template_(
			static_cast<const unsigned char *>(in), height, width, channels, x1, y1, x2, y2, 1.0 / cl_brush_smudge_circle.get_i32_subpixel_divide(), cl_brush_smudge_circle.get_dp_pixel_image());
	}
}
#include <math.h> /* floor() */

#include "igs_line_blur.h" // "brush_smudge_circle.h"

template <class T>
void igs_line_blur_brush_smudge_put_image_template_(
	double x1, double y1, double x2, double y2, double d_subsize, double *dp_image, int height, int width, int channels, T *out)
{
	/* 保存(復元)位置 */
	int x1p = (int)floor(x1 + d_subsize / 2.0);
	int x2p = (int)floor(x2 - d_subsize / 2.0);
	int y1p = (int)floor(y1 + d_subsize / 2.0);
	int y2p = (int)floor(y2 - d_subsize / 2.0);

	for (int yy = y1p; yy <= y2p; ++yy) {
		for (int xx = x1p; xx <= x2p; ++xx, dp_image += 5) {
			if ((0.0 < dp_image[4]) && (0 <= xx) && (xx < width) && (0 <= yy) && (yy < height)) {
				for (int zz = 0; zz < channels; ++zz) {
					out[yy * channels * width +
						xx * channels +
						zz] = static_cast<T>((double)out[yy * channels * width +
														 xx * channels +
														 zz] *
												 (1.0 - dp_image[4]) +
											 dp_image[zz]);
				}
			}
		}
	}
}

void igs_line_blur_brush_smudge_put_image_(
	brush_smudge_circle &cl_brush_smudge_circle, double d_xp, double d_yp, const int height // no_margin
	,
	const int width // no_margin
	,
	const int channels, const int bits, void *out // no_margin
	)
{
	/* 画像上に置いたブラシの範囲 */
	double x1, y1, x2, y2;
	cl_brush_smudge_circle.get_dp_area(
		d_xp, d_yp, &x1, &y1, &x2, &y2);

	if (16 == bits) {
		igs_line_blur_brush_smudge_put_image_template_(
			x1, y1, x2, y2, 1.0 / cl_brush_smudge_circle.get_i32_subpixel_divide(), cl_brush_smudge_circle.get_dp_pixel_image(), height, width, channels, static_cast<unsigned short *>(out));
	} else if (8 == bits) {
		igs_line_blur_brush_smudge_put_image_template_(
			x1, y1, x2, y2, 1.0 / cl_brush_smudge_circle.get_i32_subpixel_divide(), cl_brush_smudge_circle.get_dp_pixel_image(), height, width, channels, static_cast<unsigned char *>(out));
	}
}
#include <stdexcept>

#include "igs_line_blur.h" // "brush_smudge_circle.h"

void igs_line_blur_brush_smudge_line_(
	brush_smudge_circle &cl_brush_smudge_circle, const void *in // no_margin
	,
	const int height // no_margin
	,
	const int width // no_margin
	,
	const int channels, const int bits, void *out // no_margin
	,
	pixel_line_node *clp_line)
{
	pixel_point_node *clp_one_expand, *clp_another_expand, *clp_crnt;
	int32_t ii;
	double d_x1, d_y1, d_x2, d_y2;

	/* 両端点と中点を得る */
	clp_one_expand = clp_line->get_clp_link_one_expand();
	clp_another_expand = clp_line->get_clp_link_another_expand();

	/* ブラシマスクを円形にセット */
	cl_brush_smudge_circle.set_brush_circle();

	/* 指先(にじみ)の始点 */
	clp_crnt = clp_line->get_clp_link_middle();
	/* サブピクセルメモリーに画像を得る */
	igs_line_blur_brush_smudge_get_image_(
		cl_brush_smudge_circle, in, height, width, channels, bits, clp_crnt->get_d_xp_tgt(), clp_crnt->get_d_yp_tgt());
	/* 画像上に置いたブラシの範囲 */
	cl_brush_smudge_circle.get_dp_area(
		clp_crnt->get_d_xp_tgt(), clp_crnt->get_d_yp_tgt(), &d_x1, &d_y1, &d_x2, &d_y2);
	/* サブピクセル画像にする */
	cl_brush_smudge_circle.to_subpixel_from_pixel(
		d_x1, d_y1, d_x2, d_y2);
	/* ブラシに画像をセット */
	cl_brush_smudge_circle.copy_to_brush_from_image();

	/* 線の後ろへ向かってこする */
	for (clp_crnt = clp_crnt->get_clp_next_point(), ii = 0;
		 NULL != clp_crnt;
		 clp_crnt = clp_crnt->get_clp_next_point(), ++ii) {
		/* 偽の場合、たぶん無限ループ */
		if (clp_line->get_i32_point_count() <= ii) {
			throw std::domain_error(
				"Error : over clp_line->get_i32_point_count()");
		}
		/* 画像上に置いたブラシの範囲 */
		cl_brush_smudge_circle.get_dp_area(
			clp_crnt->get_d_xp_tgt(),
			clp_crnt->get_d_yp_tgt(),
			&d_x1, &d_y1, &d_x2, &d_y2);

		/* ブラシをかける */
		if ((0.0 <= d_x2) && (d_x1 < width) && (0.0 <= d_y2) && (d_y1 < height)) {
			/* 画像を得る */
			igs_line_blur_brush_smudge_get_image_(
				cl_brush_smudge_circle, in, height, width, channels, bits, clp_crnt->get_d_xp_tgt(), clp_crnt->get_d_yp_tgt());
			/* サブピクセル画像にする */
			cl_brush_smudge_circle.to_subpixel_from_pixel(
				d_x1, d_y1, d_x2, d_y2);
			/* かすれ処理 */
			cl_brush_smudge_circle.exec();
			/* サブピクセル画像を元ピクセルサイズにする */
			cl_brush_smudge_circle.to_pixel_from_subpixel(
				d_x1, d_y1, d_x2, d_y2);
			/* 画像を置く */
			igs_line_blur_brush_smudge_put_image_(
				cl_brush_smudge_circle, clp_crnt->get_d_xp_tgt(), clp_crnt->get_d_yp_tgt(), height, width, channels, bits, out);
		}
		/* 終点 */
		// if (clp_another_expand == clp_crnt) break;
	}

	/* 指先(にじみ)の始点 */
	clp_crnt = clp_line->get_clp_link_middle();
	/* サブピクセルメモリーに画像を得る */
	igs_line_blur_brush_smudge_get_image_(
		cl_brush_smudge_circle, in, height, width, channels, bits, clp_crnt->get_d_xp_tgt(), clp_crnt->get_d_yp_tgt());
	/* 画像上に置いたブラシの範囲 */
	cl_brush_smudge_circle.get_dp_area(
		clp_crnt->get_d_xp_tgt(),
		clp_crnt->get_d_yp_tgt(),
		&d_x1, &d_y1, &d_x2, &d_y2);
	/* サブピクセル画像にする */
	cl_brush_smudge_circle.to_subpixel_from_pixel(
		d_x1, d_y1, d_x2, d_y2);
	/* ブラシに画像をセット */
	cl_brush_smudge_circle.copy_to_brush_from_image();

	/* 線の前へ向かってこする */
	for (clp_crnt = clp_crnt->get_clp_previous_point(), ii = 0;
		 NULL != clp_crnt;
		 clp_crnt = clp_crnt->get_clp_previous_point(), ++ii) {
		/* 偽の場合、たぶん無限ループ */
		if (clp_line->get_i32_point_count() <= ii) {
			throw std::domain_error(
				"Error : over clp_line->get_i32_point_count() going front");
		}

		/* 画像上に置いたブラシの範囲 */
		cl_brush_smudge_circle.get_dp_area(
			clp_crnt->get_d_xp_tgt(),
			clp_crnt->get_d_yp_tgt(),
			&d_x1, &d_y1, &d_x2, &d_y2);

		/* ブラシをかける */
		if ((0.0 <= d_x2) && (d_x1 < width) && (0.0 <= d_y2) && (d_y1 < height)) {
			/* 画像を得る */
			igs_line_blur_brush_smudge_get_image_(
				cl_brush_smudge_circle, in, height, width, channels, bits, clp_crnt->get_d_xp_tgt(), clp_crnt->get_d_yp_tgt());
			/* サブピクセル画像にする */
			cl_brush_smudge_circle.to_subpixel_from_pixel(
				d_x1, d_y1, d_x2, d_y2);
			/* かすれ処理 */
			cl_brush_smudge_circle.exec();
			/* サブピクセル画像を元ピクセルサイズにする */
			cl_brush_smudge_circle.to_pixel_from_subpixel(
				d_x1, d_y1, d_x2, d_y2);
			/* 画像を置く */
			igs_line_blur_brush_smudge_put_image_(
				cl_brush_smudge_circle, clp_crnt->get_d_xp_tgt(), clp_crnt->get_d_yp_tgt(), height, width, channels, bits, out);
		}
		/* 終点 */
		// if (clp_one_expand == clp_crnt) break;
	}
}

#include <iostream>
#include <string.h> /* memcpy() */

#include "igs_line_blur.h" // "pri.h"

void igs_line_blur_brush_smudge_all_(
	int mv_sw, int pv_sw, int cv_sw, brush_smudge_circle &cl_brush_smudge_circle, pixel_line_root &cl_pixel_line_root, const void *in // no_margin
	,
	const int height // no_margin
	,
	const int width // no_margin
	,
	const int channels, const int bits, void *out // no_margin
	)
{
	/* 処理ごとのメッセージ */
	if (ON == mv_sw) {
		std::cout << "igs::line_expand::_brush_smudge_all()" << std::endl;
	}

	if (ON == pv_sw) {
		std::cout
			<< " smudge ratio "
			<< cl_brush_smudge_circle.get_d_ratio() << std::endl
			<< " smudge brush size by pixel "
			<< cl_brush_smudge_circle.get_i32_size_by_pixel() << std::endl
			<< " smudge subpixel divide "
			<< cl_brush_smudge_circle.get_i32_subpixel_divide() << std::endl;
	}

	/* ブラシメモリの確保 */
	if (OK != cl_brush_smudge_circle.mem_alloc()) {
		throw std::domain_error(
			"Error : cl_brush_smudge_circle.mem_alloc() returns NG");
	}

	/* 画像をinからoutへコピーしておく */
	(void)memcpy(out, in, height * width * channels * ((16 == bits) ? 2 : 1));

	/* カウントダウン表示始め */
	if (ON == cv_sw) {
		pri_funct_cv_start(cl_pixel_line_root.get_i32_count());
	}

	/* 汚れ線描画 */
	pixel_line_node *clp_line =
		(pixel_line_node *)cl_pixel_line_root.get_clp_first();
	for (int ii = 0; NULL != clp_line; clp_line = (pixel_line_node *)clp_line->get_clp_next(), ++ii) {
		if (cl_pixel_line_root.get_i32_count() <= ii) {
			throw std::domain_error(
				"Error : over cl_pixel_line_root.get_i32_count()");
		}

		/* カウントダウン表示中 */
		if (ON == cv_sw) {
			pri_funct_cv_run(ii);
		}

		igs_line_blur_brush_smudge_line_(
			cl_brush_smudge_circle, in, height, width, channels, bits, out, clp_line);
	}
	/* カウントダウン表示終了 */
	if (ON == cv_sw) {
		pri_funct_cv_end();
	}
}

#include <iostream>
#include <stdexcept>

#include "igs_line_blur.h" // "pri.h"

void igs_line_blur_image_get_(
	const int mv_sw, const int cv_sw, const long reference_channel, thinnest_ui16_image &cl_thinnest_ui16_image

	,
	const void *in // no_margin
	,
	const int height // no_margin
	,
	const int width // no_margin
	,
	const int channels, const int bits)
{
	/* 処理ごとのメッセージ */
	if (ON == mv_sw) {
		std::cout
			<< "igs_line_blur_image_get_()"
			<< std::endl
			<< "com : reference channel " << reference_channel
			<< std::endl;
	}

	int i32_xs = cl_thinnest_ui16_image.get_i32_xs();
	int i32_ys = cl_thinnest_ui16_image.get_i32_ys();
	unsigned short *out_incr = cl_thinnest_ui16_image.get_ui16p_src_channel();

	/* カウントダウン表示始め */
	if (ON == cv_sw) {
		pri_funct_cv_start(i32_ys);
	}

	if (8 == bits) {
		const unsigned char *in_incr =
			static_cast<const unsigned char *>(in);
		in_incr += reference_channel;
		for (int yy = 0; yy < i32_ys; ++yy) {
			/* カウントダウン表示中 */
			if (ON == cv_sw) {
				pri_funct_cv_run(yy);
			}

			for (int xx = 0; xx < i32_xs; ++xx) {
				/* 8bits -> 16bits変換して格納 */
				*out_incr =
					(((unsigned short)(*in_incr)) << 8) +
					(unsigned short)(*in_incr);
				/* 参照位置移動 */
				in_incr += channels;
				++out_incr;
			}
		}
	} else if (16 == bits) {
		const unsigned short *in_incr =
			static_cast<const unsigned short *>(in);
		in_incr += reference_channel;

		for (int yy = 0; yy < i32_ys; ++yy) {
			/* カウントダウン表示中 */
			if (ON == cv_sw) {
				pri_funct_cv_run(yy);
			}

			for (int xx = 0; xx < i32_xs; ++xx) {
				*out_incr = *in_incr;
				/* 参照位置移動 */
				in_incr += channels;
				++out_incr;
			}
		}
	} else {
		throw std::domain_error("Error : bits is not 8 or 16");
	}

	/* カウントダウン表示終了 */
	if (ON == cv_sw) {
		pri_funct_cv_end();
	}
}
}

#include <iostream>
#include <stdexcept>
#include "igs_line_blur.h" // "thinnest_ui16_image.h" "pixel_point_root.h" "pixel_line_root.h" "brush_curve_blur.h" "brush_smudge_circle.h" "pixel_select_same_way.h" "pixel_select_curve_blur.h" "igs_line_blur.h"

void igs::line_blur::convert(
	/* 入出力画像 */
	const void *in // no_margin
	,
	void *out // no_margin

	,
	const int height // no_margin
	,
	const int width // no_margin
	,
	const int channels, const int bits

	/* Action */
	,
	const int b_blur_count /* min=1   def=51   incr=1   max=100  */
	,
	const double b_blur_power /* min=0.1 def=1    incr=0.1 max=10.0 */
	,
	const int b_subpixel /* min=1   def=1    incr=1   max=8	 */
	,
	const int b_blur_near_ref /* min=1   def=5    incr=1   max=100	 */
	,
	const int b_blur_near_len /* min=1   def=160  incr=1   max=1000 */

	,
	const int b_smudge_thick /* min=1   def=7    incr=1   max=100	 */
	,
	const double b_smudge_remain
	/* min=0.0 def=0.85 incr=0.001 max=1.0*/

	,
	const int v_smooth_retry /* min=0   def=100  incr=1   max=1000 */
	,
	const int v_near_ref /* min=0   def=4    incr=1   max=100	 */
	,
	const int v_near_len /* min=2   def=160  incr=1   max=1000 */

	,
	const int mv_sw /* OFF */
	,
	const int pv_sw /* OFF */
	,
	const int cv_sw /* OFF */
	,
	const long reference_channel /* 3	=Alpha:RGBA orBGRA */
	,
	const int debug_save_sw /* OFF */
	,
	const int brush_action /* 0 =Curve Blur ,1=Smudge Brush */
	)
{
	/* --- 動作クラスコンストラクション --- */
	thinnest_ui16_image cl_thinnest_ui16_image;
	pixel_point_root cl_pixel_point_root;
	pixel_line_root cl_pixel_line_root;
	pixel_select_same_way_root cl_pixel_select_same_way_root;
	pixel_select_curve_blur_root cl_pixel_select_curve_blur_root;
	brush_smudge_circle cl_brush_smudge_circle;
	brush_curve_blur cl_brush_curve_blur;

	/* --- 標準出力への動作表示スイッチ --- */
	cl_thinnest_ui16_image.set_i_mv_sw(mv_sw);
	cl_pixel_point_root.set_i_mv_sw(mv_sw);
	cl_pixel_line_root.set_i_mv_sw(mv_sw);
	cl_brush_smudge_circle.set_i_mv_sw(mv_sw);
	cl_brush_curve_blur.set_i_mv_sw(mv_sw);

	cl_thinnest_ui16_image.set_i_pv_sw(pv_sw);
	cl_pixel_point_root.set_i_pv_sw(pv_sw);
	cl_pixel_line_root.set_i_pv_sw(pv_sw);
	cl_brush_smudge_circle.set_i_pv_sw(pv_sw);
	cl_brush_curve_blur.set_i_pv_sw(pv_sw);

	cl_thinnest_ui16_image.set_i_cv_sw(cv_sw);
	cl_pixel_point_root.set_i_cv_sw(cv_sw);
	cl_pixel_line_root.set_i_cv_sw(cv_sw);
	cl_brush_smudge_circle.set_i_cv_sw(cv_sw);
	cl_brush_curve_blur.set_i_cv_sw(cv_sw);

	/* --- 処理ごとのメッセージ --- */
	if (ON == mv_sw) {
		std::cout << "igs::line_blur::convert()" << std::endl;
	}

	/* --- 引数値をセット --- */

	/* 細線化用メモリの大きさを設定 */
	cl_thinnest_ui16_image.set_i32_xs(width);
	cl_thinnest_ui16_image.set_i32_ys(height);

	/* 0: $type : Wheel (min=0, default=0, increment=1, max=1)              */
	/******	if (0 == b_action_mode) {
		this->set_e_brush_smudge_action();
	} else {
		this->set_e_brush_curve_blur_action();
 		cl_pixel_line_root.set_i_same_way_exec_sw(OFF);
	}******/

	/* 1: $b_blur_count   :Wheel(min=1, default=51, increment=1, max=100)   */
	cl_brush_curve_blur.set_i32_count(b_blur_count);
	cl_brush_curve_blur.set_d_effect_area_radius(
		(double)(b_blur_count / 2));

	/* 2: $b_blur_power   :Wheel(min=0.1,default=1,increment=0.1,max=hh0)   */
	cl_brush_curve_blur.set_d_power(b_blur_power);

	/* 3. $b_subpixel : Wheel (min=1, default=1, increment=1, max=8)	*/
	cl_brush_curve_blur.set_i32_subpixel_divide(b_subpixel);

	/* 4: $b_blur_near_ref:Wheel(min=1, default=5, increment=1, max=100)    */
	cl_pixel_select_curve_blur_root.set_i32_count_max(
		b_blur_near_ref);

	/* 5: $b_blur_near_len:Wheel(min=1, default=160, increment=1, max=1000) */
	cl_pixel_select_curve_blur_root.set_d_length_max(
		b_blur_near_len);

	/* 6: $b_smudge_thick :Wheel(min=1, default=7, increment=1, max=100)    */
	cl_brush_smudge_circle.set_i32_size_by_pixel(b_smudge_thick);

	/* 7: $b_smudge_remain:Wheel(min=0.0,default=0.85,increment=0.001,max=1.0)*/
	cl_brush_smudge_circle.set_d_ratio(b_smudge_remain);

	/* 8: $v_smooth_retry :Wheel(min=0, default=100, increment=1, max=1000) */
	cl_pixel_line_root.set_i32_smooth_retry(v_smooth_retry);

	/* 9: $v_near_ref     :Wheel(min=0, default=4, increment=1, max=100)    */
	cl_pixel_select_same_way_root.set_i32_count_max(v_near_ref);

	/*10: $v_near_len     :Wheel(min=2, default=160, increment=1, max=1000) */
	cl_pixel_select_same_way_root.set_d_length_max(v_near_len);

	/* 細線化用メモリ確保 */
	if (OK != cl_thinnest_ui16_image.mem_alloc()) {
		throw std::domain_error(
			"Error : cl_thinnest_ui16_image.mem_alloc() returns NG");
	}

	/* 画像情報を細線化用メモリに移す */
	igs_line_blur_image_get_(
		mv_sw, cv_sw, reference_channel, cl_thinnest_ui16_image, in, height, width, channels, bits);

	/****** 細線化処理 start ******/

	/* 元画像の大きさの時に線分の点ノイズを消す */
	cl_thinnest_ui16_image.exec01_fill_noise_pixel();

	/* 精度をあげるための画像の拡大のための外周１ピクセル付加 */
	cl_thinnest_ui16_image.exec02_scale_add_edge_pixel();

	/* 精度をあげるための画像の拡大 */
	cl_thinnest_ui16_image.exec03_scale_liner();

	/* 白黒２値化 */
	cl_thinnest_ui16_image.exec04_bw();

	/* 細線化処理 */
	if (OK != cl_thinnest_ui16_image.exec05_thin()) {
		throw std::domain_error(
			"Error : cl_thinnest_ui16_image.exec05_thin() returns NG");
	}

	/****** 細線化処理 end ******/

	/****** ベクトルリスト処理 start ******/

	/* 細線化した画像をリストにする */
	if (OK != cl_pixel_point_root.alloc_mem_and_list_node(cl_thinnest_ui16_image.get_i32_xs(), cl_thinnest_ui16_image.get_i32_ys(), cl_thinnest_ui16_image.get_ui16p_src_channel())) {
		throw std::domain_error(
			"Error : cl_pixel_point_root.alloc_mem_and_list_node() returns NG");
	}

	/* 単なるポイントリストを線分リストにする */
	if (OK != cl_pixel_line_root.exec01020304(&(cl_pixel_point_root))) {
		throw std::domain_error(
			"Error : cl_pixel_line_root.exec01020304() returns NG");
	}
	if (ON == debug_save_sw) {
		if (OK != cl_pixel_line_root.save_lines("tmp08_jaggy_lines.txt")) {
			throw std::domain_error(
				"Error : cl_pixel_line_root.save_lines(-) returns NG");
		}
		if (OK != cl_pixel_line_root.save_one_point("tmp09_one_point.txt")) {
			throw std::domain_error(
				"Error : cl_pixel_line_root.save_one_point(-) returns NG");
		}
		if (OK != cl_pixel_line_root.save_another_point("tmp10_another_point.txt")) {
			throw std::domain_error(
				"Error : cl_pixel_line_root.save_another_point(-) returns NG");
		}
		if (OK != cl_pixel_line_root.save_not_include(&(cl_pixel_point_root), "tmp11_not_include.txt")) {
			throw std::domain_error(
				"Error : cl_pixel_line_root.save_not_include(-) returns NG");
		}
	}

	/* 中点をセットする */
	cl_pixel_line_root.exec05_set_middle();
	if (ON == debug_save_sw) {
		if (OK != cl_pixel_line_root.save_middle_point("tmp12_middle_point.txt")) {
			throw std::domain_error(
				"Error : cl_pixel_line_root.save_middle_point(-) returns NG");
		}
	}

	/* ベクトル化した線分の座標値を整数からdoubleに */
	cl_pixel_line_root.exec06_int2double_body();

	/* ベクトル化した線分をスムースに */
	cl_pixel_line_root.exec07_smooth_body();

	/* 線分を伸長する */
	if (OK != cl_pixel_line_root.exec08_expand_lines(&(cl_pixel_point_root))) {
		throw std::domain_error(
			"Error : cl_pixel_line_root.exec08_expand_lines(-) returns NG");
	}

	/* 伸ばした線分の方向をそろえる */
	if (ON == debug_save_sw) {
		if (OK != cl_pixel_line_root.save_expand_vector("tmp13_expand_vector.txt")) {
			throw std::domain_error(
				"Error : cl_pixel_line_root.save_expand_vector(-) returns NG");
		}
	}
	cl_pixel_line_root.exec09_same_way_expand(
		&(cl_pixel_select_same_way_root));
	if (ON == debug_save_sw) {
		if (OK != cl_pixel_line_root.save_expand_vector("tmp14_same_way_vector.txt")) {
			throw std::domain_error(
				"Error : cl_pixel_line_root.save_expand_vector(-) returns NG");
		}
	}

	/* 伸ばした線分を再度スムースに */
	cl_pixel_line_root.exec10_smooth_expand();

	if (ON == debug_save_sw) {
		if (OK != cl_pixel_line_root.save_expand_lines("tmp15_expand_lines.txt")) {
			throw std::domain_error(
				"Error : cl_pixel_line_root.save_expand_lines(-) returns NG");
		}
		if (OK != cl_pixel_line_root.save_one_expand_point("tmp16_one_expand_point.txt")) {
			throw std::domain_error(
				"Error : cl_pixel_line_root.save_one_expand_point(-) returns NG");
		}
		if (OK != cl_pixel_line_root.save_another_expand_point("tmp17_another_expand_point.txt")) {
			throw std::domain_error(
				"Error : cl_pixel_line_root.save_another_expand_point(-) returns NG");
		}
	}

	/* 各ラインと全体のバウンダリーボックスをセット */
	cl_pixel_line_root.exec11_set_bbox();

	/****** ベクトルリスト処理 end ******/

	/* 画像の加工 */
	if (0 == brush_action) {
		/* 線分情報からだいたいその方向にぼかす */
		igs_line_blur_brush_curve_blur_all_(
			mv_sw, pv_sw, cv_sw, cl_brush_curve_blur, cl_pixel_select_curve_blur_root, cl_pixel_line_root, in, height, width, channels, bits, out);
	} else if (1 == brush_action) {
		/* 画像をコピーしてから、指先ツールのようにこする */
		igs_line_blur_brush_smudge_all_(
			mv_sw, pv_sw, cv_sw, cl_brush_smudge_circle, cl_pixel_line_root, in, height, width, channels, bits, out);
	}

	/* 指先ツール用メモリ開放 */
	cl_brush_smudge_circle.mem_free();

	/* 線ぼかしツール用メモリ開放 */
	cl_brush_curve_blur.mem_free();

	/* 線ぼかし方向参照リスト開放 */
	cl_pixel_select_curve_blur_root.mem_free();

	/* 同方向線曲げ方向参照リスト開放 */
	cl_pixel_select_same_way_root.mem_free();

	/* ラインリストメモリ開放 */
	cl_pixel_line_root.mem_free();

	/* ポイントリストメモリ開放 */
	cl_pixel_point_root.mem_free();

	/* 細線化用メモリ開放 */
	cl_thinnest_ui16_image.mem_free();
}