osse/utils/cart2cyl.m

function [z] = cart2cyl(c,thetai,rhoi,varargin)
% z=cart2cyl(c,xi,yi);
%

NN=size(c);
[x y nz]=size(c);
stepx=(x+1)/(x);
stepy=(y+1)/(y);
[Y,X] = meshgrid(-(x)/2:stepx:(x)/2,-(y)/2:stepy:(y)/2);
[theta,rho] = cart2pol(X,Y);
theta=(theta./(pi*2)+0.5).*(size(thetai,2));
rho=max(rhoi)*rho/((x)/2);
ntheta=x;nrho=y;

THETA=reshape(theta,[1 ntheta*nrho]);
RHO=reshape(rho,[1 ntheta*nrho]);
del=griddata_preprocess(RHO,THETA,rhoi,thetai',varargin{:});

for k=1:nz;
 C=reshape(c(:,:,k),[1 ntheta*nrho]);
z(:,:,k)=griddata(RHO,THETA,C,rhoi,thetai',varargin{:});
end % k

% Split vertical and time dimensions
if size(NN,2)>2
z=reshape(z,[size(z,1) size(z,2) NN(3:end)]);
end
1	afe	1.1	function [z] = cart2cyl(c,thetai,rhoi,varargin)
2			% z=cart2cyl(c,xi,yi);
3			%
4
5			NN=size(c);
6			[x y nz]=size(c);
7	afe	1.3	stepx=(x+1)/(x);
8			stepy=(y+1)/(y);
9			[Y,X] = meshgrid(-(x)/2:stepx:(x)/2,-(y)/2:stepy:(y)/2);
10	afe	1.1	[theta,rho] = cart2pol(X,Y);
11	afe	1.3	theta=(theta./(pi2)+0.5).(size(thetai,2));
12			rho=max(rhoi)*rho/((x)/2);
13	afe	1.1	ntheta=x;nrho=y;
14
15			THETA=reshape(theta,[1 ntheta*nrho]);
16			RHO=reshape(rho,[1 ntheta*nrho]);
17			del=griddata_preprocess(RHO,THETA,rhoi,thetai',varargin{:});
18
19			for k=1:nz;
20			C=reshape(c(:,:,k),[1 ntheta*nrho]);
21			z(:,:,k)=griddata(RHO,THETA,C,rhoi,thetai',varargin{:});
22			end % k
23
24			% Split vertical and time dimensions
25			if size(NN,2)>2
26			z=reshape(z,[size(z,1) size(z,2) NN(3:end)]);
27			end