画像処理プログラム

質問

以下の問題でプログラムを改良したのですが合っていますか？間違ってたら教えてください

図3 は、画像中の全ての画素に対する、濃度値の頻度分布を示すヒストグラムの例である。サンプルプロググラムhistogram.c をもとに、入力画像の濃度値のヒストグラムを作成するプログラムを完成させる。
1. Histogram.c をそのまま実行すると、ダミーデータを用いた結果画像が出力される。
$ gcc □ Histogram.c -o □ Histogram ↓
$ ./Histogram ○○○.pgm □ Histogram.pgm ↓
2. Histogram.c を次のように変更し、プログラムを完成させる。
（a）入力画像のヒストグラムを求める。
各画素を参照し、その画素値img[i][j] に対応する番号の配列histogram[0255] の値を1 増やす。
（b）頻度の最大値を求める。
濃度値i=0255 について調べ、もし頻度histogram[i] の値がそれまでの最大値max frequency より大きければ、最大値max frequency を更新する。
（c）ヒストグラム画像を作成する。
最大頻度を1 としたときの割合（histogram[i]/max frequency）を求め、それに画像の大きさH_height を掛けることで各濃度値のグラフの長さlength を求める。
length = H height * histogram[i] / max frequency

#include <stdio.h>
#include <stdlib.h>

#define V_width 320
#define V_height 240
#define amp 1.0

#define GLEVEL 256 /* 最大階調 */
#define MAX_BRIGHTNESS 255

#define H_width GLEVEL
#define H_height 256

int d[9];
int i, j, dat;
float z, zz;

unsigned char img [V_height][V_width]; /* input image */
unsigned char work[H_height][H_width]={0}; /* work space */

long int histogram[GLEVEL]={0}; /* ヒストグラム */
long int max_frequency=0; /* 頻度の最大値 */
float length; /* 頻度を表すグラフの長さ */

int main(int argc, char* argv[])
{
FILE *fpi, *fpo;
unsigned char data;

char str[ 256];
int width, height;

/* check arg number */
if (argc != 3) {
fprintf(stderr, "Usage: %s [input] [output]\n", argv[0]);
exit(1);
}

/* open input file */
if ((fpi = fopen(argv[1], "rb")) == NULL) {
fprintf(stderr, "input file open error!\n");
exit(1);
}

/* open output file */
if ((fpo = fopen(argv[2], "wb")) == NULL) {
fprintf(stderr, "output file open error!\n");
exit(1);
}

/* read PGM header */
while(1){
fgets(str, 256, fpi);
if(str[0] == '#'){
fprintf( stderr, "%s", str);
}
else if( !strncmp( str, "P5", 2)){
fprintf( stderr, "This file is a PGM file.\n");
}
else if( !strncmp( str, "255", 3)){
fprintf( stderr, "The file is opened.\n");
break;
}
else{
sscanf( str, "%d %d", &width, &height);
fprintf( stderr, "SIZE: %3d x %3d\n", width, height);
}
}

/* read data */
for (i=0; i<height; i++) {
for (j=0; j<width; j++) {
fread( &data, 1, 1, fpi);
img[i][j] = data;
}
}

/* (a) 入力画像のヒストグラムを求める */
for(i=0; i<height; i++){
for(j=0; j<width; j++){
// histogram([i][j])++;
}
}

/* (b) 頻度の最大値を求める */
for ( i=0; i<GLEVEL; i++) {
// if (histogram[i]>max_frequency) max_frequency = histogram[i];
}

/* (c) ヒストグラム画像を作成する */
for (j=0; j<H_width; j++) {
length = (float) H_height * j / H_width; //仮のデータ
for (i=0; i<length; i++) {
work[H_height-1-i][j] = MAX_BRIGHTNESS;
length = H_height * histogram[i] / max_frequency
}

}

/* write data */
fprintf(fpo,"P5\n %d %d\n 255\n", H_width, H_height);
for (i= 0; i<H_height; i++) {
for (j= 0; j<H_width; j++) {
fwrite( &work[i][j], 1, 1, fpo);
}
}

fclose(fpi);
fclose(fpo);
return 0;
}

回答No.1

noname#208507 · ベストアンサー

「間違ってたら」も何も、まともに実行できないですよ。
コメント外してませんし。
それに「図3 は」と言われても。

とりあえず、濃度値の頻度の意味を理解しましょう。
(a)は、histogram[img[i][j]]++; です。

(b)は（コメントを外せば）よいとして。

(c)は、下記をヒントに。

length = H_height * histogram[j] / max_frequency;
for (i=0; i<H_height; i++) {
if (i <= length) work[H_height-1-i][j] = 255;
else work[H_height-1-i][j] = 0;
}