Shading error when using surf to plot symmetric data set

Question

Lucas 2015-5-3

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/214873-shading-error-when-using-surf-to-plot-symmetric-data-set

评论： Nalini Vishnoi 2015-5-5

I'm trying to plot a temperature distribution using the surf command, but I've run into an error where my graph ends up shading the domain differently in the upper right quadrant despite the fact that the data set is symmetric. Here's a simple example to illustrate my problem:

x=[-1 0 1]

y=[-1 0 1]

z=

   0     0
   1     0
   0     0

surf(x,y,z)

This results in the following plot:

No matter how many additional data points I add, I still run into this problem. Why is this happening? Is there something about the surf command that I don't understand? Thanks in advance for the help.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Nalini Vishnoi 2015-5-5

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/214873-shading-error-when-using-surf-to-plot-symmetric-data-set#answer_177805

编辑：Nalini Vishnoi 2015-5-5

在 MATLAB Online 中打开

Hi Lucas,

You are correct, this behavior is reproducible and can be explained as follows: The 'Facecolor' property of surf (patch object) can be set to 'none', 'flat' and 'interp': link to the documentation. The default value of the 'FaceColor' is 'flat' which uses uniform face colors. The color data at the first vertex determines the color for the entire face (and hence one quadrant was red since its first vertex had value '1'). You cannot use this value when the 'FaceAlpha' property is set to 'interp. Find below the same surf plot with 'FaceColor' property being set to 'interp':

>> s = surf(x,y,z);
>> set(s,'FaceColor', 'interp');

Now if we rotate the plot appropriately we can see that each quadrant is made up of two triangles (since all the 4 points that make up the patch are not co-planar, so the patch needs to be drawn as two triangles). You can also notice how the triangles are arranged/joined differently in the lower left and lower right quadrants (the demarcation is shown by red and green lines respectively):

>> l1 = line([0 1],[-1 0],[0 0],'Color','green');
>> l2 = line([-1 -0],[-1 0],[0 1],'Color','red');

The order and the way in which the triangles are created is decided by the 'surf' command. We can explicitly define the order of the faces/triangles, for example in the following way (though it would be a little tedious):

>> [X,Y] = meshgrid(-1:1,-1:1);
>> Z = [0 0 0; 0 1 0; 0 0 0];
>> V = [X(:) Y(:) Z(:)]; % vertices in 3D 
>> F = [2 5 3; 3 5 6; 6 5 9; 9 5 8; 8 5 7; 7 5 4; 4 5 1; 1 5 2]; % order of the vertices making up the patch
>> p = patch('Faces',F,'Vertices',V); 
>> set(p,'FaceVertexCData',Z(:)); 
>> set(p,'FaceColor','interp');