編集履歴

回答編集履歴

説明の誤字を修正

2020/06/28 11:36

投稿

スコア84

test CHANGED Viewed

@@ -1,9 +1,17 @@
 回答が長くなってしまったので、「詰め直し」は別に回答します。
-こちらは府全般的な参考に、ということで
+こちらは全般的な参考に、ということで。
+なお、コンパイル前に先頭付近のマクロ NEIGHBOR により、4連結、8連結のどちらかを選んで実行できます。
 ```C
+#define NEIGHBOR 4		/*連結数:4か8*/
+```
+```C
 /*
  ・マクロの使い方など、かなり行儀が悪い部分がありますがご了承ください。

コードを少々修正し、説明を別枠で行うことを追記

2020/06/28 11:35

投稿

e-watt

スコア84

test CHANGED Viewed

@@ -1,8 +1,6 @@
-肝心の「ルックアップテーブルの詰め直し」で「なぜ」「何を」やっているかの説明が未記入ですがとりあえずコードだけ這っておきます。
-説明は追って記します。
+回答が長くなってしまったので、「詰め直し」は別に回答します。
+こちらは府全般的な参考に、ということで
 ```C
@@ -10,10 +8,6 @@
  ・マクロの使い方など、かなり行儀が悪い部分がありますがご了承ください。
- ・変数名などかなりいじってしまいました。悪しからず
- ・インデントなども雑です。すみません
  */
 #include<stdio.h>
