Utilisation de la fonction blockproc

Bonjour,

Je ne parviens pas à utiliser la fonction blockproc pour faire un calcul par blocs 8x8 sur une image 1920x1080.
J'ai une fonction simple appelée vqm2 :

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

function [dy]=vqm2(im)   % Color Tranfform %% r = im(:,:,1); g = im(:,:,2); b = im(:,:,3); [y,u,v]=rgb2yuv(r,g,b);   %% DCT %% dy = dct2(y);   %% Local Contrast %% DC = dy(1,1); LC = (dy * power( DC/1024, 0.65 )) / DC;   %% SCSF %% Q = [8 16 19 22 26 27 29 34; 16 16 22 24 27 29 34 37; 19 22 26 27 29 34 34 38; 22 22 26 27 29 34 37 40; 22 26 27 29 32 35 40 48; 26 27 29 32 35 40 48 58; 26 27 29 34 38 46 56 69; 27 29 35 38 46 56 69 83];   dyq = dy.*Q;

---

J'utilise ensuite la fonction blockproc de la façon suivante :

Code:

---

B1 = blockproc(im1,[8 8],vqm2(im));

---

Matlab me donne l'erreur suivante :

Code:

---

1 2 3

??? Error using ==> BLOCKPROC Argument 'Function' failed validation with error: Invalid block function. BLOCKPROC expects the user function, FUN, to be a valid function handle.

---

Il semble donc que la fonction vqm2 ne fonctionne pas avec blockproc, mais je ne comprend pas pourquoi. J'ai regardé ce qu'est une fonction handle sans plus de résultat...

Merci d'avance pour votre aide

Tu es sur de l'écriture

Code:

---

>> doc blkproc

---

Sinon, essayes ceci :

Code:

---

B1 = blockproc(im1,[8 8],vqm2);

---

En plus blkproc demande la déclaration d'un handle sur une fonction, du genre :

Code:

---

1 2 3 4

I = imread('liftingbody.png'); fun = @(x) std2(x)*ones(size(x)); I2 = blkproc(I,[32 32],fun); imshow(I), figure, imshow(I2,'DisplayRange',[])

---

Merci pour ta réponse,

D'après ce que j'ai compris, blockproc est la nouvelle version de blkproc.

J'ai déjà essayé de taper :

Code:

---

B1 = blockproc(im1,[8 8],vqm2);

---

Le problème est qu'il ne sait ce qu'est la variable im dans la fonction vqm2.

J'ai essayer d'utiliser un handle comme indiqué dans la doc :

Code:

---

1 2 3

fun = @vqm2; B = blockproc(im1,[8 8],fun);

---

Mais dans ce cas j'ai une erreur de dimension de matrice :

Code:

---

1 2 3

Index exceeds matrix dimensions. Error in ==> blockproc at 214 output_block = userfunDispatcher(fun,input_struct,trim_border);

---

Cependant je suis certaine que la fonction vqm2 fonctionne puisqu'il n'y a pas d'erreur quand je l'exécute avec une image 8x8 en entrée...

Je viens de résoudre mon problème, je poste la réponse au cas où ça intéresse quelqu'un...

Il y avait deux erreurs dans mon code:
Premièrement, blockproc n'accepte pas de matrices 3D (en tout cas je n'ai pas réussi à le faire fonctionner), dons je fait le changement d'espace colorimétrique avant d'utiliser blockproc et je lui donne en entrée seulement la luminance.

Deuxièmement, la syntaxe de la fonction handle (fun) était incomplète (cf ci-dessous)

Code:

---

1 2 3 4 5 6 7 8

r = im(:,:,1); g = im(:,:,2); b = im(:,:,3); [y,u,v]=rgb2yuv(r,g,b);   fun = @ (block_struct)vqm2(block_struct.data); B1 = blockproc(y,[8 8],fun);

---

et la fonction utilisée par blockproc :

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14

function [dyq]=vqm2(y) dy = dct2(y); DC = dy(1,1); LC = (dy * power( DC/1024, 0.65 )) / DC; Q = [8 16 19 22 26 27 29 34; 16 16 22 24 27 29 34 37; 19 22 26 27 29 34 34 38; 22 22 26 27 29 34 37 40; 22 26 27 29 32 35 40 48; 26 27 29 32 35 40 48 58; 26 27 29 34 38 46 56 69; 27 29 35 38 46 56 69 83]; dyq = dy.*Q;

---

Voila! en espérant que ça puisse être utile