Draw contourslice on patch surface

4 次查看（过去 30 天）

显示更早的评论

Nathaniel H Werner 2019-8-2

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/474609-draw-contourslice-on-patch-surface

评论： Nathaniel H Werner 2019-8-6

采纳的回答： darova

Data_8519.mat

在 MATLAB Online 中打开

I am trying to get the XData and YData on a contourslice in matlab. I want this contourline to be on a patch that I have already made in 3D. The patch corresponds to a rectangular plane orthogonal to the -x axis rotated about the y axis from -x by some angle dphi. I have provided some of my code below with images.

AoAd = 45;
Ty = 1/2*cosd(AoAd);   % wing position in y
Tz = 1/2*sind(AoAd);   % wing position in z
yslice = [Ty, Ty, -Ty, -Ty]+[1/2, 1/2, -1/2, -1/2];
zslice = [-Tz, Tz, Tz, -Tz]+[-1/2, 1/2, 1/2, -1/2];
xslice = -.15;
S = [xslice*ones(size(yslice));
yslice;
zslice];
dphi = 7.5;    % degrees
MR = [[cosd(dphi) 0 -sind(dphi) 0];...
        [0 1 0 0];...
        [sind(dphi) 0 cosd(dphi) 0];...
        [0 0 0 1]];
S_prime = MR*[S,[0,0,0]';
[0,0,0,0],1];
x = S_prime(1,1:4);
y = S_prime(2,1:4);
z = S_prime(3,1:4);
figure
hold on
patch(x,y,z,'r')
cc = contourslice(Xw,Yw,Zw,D_ta,z,y,x,[-3,-3],'linear');

Now, I expect the there two be a single iso-line on the red patch. But instead I get several iso-lines that look orthogonal to the red patch. Can I get some help figuring out what I'm doing wrong here.

It seems to me that the smaller region is at least on a parallel plane but I don't know why the others aren't.

Addendum:

I also tried using the XYZ data from patch and this was the result. Still not on the plane, and there are some slices on different planes.

figure
hold on
p = patch(x,y,z,'r');
xp = p.XData;
yp = p.YData;
zp = p.ZData;
cc = contourslice(Xw,Yw,Zw,D_ta,xp,yp,zp,[-3,-3],'linear');

2 个评论
显示无隐藏无

darova 2019-8-3

Just to clarify: you want to project your data onto red plane?

Nathaniel H Werner 2019-8-3

编辑：Nathaniel H Werner 2019-8-3

Yes, I am trying to get a volume slice of the data on the red patch. If you can't run the code try using the built in flow or peaks.

请先登录，再进行评论。

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

采纳的回答

darova 2019-8-3

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/474609-draw-contourslice-on-patch-surface#answer_386041

Untitled.m

Try to refine mesh of your plane

I attached a simple example (help example modified)

Can you attach your data

16 个评论
显示 14更早的评论隐藏 14更早的评论

Nathaniel H Werner 2019-8-6

编辑：Nathaniel H Werner 2019-8-6

在 MATLAB Online 中打开

I'm sorry I was mistaken when I previously accepted the answer.

There is some misunderstanding about how to rotate the plane.

Take this section of code.

yw = linspace(-5,5,50);
zw = linspace(-5,5,50);
figure
RoLoc = linspace(0.1:3,50)
for r=1:length(RoLoc)
    R = [RoLoc(r)*ones(50,1),yw(:),zw(:),ones(50,1)]';
    MR = [cosd(7.5) 0 -sind(7.5) 0;...
        0 1 0 0;...
        sind(7.5) 0 cosd(7.5) 0;
        0 0 0 1];
    R_prime = MR*R;
    xw = R_prime(1,:);
    zw = R_prime(3,:);
    
    [xs,ys,zs] = meshgrid(xw,yw,zw);
    
    hold on
    surf(xs,ys,zs);
    contourslice(Xw,Yw,Zw,D_ta,xs,ys,zs,[-3,-3])
end

What I need is for the contourslice to be on the plane xs, ys, zs. I am unable to create a surf with that 3D meshgrid data. Running the above code gives me this error.

Error using matlab.graphics.chart.primitive.Surface/set

Value must be a vector or 2D array of numeric type

Error in matlab.graphics.chart.internal.ctorHelper (line 8)

set(obj, pvpairs{:});

Error in matlab.graphics.chart.primitive.Surface

Error in surf (line 144)

hh = matlab.graphics.chart.primitive.Surface(allargs{:});

To summarize the process.

The first plane should start parallel to the yz plane (normal aligned with -x).
Rotate the plane about the y-axis by an angle dphi = 7.5 deg measured from the -x axis.
Extract the rotated x, y, and z data.
Make a contourslice on this new rotated surface.

If this is corrected I can accept the answer. I will continue trying to figure out why I am getting this error.

Additionally

If you run the portion of the code just with contourslice you get this error.

Error using contourc

Input arguments for contourc must have at most 2 dimensions.

Error in contourslice (line 239)

c = contourc(vi, cvals);

darova 2019-8-6

Maybe you wanted slice of XZ plane? I like it more

Nathaniel H Werner 2019-8-6

The right image looks exactly like what I need, the red planes are XY or YZ planes rotated around the Y axis. Does the new Untitled3 script make that figure?

请先登录，再进行评论。

类别

MATLAB Graphics 2-D and 3-D Plots

在 Help Center 和 File Exchange 中查找有关 2-D and 3-D Plots 的更多信息

产品

MATLAB

版本

R2018b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by

Draw contourslice on patch surface

2 个评论
显示无隐藏无

采纳的回答

16 个评论
显示 14更早的评论隐藏 14更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

Draw contourslice on patch surface

2 个评论 显示 无隐藏 无

采纳的回答

16 个评论 显示 14更早的评论隐藏 14更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

2 个评论
显示无隐藏无

16 个评论
显示 14更早的评论隐藏 14更早的评论