Filtre gauss (flou)

Version imprimable

17/06/2010, 09h43
Opencl_jedi

Filtre gauss (flou)

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?

Bonjour,

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

Code:

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.;
    }
}

17/06/2010, 11h58
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.
17/06/2010, 12h12
Invité(e)

Peux tu nous montrer ton code ??

OK voilà mon code

Code:

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 };