Discussion :

[cindy63]

Bonjour,
Je souhaiterais utiliser les moments géométriques dans une méthode de segmentation. J'ai trouvé l'article en pièce jointe pour effectuer le calcul mais je n'arrive pas à retrouver les mêmes résultats qu'eux (page 5 et 6 du pdf).
Voici l'algorithme que j'ai utilisé:

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
 
pour k = -n/2 à n/2 faire
  pour j = -n/2 à n/2 faire
    pour i = -n/2 à n/2 faire
      coeff = 0
      pour k’ = k–m/2 à k+m/2 faire
        pour j’ = j–m/2 à j+m/2 faire
          pour i’ = i–m/2 à i+m/2 faire
            coeff += (i’)^p + (j’)^q + (k’)^r
          finpour
        finpour
      finpour
      C(i+n/2,j+n/2,k+n/2) = coeff 
    finpour
  finpour
finpour

Pensez-vous que ce soit correcte? Où est l'erreur?
Merci de vos réponses.

[ToTo13]

Bonjour,

merci d'utiliser les balises et d'indenter votre code afin de le rendre lisible.
Cela nous permettra de lire et donc de comprendre plus facilement votre algo.

[pseudocode]

Pensez-vous que ce soit correcte? Où est l'erreur?
Merci de vos réponses.

C'est la formule (20) ou la (19b) que tu veux calculer ?

Parceque dans la (19b) il y a un facteur multiplicatif devant la somme.

[cindy63]

La (20).

[pseudocode]

Each voxel in a nxnxn cubic window is subdivided into mxmxm subvoxels

(i',j',k') are integers on a different scale from (i,j,k).

Ce que je comprend c'est que chaque voxel est découpé en m^3 subvoxels, et que 2 voxels voisins n'ont AUCUN subvoxels commun. Donc je verrai plus un code comme ca:

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
 
FOR k = -(n*m)/2 TO (n*m)/2 STEP m DO
  FOR j = -(n*m)/2 TO (n*m)/2 STEP m DO
    FOR i = -(n*m)/2 TO (n*m)/2 STEP m DO
 
    coeff = 0
    FOR k'= k-m/2 TO k+m/2 STEP 1 DO
      FOR j'= j-m/2 TO j+m/2 STEP 1 DO
        FOR i'= i-m/2 TO i+m/2 STEP 1 DO
          coeff += (i’)^p + (j’)^q + (k’)^r
        NEXT i'
      NEXT j'
    NEXT k'
 
    C(i+n/2,j+n/2,k+n/2) = coeff
 
    NEXT i
  NEXT j
NEXT k

[cindy63]

J'ai essayé ton code mais je n'ai toujours pas les bons résultats.

[pseudocode]

J'ai essayé ton code mais je n'ai toujours pas les bons résultats.

Bon, bah je ne sais pas.

Quelles sont les valeurs que tu trouves et quelles sont celles que tu devraient trouver ?