Plotting multi-radius spherical mesh

2014 5 16
2 个回答
回答已采纳
更新时间：2014 5 19
4 次查看（30 天）

1 个投票

I'm attempting to plot a mesh which represents a sphere sectioned into 12 logitudinal and 22 latitudinal segments. The difficultly is that each of the segments has a different radius to show it's importance. I have used the code below to draw the sphere but can not work out how to impose the importance onto the mesh points. The commented code if used with phi and theta of lengths 12, and 22 respectively gives me points, I however, wanted the whole segment surface.
out=double(int32(rand(22,12)*9));
% phi=linspace(0,pi,12);
% theta=linspace(0,2*pi,22);
phi=linspace(0,pi,13);
theta=linspace(0,2*pi,23);
[phi,theta]=meshgrid(phi,theta);
for i = 1:22
for j=1:12
% x=out.*sin(phi).*cos(theta);
% y=out.*sin(phi).*sin(theta);
% z=out.*cos(phi);
x=out.*sin(phi).*cos(theta);
y=out.*sin(phi).*sin(theta);
z=out.*cos(phi);
end
end
mesh(x,y,z)

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

 采纳的回答

Daniel
Daniel 2014-5-19

3 个投票

For any one that is interested I have used the following code to anwser the question;
out=double(int32(rand(22,12)*9));
phi=linspace(0,pi,13);
theta=linspace(0,2*pi,23);
phi=[phi; phi];
theta=[theta; theta];
phi=reshape(phi,1,2*length(phi));
theta=reshape(theta,1,2*length(theta));
[phi,theta]=meshgrid(phi,theta);
r = zeros(size(phi,1),size(phi,2));
for i = 1:12
for j=1:22
r(j*2:j*2+1,i*2:i*2+1) = out(j,i);
end
end
x=r.*sin(phi).*cos(theta);
y=r.*sin(phi).*sin(theta);
z=r.*cos(phi);
mesh(x,y,z)
axis square
axis off
hidden off

更多回答（1 个）

Daniel
Daniel 2014-5-16
编辑：Daniel 2014-5-16

2 个投票

I do belive I have got much closer to the answer, see code below. So the question is how do I manipulate out a [22,12] matrix to scale the X,Y,Z values which are [46,26] matrix of the vertics.
out=double(int32(rand(22,12)*9));
phi2=linspace(0,pi,13);
theta2=linspace(0,2*pi,23);
phi3=[phi2; phi2];
theta3=[theta2; theta2];
phi4=reshape(phi3,1,2*length(phi3));
theta4=reshape(theta3,1,2*length(theta3));
[phi5,theta5]=meshgrid(phi4,theta4);
x=1*sin(phi5).*cos(theta5);
y=1*sin(phi5).*sin(theta5);
z=1*cos(phi5);
mesh(x,y,z)out=double(int32(rand(22,12)*9));

类别

在 帮助中心 和 File Exchange 中查找有关 Surface and Mesh Plots 的更多信息

产品

标签

Community Treasure Hunt

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

Start Hunting!

Translated by