@@ -24,11 +18,11 @@
-#define NEIGHBOR 4		/*連結数を4か8で指定*/
+#define NEIGHBOR 4		/*連結数:4か8*/
-/* デバック用：画像っぽいものをなんとなく表示 */
+/* デバッグ用 */
 static void showIntMatrix(const char* name, const int* array, int width, int height, int line) {
@@ -66,11 +60,11 @@
-/*デバッグ用：配列の要素を横一列に表示*/
+/*デバッグ用*/
 static void showIntArray(const char* caption, const char* name, const int* array, int N) {
-	printf("%s:\n", caption);
+	printf("%s:\n%s ", caption, name);
 	for (int i=1; i<=N; i++){
@@ -82,11 +76,11 @@
 }
-#define	showDST(expr)		showIntArray(expr, "dst", dst, ((w*h)+1)/2)
+#define	showDST(expr)		showIntArray(expr, "Dst", dst, ((t)+1)/2)
 #define	PartN (14)
-#define	showPartOfDST(expr)	showIntArray(expr, "dst", dst, PartN)	/*全部出すと冗長なので先頭付近のみ表示*/
+#define	showPartOfDST(expr)	showIntArray(expr, "Dst", dst, PartN)	/*全部出すと冗長なので先頭付近のみ表示*/
@@ -128,7 +122,7 @@
   fscanf(fp, "%d", &m); //画像の輝度の階調
-  fgets(junk,10,fp); // 改行を読み飛ばす
+  fgets(junk,10,fp); //改行を読み飛ばす
 	printf("J[%s]len(%d)\n", junk, strlen(junk));
@@ -252,11 +246,11 @@
   fp = fopen("test.pgm", "rb");
-  fgets(junk,10,fp); // ヘッダを読み飛ばす
+  fgets(junk,10,fp); //ヘッダを読み飛ばす
-  fgets(junk,10,fp); // ヘッダを読み飛ばす
+  fgets(junk,10,fp); //ヘッダを読み飛ばす
-  fgets(junk,10,fp); // ヘッダを読み飛ばす
+  fgets(junk,10,fp); //ヘッダを読み飛ばす
@@ -332,7 +326,7 @@
   int Label=0;	//ラベル番号
-  int NL = ((w*h)+1)/2;		// ラベルの最大数
+  int NL = ((w*h)+1)/2;		//ラベルの最大数
   int dst[NL+1];//ラベリング配列 ; Srcは書き換えないし連番だから、Dstしか用意しないぜ
@@ -340,11 +334,11 @@
   for (i=1;i<=NL;i++){
-  	dst[i] = i;	// [1]から[((w*h)+1)/2]にその値を入れるぜ;
+  	dst[i] = i;	//[1]から[((w*h)+1)/2]にその値を入れるぜ;
-  					// ただし質問者様は要素数と末尾の処理を間違ってるっぽい
+  					//ただし質問者様は要素数と末尾の処理を間違ってるっぽい
-  					// + 領域数の最大値は、市松模様に塗り分けられた時の ceil(W*H/2)  3x3マスだと5個まであり
+  					//+ 領域数の最大値は、市松模様に塗り分けられた時の ceil(W*H/2)  3x3マスだと5個まであり
   }
@@ -370,7 +364,7 @@
-/* (x,y)のラベル番号を与える。画像の範囲外に対してをゼロを与える */
+/* (x,y)のラベル番号を与える。画像の範囲外に対してはゼロを与える */
 #define getLabel(x,y)	((0<=(x)&&(x)<w && 0<=(y)) ? labeled_image[x][y] : 0)
@@ -388,13 +382,13 @@
       if(image[j][i]==1){	//対象画素が1
-			int	up_m = getLabel(j, i-1);	// 上	4／8連結で共通
+			int	up_m = getLabel(j, i-1);	//上	4／8連結で共通
-			int left = getLabel(j-1, i);	// 左	4／8連結で共通
+			int left = getLabel(j-1, i);	//左	4／8連結で共通
 #if NEIGHBOR == 4 //4連結
-        	if (up_m == 0 && left == 0) {	// コメントは8連結側を参照されたし
+        	if (up_m == 0 && left == 0) {	//コメントは8連結側を参照されたし
         		labeled_image[j][i] = ++Label;
@@ -410,7 +404,7 @@
         		labeled_image[j][i] = min;
-        		dst[up_m] = nonzeroMin(dst[up_m], min);		// 双方ゼロではないので実は普通のmin(a,b)でよい
+        		dst[up_m] = nonzeroMin(dst[up_m], min);		//双方ゼロではないので実は普通のmin(a,b)でよい
         		dst[left] = nonzeroMin(dst[left], min);
@@ -420,19 +414,19 @@
 #elif NEIGHBOR == 8
-			int	up_l = getLabel(j-1, i-1);	// 左上
+			int	up_l = getLabel(j-1, i-1);	//左上
-			int up_r = getLabel(j+1, i-1);	// 右上
+			int up_r = getLabel(j+1, i-1);	//右上
-      		// 『』内は参照ページの記述の引用；なるべく書いてある通りに作るぜ
+      		//『』内は参照ページの記述の引用；なるべく書いてある通りに作るぜ
-        	// 『左上、上、右上、左の画素のラベル番号を参照し、全て０（ゼロ）の場合は』
+        	//『左上、上、右上、左の画素のラベル番号を参照し、全て０（ゼロ）の場合は』
         	if (up_l == 0 && up_m == 0 && up_r == 0 && left == 0) {
-        		labeled_image[j][i] = ++Label;	// 『最後に割り振った番号＋１のラベル番号を割り振ります。』
+        		labeled_image[j][i] = ++Label;	//『最後に割り振った番号＋１のラベル番号を割り振ります。』
 _TL_
@@ -454,11 +448,11 @@
         		min = nonzeroMin(left, min);
-        		labeled_image[j][i] = min;		// （0以外の）『最小の番号を割り振ります。』
+        		labeled_image[j][i] = min;		//（0以外の）『最小の番号を割り振ります。』
 _Tm_
-        		// 『用いなかったラベル番号のルックアップテーブルの番号を最小の番号に書き換えます。』
+        		//『用いなかったラベル番号のルックアップテーブルの番号を最小の番号に書き換えます。』
         		/* ループ内でdst[0] は未使用だ（参照しない）からチェックせずにガンガン上書きするぜ*/
@@ -484,7 +478,7 @@
     }
-  	showLabeledImage();	// 一行操作するごとに表示する…と、解説ページと似た感じの出力に
+  	showLabeledImage();	//一行操作するごとに表示する…と、解説ページと似た感じの出力に
   }
@@ -500,15 +494,15 @@
-/* ルックアップテーブルの詰め替え；思いのほか面倒；もっといいやり方はありそう */
+/* ルックアップテーブルの詰め替え */
-  int conv[NL+1];	//
+  int conv[NL+1];	//
   for (i=1;i<=NL;i++){
-  	conv[i] = -1;	// すべて「未定」にする
+  	conv[i] = -1;	//すべて「未定」にする
   }
@@ -518,56 +512,54 @@
-  //再帰的結合  ;
+	//再帰的結合; ←これは回答者の勝手な命名で、専門用語ではない
-  for (i=1;i<=NL;i++){
-	  for (j=1;j<=NL;j++){
+	for (j=1;j<=NL;j++){
-	  	while (dst[j] != dst[ dst[j] ])
+		while (dst[j] != dst[ dst[j] ])
-	  		dst[j] = dst[ dst[j] ];
+			dst[j] = dst[ dst[j] ];
-	  }
+	}
+  showPartOfDST("recursed");
+	//番号を圧縮; 説明しないとわからんと思う
+	int seqn = 1;
+  for (i=1;i<=13;i++){
+  	int dstQ = dst[i];
+  	printf("[%2d] dstQ=%2d conv[dstQ]=%2d...",i, dstQ, conv[dstQ]);
+  	if (conv[ dstQ ] == -1) {
+  		conv[ dstQ ] = seqn;
+  		printf("conv[%2d] gets-to %2d", dstQ, seqn);
+  		++seqn;
+  	}else{
+  		printf("%19s", "");
+  	}
+  	printf(" dst[%2d]<<%2d, dstQ=%2\n", i, conv[dstQ], dstQ);
+  	dst[i] = conv[dstQ];
   }
-  showPartOfDST("recursed");
-	// 番号を圧縮; 説明しないとわからんと思う
-	int seqn = 1;
-  for (i=1;i<=13;i++){
-  	int dstQ = dst[i];
-  	printf("[%2d] dstQ=%2d conv[dstQ]=%2d...",i, dstQ, conv[dstQ]);
-  	if (conv[ dstQ ] == -1) {
-  		conv[ dstQ ] = seqn;
-  		printf("conv[%2d] gets-to %2d", dstQ, seqn);
-  		++seqn;
-  	}else{
-  		printf("%19s", "");
-  	}
-  	printf(" dst[%2d]<<%2d, dstQ=%2\n", i, conv[dstQ], dstQ);
-  	dst[i] = conv[dstQ];
-  }
 	showIntArray("Conv:", "conv", conv, PartN);
@@ -578,9 +570,7 @@
-  // ルックアップテーブルに従い結果を書き換え
+  //ルックアップテーブルに従い結果を書き換え
-	// 上でさんざんdst[0]に好き放題したが、ここでは参照しうるので本来の値を入れる
 	dst[0] = 0;
@@ -594,13 +584,13 @@
 	}
-  // 最終結果を表示
+  //最終結果を表示
   	showLabeledImage();
-	printf("%d連結の連結成分の個数=%d\n", NEIGHBOR, dst[Label]);	/*使われたラベルの値（最大値）の圧縮結果が領域の数*/
+	printf("%d連結の連結成分の個数=%d?\n", NEIGHBOR, dst[Label]);	/*使われたラベルの値（最大値）の圧縮結果が領域の数…とは限らなそう*/