Discussion :

[PaM...]

Bonjour,

J'ai 2 petites questions.
La première:
J'ai un bout de code dont je me demandais si y avait pas moyen d'optimiser, le voici, dites moi ce que vous en pensez.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
l=size(matr);
a=[0 0 0]; 
 
for i=1:l(1)
    for j=1:l(2)
            if matr(i,j)~=0
                a=[a; i j ];
            end
    end
end
a=a(2:length(a),:);

La seconde:
J'aimerai calculer la surface occupé par mon nuage de point formé par les ligne de a (cette surface n'est pas proportionnel à la taille de a car la répartition des points n'est pas homogène). Existe t-il une fonction matlab qui fait ca ?

Merci d'avance

[Caro-Line]

1ère question :
note : ton code est faux il généère une erreur car tu concatènes un vecteur de longueur 3 avec des vecteurs de longueur 2.

Objetcif : récupérer les indices i,j tels que mat(i,j) différent de 0 :

Code :

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

1
2
[idl,idc]=find(matr~=0);
a=[idl,idc];

Ceci dit quel est l'intérêt car il est souvent inutile de passer par ces indices et de préférer l'indexage logique.

2ème question :
Je sais pas.
Avant de se demander s'il existe une fonction MATLAB toute faite : existe-t-il un algorithme pour faire cela, un calcul ?

[Jerome Briot]

Il faudrait que tu nous montres le nuages de points et le polygone représentant les bord de la surface à calculer... pour préciser un peu plus ce que tu entends par "surface occupé par mon nuage de point"

[PaM...]

Bon je pense qu'il va falloir que j'explique tout:

J'ai un volume représenté par des valeurs dans un matrice 3d (matr), je dois trouver l'angle de rotation du volume qui minimise la surface de la projection de ce volume sur le plan 0xy.

Voici mon programe:

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
l=size(matr);
a=[0 0 0]; %initialisation
 
for i=1:l(1)
    for j=1:l(2)
        for k=1:l(3)
            if matr(i,j,k)~=0
                a=[a; i j k];
            end
        end
    end
end
a=a(2:length(a),:);  %on vire le premier vecteur qui a servi a initialiser a
a=transpose(a);      % "a" contient les coordonées de mon volume
 
 
for alpha=0:pi/2:2*pi
            rot=rotx(alpha);   % rotx fabrique ma matrice de rotation
            b=rot*a;             
            b=b(1:2,:);         % je supprime les coordonées en z
            figure
            plot(b(1,:),b(2,:),'x')
 
            % il faudrait maintenant que je puisse calculer la surface occupé
            % par ces points (b), pour trouver le alpha qui va bien
end

EDIT: j'ai trouvé comment calculé cette surface:
j'ai fait un changement de repere des coordonnée dans b pour quel ne soit que positif puis je les ai arrondi pour quelles soient entiere , puis j'ai pris une matrice de 0 assez grande et mis des "1" au coordonnée . Enfin je n'ai plus qu'a compter mon nombre de 1. Je sais pas si c'est tres clair, si ca interesse quelqun je peux donner le code.
Il reste néamoin une questions:

par contre caro-line
je comprend pas, j'ai remplacé le début de mon code par:

Code :

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

1
2
3
[idl,idc,idcf]=find(matr~=0);
a=[idl,idc,idcf];
a=transpose(a);

mais je n'obtient pas le même résultat, j'aimerai pouvoir quand même l'utiliser parcque c'est nettement plus rapide.

[Jerome Briot]

mais je n'obtient pas le même résultat, j'aimerai pouvoir quand même l'utiliser parcque c'est nettement plus rapide.

Simplement parce que FIND ne fonctionne pas en 3D. La 3ème valeur retournée (ici idcf) ne correspond aux indices selon la 3ème dimension de matr mais :

[row,col,v] = find(X, ...) returns a column or row vector v of the nonzero entries in X

Voir les exemples de la documentation pour plus d'informations

Il faut faire comme ceci :

Code :

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

1
2
3
4
idx = find(matr~=0);
[idl,idc,idcf]=ind2sub(size(matr),idx);
a=[idl,idc,idcf];
a=transpose(a);

Ce qui se simplifie grandement en :

Code :

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

1
2
3
4
 
idx = find(matr~=0);
% Si a existe déjà, faire a = []; avant la ligne suivante
[a(1,:),a(2,:),a(3,:)] = ind2sub(size(matr),idx);

Les valeurs sont les même mais n'apparaissent pas de le même ordre qu'avec la boucle FOR-END. Pour preuve :

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
%% METHODE 1
l=size(matr);
a=[0 0 0]; %initialisation
 
for i=1:l(1)
    for j=1:l(2)
        for k=1:l(3)
            if matr(i,j,k)~=0
                a=[a; i j k];
            end
        end
    end
end
a=a(2:length(a),:);  %on vire le premier vecteur qui a servi a initialiser a
a=transpose(a);      % "a" contient les coordonées de mon volume
 
%% METHODE 2
idx = find(matr~=0);
[b(1,:),b(2,:),b(3,:)] = ind2sub(size(matr),idx);
 
%% VERIFICATION
a = sortrows(a.');
b = sortrows(b.');
all(a(:)==b(:))

[Jerome Briot]

Pour ton problème de calcul de surface, es-tu sûr de la méthode choisie ?

Dans le cas contraire, et si tu veux d'autres avis, je te conseille de pauser une question sur le choix de la méthode ici : Algorithmes > Mathématiques (en donnant le lien de cette discussion )

[PaM...]

Merci pour cette réponse

J'ai testé et ton code marche pour une matrice de test de 3*2*2, mais des que je passe avec mes données (matrice de 512*512*128), j'obtiens ce message d'erreur:

Code :

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

1
2
???  In an assignment A(matrix,:) = B, the number of elements in the subscript of A and
the number of rows in B must be the same.

Apparement la fonction ind2sub(size(matr),find(matr~=0)) retourne le même résultat que find(matr~=0).
Je comprend pas ...

Et pour le calcul de la surface ca à l'air de fonctionner pour les premières données que j'ai passé, tant que je rencontre pas d'erreur je creuse pas plus

[Jerome Briot]

Code :

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

1
2
3
4
idx = find(matr~=0);
% Si a existe déjà, faire a = []; avant la ligne suivante
[a(1,:),a(2,:),a(3,:)] = ind2sub(size(matr),idx);

[PaM...]

Non j'ai essayé c'est pas ca


m'enfin c'est pas trop grave c'est du détail
comment on fait pour mettre le tag résolu ?

[nahouto]

comment on fait pour mettre le tag résolu ?

En bas de cette page, il y a au centre un bouton "Résolu"