Matにおける同一画素値とその座標の取得

Question

前提・実現したいこと

CV_32SC2型のMat(以下img)を取り扱うプログラムになります。
imgの１チャネル目には0~PROJ_WIDTHPROJ_HEIGHTの画素値が重複ありで格納されています。(PROJ_WIDTHPROJ_HEIGHTは200万ほどです。)
ただし、画素値が重複する画素同士は必ず隣接しています。
このようなsrcにおいて画素値が同じ領域の重心座標を計算し、２チャネル目に重心座標に相当する画素に目印となる値を格納するようなプログラムを考えています。
愚直に記述すると下記のような形になると思うのですが、for文が3重になっていたりと複雑なコードになってしまいます。
もっとシンプルな形で記述をしたいためアドバイスをお願いします。
初心者のため読みにくいコードだと思いますがよろしくお願いいたします。
###試したこと
for (int y = 0;y < CAMERA_HEIGHT;y++)
for (int x = 0;x < CAMERA_WIDTH;x++)
の2重ループの部分を簡略化するために、「指定した画素値をもつピクセルの総数を返す関数」、「指定した画素値をもつピクセルの座標を返す関数」のようなものがあれば簡略化できると思い探してみたのですが見つけることができませんでした。もし、そのような関数があるのでしたら教えていただきたいです。

該当のソースコード

読みやすいように一部省略しています。

//Mat imgはあらかじめ用意してあります。
	for (int proj_index = 0;proj_index < PROJ_HEIGHT * PROJ_WIDTH;proj_index++)
	{
		double sum_x = 0;
		double sum_y = 0;
		double ave_x = 0;
		double ave_y = 0;
		int count = 0;

		for (int y = 0;y < CAMERA_HEIGHT;y++)
		{
			for (int x = 0;x < CAMERA_WIDTH;x++)
			{
				if (img.data[y*CAMERA_WIDTH+x*img.channels()] == proj_index)
				{
					sum_x += x;
					sum_y += y;
					count++;
				}
				ave_x = round(sum_x / count);
				ave_y = round(sum_y / count);
				img.at<int[2]>(ave_y, ave_x)[1] = proj_index;
			}
		}
	}

Answer 1

画素値が重複する画素同士は必ず隣接しています。

とのことであれば，まずはラベリングを行い，次にその結果の領域毎に重心を計算すれば良いのではないでしょうか．

Comments

[kpkq307]

fanaさん
ご回答ありがとうございます。
ラベリング処理についてなのですがCV::conectedComponents関数をつかうのでしょうか？
その場合だと入力は２値画像になるため、画素値が0と1以外の値も持つ今回の場合は適用できないのではないでしょうか？
初歩的な質問で申し訳ありませんが、ご返信いただけると幸いです。

[fana]

私が知らないだけなのかもしれませんが，OpenCVには多値画像に対する扱いやすいラベリングの実装が無いように見えます．（逆に，見つけたら教えていただきたい）

PROJ_WIDTH*PROJ_HEIGHT　なる値がかなり小さい（数パターンとか）ならば，
「"ある特定の画素値を持つか否か"という判定による二値化を行ってはconnectedComponents()を使う」
ということを繰り返す，という手はあるでしょうが……

PROJ_WIDTH*PROJ_HEIGHT がある程度大きいようであれば，ラベリング処理を自前で実装することを考えた方が良いかもしれません．

[fana]

（状況次第では floodfill() あたりを利用できる可能性もあるかも？）