Discussion :

[cedji]

Bonjour ,

Je voudrais superposer deux images pour les comparer. Le problème est que l'une des images est codée en "int16" et l'autre en "int32" donc évidemment, lorsque je fais la commande

Code :

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

1
2
3
4
5
6
7
8
imsup=zeros(height_temp,width_temp,NbSlices_temp);
for i=1:height_temp
    for j=1:width_temp
        for k=1:NbSlices_temp
    im_sup(i,j,k)=im_temp(i,j,k)+im2(i,j,k);
        end
    end
end

je me fais insultée par Matlab...

J'ai essayé les deux solutions
Tout d'abord, j'ai converti l'image im2 en 32bits (en utilisant la fonction int32) mais les valeurs des pixels sont alors trop faibles par rapport a ceux de l'image im_temp. Donc la somme des deux images donnent l'image im_temp
J'ai ensuite converti l'image im_temp en 16 bits (en utilisant la fonction int16) mais les valeurs des pixels sont alors toutes égales à 2^15..

Est-ce que qqn aurait une solution pour m'aider???
merci!!

[Jean Dumoncel]

Bonjour,

Tout d'abord une remarque : le format int16 représente des nombres compris entre -2^15 et 2^15-1 et le format int32 représente des nombres compris entre -2^31 et 2^31-1. Voila pourquoi si tes int32 sont tous supérieurs à 2^15-1, la conversion en int16 les convertit à 2^15-1 :

Code :

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

1
2
3
4
5
>> int16(int32(32768))
 
ans =
 
  32767

int16 et int32 change le format de tes données, toi tu veux changer le format et les valeurs, regarde cette fonction :

Code :

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

doc typecast

Une autre remarque : pourquoi veux-tu faire la somme de tes images pour les comparer? Quel est le but?
Car ta prochaine problématique sera certainement le dépassement à cause de ces additions :

Code :

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

1
2
3
4
5
>> int16(16700)+int16(26700)
 
ans =
 
  32767

par exemple

[cedji]

Merci pour ta réponse Magelan.
Je vais mieux expliquer ce que je fais. En fait je dispose d'images acquise par IRM d'une part, et d'un template déjà recalé sur un atlas d'autre part. Ce que je veux faire est de recaler les images acquises de l'IRM sur les images du template. Pour cela j'ai d'abord fait une normalisation en utilisant le logiciel SPM; Maintenant je veux vérifier la qualité de la normalisation. C'est pourquoi je veux faire la superposition des images (ou la difference) pour voir les differences restantes et donc voir si je dois améliorer le recalage fait avec SPM.
Depuis tout à l'heure j'ai trouvé une "pseudo solution" à mon problème. J'ai converti les images acquises en 32 bits et j'ai changé la dynamique en multipliant les valeurs de chaque pixel par (2^31)/(2^15). AInsi les valeurs des deux images sont comparables mais le contraste différe encore trop pour pouvoir vraiment faire des comparaisons