Discussion :

[LordPeterPan2]

Bonjour,

Je voudrais trouver une manière assez rapide pour parcourir une matrice n x m de la façon suivant (dans l'exemple n = 6 et m=4) :

Code :

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

1
2
3
4
5
6
7
 
      |   1   2   6    7 |
      |   3   5   8   14 |
      |   4   9  13   15 |
      |  10  12  16   21 |
      |  11  17  20   22 |
      |  18  19  23   24 |

Ma première idée est de créer la matrice ci dessous pour ensuite me servir de l'indexage des matrices, mais là encore ce n'est pas évident (en ne faisant pas de boucle for).

Pour info je cherche à refaire l'algorithme de codage JPEG.

[Jerome Briot]

=> FEX : Zig-zag and De zig-zag scan

[Jerome Briot]

Sinon en fouillant dans mes archives personnelles

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
function [r,c]=zigzag(n,m)
% ZIGZAG...(Complete H1 line)
%   ZIGZAG...(Complete file description)
%
 
%   Author: Jérôme Briot, Matlab 6.1.0.450 (R12.1)
%   Contact: dutmatlab@yahoo.fr
%   Revision: 1.0 (15-Feb-2007)
%   Comments:
%
 
% Check the number of input arguments
error(nargchk(1,2, nargin))
 
if nargin==1
    m=n;
end
 
 
idx=reshape(1:n*m,n,m);
idx=idx(:,end:-1:1);
 
r=zeros(1,numel(idx));
 
flag=0;
x=1;
j=1;
 
if n>m
 
    idx=rot90(idx,3);
    idx=idx(:,end:-1:1);
    [n,m]=size(idx);
end
 
if n==m
    dim=min([1:m-1 m m-1:-1:1],n);
else
    dim=min([1:n repmat(n,1,m-n-1) n:-1:1],n);
end
 
 
for k=m-1:-1:1-m
 
	if isempty(diag(idx,k))
 
        break
	end
 
	if ~flag
        r(x:x+dim(j)-1)=diag(idx,k);
	else
        r(x+dim(j)-1:-1:x)=diag(idx,k);
	end
 
	flag=~flag;
	x=x+dim(j);
    j=j+1;
 
end
 
if n>m
    r=r(end:-1:1);
end
 
if nargout==2
    [r,c]=ind2sub(size(idx),r);
end

A utiliser comme ceci :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
M = [1   2   6    7
     3   5   8   14
     4   9  13   15
    10  12  16   21
    11  17  20   22
    18  19  23   24];
 
M = M.';
 
idx = zigzag(size(M,1),size(M,2));
 
M(idx)

Le code n'est pas vraiment achevé mais c'est l'idée...

[LordPeterPan2]

Merki bicoup,

Je préfère (à priori) la seconde solution qui ne se limite pas aux matrices carrées

[LordPeterPan2]

En faite çà fait un zigzag dans l'autre sens :

Avec ta macro la matrice (pour une 3 x 5) est lu comme suit :

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
 
>> A=reshape(1:15,3,5)
 
A =
 
     1     4     7    10    13
     2     5     8    11    14
     3     6     9    12    15
 
>> zigzag(3,5)
 
ans =
 
  Columns 1 through 14
 
     1     2     4     7     5     3     6     8    10    13    11     9    12    14
 
  Column 15
 
    15

Alors que je voudrai :

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
 
>> A=reshape(1:15,3,5)
 
A =
 
     1     4     7    10    13
     2     5     8    11    14
     3     6     9    12    15
 
>> zigzag(3,5)
 
ans =
 
  Columns 1 through 14
 
     1     4     2     3     5     7    10     8     6     9    11    13    14    12
 
  Column 15
 
    15

[LordPeterPan2]

Merci Dut, en modifiant ton code (y avait pas de copyright ) j'ai réussi à faire ce que je voulais, je met le code aux cas où çà intéresse quelqu'un d'autre :

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
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
function [vect_elt]=zigzag(A);

[n,m] = size(A) ;
vect_elt = A(zigzag_aux(n,m));

function [r,c]=zigzag_aux(n,m)
% ZIGZAG...(Complete H1 line)
%   ZIGZAG...(Complete file description)
%

%   Author: Jérôme Briot, Matlab 6.1.0.450 (R12.1)
%   Contact: dutmatlab@yahoo.fr
%   Revision: 1.0 (15-Feb-2007)
%   Comments:
%

% Check the number of input arguments
error(nargchk(1,2, nargin))

if nargin==1
    m=n;
end

idx=reshape(1:n*m,n,m);
idx=idx(:,end:-1:1);

r=zeros(1,numel(idx));

flag=0;
x=1;
j=1;
flip = 0 ;

if n>m
    idx=rot90(idx,3);
    idx=idx(:,end:-1:1);
    [n,m]=size(idx);
    flip = 1;
end

if n==m
    dim=min([1:m-1 m m-1:-1:1],n);
else
    dim=min([1:n repmat(n,1,m-n-1) n:-1:1],n);
end


for k=m-1:-1:1-m

	if isempty(diag(idx,k))
        
        break
	end

    % 	if ~flag
	if flag
        r(x:x+dim(j)-1)=diag(idx,k);
	else
        r(x+dim(j)-1:-1:x)=diag(idx,k);
	end
    
	flag=~flag;
	x=x+dim(j);
    j=j+1;
    
end

if flip
    r=r(end:-1:1);
end

if nargout==2
    [r,c]=ind2sub(size(idx),r);
endfunction [vect_elt]=zigzag(A);

[n,m] = size(A) ;
vect_elt = A(zigzag_aux(n,m));

function [r,c]=zigzag_aux(n,m)
% ZIGZAG...(Complete H1 line)
%   ZIGZAG...(Complete file description)
%

%   Author: Jérôme Briot, Matlab 6.1.0.450 (R12.1)
%   Contact: dutmatlab@yahoo.fr
%   Revision: 1.0 (15-Feb-2007)
%   Comments:
%

% Check the number of input arguments
error(nargchk(1,2, nargin))

if nargin==1
    m=n;
end

idx=reshape(1:n*m,n,m);
idx=idx(:,end:-1:1);

r=zeros(1,numel(idx));

flag=0;
x=1;
j=1;
flip = 0 ;

if n>m
    idx=rot90(idx,3);
    idx=idx(:,end:-1:1);
    [n,m]=size(idx);
    flip = 1;
end

if n==m
    dim=min([1:m-1 m m-1:-1:1],n);
else
    dim=min([1:n repmat(n,1,m-n-1) n:-1:1],n);
end


for k=m-1:-1:1-m

	if isempty(diag(idx,k))
        
        break
	end

 	if flag
        r(x:x+dim(j)-1)=diag(idx,k);
	else
        r(x+dim(j)-1:-1:x)=diag(idx,k);
	end
    
	flag=~flag;
	x=x+dim(j);
    j=j+1;
    
end

if flip
    r=r(end:-1:1);
end

if nargout==2
    [r,c]=ind2sub(size(idx),r);
end