Message Boards Message Boards

[GIF] Around and Around (Stereographic projection of rotating circles)


Stereographic projection of rotating circles

Around and Around

This is the same basic idea as Move Along; quoting from my description there:

Consider a circle of latitude at height $z=-3/4$ on the unit sphere. The orthogonal complements of unit vectors on the circle are planes through the origin in $\mathbb{R}^3$, which we can also interpret as lines in $\mathbb{RP}^2$, the real projective plane.

First of all, the more obvious identification of a point on the sphere with a point in $\mathbb{RP}^2$ is simply that a point on the sphere determines a unique line through the origin. That being said, the point of Move Along was really that we can interpret (most) points on the sphere as lines (not through the origin) in the plane.

Anyway, for this GIF, instead of thinking of the orthogonal complement of the unit vector (or line) as a plane, take the intersection of that plane with the unit sphere, giving a great circle on the sphere (equivalently, take all points on the sphere at a distance of $\pi/2$ from the original point).

Since I'm considering a circle of points on the sphere, this gives a whole circle's worth of great circles on the sphere. In the animation, I stereographically project several of these great circles to the plane while rotating the original circle around the $y$-axis. To be slightly glib, this shows a representative selection from a circle of circles rotating in a circle.

Stereo[{x_, y_, z_}] := {x/(1 + z), y/(1 + z)};

m[t_, ϕ_, θ_] := RotationMatrix[t, {Cos[θ] Sin[ϕ], Sin[θ] Sin[ϕ], Cos[ϕ]}];
p[r_, s_, t_, ϕ_, θ_] := m[t, ϕ, θ].{Sqrt[1 - r^2] Cos[s], Sqrt[1 - r^2] Sin[s], r};
basis[p_] := Orthogonalize[NullSpace[{p}]];

circle3[p_, ψ_] := Cos[ψ] basis[p][[1]] + Sin[ψ] basis[p][[2]];
stereocircle3[p_, ψ_] := Stereo[circle3[p, ψ]];

With[{θ = π/2., t = π/1., r = -0.75, cols = RGBColor /@ {"#FC85AE", "#303A52"}},
     Evaluate@Table[stereocircle3[p[r, s, t, ϕ, θ], ψ], {s, 0., 2 π - 2 π/12, 2 π/12}],
     {ψ, 0, 2 π}, 
     PlotStyle -> Directive[Thickness[.01 + .004 i], cols[[i]]], Axes -> False],
    {i, {1, -1}}],
   PlotRange -> 2.4, ImageSize -> 540, Background -> cols[[-1]]],
  {ϕ, 0, π/2}]
POSTED BY: Clayton Shonkwiler
5 months ago

enter image description here - Congratulations! This post is now Staff Pick! Thank you for your wonderful contributions. Please, keep them coming!

POSTED BY: Moderation Team
5 months ago

Group Abstract Group Abstract