Showing posts with label Stereographic projection. Show all posts
Showing posts with label Stereographic projection. Show all posts

Tuesday, April 29, 2025

Lie Sphere Geometry Part 12: Preimages of planes

Oriented planes in R3 are images by stereographic projection of those oriented spheres in S3 that contain the point -e0. In this post we use the same method as we did it in Part 11 when looking for preimages of oriented spheres. At the end we will have to decide again about the value of  ε to obtain the correct orientation. While for spheres it required some work with somewhat complicated algebraic expression, here it will be much easier - a merry and fast slide downhill.

Easy epsilon

With n a unit vector in R3, and h a real number, consider the oriented plane Πh(n) in R3:

Πh(n) = {yR3: y·n = h}.   (1)

We expect it to be an image, by the stereographic projection, of a certain oriented sphere

St(m) =  {xS3: x·m = cos(t)},     0<t<π ,     (2)

that contains the point -e0=(-1,0,0,0) - the origin of the stereographic projection. Our aim is to find an explicit expression of (t,m) through (n,h). Let us stress that the formulas (1) and (2) should be understood having in mind the fact Πh(n) and St(m)are more than just "sets". They include the orientations, and orientations are contained in the pairs (n,h) and (m,t). Pairs (n,h) and (-n,-h) (resp. (m,t) and (-m,π-t)) define the same sets, but correspond to opposite "orientations".

We will proceed the same way as we have done it with spheres in Part 11. This time it will be even simpler, because the equation (1) defining Πh(n) is linear. Thus, with x=(x0,...,x3)∈R4, ||x||=1, we substitute the stereographic projection formula

yi  =  xi/(1+x0),         (i=1,2,3)        (3)

into (1) to obtain

xn = (1+x0)h = h+x0h,                (4)

where x'=(x1,x2,x3).

Then we rewrite (4) as


- x0h + xn = h,  

which suggests the candidate for m, namely (-h,n). But m should have the norm 1, while (-h,n) has the norm squared h2+n2=1+h2. We must  also remember that (n,-h) and (-n,h) lead to the same sphere as a set, so our general solution is

m = (-εh, εn )/√(1+h2),

t = arccos(εh/√(1+h2)),       0<t<π


where ε=±1, and we have to decide the sign taking into account the orientations. To this end we follow the same path as we have taken in Part 11 with the spheres. From Part 10, Proposition 1 therein, we know that to have the correct orientation we should have h=cot(t).  Now, cos(arccos( εh/√(1+h2) ) )=εh/√(1+h2), while sin( arccos(εh/√(1+h2)) )=(1-(εh/√(1+h2))2 )1/2= 1/√(1+h2). Thus cot(t) = εh, therefore ε=+1. 

Et voilà: ε=+1!

This way we have arrived at the following Proposition:

Proposition 1. With n a unit vector in R3, and h a real number, the pre-image of the oriented plane Πh(n)with respect to the stereographic projection is the oriented sphere Sarccos(h/√(1+h2))( (-h, n )/√(1+h2) ).

Exercise 1. Verify that m0+cos(t)=0 as it should be. (Why?)

In the next post we will return to the 6D universe and we will place our 3D world of  spheres and planes in there.


Thursday, April 24, 2025

Lie Sphere Geometry Part 10: Planes

 Spheres of infinite radius, those that cross the infinity point,  are just planes. Spheres are mysterious, planes are mysterious too, perhaps even more so.

Even more so...


This post is a straightforward continuation of Lie Sphere Geometry: Part 9: Spheres of negative radius. The stereographic projection of the oriented sphere

St(m) = {xS3: x·m = cos(t)},     0<t<π ,      (0)

leads to the following equation (Eq. (4) of Part 9) for its image in R3:


y2(m0+cos(t)) - 2m'·ym0 - cos(t),                (1)

where

m '= m1e1+m2e2+m3e3.                (2)

We have already considered the case of m0+cos(t) ≠ 0. We consider now the case of m0+cos(t) = 0, i.e. m0 = -cos(t). Eq (1) then becomes

m'·y = cos(t).                 (3)     

Definition 1. With nS2R3 and h∈R, the oriented plane Πh(n) with unit normal n and signed height h, is a pair (n,h) and  associated with it the set [Πh(n)] defined as:

