Discussion :

[atoto]

Bonjour,

Je dispose de deux matrices composés d'index comme présenté ci-dessous :

Code :

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

1
2
index1=[2;1;3]
index2=[4;2;4]

En posant C la matrice C = [0 0 0 0;0 0 0 0; 0 0 0 0];, je souhaiterai mettre la valeur 1 à tous les éléments compris entre index1 et index2 pour chaque ligne.

Le résultat serait ainsi [0 1 1 1;1 0 0 0;0 0 1 1].
Je souhaite vectorialiser ça car en réalité je dispose de matrice de 100 000 lignes
Je pensais réaliser quelque chose comme ça : C(:,index1(:):index2(:))=1 mais cela ne fonctionne pas.

Merci par avance

Anthony

[Gooby]

Bonjour

Une solution pas optimale en terme de temps mais qui fonctionne:

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
>> C = zeros(3,4)
index1 = [2;1;3]
index2 = [4;2;4]
index = cat(2,index1,index2)
 
for i = 1:size(C,1)
    C(i,index(i,1):index(i,2)) = 1;
end
C
 
C =
 
     0     0     0     0
     0     0     0     0
     0     0     0     0
 
 
index1 =
 
     2
     1
     3
 
 
index2 =
 
     4
     2
     4
 
 
index =
 
     2     4
     1     2
     3     4
 
 
C =
 
     0     1     1     1
     1     1     0     0
     0     0     1     1

[shaiHulud]

Bonjour,

Cette vectorisation ne peut se faire qu'en indicage linéaire. row_idx et col_idx sont les indices de ligne et de colonne de tous éléments à assigner dans C.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
 
row_idx= arrayfun(@(d,f) f-d+1  ,index1,index2); % longueurs de la plage sur chaque ligne
row_idx= arrayfun(@(k,l) k*ones(l,1) ,(1:length(index1))',row_idx(:),'uni',0)
row_idx= cat(1,row_idx{:}); % num. de ligne du kème 1
 
col_idx= arrayfun(@(d,f) (d:f)'  ,index1,index2,'uni',0);
col_idx= cat(1,col_idx{:});  % num. de colonne du kème 1
 
assert(length(col_idx)==length(row_idx));% check dims
 
linearInd =sub2ind(size(C), row_idx, col_idx);
 
C(linearInd) = 1; % <-- Assignation vectorisée