Problem with reshape and Surf plot

Question

ragnor 2021-7-26

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/886039-problem-with-reshape-and-surf-plot

评论： Star Strider 2021-7-28

采纳的回答： Star Strider

Data123.xlsx

在 MATLAB Online 中打开

Hello MATLAB Community,

I have the following script and data where i would like to plot a 3d Surface and it works well for some data and it doesn't work because of the dimensions in some cases. Can someone please suggest where is the mistake in the script and how to fix it?

opts = detectImportOptions('Data123.xlsx');
opts = opts.setvartype('Y','double');
opts = opts.setvartype('Z','double');
% Load data
data = readtable('Data123.xlsx', opts);
data.Properties.VariableNames = ["X","F_sw","Y","Z"];
% order the data
data = sortrows(data,["X","Y"]);
X = unique(data.X);
Y = unique(data.Y);
Z = reshape(data.Z,[length(Y),length(X)]);
% View
figure(1)
scatter3(data.X,data.Y,data.Z,12,[1 0 0],"filled","MarkerEdgeColor",[0 0 0])
hold on
s = surf(X,Y,Z,"FaceColor",'interp');
hold off

6 个评论
显示 4更早的评论隐藏 4更早的评论

ragnor 2021-7-26

Dear @Star Strider,

If i don't use unique, is there any other way, i can create a surface plot with the data.

I just want to make a 3d surface plot to display the results.

Thanks in advance.

ragnor 2021-7-26

untitled.fig

Dear @Nikhil Sapre,

I just want to make a 3D surface plot as displayed in the attached figure.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Star Strider 2021-7-27

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/886039-problem-with-reshape-and-surf-plot#answer_754449

在 MATLAB Online 中打开

The only solution I can propose is to use either griddata or scatteredInterpolant .

Try this —

opts = detectImportOptions('https://www.mathworks.com/matlabcentral/answers/uploaded_files/695269/Data123.xlsx');

opts = opts.setvartype('Y','double');

opts = opts.setvartype('Z','double');

% Load data

data = readtable('https://www.mathworks.com/matlabcentral/answers/uploaded_files/695269/Data123.xlsx', opts)

data = 184×4 table

X F_sw Y Z _____ ____ ___ ______ 10000 350 0.5 2.9008 10000 350 1 3.0802 10000 350 1.5 3.2261 10000 350 2 3.3174 10000 350 2.5 3.3299 10000 350 3 3.3618 10000 350 3.5 3.4404 10000 350 4 3.5091 10000 350 4.5 3.5816 10000 350 5 3.6302 10000 350 5.5 3.6945 10000 350 6 3.7262 10000 350 6.5 3.7544 10000 350 7 3.7761 10000 350 7.5 3.7944 10000 350 8 3.8067

data.Properties.VariableNames = ["X","F_sw","Y","Z"];

% order the data

data = sortrows(data,["X","Y"]);

% X = unique(data.X);

% Y = unique(data.Y);

X = data.X;

Y = data.Y;

Z = data.Z;

N = 50;

xv = linspace(min(X), max(X), N);

yv = linspace(min(Y), max(Y), N);

[Xm,Ym] = ndgrid(xv, yv);

Zm = griddata(X,Y,Z,Xm,Ym);

% Z = reshape(data.Z,[length(Y),length(X)]);

% View

figure(1)

scatter3(data.X,data.Y,data.Z,12,[1 0 0],"filled","MarkerEdgeColor",[0 0 0])

hold on

s = surf(Xm,Ym,Zm,"FaceColor",'interp', 'EdgeColor','none');

hold off

If it is not the result you want, I will delete my Answer.

.

4 个评论
显示 2更早的评论隐藏 2更早的评论

Star Strider 2021-7-27

在 MATLAB Online 中打开

As always, my pleasure!

Joining the red dots rerquires switching from scatter3 to plot3 —

opts = detectImportOptions('https://www.mathworks.com/matlabcentral/answers/uploaded_files/695269/Data123.xlsx');

opts = opts.setvartype('Y','double');

opts = opts.setvartype('Z','double');

% Load data

data = readtable('https://www.mathworks.com/matlabcentral/answers/uploaded_files/695269/Data123.xlsx', opts)

data = 184×4 table

X F_sw Y Z _____ ____ ___ ______ 10000 350 0.5 2.9008 10000 350 1 3.0802 10000 350 1.5 3.2261 10000 350 2 3.3174 10000 350 2.5 3.3299 10000 350 3 3.3618 10000 350 3.5 3.4404 10000 350 4 3.5091 10000 350 4.5 3.5816 10000 350 5 3.6302 10000 350 5.5 3.6945 10000 350 6 3.7262 10000 350 6.5 3.7544 10000 350 7 3.7761 10000 350 7.5 3.7944 10000 350 8 3.8067

data.Properties.VariableNames = ["X","F_sw","Y","Z"];

% order the data

data = sortrows(data,["X","Y"]);

X = data.X;

Y = data.Y;

Z = data.Z;

N = 50;

xv = linspace(min(X), max(X), N);

yv = linspace(min(Y), max(Y), N);

[Xm,Ym] = ndgrid(xv, yv);

Zm = griddata(X,Y,Z,Xm,Ym);

% Z = reshape(data.Z,[length(Y),length(X)]);

% View

[Xu,Xix] = unique(data.X,'stable');

Xix = [Xix; numel(X)];

for k = 1:numel(Xix)-1

C{k} = data{(Xix(k):Xix(k+1)-1),[1 3 4]};

end

% C

figure(1)

hold on

% scatter3(data.X,data.Y,data.Z,12,[1 0 0],"filled","MarkerEdgeColor",[0 0 0])

for k = 1:numel(C)

plot3(C{k}(:,1),C{k}(:,2),C{k}(:,3),'o-r', 'MarkerSize',6, 'MarkerFaceColor',[1 0 0],"MarkerEdgeColor",[0 0 0], 'LineWidth',2.5)

end

s = surf(Xm,Ym,Zm,"FaceColor",'interp', 'EdgeColor','none', 'FaceAlpha',0.8);

hold off

xlabel('X')

ylabel('Y')

zlabel('Z')

grid on

view(-30,30)

However this did not turn out to be straightforward, since there were connecting lines between the plotted line segments. Sectioning them in the ‘C’ cell array and then plotting them separately in the loop solved this. (An alternative approach would have been to section them, concatenating ‘NaN(1,3)’ to the end of each section, then recombining them into a single matrix. Thie approach I took here is easier.)

Experiment with the various line properties to get the result you want.

.

ragnor 2021-7-28

@Star Strider: Thank you very much. It was indeed very uselful and it is something i wanted to plot.

Thanks again.

Star Strider 2021-7-28

As always, my pleasure!

请先登录，再进行评论。

Problem with reshape and Surf plot

6 个评论
显示 4更早的评论隐藏 4更早的评论

采纳的回答

4 个评论
显示 2更早的评论隐藏 2更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

Problem with reshape and Surf plot

6 个评论 显示 4更早的评论隐藏 4更早的评论

采纳的回答

4 个评论 显示 2更早的评论隐藏 2更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

6 个评论
显示 4更早的评论隐藏 4更早的评论

4 个评论
显示 2更早的评论隐藏 2更早的评论