Draw contourslice on patch surface

Question

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

Data_8519.mat

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.

请先登录，再进行评论。

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

Answer 1

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-3

在 MATLAB Online 中打开

ForDrawingOnPlane.fig

Well it works, and it doesn't work. For the planes you made it works it just fine with any 3D data. But I need to use the rotated plane described something like this.

[ys,zs] = meshgrid(-2:0.2:2);
xs = -.1*ones(size(ys));
xv = xs(:)';
yv = ys(:)';
zv = zs(:)';
M = [cosd(angle) 0 -sind(angle) 0;...
    0 1 0 0;
    sind(angle) 0 cosd(angle) 0;
    0 0 0 1];
v = M*[xv;
    yv;
    zv;
    ones(size(xv))];
xs_prime = reshape(v(1,:),[21,21]);
ys_prime = reshape(v(2,:),[21,21]);
zs_prime = reshape(v(3,:),[21,21]);

You can see that if I make a patch with this data the contourslice is not on the rotated plane.

x = [max(xs_prime(:)) max(xs_prime(:)) min(xs_prime(:)) min(xs_prime(:))];
y = [max(ys_prime(:)) max(ys_prime(:)) min(ys_prime(:)) min(ys_prime(:))];
z = [min(zs_prime(:)) max(zs_prime(:)) max(zs_prime(:)) min(zs_prime(:))];
hold on
% mesh(xs,ys,zs)
h = patch(x,y,z,'r');
alpha(h,0.5)
hold off
view(3)
grid on
xlabel('X axis')
ylabel('Y axis')
axis equal

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?

请先登录，再进行评论。

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更早的评论