Discussion :

[totoc1001]

Bonsoir a tous,
J'aurais besoin d'un tout petit peu d'aide svp.
J'ai 2 boucles For que j'aimerai simplifier, car elle diminuent grandement le temps d'exection de mon soft.

Les voici:
La première:

Code :

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

1
2
3
4
5
6
 
        for x=1:100
            for y=1:100
                RotatedTest(y,x)=0;
            end
        end

Ceci ne correspondant qu'a une seule partie de l'image, je chercher a mettre a 0 les 100*100 premiers pixels de mon image, je ne sais pas comment utiliser la fonction.

La seconde je pense qu'on doit pouvoir untiliser la fonction find non? Mais comment assigner la contition?

Code :

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

1
2
3
4
5
6
7
8
9
10
 
       [Y,X]=size(MergedArray);
       %find color grayscale greater than 245 
         for y=1:Y
             for x=1:X
                 if(MergedArray(y,x)>245)
                    MergedArray(y,x)=0;
                 end
             end
         end

Merci beaucoup
Thomas

[Caro-Line]

Première boucle, si RotatedTest n'a pas été défini avant :

Code :

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

1
2
 
RotatedTest = zeros(100,100);

Sympathique, non ?

La 2ème :

Code :

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

MergedArray(MergedArray>245) = 0

Je n'ai pas Matlab sous la main pour vérifier, j'ai quand même une crainte c'est que cette méthode risque de te redimensionner ta matrice en vecteur...
Après il faut utiliser reshape si c'est le cas.

[totoc1001]

Première boucle, si RotatedTest n'a pas été défini avant :

Code :

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

1
2
 
RotatedTest = zeros(100,100);

Sympathique, non ?

Merci beaucoup pour la réponse, mais en fait, ca ca me redimentionne toute mon image comme étant une image de 100*100 avec tous les caracteres a 0.
ALors que ce que j'essaye de faire c'est ne mettre a 0 qu'une partie de mon image.
En effet, ceci correspond uniquement au coint supérieur gauche de mon image qui elle est 10fois plus grande. Et je ne veux pas toucher au reste de l'image.

Merci pour votre aide;
Thomas

[rostomus]

Salut,

Par indexage:

Code :

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

RotatedTest(1:100,1:100)=0;

[totoc1001]

La manière suivante fonctionne:

Code :

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

1
2
3
x = 1:10;
y = 1:10;
RotatedTest(y,x)=0;

Merci Rostomus, je viens juste de voir ta réponse, je vais tester.

Bye
Thomas

[totoc1001]

Désolé, j'ai vraiment du mal avec le passage par indexation pour supprimer les boucles for.
J'ai encore un cas dans lequel j'ai 1 image, que je veux redimensionner uniquement sur la zone contenant des données positives dedans.
Avec 2 boucles For, ca donne:

Code :

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

1
2
3
4
5
6
7
8
9
10
 
    [yMax,xMax]=size(MergedArray);
    for y=1:yMax
        for x=1:xMax
            if(MergedArray(y,x)>0)
                NewMergeArray(y,x)=MergedArray(y,x);
            end
        end
    end
    MergedArray=NewMergeArray;

Mais une fois de plus l'image étant très grosse, ca prend énormément de temps.
Est ce que quelqu'un pourrait m'aider a résoudre la condition if a l'intérieur des 2 boucles, et par la même occasion est-il possible de supprimer aussi les boucles?

Merci beaucoup
Thomas

[Jerome Briot]

Je ne suis pas sur d'avoir tout suivi mais :

Code :

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

1
2
M=rand(5)
M(1:2,1:2)=M(1:2,1:2).*(M(1:2,1:2)>.5)

Est-ce que c'est ce que tu cherches à faire ? Ce qui se traduirait dans ton problème :

Code :

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

NewMergeArray(1:yMax,1:xMax)=MergeArray(1:yMax,1:xMax).*(MergeArray(1:yMax,1:xMax)>0)

[totoc1001]

En fait, j'ai un peu de mal a comprendre ce que tu fais avec cette opération.
J'ai essayé de mettre ton code tel kel, mais je n'obtiens pas le bon résultat.
En fait j'essaye de redimensionner mon image en ne copiant dans une nouvelle image, quel la partie ou les valeur a l'intérieur de mon image sont non nulle, (on enlève les parties noires sur les bords de mon image).
Ton code ne redimensionne pas mon image, il contiens toujours les bandes noires sur les bords.

Merci pour ton aide

[rostomus]

Salut,

c'est un peu ambigu je pense, un exemple d'une matrice:

Code :

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

1
2
3
4
5
6
7
 
x=
 
1 2 3 4 0 3
3 1 3 0 0 0
3 2 1 3 0 2
0 0 2 3 4 5

comment tu veux faire?

[totoc1001]

Voici l'image que j'essaye de modifier:



LA partie en Rouge est normalement costituée de 0, et c'est la partie que j'essaye d'enlever.
Avec les boucles for, nous obtenons l'image sans la partie Rouge, ce qui est ce que je veux.
Pour cela, je test simplement quand la valeur des pixels de mon image est positive, je recopie ces pixels ds ma nouvelle image. La partie noire a gauche n'est pas importante.
Mais avec les boucles, cela demande beaucoup plus de temps au systeme.
J'aurais aimé faire la meme chose, mais en passant par l'indéxation et la vectorisation, car cette methode est beaucoup plus rapide.
Je pense avoir compris la methode de DutMatlab, et ca me parait juste, mais le résultat obtenu est le meme que l'image avec les bandes rouge.

