extraire une sous matrice d'une matrice sous condition

**jeuneprincesse** · 30/04/2016, 23h46

Bonsoir à tous,
j'ai une matrice 2000x2000 contenant des 0 et des 1 je veux faire un algorithme qui me permet de vérifier l'existence d'un bloc 50x50 de 0 et de l'extraire (cà dire extraire l'intervalle de ligne et de colonne)
exemple:
si le bloc 1:50x1:50 contient des zeros seulement pas des 1 le programme me donne oui le bloc est 1:50x1:50

Merci infiniment d'avance pour toute réponse

**Jerome Briot** · 01/05/2016, 21h37

Utilise la fonction any pour tester la présence de 1 dans les sous-matrices.

**jeuneprincesse** · 01/05/2016, 23h15

Envoyé par Jerome Briot

Utilise la fonction any pour tester la présence de 1 dans les sous-matrices.

merci pour la réponse en fait je veux un algo qui me fait sortir les bloc contenant de 0 càd s'il existe une sous matrice de dimension 50x50 qui ne contient que des 0 alors l'algo m'affiche il existe le bloc (1:50,1:50) qui contient des 0 purs et si un autre bloc contient des 0 il me le donne aussi j'espère que vous avez compris ce que je viens de dire
merci beaucoup encore une fois pour tous ceux qui puissent m'aider

**vbbarent** · 02/05/2016, 15h28

Bonjour,

C'est à toi d'écrire une telle fonction (c'est pas si compliquer), Il suffit de découper ta matrice en blocs (une boucle for) et de tester chaque bloc, sachant tu connais les indices des blocs que tu traites.

@++

**jeuneprincesse** · 02/05/2016, 23h23

Envoyé par vbbarent

Bonjour,

C'est à toi d'écrire une telle fonction (c'est pas si compliquer), Il suffit de découper ta matrice en blocs (une boucle for) et de tester chaque bloc, sachant tu connais les indices des blocs que tu traites.

@++

en fait je suis débutante en matlab et pour les blocs c l'algo qui va les sortir en fait le bloc de 50 peut être 30:79 x 50:99
un peu plus d'éclaircissement SVP si c possible
merci d'avance

**vbbarent** · 03/05/2016, 21h06

Envoyé par jeuneprincesse

en fait je suis débutante en matlab et pour les blocs c l'algo qui va les sortir en fait le bloc de 50 peut être 30:79 x 50:99
un peu plus d'éclaircissement SVP si c possible
merci d'avance

Essaie ceci

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
 
idxbloc=1;
for l=1:1950 % on parcourt toutes les lignes
    idl = l:l+50; % on prend les indices des lignes par groupe de 50 en decalant d'une unité
        for c=1:1950% on parcourt toutes les colonnes
        idc = c:c+50; % idem que pour les lignes
        if isempty(find(im(idl,idc))~=0) % on teste le bloc pour voir si y a des 1
          bloc{dxbloc}= [idl;idc]; % on garde les les indices correspondants dans une structure
          idxbloc = idxbloc+1; % on incremente la taille de la structure
        end
    end
end

**jeuneprincesse** · 03/05/2016, 23h34

merci pour la réponse
mais il me donne ce message d'erreur

Undefined function or variable 'dxbloc'.

que puisse je faire s'il vous plait ?

**vbbarent** · 04/05/2016, 09h30

Envoyé par jeuneprincesse

merci pour la réponse
mais il me donne ce message d'erreur

Undefined function or variable 'dxbloc'.

que puisse je faire s'il vous plait ?

C'est une erreur de frappe; Désolé! c'est plutôt

Code :

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

1
2
 
idxbloc

**jeuneprincesse** · 05/05/2016, 01h09

Envoyé par vbbarent

C'est une erreur de frappe; Désolé! c'est plutôt

Code :

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

1
2
 
idxbloc

merci
et voilà le résultat que j'ai obtiens en fait le prog ne montre pas tous les blocs d'une et d'autre par le bloc qu'il a montré dépasse 1950 en ligne et colonne ?

[idl;idc]

ans =

Columns 1 through 10

1950 1951 1952 1953 1954 1955 1956 1957 1958 1959
1950 1951 1952 1953 1954 1955 1956 1957 1958 1959

Columns 11 through 20

1960 1961 1962 1963 1964 1965 1966 1967 1968 1969
1960 1961 1962 1963 1964 1965 1966 1967 1968 1969

Columns 21 through 30

1970 1971 1972 1973 1974 1975 1976 1977 1978 1979
1970 1971 1972 1973 1974 1975 1976 1977 1978 1979

Columns 31 through 40

1980 1981 1982 1983 1984 1985 1986 1987 1988 1989
1980 1981 1982 1983 1984 1985 1986 1987 1988 1989

Columns 41 through 50

1990 1991 1992 1993 1994 1995 1996 1997 1998 1999
1990 1991 1992 1993 1994 1995 1996 1997 1998 1999

Column 51

2000
2000

**phryte** · 05/05/2016, 10h06

Bonjour,

Un essai pour répondre à ta première question sur la détection matrices zéros (limité de 1x1 à 4x4) :

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
23
24
clear
M =[0     1     0     0     0     0     0     0     1     0
    1     0     0     0     0     0     1     0     1     0
    0     1     0     0     0     0     1     0     1     0
    1     1     0     0     0     0     0     1     1     1
    0     0     1     0     1     1     1     0     1     0
    0     0     1     1     1     1     0     0     0     0
    1     1     0     0     1     0     0     0     0     0
    1     0     0     1     0     1     0     0     0     0
    1     1     0     1     0     1     0     0     0     0
    1     1     0     1     1     1     1     0     1     0];
MM=num2str(M);% Conversion en caractères
dx=4;dy=4;% Dimensions de la matrice de zéros recherchée
for n=1:length(M)
    r=strfind(MM(n,:), num2str(zeros(dx,1)'));% Recherche des chaînes de caractères candidates
    T(n,1:length(r))=r;% Rang des colonnes avec chaîne recherchée
end
resu=inters(M,T,dx,dy);% Fonction de calcul des zones
disp(['Il y a : '  num2str(nnz(resu)) ' zone(s) de zéros de cotés : ' num2str(dx) 'x' num2str(dy)])
[I,J,V]=find(resu);% calcul des indices
for k=1:length(I)
    disp(['Matrice : ' num2str(k)]);
    disp(MM(I(k):I(k)+dy-1,V(k):V(k)+2*dx+1))
end

La fonction (je cherche à la réduire avec la récursivité) :

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
23
function res=inters(M,T,dx,dy)
res=zeros(size(T));
switch dy
    case 1
        res=T;
    case 2
        for n=1:length(T)-dy
            v=intersect(T(n,:),T(n+1,:));
            res(n,1:length(v))=v;
        end
    case 3
 
        for n=1:length(T)-dy
            v=intersect(T(n,:),intersect(T(n+1,:),T(n+2,:)));
            res(n,1:length(v))=v;
        end
    case 4
        for n=1:length(T)-dy
            v=intersect(T(n,:),intersect(T(n+1,:),intersect(T(n+2,:),T(n+3,:))));
            res(n,1:length(v))=v;
        end
end
end

Dans cet exemple :
Il y a : 2 zone(s) de zéros de cotés : 4x4
Matrice : 1
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
Matrice : 2
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0

**jeuneprincesse** · 05/05/2016, 11h14

Envoyé par phryte

Bonjour,

Un essai pour répondre à ta première question sur la détection matrices zéros (limité de 1x1 à 4x4) :

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
23
24
clear
M =[0     1     0     0     0     0     0     0     1     0
    1     0     0     0     0     0     1     0     1     0
    0     1     0     0     0     0     1     0     1     0
    1     1     0     0     0     0     0     1     1     1
    0     0     1     0     1     1     1     0     1     0
    0     0     1     1     1     1     0     0     0     0
    1     1     0     0     1     0     0     0     0     0
    1     0     0     1     0     1     0     0     0     0
    1     1     0     1     0     1     0     0     0     0
    1     1     0     1     1     1     1     0     1     0];
MM=num2str(M);% Conversion en caractères
dx=4;dy=4;% Dimensions de la matrice de zéros recherchée
for n=1:length(M)
    r=strfind(MM(n,:), num2str(zeros(dx,1)'));% Recherche des chaînes de caractères candidates
    T(n,1:length(r))=r;% Rang des colonnes avec chaîne recherchée
end
resu=inters(M,T,dx,dy);% Fonction de calcul des zones
disp(['Il y a : '  num2str(nnz(resu)) ' zone(s) de zéros de cotés : ' num2str(dx) 'x' num2str(dy)])
[I,J,V]=find(resu);% calcul des indices
for k=1:length(I)
    disp(['Matrice : ' num2str(k)]);
    disp(MM(I(k):I(k)+dy-1,V(k):V(k)+2*dx+1))
end

La fonction (je cherche à la réduire avec la récursivité) :

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
23
function res=inters(M,T,dx,dy)
res=zeros(size(T));
switch dy
    case 1
        res=T;
    case 2
        for n=1:length(T)-dy
            v=intersect(T(n,:),T(n+1,:));
            res(n,1:length(v))=v;
        end
    case 3
 
        for n=1:length(T)-dy
            v=intersect(T(n,:),intersect(T(n+1,:),T(n+2,:)));
            res(n,1:length(v))=v;
        end
    case 4
        for n=1:length(T)-dy
            v=intersect(T(n,:),intersect(T(n+1,:),intersect(T(n+2,:),T(n+3,:))));
            res(n,1:length(v))=v;
        end
end
end

Dans cet exemple :
Il y a : 2 zone(s) de zéros de cotés : 4x4
Matrice : 1
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
Matrice : 2
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0

Bonjour et merci pour la réponse j'ai essayé le code et il m'a donné ce message d'erreur
Undefined function or variable 'inters'. ?

**phryte** · 05/05/2016, 11h30

Bonjour
Il faut mettre les deux programmes dans le même dossier.

**jeuneprincesse** · 05/05/2016, 18h33

Envoyé par phryte

Bonjour
Il faut mettre les deux programmes dans le même dossier.

quel dossier ? j'ai copié les deux programmes et les coller dans command window

**phryte** · 05/05/2016, 18h38

Bonjour,

Tu copies chaque programme dans un new script que tu sauvegardes dans un même dossier.

**jeuneprincesse** · 05/05/2016, 19h19

Envoyé par phryte

Bonjour,

Tu copies chaque programme dans un new script que tu sauvegardes dans un même dossier.

bonjour et merci beaucoup maitenant le programme s'exécute mais comme même il y a une autre erreur que puisse je faire ?

**phryte** · 06/05/2016, 07h55

Bonjour,

"inters" est un sous programme "function", en général il ne s'exécute pas, il est appelé par d'autres programmes avec les paramètres initialisés.
l'autre script est le programme principal (clear; M=.....). C'est lui qu'il faut exécuter.

Ce programme n'a pas de limitations en dx par contre dy est limité à 4. (il faut améliorer la fonction "inters" par de la récursivité.

**jeuneprincesse** · 06/05/2016, 11h56

Envoyé par phryte

Bonjour,

"inters" est un sous programme "function", en général il ne s'exécute pas, il est appelé par d'autres programmes avec les paramètres initialisés.
l'autre script est le programme principal (clear; M=.....). C'est lui qu'il faut exécuter.

Ce programme n'a pas de limitations en dx par contre dy est limité à 4. (il faut améliorer la fonction "inters" par de la récursivité.

Bonjour,
donc j'ajoute d'autre cases jusqu'à 50 ?
disp(MM(I(k):I(k)+dy-1,V(k):V(k)+2*dx+1)) cette commande ça sert à quoi ?

je veux aussi que le message suivant s'affiche c à dire insérer les intervalles des matrice

Code :

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

1
2
3
4
5
6
7
8
9
10
11
Il y a : 2 zone(s) de zéros de cotés : 4x4
Matrice : 1 (1:3,3:6)
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
Matrice : 2(5:8,7:10)
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0

merci infiniment d'avance pour toute réponse

**jeuneprincesse** · 06/05/2016, 19h32

en plus quand j'ai testé le programme sur ma matrice qui est de type logique
du coup j'ai fait la modification suivante au code:

Code :

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

1
2
3
M =matrice;
s=double(M);
MM=num2str(s)

et ça n'a pas marché j'ai eu ce message d'erreur

Code :

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

1
2
3
4
5
6
7
Out of memory. Type HELP MEMORY for your options.
 
Error in num2str (line 66)
        xmax = double(max(abs(widthCopy( : ))));
 
Error in fonction1 (line 3)
MM=num2str(s);% Conversion en caractères

:'(
que puisse je faire s'il vous plait ?

**jeuneprincesse** · 07/05/2016, 00h10

[QUOTE=phryte;8622181]Bonjour,

Un essai pour répondre à ta première question sur la détection matrices zéros (limité de 1x1 à 4x4) :

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
23
24
clear
M =[0     1     0     0     0     0     0     0     1     0
    1     0     0     0     0     0     1     0     1     0
    0     1     0     0     0     0     1     0     1     0
    1     1     0     0     0     0     0     1     1     1
    0     0     1     0     1     1     1     0     1     0
    0     0     1     1     1     1     0     0     0     0
    1     1     0     0     1     0     0     0     0     0
    1     0     0     1     0     1     0     0     0     0
    1     1     0     1     0     1     0     0     0     0
    1     1     0     1     1     1     1     0     1     0];
MM=num2str(M);% Conversion en caractères
dx=4;dy=4;% Dimensions de la matrice de zéros recherchée
for n=1:length(M)
    r=strfind(MM(n,:), num2str(zeros(dx,1)'));% Recherche des chaînes de caractères candidates
    T(n,1:length(r))=r;% Rang des colonnes avec chaîne recherchée
end
resu=inters(M,T,dx,dy);% Fonction de calcul des zones
disp(['Il y a : '  num2str(nnz(resu)) ' zone(s) de zéros de cotés : ' num2str(dx) 'x' num2str(dy)])
[I,J,V]=find(resu);% calcul des indices
for k=1:length(I)
    disp(['Matrice : ' num2str(k)]);
    disp(MM(I(k):I(k)+dy-1,V(k):V(k)+2*dx+1))
end

Bonsoir,

C'est quoi nnz(resu) c'est une fonction de matlab ?
j'ai étalé le programme à 14 bien sûr pour une matrice plus grande ça marche mais il m'affiche une matrice 10x14 et non 14 x 14 le problème est dû à dx ?
voilà le code 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
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
function res=inters(M,T,dx,dy)
res=zeros(size(T));
switch dy
    case 1
        res=T;
    case 2
        for n=1:length(T)-dy
            v=intersect(T(n,:),T(n+1,:));
            res(n,1:length(v))=v;
        end
    case 3
 
        for n=1:length(T)-dy
            v=intersect(T(n,:),intersect(T(n+1,:),T(n+2,:)));
            res(n,1:length(v))=v;
        end
    case 4
        for n=1:length(T)-dy
            v=intersect(T(n,:),intersect(T(n+1,:),intersect(T(n+2,:),T(n+3,:))));
            res(n,1:length(v))=v;
        end
        case 5
        for n=1:length(T)-dy
            v=intersect(T(n,:),intersect(T(n+1,:),intersect(T(n+2,:),intersect(T(n+3,:),T(n+4,:)))));
            res(n,1:length(v))=v;
        end
         case 6
        for n=1:length(T)-dy
            v=intersect(T(n,:),intersect(T(n+1,:),intersect(T(n+2,:),intersect(T(n+3,:),intersect(T(n+4,:),T(n+5,:))))));
            res(n,1:length(v))=v;
        end
          case 7
        for n=1:length(T)-dy
            v=intersect(T(n,:),intersect(T(n+1,:),intersect(T(n+2,:),intersect(T(n+3,:),intersect(T(n+4,:),intersect(T(n+5,:),T(n+6,:)))))));
            res(n,1:length(v))=v;
        end
         case 8
        for n=1:length(T)-dy
            v=intersect(T(n,:),intersect(T(n+1,:),intersect(T(n+2,:),intersect(T(n+3,:),intersect(T(n+4,:),intersect(T(n+5,:),intersect(T(n+6,:),T(n+7,:))))))));
            res(n,1:length(v))=v;
        end
        case 9
        for n=1:length(T)-dy
            v=intersect(T(n,:),intersect(T(n+1,:),intersect(T(n+2,:),intersect(T(n+3,:),intersect(T(n+4,:),intersect(T(n+5,:),intersect(T(n+6,:),intersect(T(n+7,:),T(n+8,:)))))))));
            res(n,1:length(v))=v;
        end
        case 10
        for n=1:length(T)-dy
            v=intersect(T(n,:),intersect(T(n+1,:),intersect(T(n+2,:),intersect(T(n+3,:),intersect(T(n+4,:),intersect(T(n+5,:),intersect(T(n+6,:),intersect(T(n+7,:),intersect(T(n+8,:),T(n+9,:))))))))));
            res(n,1:length(v))=v;
        end
         case 11
        for n=1:length(T)-dy
            v=intersect(T(n,:),intersect(T(n+1,:),intersect(T(n+2,:),intersect(T(n+3,:),intersect(T(n+4,:),intersect(T(n+5,:),intersect(T(n+6,:),intersect(T(n+7,:),intersect(T(n+8,:),intersect(T(n+9,:),T(n+10,:)))))))))));
            res(n,1:length(v))=v;
        end
 
         case 12
        for n=1:length(T)-dy
            v=intersect(T(n,:),intersect(T(n+1,:),intersect(T(n+2,:),intersect(T(n+3,:),intersect(T(n+4,:),intersect(T(n+5,:),intersect(T(n+6,:),intersect(T(n+7,:),intersect(T(n+8,:),intersect(T(n+9,:),intersect(T(n+10,:),T(n+11,:))))))))))));
            res(n,1:length(v))=v;
        end
        case 13
        for n=1:length(T)-dy
            v=intersect(T(n,:),intersect(T(n+1,:),intersect(T(n+2,:),intersect(T(n+3,:),intersect(T(n+4,:),intersect(T(n+5,:),intersect(T(n+6,:),intersect(T(n+7,:),intersect(T(n+8,:),intersect(T(n+9,:),intersect(T(n+10,:),intersect(T(n+11,:),T(n+12,:)))))))))))));
            res(n,1:length(v))=v;
        end
        case 14
        for n=1:length(T)-dy
            v=intersect(T(n,:),intersect(T(n+1,:),intersect(T(n+2,:),intersect(T(n+3,:),intersect(T(n+4,:),intersect(T(n+5,:),intersect(T(n+6,:),intersect(T(n+7,:),intersect(T(n+8,:),intersect(T(n+9,:),intersect(T(n+10,:),intersect(T(n+11,:),intersect(T(n+12,:),T(n+13,:))))))))))))));
            res(n,1:length(v))=v;
        end
end
end

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
M=mat;
MM=num2str(M);% Conversion en caractères
dx=14;dy=14;% Dimensions de la matrice de zéros recherchée
for n=1:length(M)
    r=strfind(MM(n,:), num2str(zeros(dx,1)'));% Recherche des chaînes de caractères candidates
    T(n,1:length(r))=r;% Rang des colonnes avec chaîne recherchée
end
resu=interstest(M,T,dx,dy);% Fonction de calcul des zones
disp(['Il y a : '  num2str(nnz(resu)) ' zone(s) de zéros de cotés : ' num2str(dx) 'x' num2str(dy)])
[I,J,V]=find(resu);% calcul des indices
for k=1:length(I)
    disp(['Matrice : ' num2str(k)]);
    disp(MM(I(k):I(k)+dy-1,V(k):V(k)+2*dx+1))
end

et encore il y a des zones qui présentent des 1 ?
merci d'avance pour toute aide

**jeuneprincesse** · 07/05/2016, 00h16

Envoyé par jeuneprincesse

en plus quand j'ai testé le programme sur ma matrice qui est de type logique
du coup j'ai fait la modification suivante au code:

Code :

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

1
2
3
M =matrice;
s=double(M);
MM=num2str(s)

et ça n'a pas marché j'ai eu ce message d'erreur

Code :

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

1
2
3
4
5
6
7
Out of memory. Type HELP MEMORY for your options.
 
Error in num2str (line 66)
        xmax = double(max(abs(widthCopy( : ))));
 
Error in fonction1 (line 3)
MM=num2str(s);% Conversion en caractères

:'(
que puisse je faire s'il vous plait ?

j'ai fait un clear et j'ai travaillé sur une partie de la matrice sans passer du logique au numérique j'ai appliqué le programme il m'affiche des résultats mais d'une part il ya des zones qui présentent des 1 et la matrice n'est pas carré càd si je lui demande matrice 14x14 il me donne une 10x14
merci d'avance pour toute aide

extraire une sous matrice d'une matrice sous condition

MATLAB

Discussions similaires

Partager

Partager