Diagnostics/DiagUtility/tiles.m

function [a] = tile(b,varargin)
% 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 nargin==1
 n=1:6;
else
 n=varargin{1};
end

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,varargin)
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	if nargin==1
11	n=1:6;
12	else
13	n=varargin{1};
14	end
15
16	if min(n)<1 \| max(n)>6
17	disp(sprintf('n=',n));
18	error('tile: second argument n is out of range');
19	end
20
21	if size(b,ndims(b))==6
22	switch ndims(b)
23	case 3,
24	a=b(:,:,n);
25	case 4,
26	a=b(:,:,:,n);
27	otherwise
28	error('tile: it seems that b has too many dimensions');
29	end
30	elseif size(b,2)==6
31	m=size(b,1);
32	k=1;
33	for N=n;
34	switch ndims(b)
35	case 3,
36	a(:,:,k)=squeeze(b(:,N,:));
37	case 4,
38	a(:,:,:,k)=squeeze(b(:,N,:,:));
39	otherwise
40	error('tile: it seems that b has too many dimensions');
41	end
42	k=k+1;
43	end
44	elseif size(b,1)==size(b,2)*6
45	m=size(b,2);
46	k=1;
47	for N=n;
48	switch ndims(b)
49	case 2,
50	a(:,:,k)=b((N-1)m+1:Nm,:);
51	case 3,
52	a(:,:,:,k)=b((N-1)m+1:Nm,:,:);
53	otherwise
54	error('tile: it seems that b has too many dimensions');
55	end
56	k=k+1;
57	end
58	else
59	disp(sprintf('Size(b) = %i %i %i %i %i %i',size(b)));
60	error('tile: Size of first argument is not consistent with cubed array');
61	end