Code source pour appliquer une DCT

**lemenyau** · 30/05/2006, 18h39

Je suis en train de programmer un compresseur d'image et j'ai un problème avec ma DCT

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
 
struct s_dct dct(struct s_image bloc_in[],int k)//DCT
{
    int x,y;
    struct s_dct bloc_out;
 
    bloc_out.bleu = bloc_out.rouge = bloc_out.vert = 0;
    for(x=0;x<8;x++)
    {
 
                    for(y=0;y<8;y++)
                    {
                                    bloc_out.bleu=bloc_out.bleu + (bloc_in[x*y].bleu*cos((((2*x)+1)*(bloc_in[k].x)*M_PI)/(16))*cos((((2*x)+1)*(bloc_in[k].y)*M_PI)/(16)));
                                    bloc_out.vert=bloc_out.vert + (bloc_in[x*y].vert*cos((((2*x)+1)*(bloc_in[k].x)*M_PI)/(16))*cos((((2*x)+1)*(bloc_in[k].y)*M_PI)/(16)));
                                    bloc_out.rouge=bloc_out.rouge + (bloc_in[x*y].rouge*cos((((2*x)+1)*(bloc_in[k].x)*M_PI)/(16))*cos((((2*x)+1)*(bloc_in[k].y)*M_PI)/(16)));
 
                    }
    }
    bloc_out.bleu=0.25*C(bloc_in[k].x,bloc_in[k].y)*bloc_out.bleu;;
    bloc_out.vert=0.25*C(bloc_in[k].x,bloc_in[k].y)*bloc_out.vert;
    bloc_out.rouge=0.25*C(bloc_in[k].x,bloc_in[k].y)*bloc_out.rouge;
    return (bloc_out);
}

bloc_in[] est une matrice de 8*8 et la fonction retourne un entier à placer dans un tableau à l'indice k.

Merci d'avance.

EDIT : la formule de la DCT ce trouve là : http://fr.wikipedia.org/wiki/Compression_JPEG

Les structure s_dct et s_imagecontiennent les coordonnées (x,y) d'un point dans une matrice "virtuelle" ainsi que la valeur des composante de couleur du point.

En fait j'envoie dans ma fonction un bloc de 8*8 "pixels" (bloc_in[]) et je voudrai qu'elle me retourne un entier que je place dans un tableau à l'indice k ca pour les trois couleur de la matrice avec derrier une boucle for pour remplire un tableau de la meme structure contenant les résultat de la DCT

**granquet** · 30/05/2006, 19h58

avant qu'on verifie l'implementation de l'algo, ce serait bien que tu nous explique ce que tu attend de ta fonction, et ce qu'il se passe.

Code :

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

    bloc_out.bleu = bloc_out.rouge = bloc_out.vert = 0;

c'est peut etre autorisé, mais je trouve ça moche quand meme.

**miron** · 30/05/2006, 20h07

Envoyé par Dark_Ebola

Code :

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

    bloc_out.bleu = bloc_out.rouge = bloc_out.vert = 0;

Je savais meme pas que c'aitais posible ...

**Jack_serious** · 30/05/2006, 20h12

Envoyé par miron

Je savais meme pas que c'aitais posible ...

Possible mais extremement deconseille, car tres nuisible a la lisibilite du code.

(caches toi je crois que tu vas te faire allumer si tonton emdel passe par la...)

**lemenyau** · 30/05/2006, 20h21

Envoyé par Jack_serious

Possible mais extremement deconseille, car tres nuisible a la lisibilite du code.

Désolé ca fait un mois que je fais du C

**Jack_serious** · 30/05/2006, 20h25

Envoyé par lemenyau

Désolé ca fait un mois que je fais du C

Pas de probleme. C'est juste un conseil, pas une attaque.

**lemenyau** · 30/05/2006, 20h28

Envoyé par Jack_serious

Pas de probleme. C'est juste un conseil, pas une attaque.

Pour moi non plus mais c'est vrai que ma réponse pouvait paraitre un peu sêche.

**miron** · 30/05/2006, 21h36

Envoyé par Jack_serious

Possible mais extremement deconseille, car tres nuisible a la lisibilite du code.

(caches toi je crois que tu vas te faire allumer si tonton emdel passe par la...)

MDR
Il est 21h35 , je vais pas bosser t tla journée!!!

Non mais il est super simpa EMDEL.
Un peu sec par moment, si on lui donne pas ce qu'il fautJustifié souventMais bon aller bonne nuit a tous!!!

