Masquer plusieurs polygones (ROI) simultanement [Résolu]

[KrokroAus]

Bonjour à tous,

Je souhaiterai masquer simultanément plusieurs polygones (ROI) dont j'ai les coordonnées géographiques sur une image georéférencée. En gros, mon image est une image satellite d'une partie de côte avec un peu de terre, un peu de mer et des iles. Je souhaiterai garder seulement les valeurs satellites correspondant au domaine marin et enlever les valeurs correspondants aux iles et la terre.

Normalement, j'utilise la fonction inpolygon pour supprimer le domaine terrestre:

Code :

1
2
3
4
5
6
 
Mask=load([rep '../Coastline.txt']);   
xv=Mask(:,2);
yv=Mask(:,1);
IN = inpolygon(Xp,Yp,xv,yv);  
ValeurSatellites(IN)=0;

Avec le fichier coastline.txt qui contient les coordonnées du trait de côte:

Code :

1
2
3
4
5
6
-18.247528	146.224556
-18.255861	146.243750
-18.230028	146.213750
-18.193722	146.219167
         ....
-18.247528	146.224556

Je voudrais supprimer également en même temps les iles aux alentours, c'est à dire lire un fichier txt unique contenant les coordonnées du trait de cote + les coordonnées du trait de cote des différentes iles.
J'ai essayé un fichier txt du genre:

Code :

1
2
3
4
5
6
7
8
9
-18.247528	146.224556
-18.255861	146.243750
-18.230028	146.213750
-18.193722	146.219167
-18.247528	146.224556
nan	nan
-18.750833	146.566250
-18.749583	146.567500
-18.750833	146.566250

Mais cela ne marche pas: les masques se superposent. D’où les questions suivantes:
- peut on utiliser la fonction inpolygone pour masquer plusieurs polygones dont les coordonnées géographiques seraient indiquées sur le même fichier txt?
- Ou bien faut-il que je créer un fichier txt pour chacune des iles....(sachant qu'il y en a plus de 100...)
- ou existe il une autre fonction matlab plus adaptée?

Merci d'avance pour votre aide! :-)

[Dut]

Dans le cas où tu insérerais des NaN entre chaque ROI dans ton fichier texte, essaie ceci :

Code :

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
idx = find(isnan(Mask(:,1)));

xv = Mask(1:idx(1)-1,2);
yv = Mask(1:idx(1)-1,1);

IN = inpolygon(Xp,Yp,xv,yv);

for n = 1:numel(idx)-1
    
    xv = Mask(idx(n)+1:idx(n+1)-1,2);
    yv = Mask(idx(n)+1:idx(n+1)-1,1);

    IN = IN | inpolygon(Xp,Yp,xv,yv);
    
end

xv = Mask(idx(end)+1:end,2);
yv = Mask(idx(end)+1:end,1);

IN = IN | inpolygon(Xp,Yp,xv,yv);

ValeurSatellites(IN) = 0;

[KrokroAus]

C'est parfait, ça marche très bien! Merci beaucoup pour l'aide.