h(n)] = {yR3: y·n = h}.                (4)

Remark. To see the meaning of h, suppose we orient the Cartesian axes of R3 so that n is the unit vector in the z-direction. Then h is the value of the z-coordinate at which the plane intersects the z-axis. Thus the name "signed height". Notice that while Πh(n)  and Π-h(-n) define the same set, these are two different oriented planes, as they have opposite normals.

Proposition 1. The image of the oriented sphere St(m)S3 , containing -e0, under the stereographic projection xy  is the oriented Πh(n') , where n' = m'/sin(t), h = cot(t).

Proof.  Comparing Eq. (3) with (4) we get instantly that n' has to be ε=±1 times the normalized vector m'. Now, (m0)2+m'2=1, and m0 = -cos(t). It follows that m'2 = 1- cos2(t) = sin(t)2. Since 0<t<π, we deduce that ||m'|| = sin(t).
Thus n'=ε m'/sin(t), h = ε cot(t). In Part 9 we have calculated (Eq. (11)) the image n'' of n:

n''(x)i = (1+x0)-1(mi/sin(t) -xi(m0+cos(t))/(sin(t)(1+x0)) ,

which now reduces to

n''(x)i = (1+x0)-1(mi/sin(t)).

Since 1+x0 > 0, we deduce that ε = 1.



In the next post we will deal with the inverse stereographic projection of oriented spheres and planes.

Monday, April 21, 2025

Lie Sphere Geometry: Part 9: Spheres of negative radius


The Fourth Way – P D Ouspensky [from the oral teachings of G I Gurdjieff]

"The first step is to try not to express  negative emotions;  - fear, anger, jealousy, possessiveness, pride and so on; the second step is the study of these negative emotions themselves, making lists of them, finding their connections ….. and trying to understand that they are quite useless.

Question:  In some cases the negative emotion of fear seems useful, otherwise people would cross the road at any time without looking

You speak about instinctive fear, emotional fear is different, it is based on our imaginings of what might happen

Question: are there no negative emotions that have a use?

It sounds strange, but it is very important to understand that all negative emotions are all absolutely useless; they do not serve any useful purpose; they do not make us acquainted with new things or bring us nearer to new things; they do not give us energy; they only waste energy and create unpleasant illusions."

Well, this post is not about negative emotions. Instead, it introduces spheres of negative radius. Which concept (thanks to Anna for letting me know) has been  discussed by another Russian thinker, Pavel Florensky.

"Referring to Dante's Divine Comedy, Florensky opposes Copernicus' heliocentric system. Interprets the Michelson-Morley experience as proof of the immobility of the Earth. Declares “the notorious Foucault's experience” fundamentally unproven. Commenting on Einstein's special theory of relativity, Florensky argues that beyond the limit of the speed of light begins non-physical “that light”. This otherworld of imaginary magnitudes provides a description of the ultimate eternal reality. Based on the geocentric system, Florensky calculates the distance to this world as the distance at which a body orbiting the Earth in one day would travel at the speed of light."

This otherworld of imaginary magnitudes provides a description of the ultimate eternal reality.


We denote vectors in R4 by bold letters. With m in S3R4, consider a sphere St(m) on S3 given by (cf. Part 4, Eq. (2))

St(m) = {xS3: x·m = cos(t)} .               (0)

For t in [0,2π), t∉{0,π},  we orient it by the normal vector
(cf. Part 4, Eq. (3))

n(x) = ( m-cos(t)x )/sin(t).               (1)

In the previous post we have derived the formulas for stereographic projection xyS3R3 (S3 is taken without its South Pole!), and our aim now is to identify the surface obtained by this projection from St(m). Let us recall the formulas for stereographic projection and its inverse (cf. Part 8, Eq. (2),(3a),(3b)). With i=1,2,3, we have:

yi = xi/(1+x0).            (2)

x0 = (1-y2)/(1+y2).               (3a)

xi = 2yi/(1+y2).                (3b)

Substituting (3a) and (3b) in (0), multiplying both sides by (1+y2), and collecting y2 terms,  we obtain:

y2(m0+cos(t)) - 2m'·ym0 - cos(t),                (4)

where

m '= m1e1+m2e2+m3e3.

There will be now two cases. The first case is when m0+cos(t) ≠ 0. When this happens?

Proposition 1. m0+cos(t)) = 0 if and only if the origin -e0 of the stereographic projection is on the sphere St(m).

Proof. Using Eq. (0), with m=m0e0+...+m3e3, we see that x=-e0 is on St(m) if and only if m0+cos(t) = 0.

Suppose now that  -e0 is not on the sphere. We divide both sides of (4) by (m0+cos(t)):

y2 - 2m'·y/(m0+cos(t))  = (m0 - cos(t))/(m0+cos(t)) .

We now add m'2/(m0+cos(t))2 to both sides to obtain:

(y-m'/(m0+cos(t)))2 = (m0 - cos(t))/(m0+cos(t)) + m'2/(m0+cos(t))2.                (5)
Since m2=1, we get m'2=1-(m0)2, and the right hand side of (5) simplifies to


(y-m'/(m0+cos(t)))2= sin(t)2/(m0+cos(t))2.                 (6)

From (6) we conclude that the image of St(m) is a sphere in R3 with the center at

c = m'/(m0+cos(t)),                (7)

and radius ρ

ρ = ε sin(t)/(m0+cos(t)),                 (8)

where ε = ±1.
For t>0, t<π, sin(t)>0, but what about the sign of (m0+cos(t))? It can be positive, but it can also be negative! This leads us to the concept of spheres with negative radius.

Oriented spheres in R3.

In the definition below the term "signed" means that ρ can have positive or negative value.

Definition 1. The oriented sphere with center p in R3 and signed radius ρ, 0≠ρ∈R, is

Sρ(p) = {yR3: (y-p)2=ρ2}                (9)

with unit normal vector field

n'(y) = (p-y)/ρ.                (10)

To understand what is going on, assume first that the center of the sphere is at the origin p=0. Then n'(y) = -y/ρ. Thus the normal n'(y) is inward for ρ>0, and outward for ρ<0. If p is arbitrary, not necessarily the origin of R3, the situation is the same, just translated by p.

Proposition 2. The image of the oriented sphere St(m)S3, not containing -e0, under the stereographic projection xy given by (3a,3b), is the oriented sphere Sρ(c)R3, where

c = m'/(m0+cos(t)),
ρ
= sin(t)/(m0+cos(t)).

Proof. We have already calculated the center c of the projected sphere. It remains to be checked that the direction of the image of the normal n(x) given by (1) is the same as the direction of n'(y(x)) = (c-y(x))/ρ. From (2),(7),(8),(10) we get

n'(y(x))i =mi/(ε sin(t)) - xi(m0+cos(t))/(ε sin(t) (1+x0)).      

Now we need to find the image, let us call it  n''(x) of the normal unit vector n(x), Eq. (1),  under the stereographic projection. While coordinates transform using the formulas (2): yi=yi(xμ), i=1,2,3; μ=0,1,2,3, tangent vectors transform using the matrix of partial derivatives (Why?), thus

n''(x)i = (∂yi/∂xμ) nμ(x).

The result of this calculation is (you can use any symbolic algebra or, nowadays, probably also  online AI, software to do it for you):

n''(x)i = (1+x0)-1(mi/sin(t) -xi(m0+cos(t))/(sin(t)(1+x0)) .        (11)

Since (1+x0) is positive (Why?), for n' and n'' to have the same direction we need to take ε=+1.
Notice that while stereographic projection preserves angles between tangent vectors, it does not preserve their lengths. That is why n'' is only proportional to n'. To become a unit normal it would have to be normalized. Then it should coincide with n'.

Exercise 1. Verify the last sentence.

Stereographic projection can indeed exchange the "inside" with the "outside". Here is an illustration for circles on the sphere, projected onto the plane:





Proposition 2 takes care of this problem.

The case when the origin -e0 of the stereographic projection, the ∞ point,  is on St(m), i.e. when 
m0+cos(t) = 0, will be discussed in the next note.


P.S. 23-04-25 19:21 I was busy writing another review






Biolocation

  On Tuesday, December 23, Vlad Zhigalov (see e.g. here ) had a talk at the " Temporology " seminar hosted at Omsk.  He spoke abo...