**Pouic** · 30/05/2006, 21h38

Hmm, je comprends pas, j'ai du rater un fil aujourd'hui : pourquoi Emmanuel te poursuit-il aujourd'hui au travers de tous les messages à propos de tes révisions ?? Il y a eu des messages spécifiques que j'ai manqués ?

J'aimerai comprendre?? ^__^

**lemenyau** · 30/05/2006, 21h58

Pour les deux posteur précédent c'est gentil de répondre mais ce serait encore mieu si ca avait un rapport avec le sujet.

**granquet** · 30/05/2006, 22h53

petite question: ton tableau bloc_in, tu l'as remplis comment?

et: bloc_in[x*y]
sans vouloir m'avancer trop, ca me semble etre le probleme ... 8*8 ca fait bien 64 ...
sans faire de math, il me semble que ca fait pas 64 valeurs uniques.
comment tu differencie 3*8 et 8*3 ?

moi j'aurai fait bloc_in[8*x + y]

je me trompe peut etre. je dis ça "comme ça", parce que ca me parait "bizzare".

**lemenyau** · 30/05/2006, 23h56

Le remplissage de bloc_in ce fait normalement sans problème avec des boucle for à partir d'un tableau de structure contenant les données concernant les pixels de l'image (coordonées, couleurs). Sinon pour le reste en y réflechissant, tu doit avoir raison. Faut que je test pourtant il me sempblait l'avoir déjà fait.

EDIT:

Envoyé par Dark_Ebola

et: bloc_in[x*y]
sans vouloir m'avancer trop, ca me semble etre le probleme ... 8*8 ca fait bien 64 ...
sans faire de math, il me semble que ca fait pas 64 valeurs uniques.
comment tu differencie 3*8 et 8*3 ?

moi j'aurai fait bloc_in[8*x + y]

Sur ce point je pense que ta as raison, mais ça ne resout pas mon problème. Pour tester j'ai utiliser une image avec une seul couleur donc en theorie je devrais après la DCT obtenir seulement la première valeur non nul (enfin il me semble) alors que là j'obtient dans ma matrice une diagonale de valeur assez importantes et je n'arrive pas à voir d'ou cela vient.

**miron** · 31/05/2006, 10h55

Envoyé par Pouic

Hmm, je comprends pas, j'ai du rater un fil aujourd'hui : pourquoi Emmanuel te poursuit-il aujourd'hui au travers de tous les messages à propos de tes révisions ?? Il y a eu des messages spécifiques que j'ai manqués ?

J'aimerai comprendre?? ^__^

http://www.developpez.net/forums/sho...49&postcount=8
http://www.developpez.net/forums/new...eply&p=1005749
Voila ca peut peut etre vous aider...

**Pouic** · 31/05/2006, 16h12

Aahh, d'accord ^^
Bon, ben retourne réviser alors

Hop hop hop !

**lemenyau** · 31/05/2006, 18h38

ça y 'es j'ai réussi

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
struct s_dct dct(struct s_image bloc_in[],int k)//DCT
{
    int x,y;
    float arg1, arg2;
    struct s_dct bloc_out;
    bloc_out.bleu = bloc_out.rouge = bloc_out.vert = 0;
    for(x=0;x<8;x++)
    {
 
                    for(y=0;y<8;y++)
                    {
                                    arg1 = (2*x+1) * bloc_in[k].x*M_PI/16;
                                    arg2 = (2*y+1) * bloc_in[k].y*M_PI/16;
                                    bloc_out.bleu += bloc_in[x*8+y].bleu*cos(arg1)*cos(arg2);
                                    bloc_out.vert += bloc_in[x*8+y].vert*cos(arg1)*cos(arg2);
                                    bloc_out.rouge += bloc_in[x*8+y].rouge*cos(arg1)*cos(arg2);
                    }
    }
    bloc_out.bleu=0.25*C(bloc_in[k].x,bloc_in[k].y)*bloc_out.bleu;;
    bloc_out.vert=0.25*C(bloc_in[k].x,bloc_in[k].y)*bloc_out.vert;
    bloc_out.rouge=0.25*C(bloc_in[k].x,bloc_in[k].y)*bloc_out.rouge;
    return (bloc_out);
}

Problème résolu après quelques jours en mettant un peu d'ordre d'un le code.

pour arg2 c'était (2*y+1) et pas (2*x +1)

Merci quand même

Code source pour appliquer une DCT

C

Discussions similaires

Partager

Partager