utils/matlab/tiles.m

function [a] = tile(b,n)
% a=tile(b,n);
%
% Extract single tile from cubed array
%
% b can have dimensions (M*6,M,Nr) or (M,M,Nr,6)
%
% n can be vector of integers between 1 and 6


if min(n)<1 | max(n)>6
 disp(sprintf('n=',n));
 error('tile: second argument n is out of range');
end

if size(b,ndims(b))==6
 switch ndims(b)
  case 3,
   a=b(:,:,n);
  case 4,
   a=b(:,:,:,n);
  otherwise
   error('tile: it seems that b has too many dimensions');
 end
elseif size(b,2)==6
 m=size(b,1);
 k=1;
 for N=n;
 switch ndims(b)
  case 3,
   a(:,:,k)=squeeze(b(:,N,:));
  case 4,
   a(:,:,:,k)=squeeze(b(:,N,:,:));
  otherwise
   error('tile: it seems that b has too many dimensions');
 end
 k=k+1;
 end
elseif size(b,1)==size(b,2)*6
 m=size(b,2);
 k=1;
 for N=n;
 switch ndims(b)
  case 2,
   a(:,:,k)=b((N-1)*m+1:N*m,:);
  case 3,
   a(:,:,:,k)=b((N-1)*m+1:N*m,:,:);
  otherwise
   error('tile: it seems that b has too many dimensions');
 end
 k=k+1;
 end
else
 disp(sprintf('Size(b) = %i %i %i %i %i %i',size(b)));
 error('tile: Size of first argument is not consistent with cubed array');
end
1	function [a] = tile(b,n)
2	% a=tile(b,n);
3	%
4	% Extract single tile from cubed array
5	%
6	% b can have dimensions (M*6,M,Nr) or (M,M,Nr,6)
7	%
8	% n can be vector of integers between 1 and 6
9
10
11	if min(n)<1 \| max(n)>6
12	disp(sprintf('n=',n));
13	error('tile: second argument n is out of range');
14	end
15
16	if size(b,ndims(b))==6
17	switch ndims(b)
18	case 3,
19	a=b(:,:,n);
20	case 4,
21	a=b(:,:,:,n);
22	otherwise
23	error('tile: it seems that b has too many dimensions');
24	end
25	elseif size(b,2)==6
26	m=size(b,1);
27	k=1;
28	for N=n;
29	switch ndims(b)
30	case 3,
31	a(:,:,k)=squeeze(b(:,N,:));
32	case 4,
33	a(:,:,:,k)=squeeze(b(:,N,:,:));
34	otherwise
35	error('tile: it seems that b has too many dimensions');
36	end
37	k=k+1;
38	end
39	elseif size(b,1)==size(b,2)*6
40	m=size(b,2);
41	k=1;
42	for N=n;
43	switch ndims(b)
44	case 2,
45	a(:,:,k)=b((N-1)m+1:Nm,:);
46	case 3,
47	a(:,:,:,k)=b((N-1)m+1:Nm,:,:);
48	otherwise
49	error('tile: it seems that b has too many dimensions');
50	end
51	k=k+1;
52	end
53	else
54	disp(sprintf('Size(b) = %i %i %i %i %i %i',size(b)));
55	error('tile: Size of first argument is not consistent with cubed array');
56	end