Discussion :

[Opencl_jedi]

Bonjour,

j'essaie de faire un flou avec un filtre gauss sur une image en C,
mais tous ce que je trouve c'est une implémentation en 2 passage sur l'image (horizontalement puis verticalement), je voudrais faire un filtre 2-dim en un seul passage.

avez-vous une idée sur comment faire?

[Invité(e)]

Bonjour,

Il suffit d'imbriquer plein de boucles for ensemble. Voilà à quoi ca peut ressembler :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
 
double noyau[5][5] = {
 1,  4,  7,  4,  1,
 4, 16, 26, 16,  4,
 7, 26, 41, 26,  7, 
 4, 16, 26, 16,  4,
 1,  4,  7,  4,  1 };
 
double **image_source = ...;
double **image_floue = ...;
 
/* On commence à 2 au lieu de 0 et on finit à hauteur -2 au lieu de hauteur 
    Pour ne pas s'embêter avec la gestion des bords */
for(j = 2; j < hauteur - 2; ++j) {
    /* idem */
    for(i = 2; i < largeur - 2; ++i) {
        /* initialisation du pixel de destination */
        image_floue[j][i] = 0;
 
        /* application du noyau */
        for(ii = 0; ii < 5; ++ii) {
            for(jj = 0; jj < 5; ++jj) {
                /* application de chaque élément du noyau sur l'image source
                    dans l'image finale */
                image_floue[j][i] += noyau[jj][ii] * image_source[i+ii-2][j+jj-2];
            }
        }
        /* normalisation : somme des éléments du noyau = 273. */
        image_floue[j][i] /= 273.;
    }
}

[Opencl_jedi]

Merci pour la réponse, j'ai essayé ta méthode, en 5x5 puis en 3x3 avec des coefs que j'ai trouvé sur le net, mais le résultat n'est pas bon, l'image n'est pas floutée.

[Invité(e)]

Peux tu nous montrer ton code ??

[Opencl_jedi]

OK voilà mon code

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
void gauss(
    unsigned int* uiDataIn,
    unsigned int* uiDataOut,
    float *h_Kernel,
    int kernelR,
    int Height,
    int Width
){
	for(int y = 2; y < Height - 2; y++){
             for(int x = 2; x < Width - 2; x++){
            double sum = 0;
 
		for(int ky = -kernelR; ky <= kernelR; ky++){
			for(int kx = -kernelR; kx <= kernelR; kx++){
				int dx = x + kx;
				int dy = y + ky;
 
 sum += uiDataIn[dy*Width + dx] * h_Kernel[(kernelR+ky)*5 + (kernelR+kx)];
				}   
			}
 
			uiDataOut[y * Width + x] = (float)sum/273;
		}
	}
}

kernelR = 2;
h_kernel[25] = {
1, 4, 7, 4, 1,
4, 16, 26, 16, 4,
7, 26, 41, 26, 7,
4, 16, 26, 16, 4,
1, 4, 7, 4, 1 };