Plotting concentric circles on sphere

Question

Adam Danz 2015-6-23

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/224988-plotting-concentric-circles-on-sphere

评论： darova 2019-8-19

Hello matlabers,

I've been trying to plot concentric circles on the surface of a sphere as in the left half of the attached image. Ideally I'd like to make that plot and then perform a coordinate transformation of the circles into Cartesian coordinates which would produce the image on the right of the attachment. Anyone have experience in plotting concentric circles on the surface of a sphere?

Adam

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Adam Danz 2015-6-23

concentric c.png

I forgot to attach the picture. here it is https://dl.dropboxusercontent.com/u/2388663/concentric%20c.png

请先登录，再进行评论。

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

Answer 1

Mukul Rao 2015-6-24

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/224988-plotting-concentric-circles-on-sphere#answer_183835

编辑：Mukul Rao 2015-6-24

在 MATLAB Online 中打开

Hi, please execute the code below to generate the sphere with concentric circles. There are many ways to implement this , I have just suggested one such workflow. To better understand what the terms mean, I would recommend reading a little about spherical co-ordinates to interpret the code.

https://en.wikipedia.org/wiki/Spherical_coordinate_system

%Generate custom sphere or use the sphere function
%Equation of sphere polar co-ordinates defined by theta , phi , z
R = 5; %Sphere radius
N = 200; %Number of phi-theta divisions
phi = linspace(0,pi,N);
theta = linspace(0,2*pi,N);
[PHI,THETA] = meshgrid(phi,theta);
Z = R* cos(PHI);
Rcircle = sqrt(R^2 - Z.^2);
%Map from spherical to cartesian while generating surf-plot
surf(Rcircle.*cos(THETA),Rcircle.*sin(THETA),Z,'LineStyle','none')
%Define the phi-coordinates for the circles you'd like to plot on the sphere
phiplot = linspace(0,pi,10) ; %You can replace this array to plot multiple circles of your choice
Zplot = R*cos(phiplot);
Rcircleplot = sqrt(R^2 - Zplot.^2); %Radius of the corresponding circle
hold on
plot3(cos(theta')*Rcircleplot,sin(theta')*Rcircleplot...
    , Zplot'*ones(size(theta)),'LineWidth',2)