Merci pour votre aide
Thomas

[Caro-Line]

Le mot-clé c'est que tu veux "enlever" des éléments à ta matrice.
Si je comprends bien tu souhaites supprimer les lignes et colonnes de zeros.
si c'est bien ça une méthode (il y en a sûrement d'autres...) :

Code :

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

1
2
3
4
5
6
7
8
9
10
 
%Ecriture d'une matrice avec des lignes et colonnes complètes de 0
A=rand(5,6);
A=[A;repmat(zeros(1,size(A,2)),3,1)];%ajoute 3 lignes de zeros
A=[A,repmat(zeros(size(A,1),1),1,5)]%ajoute 5 colonnes de zeros
cols=find(all(A==0)) %trouve les colonnes remplies que de 0
B=A;
B(:,cols)=[] %enleve les colonnes remplies que de 0
rows = find(all(B'==0)) %trouve les lignes remplies que de 0
B(rows,:)=[] %enleve les lignes remplies que de 0

[Jerome Briot]

D'autres méthodes :

Code :

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

1
2
3
4
5
A=rand(5,6);n=2;m=3;
A=[A zeros(size(A,1),n) ; zeros(m,size(A,2)+n)]
B=A;
B(:,~any(A,1))=[];
B(~any(A,2),:)=[]

Code :

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

1
2
3
4
A=rand(5,6);n=2;m=3;
A=[A zeros(size(A,1),n) ; zeros(m,size(A,2)+n)]
[I,J]=find(A);
B(1:I,1:J)=A(1:I,1:J)

Et la plus tricky comme on dit sur le CSSM :

Code :

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

1
2
3
4
A=rand(5,6);n=2;m=3;
A=[A zeros(size(A,1),n) ; zeros(m,size(A,2)+n)]
B=deblank(A.');
B=deblank(B.')

[totoc1001]

Merci beaucoup,

Je v essayer d'implémenter ca, mais j'ai pas mal de difficultés a comprendre les vectorisations.
Par exemple d'apres vos exemple, pour vectoriser une condition If, il faut passer par la fonction find.
Et pour éliminer les boucles for, il suffit de définir entre parenthèses les rages pour les matrices?
C'est bien ca?

J'ai essayé de trouver des docments sur internet qui explique ceci plus précisément, mais c'est assez difficile. Je pense que le seul moyen de bien comprendre c'est de pratiquer.

Si vous savez ou je peux trouver de la doc a ce sujet, je vous en serait très reconnaissant.

Merci beaucoup

J'essaye d'implémenter ca et si j'y arrive pas je vous dis ;-)

Thanks
Thomas

[Jerome Briot]

A lire : Code Vectorization Guide

[totoc1001]

Merci bien, je vais regarder ca.

THomas

[totoc1001]

Code :

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

1
2
A=rand(5,6);n=2;m=3;
A=[A zeros(size(A,1),n) ; zeros(m,size(A,2)+n)]

Bonsoir,

J'ai un peu du mal a comprendre ta méthode cidessus.
Je voudrais m'en servir pour redimensionner une matrice et l'agrandir en ajoutant des zeros a la fin.
Je voulais appliquer ton cas a ma matrice, mais je trouve des résultats, bizarres.
Pourrais-tu m'expliquer ce que fait exactement ton opératioN?

Merci beaucoup
Thomas

[Jerome Briot]

Dans le code, n est le nombre de colonnes de 0 et m est le nombre de lignes de 0 ajoutées ajoutées à la fin de la matrice.

Il est possible de décomposer la méthode comme ceci :

Code :

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

1
2
3
4
5
6
A=rand(5,6);
n=2;
m=3;
s=size(A);
A = [A zeros(s(1),n)]
A = [A ; zeros(m,s(2)+n)]

[totoc1001]

OK merci, c'est beaucoup plus clair.
Je vois que Matlab permet de beaucoup réduire le nombre de commandes. C'est vraiment bien.

Merci
Thomas

[totoc1001]

Rebonsoir a tous.

J'ai encore et toujours des questions lol.
Voici mon code:

Code :

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

1
2
3
4
[y,x]=find(MyRotatedTest(1:OrdinateSecondImage,1:AbscissaSecondImage)==256)
MyRotatedTest(y,x)=0;
MyRotRefX=x;
MyRotRefY=y;

Simple, je veux chercher le pixel qui a une intensité de 256 (il n'y en a qu'un "NORMALEMENT").
Je relève ses coordonnées que j'utilise après.

Mon problème c'est si il n'y a pas de point d'intensité 256. Comment faire pour assigner a y et x, la valeur 0, au lieu de vide, ce qui fait planter mon programme.

Merci beaucoup.
Thomas

[totoc1001]

Bonsoir,

J'aurais eu besoin de faire cette opération en utilisant les vecteurs:

CircleEquation=sqrt(((x-XRotationCenter)^2)+((y-YRotationCenter)^2));

ca ne fonctionne pas avec les vecteurs, mais si j'utilise 2 boucles for, ca fonctionne très bien?
y et x sont de taille différente.
quelqu'un peut-it m'aider?

Merci beaucoup
Thomas