Discussion :

[Maxnof]

Bonjour à tous.

Je cherche à n'afficher qu'une partie d'une image sur une image de fond (de même taille).
J'explique le problème en détail:
J'ai une succession d'images sur lesquelles je fait un filtrage par SVD pour en faire ressortie un évènement particulier.
Pour les curieux, ce sont des images d'échographie ultrasonore et j'utilise la SVD pour faire ressortir les zones de cavitation.
Après le filtrage, j'applique un seuil arbitraire sur la matrice de cavitation, pour chaque image. J'ai des donc des images avec que du 0 sauf sur les zones avec de la cavitation, où là il y a une valeur.
Je cherche à superposer l'image l'initiale (en niveau de gris) avec les zones qui dépassent le seuil sur mon image de cavitation (en couleur, si possible en transparence).

A priori ça ne me semble pas impossible à faire mais je bloque...

Merci pour votre aide !

[Maxnof]

Petite mise à jour:
- J'ai appliqué un mask pour rendre transparentes les parties inférieures au seuil, et j'applique une transparence de 0.8 aux parties au dessus du seuil.
- En revanche, je n'arrive pas à afficher le background en N&B et le foreground en couleurs...

Voilà mon code pour l'affichage:

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
seuil=valMax/10;
imagesc(piled(:,:,ii))  %BG
    colormap gray
    hold on
 
    %transparency mask:
    for xx=1:nx
        for zz=1:nz
            if s_b(xx,zz,ii)<seuil
                mask(xx,zz)=0;
            else 
                mask(xx,zz)=0.8;
            end
        end
    end
 
    im=imagesc(s_b(:,:,ii),CLIM);
 
    set(im,'AlphaData',mask);
    colormap jet

[Jerome Briot]

Il ne peut y avoir qu'une seule palette de couleurs (colormap) associée à un objet Axes (et donc à une image).

Dans ton cas, il faut convertir l'une des deux images en couleurs vraies (type RGB).

Par exemple :

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
seuil=valMax/10;
 
rgb = ind2rgb(piled(:,:,ii), gray);
image(rgb)
 
hold on
 
%transparency mask:
for xx=1:nx
    for zz=1:nz
        if s_b(xx,zz,ii)<seuil
            mask(xx,zz)=0;
        else
            mask(xx,zz)=0.8;
        end
    end
end
 
im=imagesc(s_b(:,:,ii),CLIM);
 
set(im,'AlphaData',mask);
colormap jet

[Maxnof]

C'est super !

Avec un peu de retard, un grand merci !