Поверхность в цилиндрических координатах / неравномерный массив

А=[1,2,3,4,5];
В=[1,2,3]; % только столбцы с заполненными ячейками 
 
q=[2,3,4,6,7; 3,4,5,7,8; 4,5,6,8,9]; % соответсвенно В 
 
[X,Y] = meshgrid(А,В);
[Q]=griddata(А,В,q,X,Y,'cubic');
 
[R,TH] = cart2pol(Y,X);
surf(R.*cos(TH), R.*sin(TH),Q);

q(:,1)={VarName1(1,1), VarName1(2,1), VarName1(3,1), VarName1(4,1), VarName1(5,1), VarName1(6,1), VarName1(7,1), VarName1(8,1), VarName1(9,1), VarName1(10,1), VarName1(11,1), VarName1(12,1), VarName1(13,1), VarName1(14,1), VarName1(15,1), VarName1(16,1), VarName1(17,1), VarName1(18,1), VarName1(19,1), VarName1(20,1), VarName1(21,1), VarName1(22,1), VarName1(23,1), VarName1(24,1), VarName1(25,1), VarName1(26,1), VarName1(27,1), VarName1(28,1), VarName1(29,1)};

clc; clear all; clf;
M=xlsread('Probe_10.csv')

betta=[0,2,5...180] 
[X,Y] = meshgrid(А,betta);
 
[R,TH] = cart2pol(Y,X);
 
surf(R.*cos(TH), R.*sin(TH),% матрица всех значений В%);

clc; clear all; clf;
 
M=xlsread('Probe_10.csv');
[r,c]=size(M);
A=M(1,2:c)/180*pi;
B=M(2:r,1)/180*pi;
q=M(2:r,2:c);
 
[X,Y] = meshgrid(A,B);
Q=griddata(A,B,q,X,Y,'cubic');
 
[R,TH] = cart2pol(Y,X)
surf(R.*cos(TH), R.*sin(TH),Q);
set(gca,'ZScale','log')

[R,TH] = cart2pol(Y,X);
surf(R.*cos(TH), R.*sin(TH),Q);

A=M(1,2:c)/180*pi;
B=M(2:r,1)/180*pi;

surf(R.*cos(TH/57.3), R.*sin(TH/57.3),Q);

Q=griddata(A,B,q,X,Y,'cubic');

Error using delaunayTriangulation
The coordinates of the input points must be finite values; Inf and NaN are not
permitted.
 
Error in griddata>cubic (line 208)
dt = delaunayTriangulation(x,y);
 
Error in griddata (line 127)
        vq = cubic(x,y,v,xq,yq);