/**
* 入力ラスタをHSV形式に変換して、出力ラスタへ書込みます。
* 画素形式は、コンストラクタに指定した形式に合せてください。
*/
public void doFilter(INyARRgbRaster i_input, INyARRaster i_output)
{
Debug.Assert(i_input.getSize().isEqualSize(i_output.getSize()) == true);
if (!this._do_filter_impl.isSupport(i_input, i_output))
{
this._do_filter_impl = this.createFilter(i_input, i_output);
}
this._do_filter_impl.doFilter(i_input, i_output, i_input.getSize());
}
public void doThFilter(INyARRaster i_raster, int i_l, int i_t, int i_w, int i_h, int i_th, INyARRaster o_raster)
{
Debug.Assert(i_raster.isEqualBufferType(NyARBufferType.BYTE1D_B8G8R8X8_32));
byte[] input = (byte[])i_raster.getBuffer();
int[] output = (int[])o_raster.getBuffer();
NyARIntSize s = i_raster.getSize();
int th = i_th * 3;
int skip_dst = (s.w - i_w);
int skip_src = skip_dst * 4;
int pix_count = i_w;
int pix_mod_part = pix_count - (pix_count % 8);
//左上から1行づつ走査していく
int pt_dst = (i_t * s.w + i_l);
int pt_src = pt_dst * 4;
for (int y = i_h - 1; y >= 0; y -= 1)
{
int x;
for (x = pix_count - 1; x >= pix_mod_part; x--)
{
output[pt_dst++] = ((input[pt_src + 0] & 0xff) + (input[pt_src + 1] & 0xff) + (input[pt_src + 2] & 0xff)) <= th?0:1;
pt_src += 4;
}
for (; x >= 0; x -= 8)
{
output[pt_dst++] = ((input[pt_src + 0] & 0xff) + (input[pt_src + 1] & 0xff) + (input[pt_src + 2] & 0xff)) <= th?0:1;
pt_src += 4;
output[pt_dst++] = ((input[pt_src + 0] & 0xff) + (input[pt_src + 1] & 0xff) + (input[pt_src + 2] & 0xff)) <= th?0:1;
pt_src += 4;
output[pt_dst++] = ((input[pt_src + 0] & 0xff) + (input[pt_src + 1] & 0xff) + (input[pt_src + 2] & 0xff)) <= th?0:1;
pt_src += 4;
output[pt_dst++] = ((input[pt_src + 0] & 0xff) + (input[pt_src + 1] & 0xff) + (input[pt_src + 2] & 0xff)) <= th?0:1;
pt_src += 4;
output[pt_dst++] = ((input[pt_src + 0] & 0xff) + (input[pt_src + 1] & 0xff) + (input[pt_src + 2] & 0xff)) <= th?0:1;
pt_src += 4;
output[pt_dst++] = ((input[pt_src + 0] & 0xff) + (input[pt_src + 1] & 0xff) + (input[pt_src + 2] & 0xff)) <= th?0:1;
pt_src += 4;
output[pt_dst++] = ((input[pt_src + 0] & 0xff) + (input[pt_src + 1] & 0xff) + (input[pt_src + 2] & 0xff)) <= th?0:1;
pt_src += 4;
output[pt_dst++] = ((input[pt_src + 0] & 0xff) + (input[pt_src + 1] & 0xff) + (input[pt_src + 2] & 0xff)) <= th?0:1;
pt_src += 4;
}
//スキップ
pt_src += skip_src;
pt_dst += skip_dst;
}
return;
}
/**
* o_histogramにヒストグラムを出力します。
* @param i_input
* @param o_histogram
* @return
* @throws NyARException
*/
public int analyzeRaster(INyARRaster i_input, NyARHistogram o_histogram)
{
NyARIntSize size = i_input.getSize();
//最大画像サイズの制限
Debug.Assert(size.w * size.h < 0x40000000);
Debug.Assert(o_histogram.length == 256); //現在は固定
int[] h = o_histogram.data;
//ヒストグラム初期化
for (int i = o_histogram.length - 1; i >= 0; i--)
{
h[i] = 0;
}
o_histogram.total_of_data = size.w * size.h / this._vertical_skip;
return(this._histImpl.createHistogram(i_input, size, h, this._vertical_skip));
}
public void analyzeRaster(INyARRaster i_input, NyARIntRect i_area, NyARHistogram o_histogram)
{
NyARIntSize size = i_input.getSize();
//最大画像サイズの制限
Debug.Assert(size.w * size.h < 0x40000000);
Debug.Assert(o_histogram.length == 256);//現在は固定
int[] h = o_histogram.data;
//ヒストグラム初期化
for (int i = o_histogram.length - 1; i >= 0; i--)
{
h[i] = 0;
}
o_histogram.total_of_data = i_area.w * i_area.h / this._vertical_skip;
this._histImpl.createHistogram(i_input, i_area.x, i_area.y, i_area.w, i_area.h, o_histogram.data, this._vertical_skip);
return;
}
public void createHistogram(INyARRaster i_raster, int i_l, int i_t, int i_w, int i_h, int[] o_histogram, int i_skip)
{
Debug.Assert(i_raster.isEqualBufferType(NyARBufferType.BYTE1D_X8R8G8B8_32));
byte[] input = (byte[])i_raster.getBuffer();
NyARIntSize s = i_raster.getSize();
int skip = (i_skip * s.w - i_w) * 4;
int pix_count = i_w;
int pix_mod_part = pix_count - (pix_count % 8);
//左上から1行づつ走査していく
int pt = (i_t * s.w + i_l) * 4;
for (int y = i_h - 1; y >= 0; y -= i_skip)
{
int x;
for (x = pix_count - 1; x >= pix_mod_part; x--)
{
o_histogram[((input[pt + 1] & 0xff) + (input[pt + 2] & 0xff) + (input[pt + 3] & 0xff)) / 3]++;
pt += 4;
}
for (; x >= 0; x -= 8)
{
o_histogram[((input[pt + 1] & 0xff) + (input[pt + 2] & 0xff) + (input[pt + 3] & 0xff)) / 3]++;
pt += 4;
o_histogram[((input[pt + 1] & 0xff) + (input[pt + 2] & 0xff) + (input[pt + 3] & 0xff)) / 3]++;
pt += 4;
o_histogram[((input[pt + 1] & 0xff) + (input[pt + 2] & 0xff) + (input[pt + 3] & 0xff)) / 3]++;
pt += 4;
o_histogram[((input[pt + 1] & 0xff) + (input[pt + 2] & 0xff) + (input[pt + 3] & 0xff)) / 3]++;
pt += 4;
o_histogram[((input[pt + 1] & 0xff) + (input[pt + 2] & 0xff) + (input[pt + 3] & 0xff)) / 3]++;
pt += 4;
o_histogram[((input[pt + 1] & 0xff) + (input[pt + 2] & 0xff) + (input[pt + 3] & 0xff)) / 3]++;
pt += 4;
o_histogram[((input[pt + 1] & 0xff) + (input[pt + 2] & 0xff) + (input[pt + 3] & 0xff)) / 3]++;
pt += 4;
o_histogram[((input[pt + 1] & 0xff) + (input[pt + 2] & 0xff) + (input[pt + 3] & 0xff)) / 3]++;
pt += 4;
}
//スキップ
pt += skip;
}
return;
}
public void doFilter(INyARRaster i_input, int[] o_output, int l, int t, int w, int h)
{
Debug.Assert(i_input.isEqualBufferType(NyARBufferType.BYTE1D_B8G8R8X8_32));
NyARIntSize size = i_input.getSize();
byte[] in_buf = (byte[])i_input.getBuffer();
int bp = (l + t * size.w) * 4;
int b = t + h;
int row_padding = (size.w - w) * 4;
for (int y = t; y < b; y++)
{
for (int x = 0; x < w; x++)
{
o_output[y * size.w + x + l] = ((in_buf[bp] & 0xff) + (in_buf[bp + 1] & 0xff) + (in_buf[bp + 2] & 0xff)) >> 2;
bp += 4;
}
bp += row_padding;
}
}
public void doThFilter(INyARRaster i_raster, int i_l, int i_t, int i_w, int i_h, int i_th, INyARRaster o_raster)
{
Debug.Assert(i_raster.isEqualBufferType(NyARBufferType.INT1D_X8R8G8B8_32));
int[] input = (int[])i_raster.getBuffer();
int[] output = (int[])o_raster.getBuffer();
int th = i_th * 3;
NyARIntSize s = i_raster.getSize();
int skip_src = (s.w - i_w);
int skip_dst = skip_src;
int pix_count = i_w;
int pix_mod_part = pix_count - (pix_count % 8);
//左上から1行づつ走査していく
int pt_dst = (i_t * s.w + i_l);
int pt_src = pt_dst;
for (int y = i_h - 1; y >= 0; y -= 1)
{
int x, v;
for (x = pix_count - 1; x >= pix_mod_part; x--)
{
v = input[pt_src++]; output[pt_dst++] = (((v >> 16) & 0xff) + ((v >> 8) & 0xff) + (v & 0xff)) <= th?0:1;
}
for (; x >= 0; x -= 8)
{
v = input[pt_src++]; output[pt_dst++] = (((v >> 16) & 0xff) + ((v >> 8) & 0xff) + (v & 0xff)) <= th?0:1;
v = input[pt_src++]; output[pt_dst++] = (((v >> 16) & 0xff) + ((v >> 8) & 0xff) + (v & 0xff)) <= th?0:1;
v = input[pt_src++]; output[pt_dst++] = (((v >> 16) & 0xff) + ((v >> 8) & 0xff) + (v & 0xff)) <= th?0:1;
v = input[pt_src++]; output[pt_dst++] = (((v >> 16) & 0xff) + ((v >> 8) & 0xff) + (v & 0xff)) <= th?0:1;
v = input[pt_src++]; output[pt_dst++] = (((v >> 16) & 0xff) + ((v >> 8) & 0xff) + (v & 0xff)) <= th?0:1;
v = input[pt_src++]; output[pt_dst++] = (((v >> 16) & 0xff) + ((v >> 8) & 0xff) + (v & 0xff)) <= th?0:1;
v = input[pt_src++]; output[pt_dst++] = (((v >> 16) & 0xff) + ((v >> 8) & 0xff) + (v & 0xff)) <= th?0:1;
v = input[pt_src++]; output[pt_dst++] = (((v >> 16) & 0xff) + ((v >> 8) & 0xff) + (v & 0xff)) <= th?0:1;
}
//スキップ
pt_src += skip_src;
pt_dst += skip_dst;
}
return;
}
public void doFilter(INyARRaster i_input, int[] o_output, int l, int t, int w, int h)
{
Debug.Assert(i_input.isEqualBufferType(NyARBufferType.BYTE1D_B8G8R8_24) || i_input.isEqualBufferType(NyARBufferType.BYTE1D_R8G8B8_24));
NyARIntSize size = i_input.getSize();
byte[] in_buf = (byte[])i_input.getBuffer();
int bp = (l + t * size.w) * 3;
int b = t + h;
int row_padding_dst = (size.w - w);
int row_padding_src = row_padding_dst * 3;
int pix_count = w;
int pix_mod_part = pix_count - (pix_count % 8);
int src_ptr = t * size.w + l;
for (int y = t; y < b; y++)
{
int x = 0;
for (x = pix_count - 1; x >= pix_mod_part; x--)
{
o_output[src_ptr++] = ((in_buf[bp++] & 0xff) + (in_buf[bp++] & 0xff) + (in_buf[bp++] & 0xff)) >> 2;
}
for (; x >= 0; x -= 8)
{
o_output[src_ptr++] = ((in_buf[bp++] & 0xff) + (in_buf[bp++] & 0xff) + (in_buf[bp++] & 0xff)) >> 2;
o_output[src_ptr++] = ((in_buf[bp++] & 0xff) + (in_buf[bp++] & 0xff) + (in_buf[bp++] & 0xff)) >> 2;
o_output[src_ptr++] = ((in_buf[bp++] & 0xff) + (in_buf[bp++] & 0xff) + (in_buf[bp++] & 0xff)) >> 2;
o_output[src_ptr++] = ((in_buf[bp++] & 0xff) + (in_buf[bp++] & 0xff) + (in_buf[bp++] & 0xff)) >> 2;
o_output[src_ptr++] = ((in_buf[bp++] & 0xff) + (in_buf[bp++] & 0xff) + (in_buf[bp++] & 0xff)) >> 2;
o_output[src_ptr++] = ((in_buf[bp++] & 0xff) + (in_buf[bp++] & 0xff) + (in_buf[bp++] & 0xff)) >> 2;
o_output[src_ptr++] = ((in_buf[bp++] & 0xff) + (in_buf[bp++] & 0xff) + (in_buf[bp++] & 0xff)) >> 2;
o_output[src_ptr++] = ((in_buf[bp++] & 0xff) + (in_buf[bp++] & 0xff) + (in_buf[bp++] & 0xff)) >> 2;
}
bp += row_padding_dst;
src_ptr += row_padding_src;
}
return;
}
public void doThFilter(INyARRaster i_raster, int i_l, int i_t, int i_w, int i_h, int i_th, INyARRaster o_raster)
{
Debug.Assert(i_raster.isEqualBufferType(NyARBufferType.WORD1D_R5G6B5_16LE));
short[] input = (short[])i_raster.getBuffer();
int[] output = (int[])o_raster.getBuffer();
int th = i_th * 3;
NyARIntSize s = i_raster.getSize();
int skip_dst = (s.w - i_w);
int skip_src = skip_dst;
int pix_count = i_w;
int pix_mod_part = pix_count - (pix_count % 8);
//左上から1行づつ走査していく
int pt_dst = (i_t * s.w + i_l);
int pt_src = pt_dst;
for (int y = i_h - 1; y >= 0; y -= 1)
{
int x, v;
for (x = pix_count - 1; x >= pix_mod_part; x--)
{
v = (int)input[pt_src++]; output[pt_dst++] = (((v & 0xf800) >> 8) + ((v & 0x07e0) >> 3) + ((v & 0x001f) << 3)) <= th?0:1;
}
for (; x >= 0; x -= 8)
{
v = (int)input[pt_src++]; output[pt_dst++] = (((v & 0xf800) >> 8) + ((v & 0x07e0) >> 3) + ((v & 0x001f) << 3)) <= th?0:1;
v = (int)input[pt_src++]; output[pt_dst++] = (((v & 0xf800) >> 8) + ((v & 0x07e0) >> 3) + ((v & 0x001f) << 3)) <= th?0:1;
v = (int)input[pt_src++]; output[pt_dst++] = (((v & 0xf800) >> 8) + ((v & 0x07e0) >> 3) + ((v & 0x001f) << 3)) <= th?0:1;
v = (int)input[pt_src++]; output[pt_dst++] = (((v & 0xf800) >> 8) + ((v & 0x07e0) >> 3) + ((v & 0x001f) << 3)) <= th?0:1;
v = (int)input[pt_src++]; output[pt_dst++] = (((v & 0xf800) >> 8) + ((v & 0x07e0) >> 3) + ((v & 0x001f) << 3)) <= th?0:1;
v = (int)input[pt_src++]; output[pt_dst++] = (((v & 0xf800) >> 8) + ((v & 0x07e0) >> 3) + ((v & 0x001f) << 3)) <= th?0:1;
v = (int)input[pt_src++]; output[pt_dst++] = (((v & 0xf800) >> 8) + ((v & 0x07e0) >> 3) + ((v & 0x001f) << 3)) <= th?0:1;
v = (int)input[pt_src++]; output[pt_dst++] = (((v & 0xf800) >> 8) + ((v & 0x07e0) >> 3) + ((v & 0x001f) << 3)) <= th?0:1;
}
//スキップ
pt_src += skip_src;
pt_dst += skip_dst;
}
return;
}
public void doFilter(INyARRaster i_input, int[] o_output, int l, int t, int w, int h)
{
Debug.Assert(i_input.isEqualBufferType(NyARBufferType.INT1D_X8R8G8B8_32));
NyARIntSize size = i_input.getSize();
int[] in_buf = (int[])i_input.getBuffer();
int bp = (l + t * size.w);
int v;
int b = t + h;
int row_padding_dst = (size.w - w);
int row_padding_src = row_padding_dst;
int pix_count = w;
int pix_mod_part = pix_count - (pix_count % 8);
int src_ptr = t * size.w + l;
for (int y = t; y < b; y++)
{
int x = 0;
for (x = pix_count - 1; x >= pix_mod_part; x--)
{
v = in_buf[src_ptr++]; o_output[bp++] = (((v >> 16) & 0xff) + ((v >> 8) & 0xff) + (v & 0xff)) >> 2;
}
for (; x >= 0; x -= 8)
{
v = in_buf[src_ptr++]; o_output[bp++] = (((v >> 16) & 0xff) + ((v >> 8) & 0xff) + (v & 0xff)) >> 2;
v = in_buf[src_ptr++]; o_output[bp++] = (((v >> 16) & 0xff) + ((v >> 8) & 0xff) + (v & 0xff)) >> 2;
v = in_buf[src_ptr++]; o_output[bp++] = (((v >> 16) & 0xff) + ((v >> 8) & 0xff) + (v & 0xff)) >> 2;
v = in_buf[src_ptr++]; o_output[bp++] = (((v >> 16) & 0xff) + ((v >> 8) & 0xff) + (v & 0xff)) >> 2;
v = in_buf[src_ptr++]; o_output[bp++] = (((v >> 16) & 0xff) + ((v >> 8) & 0xff) + (v & 0xff)) >> 2;
v = in_buf[src_ptr++]; o_output[bp++] = (((v >> 16) & 0xff) + ((v >> 8) & 0xff) + (v & 0xff)) >> 2;
v = in_buf[src_ptr++]; o_output[bp++] = (((v >> 16) & 0xff) + ((v >> 8) & 0xff) + (v & 0xff)) >> 2;
v = in_buf[src_ptr++]; o_output[bp++] = (((v >> 16) & 0xff) + ((v >> 8) & 0xff) + (v & 0xff)) >> 2;
}
bp += row_padding_dst;
src_ptr += row_padding_src;
}
return;
}
public void doFilter(INyARRaster i_input, int[] o_output, int l, int t, int w, int h)
{
Debug.Assert(i_input.isEqualBufferType(NyARBufferType.WORD1D_R5G6B5_16LE));
short[] input = (short[])i_input.getBuffer();
NyARIntSize s = i_input.getSize();
int skip_dst = (s.w - w);
int skip_src = skip_dst;
int pix_count = w;
int pix_mod_part = pix_count - (pix_count % 8);
//左上から1行づつ走査していく
int pt_dst = (t * s.w + l);
int pt_src = pt_dst;
for (int y = h - 1; y >= 0; y -= 1)
{
int x, v;
for (x = pix_count - 1; x >= pix_mod_part; x--)
{
v = (int)input[pt_src++]; o_output[pt_dst++] = (((v & 0xf800) >> 8) + ((v & 0x07e0) >> 3) + ((v & 0x001f) << 3)) >> 2;
}
for (; x >= 0; x -= 8)
{
v = (int)input[pt_src++]; o_output[pt_dst++] = (((v & 0xf800) >> 8) + ((v & 0x07e0) >> 3) + ((v & 0x001f) << 3)) >> 2;
v = (int)input[pt_src++]; o_output[pt_dst++] = (((v & 0xf800) >> 8) + ((v & 0x07e0) >> 3) + ((v & 0x001f) << 3)) >> 2;
v = (int)input[pt_src++]; o_output[pt_dst++] = (((v & 0xf800) >> 8) + ((v & 0x07e0) >> 3) + ((v & 0x001f) << 3)) >> 2;
v = (int)input[pt_src++]; o_output[pt_dst++] = (((v & 0xf800) >> 8) + ((v & 0x07e0) >> 3) + ((v & 0x001f) << 3)) >> 2;
v = (int)input[pt_src++]; o_output[pt_dst++] = (((v & 0xf800) >> 8) + ((v & 0x07e0) >> 3) + ((v & 0x001f) << 3)) >> 2;
v = (int)input[pt_src++]; o_output[pt_dst++] = (((v & 0xf800) >> 8) + ((v & 0x07e0) >> 3) + ((v & 0x001f) << 3)) >> 2;
v = (int)input[pt_src++]; o_output[pt_dst++] = (((v & 0xf800) >> 8) + ((v & 0x07e0) >> 3) + ((v & 0x001f) << 3)) >> 2;
v = (int)input[pt_src++]; o_output[pt_dst++] = (((v & 0xf800) >> 8) + ((v & 0x07e0) >> 3) + ((v & 0x001f) << 3)) >> 2;
}
//スキップ
pt_src += skip_src;
pt_dst += skip_dst;
}
return;
}
/**
*
* ラベリングの実体。
* @return
* ラベル数が上限に達したときはfalse
*/
private bool imple_labeling(INyARRaster i_raster, int i_th, int i_left, int i_top, int i_width, int i_height)
{
//ラスタのサイズを確認
Debug.Assert(i_raster.getSize().isEqualSize(this._raster_size));
//ラスタドライバのチェック
if (_last_input_raster != i_raster)
{
this._image_driver = (IRasterDriver)i_raster.createInterface(typeof(IRasterDriver));
}
IRasterDriver pixdrv = this._image_driver;
RleElement[] rle_prev = this._rle1;
RleElement[] rle_current = this._rle2;
// リセット処理
RleInfoStack rlestack = this._rlestack;
rlestack.clear();
//
int len_prev = 0;
int len_current = 0;
int bottom = i_top + i_height;
int id_max = 0;
int ypos = i_top;
// 初段登録
len_prev = pixdrv.xLineToRle(i_left, ypos, i_width, i_th, rle_prev);
for (int i = 0; i < len_prev; i++)
{
// フラグメントID=フラグメント初期値、POS=Y値、RELインデクス=行
if (addFragment(rle_prev[i], id_max, ypos, rlestack))
{
id_max++;
}
else
{
return(false);
}
}
NyARRleLabelFragmentInfo[] f_array = rlestack.getArray();
// 次段結合
for (int y = i_top + 1; y < bottom; y++)
{
// カレント行の読込
ypos++;
len_current = pixdrv.xLineToRle(i_left, ypos, i_width, i_th, rle_current);
int index_prev = 0;
for (int i = 0; i < len_current; i++)
{
// index_prev,len_prevの位置を調整する
int id = -1;
// チェックすべきprevがあれば確認
while (index_prev < len_prev)
{
if (rle_current[i].l - rle_prev[index_prev].r > 0) // 0なら8方位ラベリング
// prevがcurの左方にある→次のフラグメントを探索
{
index_prev++;
continue;
}
else if (rle_prev[index_prev].l - rle_current[i].r > 0) // 0なら8方位ラベリングになる
// prevがcur右方にある→独立フラグメント
{
if (addFragment(rle_current[i], id_max, y, rlestack))
{
id_max++;
}
else
{
return(false);
}
// 次のindexをしらべる
goto SCAN_CUR;
}
id = rle_prev[index_prev].fid; //ルートフラグメントid
NyARRleLabelFragmentInfo id_ptr = f_array[id];
//結合対象(初回)->prevのIDをコピーして、ルートフラグメントの情報を更新
rle_current[i].fid = id; //フラグメントIDを保存
//
int l = rle_current[i].l;
int r = rle_current[i].r;
int len = r - l;
//結合先フラグメントの情報を更新する。
id_ptr.area += len;
//tとentry_xは、結合先のを使うので更新しない。
id_ptr.clip_l = l < id_ptr.clip_l?l:id_ptr.clip_l;
id_ptr.clip_r = r > id_ptr.clip_r?r - 1:id_ptr.clip_r;
id_ptr.clip_b = y;
id_ptr.pos_x += (len * (2 * l + (len - 1))) / 2;
id_ptr.pos_y += y * len;
//多重結合の確認（２個目以降）
index_prev++;
while (index_prev < len_prev)
{
if (rle_current[i].l - rle_prev[index_prev].r > 0) // 0なら8方位ラベリング
// prevがcurの左方にある→prevはcurに連結していない。
{
goto SCAN_PREV;
}
else if (rle_prev[index_prev].l - rle_current[i].r > 0) // 0なら8方位ラベリングになる
// prevがcurの右方にある→prevはcurに連結していない。
{
index_prev--;
goto SCAN_CUR;
}
// prevとcurは連結している→ルートフラグメントの統合
//結合するルートフラグメントを取得
int prev_id = rle_prev[index_prev].fid;
NyARRleLabelFragmentInfo prev_ptr = f_array[prev_id];
if (id != prev_id)
{
//prevとcurrentのフラグメントidを書き換える。
for (int i2 = index_prev; i2 < len_prev; i2++)
{
//prevは現在のidから最後まで
if (rle_prev[i2].fid == prev_id)
{
rle_prev[i2].fid = id;
}
}
for (int i2 = 0; i2 < i; i2++)
{
//currentは0から現在-1まで
if (rle_current[i2].fid == prev_id)
{
rle_current[i2].fid = id;
}
}
//現在のルートフラグメントに情報を集約
id_ptr.area += prev_ptr.area;
id_ptr.pos_x += prev_ptr.pos_x;
id_ptr.pos_y += prev_ptr.pos_y;
//tとentry_xの決定
if (id_ptr.clip_t > prev_ptr.clip_t)
{
// 現在の方が下にある。
id_ptr.clip_t = prev_ptr.clip_t;
id_ptr.entry_x = prev_ptr.entry_x;
}
else if (id_ptr.clip_t < prev_ptr.clip_t)
{
// 現在の方が上にある。prevにフィードバック
}
else
{
// 水平方向で小さい方がエントリポイント。
if (id_ptr.entry_x > prev_ptr.entry_x)
{
id_ptr.entry_x = prev_ptr.entry_x;
}
else
{
}
}
//lの決定
if (id_ptr.clip_l > prev_ptr.clip_l)
{
id_ptr.clip_l = prev_ptr.clip_l;
}
else
{
}
//rの決定
if (id_ptr.clip_r < prev_ptr.clip_r)
{
id_ptr.clip_r = prev_ptr.clip_r;
}
else
{
}
//bの決定
//結合済のルートフラグメントを無効化する。
prev_ptr.area = 0;
}
index_prev++;
}
index_prev--;
break;
SCAN_PREV :;
}
// curにidが割り当てられたかを確認
// 右端独立フラグメントを追加
if (id < 0)
{
if (addFragment(rle_current[i], id_max, y, rlestack))
{
id_max++;
}
else
{
return(false);
}
}
SCAN_CUR :;
}
// prevとrelの交換
RleElement[] tmp = rle_prev;
rle_prev = rle_current;
len_prev = len_current;
rle_current = tmp;
}
//対象のラベルだけを追記
int max = this._max_area;
int min = this._min_area;
for (int i = id_max - 1; i >= 0; i--)
{
NyARRleLabelFragmentInfo src_info = f_array[i];
int area = src_info.area;
if (area < min || area > max) //対象外のエリア0のもminではじく
{
continue;
}
//値を相対位置に補正
src_info.clip_l += i_left;
src_info.clip_r += i_left;
src_info.entry_x += i_left;
src_info.pos_x /= area;
src_info.pos_y /= area;
//コールバック関数コール
this.OnLabelFound(src_info);
}
return(true);
}
public void doFilter(INyARRaster i_input, INyARRaster i_output)
{
this._do_filter_impl.doFilter(i_input, i_output, i_input.getSize());
}
public void doFilter(INyARRaster i_input, INyARRaster i_output)
{
Debug.Assert(i_input.getSize().isEqualSize(i_output.getSize()) == true);
this._dofilterimpl.doFilter(i_input, i_output, i_input.getSize());
}
public void switchRaster(INyARRaster i_ref_raster)
{
this._ref_buf = (int[])i_ref_raster.getBuffer();
this._ref_size = i_ref_raster.getSize();
}
public void doFilter(INyARRaster i_input, INyARRaster i_output)
{
this._do_filter_impl.doFilter(i_input,i_output,i_input.getSize());
}
public NyARVectorReader_INT1D_GRAY_8(INyARRaster i_ref_raster)
{
Debug.Assert(i_ref_raster.getBufferType() == NyARBufferType.INT1D_GRAY_8);
this._ref_buf = (int[])(i_ref_raster.getBuffer());
this._ref_size = i_ref_raster.getSize();
}
/**
* NyARRasterからパターンデータをセットします。
* この関数は、データを元に所有するデータ領域を更新します。
* @param i_buffer
*/
public void setRaster(INyARRaster i_raster)
{
//画素フォーマット、サイズ制限
Debug.Assert(i_raster.isEqualBufferType(NyARBufferType.INT1D_X8R8G8B8_32));
Debug.Assert(i_raster.getSize().isEqualSize(i_raster.getSize()));
int[] buf = (int[])i_raster.getBuffer();
//i_buffer[XRGB]→差分[R,G,B]変換
int i;
int ave;//<PV/>
int rgb;//<PV/>
int[] linput = this._data;//<PV/>
// input配列のサイズとwhも更新// input=new int[height][width][3];
int number_of_pixels = this._number_of_pixels;
int for_mod = this._optimize_for_mod;
//<平均値計算(FORの1/8展開)>
ave = 0;
for (i = number_of_pixels - 1; i >= for_mod; i--)
{
rgb = buf[i]; ave += ((rgb >> 16) & 0xff) + ((rgb >> 8) & 0xff) + (rgb & 0xff);
}
for (; i >= 0; )
{
rgb = buf[i]; ave += ((rgb >> 16) & 0xff) + ((rgb >> 8) & 0xff) + (rgb & 0xff); i--;
rgb = buf[i]; ave += ((rgb >> 16) & 0xff) + ((rgb >> 8) & 0xff) + (rgb & 0xff); i--;
rgb = buf[i]; ave += ((rgb >> 16) & 0xff) + ((rgb >> 8) & 0xff) + (rgb & 0xff); i--;
rgb = buf[i]; ave += ((rgb >> 16) & 0xff) + ((rgb >> 8) & 0xff) + (rgb & 0xff); i--;
rgb = buf[i]; ave += ((rgb >> 16) & 0xff) + ((rgb >> 8) & 0xff) + (rgb & 0xff); i--;
rgb = buf[i]; ave += ((rgb >> 16) & 0xff) + ((rgb >> 8) & 0xff) + (rgb & 0xff); i--;
rgb = buf[i]; ave += ((rgb >> 16) & 0xff) + ((rgb >> 8) & 0xff) + (rgb & 0xff); i--;
rgb = buf[i]; ave += ((rgb >> 16) & 0xff) + ((rgb >> 8) & 0xff) + (rgb & 0xff); i--;
}
//<平均値計算(FORの1/8展開)/>
ave = number_of_pixels * 255 * 3 - ave;
ave = 255 - (ave / (number_of_pixels * 3));//(255-R)-ave を分解するための事前計算
int sum = 0, w_sum;
int input_ptr = number_of_pixels * 3 - 1;
//<差分値計算(FORの1/8展開)>
for (i = number_of_pixels - 1; i >= for_mod; i--)
{
rgb = buf[i];
w_sum = (ave - (rgb & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum;//B
w_sum = (ave - ((rgb >> 8) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum;//G
w_sum = (ave - ((rgb >> 16) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum;//R
}
for (; i >= 0; )
{
rgb = buf[i]; i--;
w_sum = (ave - (rgb & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum;//B
w_sum = (ave - ((rgb >> 8) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum;//G
w_sum = (ave - ((rgb >> 16) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum;//R
rgb = buf[i]; i--;
w_sum = (ave - (rgb & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum;//B
w_sum = (ave - ((rgb >> 8) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum;//G
w_sum = (ave - ((rgb >> 16) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum;//R
rgb = buf[i]; i--;
w_sum = (ave - (rgb & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum;//B
w_sum = (ave - ((rgb >> 8) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum;//G
w_sum = (ave - ((rgb >> 16) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum;//R
rgb = buf[i]; i--;
w_sum = (ave - (rgb & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum;//B
w_sum = (ave - ((rgb >> 8) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum;//G
w_sum = (ave - ((rgb >> 16) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum;//R
rgb = buf[i]; i--;
w_sum = (ave - (rgb & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum;//B
w_sum = (ave - ((rgb >> 8) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum;//G
w_sum = (ave - ((rgb >> 16) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum;//R
rgb = buf[i]; i--;
w_sum = (ave - (rgb & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum;//B
w_sum = (ave - ((rgb >> 8) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum;//G
w_sum = (ave - ((rgb >> 16) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum;//R
rgb = buf[i]; i--;
w_sum = (ave - (rgb & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum;//B
w_sum = (ave - ((rgb >> 8) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum;//G
w_sum = (ave - ((rgb >> 16) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum;//R
rgb = buf[i]; i--;
w_sum = (ave - (rgb & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum;//B
w_sum = (ave - ((rgb >> 8) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum;//G
w_sum = (ave - ((rgb >> 16) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum;//R
}
//<差分値計算(FORの1/8展開)/>
double p = Math.Sqrt((double)sum);
this._pow = p != 0.0 ? p : 0.0000001;
return;
}
public void switchRaster(INyARRaster i_ref_raster)
{
this._ref_buf = (int[])i_ref_raster.getBuffer();
this._ref_size = i_ref_raster.getSize();
}
/**
* NyARRasterからパターンデータをセットします。
* この関数は、データを元に所有するデータ領域を更新します。
* @param i_buffer
*/
public void setRaster(INyARRaster i_raster)
{
//画素フォーマット、サイズ制限
Debug.Assert(i_raster.isEqualBufferType(NyARBufferType.INT1D_X8R8G8B8_32));
Debug.Assert(i_raster.getSize().isEqualSize(i_raster.getSize()));
int[] buf = (int[])i_raster.getBuffer();
//i_buffer[XRGB]→差分[R,G,B]変換
int i;
int ave; //<PV/>
int rgb; //<PV/>
int[] linput = this._data; //<PV/>
// input配列のサイズとwhも更新// input=new int[height][width][3];
int number_of_pixels = this._number_of_pixels;
int for_mod = this._optimize_for_mod;
//<平均値計算(FORの1/8展開)>
ave = 0;
for (i = number_of_pixels - 1; i >= for_mod; i--)
{
rgb = buf[i]; ave += ((rgb >> 16) & 0xff) + ((rgb >> 8) & 0xff) + (rgb & 0xff);
}
for (; i >= 0;)
{
rgb = buf[i]; ave += ((rgb >> 16) & 0xff) + ((rgb >> 8) & 0xff) + (rgb & 0xff); i--;
rgb = buf[i]; ave += ((rgb >> 16) & 0xff) + ((rgb >> 8) & 0xff) + (rgb & 0xff); i--;
rgb = buf[i]; ave += ((rgb >> 16) & 0xff) + ((rgb >> 8) & 0xff) + (rgb & 0xff); i--;
rgb = buf[i]; ave += ((rgb >> 16) & 0xff) + ((rgb >> 8) & 0xff) + (rgb & 0xff); i--;
rgb = buf[i]; ave += ((rgb >> 16) & 0xff) + ((rgb >> 8) & 0xff) + (rgb & 0xff); i--;
rgb = buf[i]; ave += ((rgb >> 16) & 0xff) + ((rgb >> 8) & 0xff) + (rgb & 0xff); i--;
rgb = buf[i]; ave += ((rgb >> 16) & 0xff) + ((rgb >> 8) & 0xff) + (rgb & 0xff); i--;
rgb = buf[i]; ave += ((rgb >> 16) & 0xff) + ((rgb >> 8) & 0xff) + (rgb & 0xff); i--;
}
//<平均値計算(FORの1/8展開)/>
ave = number_of_pixels * 255 * 3 - ave;
ave = 255 - (ave / (number_of_pixels * 3));//(255-R)-ave を分解するための事前計算
int sum = 0, w_sum;
int input_ptr = number_of_pixels * 3 - 1;
//<差分値計算(FORの1/8展開)>
for (i = number_of_pixels - 1; i >= for_mod; i--)
{
rgb = buf[i];
w_sum = (ave - (rgb & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum; //B
w_sum = (ave - ((rgb >> 8) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum; //G
w_sum = (ave - ((rgb >> 16) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum; //R
}
for (; i >= 0;)
{
rgb = buf[i]; i--;
w_sum = (ave - (rgb & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum; //B
w_sum = (ave - ((rgb >> 8) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum; //G
w_sum = (ave - ((rgb >> 16) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum; //R
rgb = buf[i]; i--;
w_sum = (ave - (rgb & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum; //B
w_sum = (ave - ((rgb >> 8) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum; //G
w_sum = (ave - ((rgb >> 16) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum; //R
rgb = buf[i]; i--;
w_sum = (ave - (rgb & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum; //B
w_sum = (ave - ((rgb >> 8) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum; //G
w_sum = (ave - ((rgb >> 16) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum; //R
rgb = buf[i]; i--;
w_sum = (ave - (rgb & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum; //B
w_sum = (ave - ((rgb >> 8) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum; //G
w_sum = (ave - ((rgb >> 16) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum; //R
rgb = buf[i]; i--;
w_sum = (ave - (rgb & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum; //B
w_sum = (ave - ((rgb >> 8) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum; //G
w_sum = (ave - ((rgb >> 16) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum; //R
rgb = buf[i]; i--;
w_sum = (ave - (rgb & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum; //B
w_sum = (ave - ((rgb >> 8) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum; //G
w_sum = (ave - ((rgb >> 16) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum; //R
rgb = buf[i]; i--;
w_sum = (ave - (rgb & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum; //B
w_sum = (ave - ((rgb >> 8) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum; //G
w_sum = (ave - ((rgb >> 16) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum; //R
rgb = buf[i]; i--;
w_sum = (ave - (rgb & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum; //B
w_sum = (ave - ((rgb >> 8) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum; //G
w_sum = (ave - ((rgb >> 16) & 0xff)); linput[input_ptr--] = w_sum; sum += w_sum * w_sum; //R
}
//<差分値計算(FORの1/8展開)/>
double p = Math.Sqrt((double)sum);
this._pow = p != 0.0 ? p : 0.0000001;
return;
}
public void doThFilter(INyARRaster i_raster, int i_l, int i_t, int i_w, int i_h, int i_th, INyARRaster o_raster)
{
byte[] input = (byte[])i_raster.getBuffer();
int[] output = (int[])o_raster.getBuffer();
NyARIntSize s = i_raster.getSize();
int th = i_th;
int skip_dst = (s.w - i_w);
int skip_src = skip_dst * 4;
int pix_count = i_w;
int pix_mod_part = pix_count - (pix_count % 8);
//左上から1行づつ走査していく
int pt_dst = (i_t * s.w + i_l);
int pt_src = pt_dst * 4;
for (int y = i_h - 1; y >= 0; y -= 1)
{
int x;
for (x = pix_count - 1; x >= pix_mod_part; x--)
{
output[pt_dst++] = input[pt_src + 0] <= th ? 0 : 1;
pt_src += 4;
}
for (; x >= 0; x -= 8)
{
output[pt_dst++] = input[pt_src + 0] <= th ? 0 : 1;
pt_src += 4;
output[pt_dst++] = input[pt_src + 0] <= th ? 0 : 1;
pt_src += 4;
output[pt_dst++] = input[pt_src + 0] <= th ? 0 : 1;
pt_src += 4;
output[pt_dst++] = input[pt_src + 0] <= th ? 0 : 1;
pt_src += 4;
output[pt_dst++] = input[pt_src + 0] <= th ? 0 : 1;
pt_src += 4;
output[pt_dst++] = input[pt_src + 0] <= th ? 0 : 1;
pt_src += 4;
output[pt_dst++] = input[pt_src + 0] <= th ? 0 : 1;
pt_src += 4;
output[pt_dst++] = input[pt_src + 0] <= th ? 0 : 1;
pt_src += 4;
}
//スキップ
pt_src += skip_src;
pt_dst += skip_dst;
}
return;
}
/**
* o_histogramにヒストグラムを出力します。
* @param i_input
* @param o_histogram
* @return
* @throws NyARException
*/
public int analyzeRaster(INyARRaster i_input,NyARHistogram o_histogram)
{
NyARIntSize size=i_input.getSize();
//最大画像サイズの制限
Debug.Assert(size.w*size.h<0x40000000);
Debug.Assert(o_histogram.length==256);//現在は固定
int[] h=o_histogram.data;
//ヒストグラム初期化
for (int i = o_histogram.length-1; i >=0; i--){
h[i] = 0;
}
o_histogram.total_of_data=size.w*size.h/this._vertical_skip;
return this._histImpl.createHistogram(i_input, size,h,this._vertical_skip);
}
public virtual void doFilter(INyARRaster i_input, INyARRaster i_output)
{
Debug.Assert (i_input.getSize().isEqualSize(i_output.getSize()) == true);
this._dofilterimpl.doFilter(i_input,i_output,i_input.getSize());
}
public void doFilter(INyARRaster i_input, INyARRaster i_output)
{
Debug.Assert(i_input != i_output);
this._do_filter_impl.doFilter(i_input, i_output, i_input.getSize());
}
public void doFilter(INyARRaster i_input, INyARRaster i_output)
{
Debug.Assert (i_input!=i_output);
this._do_filter_impl.doFilter(i_input,i_output,i_input.getSize());
}
public NyARVectorReader_INT1D_GRAY_8(INyARRaster i_ref_raster)
{
Debug.Assert(i_ref_raster.getBufferType() == NyARBufferType.INT1D_GRAY_8);
this._ref_buf = (int[])(i_ref_raster.getBuffer());
this._ref_size = i_ref_raster.getSize();
}
