Applied Conformal &
Quasiconformal Geometry

Pranav Haridas
Kerala School of Mathematics

Part I: Conformal Mappings
Mathematical Foundations of Contemporary Scientific Research
July 27, 2026

Computed Tomography (CT)

X-Ray Attenuation
  • Motorized Source
    Rotates around the patient, shooting narrow beams of X-rays.
  • Attenuated Radiation Measurement
    Detectors measure reduced radiation emerging through tissue.
  • Tomographic Reconstruction
    Calculates the X-ray attenuation coefficient at every 3D point.
  • Output Discrete Scalar Field
    Outputs tissue density field $I(x,y,z)$.
  • Clinical Utility
    Well-suited for distinguishing dense structures like bone.
CT Scanner
PHYSICAL MEASURAND
$I_{\text{CT}}(x,y,z) = \text{Attenuation Coefficient}$

1. Motorized X-Ray Source

Rotating Gantry

Rotational Acquisition Geometry

  • High-voltage X-ray tube mounted inside a heavy motorized circular gantry ring.
  • Rotates continuously 360° around the patient at high speed (> 120 RPM).
  • Emits narrow, highly collimated fan-shaped or cone-shaped X-ray beams.
  • Irradiates thin cross-sectional anatomical slices along the longitudinal Z-axis.
Z-Axis Slice Acquisition
PHYSICAL STEP 1
Continuous 360° Fan-Beam Emission

2. Attenuated Radiation Measurement

Beer-Lambert Law

Detector Array & Energy Loss

Arc of solid-state electronic detectors directly opposite the X-ray tube measures reduced photon intensity emerging through tissue.

$$I(\theta, s) = I_0 \exp\left( -\int_{L_{\theta,s}} \mu(x,y,z) \, ds \right)$$

Logarithmic ratio $\ln(I_0 / I)$ yields the line integral of total attenuation along ray path $L_{\theta,s}$.

PHYSICAL STEP 2
Line Integrals $p(\theta,s) = \ln(I_0 / I)$

3. Tomographic Reconstruction

Radon Inversion

MATH 2D Radon Inversion

Reconstructs the continuous 2D slice density $\mu(x,y,z_k)$ at slice position $z_k$ from 1D ray projections $p(\theta,s)$ for $\theta \in [0,\pi)$:

$$\mu(x,y) = \int_{0}^{\pi} \left[ \mathcal{R}\{\mu\}(\theta,s) * h(s) \right]_{s = x\cos\theta + y\sin\theta} d\theta$$
Geometric Purpose
Provides the reconstructed attenuation field $\mu(x,y,z)$ across all slice positions along the $Z$-axis.

ALGORITHM Filtered Back Projection

Unfiltered back-projection creates $1/r$ spatial blurring. High-pass ramp filtering $h(s)$ (Riesz transform) cancels low-frequency artifacts, recovering sharp tissue boundaries.

Scope & Course Organization

Radon Inversion, the Fourier Slice Theorem, and iterative CT reconstruction form a dedicated full lecture on inverse problems. In this course, we treat FBP as an established upstream provider of the 3D scalar field.

4. Output Discrete Scalar Field $I(x,y,z)$

Voxel Quantization

Discrete Scalar Field Function $I(x,y,z)$

The continuous attenuation field $\mu(x,y,z)$ is sampled at discrete voxel lattice locations $(x_i, y_j, z_k) \in \mathbb{Z}^3$, defining the discrete scalar field mapping:

$$I \colon \mathbb{Z}^3 \to \mathbf{R}, \quad I(i,j,k) = \text{HU}(x_i, y_j, z_k)$$

Values are calibrated into Hounsfield Units (HU) relative to water ($\text{HU}=0$) and air ($\text{HU}=-1000$):

$$\text{HU}(x,y,z) = 1000 \times \frac{\mu(x,y,z) - \mu_{\text{water}}}{\mu_{\text{water}} - \mu_{\text{air}}}$$
DISCRETE 3D GRID
$I(i,j,k) = \text{Voxel Density}$
Bone: $+1000 \text{ to } +3000 \text{ HU}$
Water: $0 \text{ HU}$
Air: $-1000 \text{ HU}$

5. Clinical Utility: Bone & Dense Structures

Anatomical Segmentation

High X-Ray Contrast Differentiation

  • High atomic number elements (Calcium in bone) cause high X-ray attenuation (+1000 to +3000 HU).
  • Provides sharp isosurface threshold boundaries $c_{\text{bone}}$ for Marching Cubes surface extraction.
  • Ideal for orthopedic modeling, skull extraction, and lung parenchyma boundary segmentation.

Hounsfield Tissue Contrast

• Air: -1000 HU
• Soft Tissue: +30 to +45 HU
• Cortical Bone: +1000 to +3000 HU

Magnetic Resonance Imaging (MRI)

Proton Relaxation
  • Magnetic Field Alignment
    Superconducting magnet aligns hydrogen nuclear spins in tissue water.
  • RF Pulse Excitation
    Radiofrequency pulses tip proton spin vectors out of alignment.
  • Signal Emission Dynamics
    Protons emit faint radio signals as they relax back into thermal equilibrium.
  • Output Discrete Scalar Field
    Reconstructs amplitude into discrete 3D field $I(x,y,z)$.
  • Clinical Utility
    Exquisite soft-tissue contrast via water content and relaxation times.
MRI Scanner
PHYSICAL MEASURAND
$I_{\text{MRI}}(x,y,z) = \text{MR Signal Intensity}$

1. Magnetic Field Alignment ($B_0$)

Larmor Precession

Superconducting Magnet Alignment

High-field superconducting magnet produces static magnetic field $\vec{B}_0$ (1.5T or 3.0T) along the longitudinal $Z$-axis, aligning nuclear magnetic moments of hydrogen protons ($^1\text{H}$) in tissue water.

Protons precess around $\vec{B}_0$ at the characteristic Larmor precession frequency:

$$\omega_0 = \gamma B_0 \quad (\gamma = 42.58 \text{ MHz/T for } ^1\text{H})$$
MRI Magnet
PHYSICAL STEP 1
Larmor Frequency $\omega_0 = \gamma B_0$

2. RF Pulse Excitation

Resonant Energy Transfer

Transverse Magnetization Tipping

Radiofrequency (RF) coils transmit an oscillating magnetic field $B_1(t)$ tuned precisely to the Larmor frequency $\omega_0$. Protons absorb resonant energy, tipping net magnetization vector $\vec{M}$ into the transverse $XY$-plane.

$$\vec{B}_1(t) = B_1 \cos(\omega_0 t)\hat{i} - B_1 \sin(\omega_0 t)\hat{j}$$

Creates a coherent, spinning transverse magnetization vector $M_{xy}$ that induces measurable RF signals in receiver coils.

PHYSICAL STEP 2
Transverse Spin Vector $M_{xy}$

3. Signal Emission & Relaxation Dynamics

Bloch Equations

MATH Bloch Differential Equation

Governs trajectory of magnetization vector $\vec{M}$ returning to thermal equilibrium:

$$\frac{d\vec{M}}{dt} = \vec{M} \times \gamma\vec{B} - \frac{M_x \hat{i} + M_y \hat{j}}{T_2} - \frac{(M_z - M_0)\hat{k}}{T_1}$$
Physical Meaning
Decouples longitudinal energy recovery ($T_1$) from transverse spin dephasing ($T_2$), generating tissue contrast.

PHYSICS $T_1$ vs. $T_2$ Relaxation Decay

Longitudinal $T_1$ (Spin-Lattice):
$M_z(t) = M_0(1 - e^{-t/T_1})$

Transverse $T_2$ (Spin-Spin):
$M_{xy}(t) = M_{xy}(0) e^{-t/T_2}$

Contrast Origin

Varying water binding in tissues creates unique $(T_1, T_2)$ parameters for gray matter, white matter, and CSF.

4. Output Discrete Scalar Field $I(x,y,z)$

k-Space IFFT Reconstruction

Discrete Scalar Field Function $I(x,y,z)$

Gradient coils $(G_x, G_y, G_z)$ encode spatial positions into RF signal frequencies in $k$-space. 3D Inverse Fast Fourier Transform (IFFT) yields the discrete voxel grid mapping:

$$I \colon \mathbb{Z}^3 \to \mathbf{R}, \quad I(i,j,k) = \text{MR Signal Amplitude}(x_i, y_j, z_k)$$

Reconstructed from 3D spatial frequency domain ($k$-space) data $S(k_x, k_y, k_z)$:

$$I_{\text{MRI}}(x,y,z) = \iiint S(k_x, k_y, k_z) e^{i 2\pi (k_x x + k_y y + k_z z)} \, dk_x dk_y dk_z$$
DISCRETE 3D GRID
$I(i,j,k) = \text{MR Signal Intensity}$
High Signal: Bright Voxels
Low Signal: Dark Voxels
Sub-millimeter Spatial Resolution

5. Clinical Utility: Soft-Tissue Contrast

Neuroimaging Contrast

Exquisite Soft-Tissue Differentiation

  • Tissues vary drastically in proton density and relaxation constants ($T_1, T_2$).
  • Differentiates cerebral gray matter, white matter, and cerebrospinal fluid (CSF).
  • Essential input for Marching Cubes surface extraction of the genus-zero ($g=0$) cortical manifold.

Brain Tissue Relaxation Times

• White Matter: $T_1 \approx 800\text{ ms}$
• Gray Matter: $T_1 \approx 1300\text{ ms}$
• CSF: $T_1 \approx 4000\text{ ms}$

Discretization of the Scalar Field $I(x,y,z)$

Voxel Lattice Sampling

1. Continuous Domain $\Omega$ to Voxel Lattice $\mathbb{Z}^3$

Anatomy is modeled as a continuous scalar function $\mu \colon \Omega \to \mathbf{R}_{\ge 0}$ over a compact domain $\Omega \subset \mathbf{R}^3$. Scanners partition $\Omega$ into closed cuboid voxel cells $V_{i,j,k} \subset \Omega$:

$$V_{i,j,k} = \left[x_i \pm \tfrac{\Delta x}{2}\right] \times \left[y_j \pm \tfrac{\Delta y}{2}\right] \times \left[z_k \pm \tfrac{\Delta z}{2}\right]$$
$$\mathcal{I} \colon \mathbb{Z}^3 \to \mathbf{R}, \quad \mathcal{I}(i,j,k) = \frac{1}{\text{vol}(V_{i,j,k})} \int_{V_{i,j,k}} \mu(x,y,z) \, dV$$
Voxel Discretization over Human Body
DISCRETE LATTICE SAMPLING
Continuous $\Omega \to$ Grid $\mathcal{I}(i,j,k)$

Clinical Level Set Extraction

Preimage Theorem

Clinical Objective: Anatomical Manifolds

A diagnosing physician cannot interpret dense 3D scalar arrays directly. To isolate and visualize distinct anatomical structures—such as the colon wall, cerebral cortex, or tumor boundaries—we must extract explicit geometric surface representations called level sets. Because the true physical field $\mu$ is unobservable, we compute the level set $\widehat{\mathcal{M}}_c$ of the continuous interpolant $\tilde{\mu}$ as the surrogate proxy.

$$\text{Discrete Grid } \mathcal{I} \xrightarrow{\text{Interpolate}} \tilde{\mu} \approx \mu \xrightarrow{\text{Level Set}} \widehat{\mathcal{M}}_c = \tilde{\mu}^{-1}(c) \approx \mathcal{M}_c$$

Regular Value Theorem Approximation

If $c$ is a regular value ($\nabla \mu(p) \neq 0, \forall p \in \mu^{-1}(c)$), the Preimage Theorem mathematically guarantees $\mu^{-1}(c)$ is a smooth 2D manifold. We computationally analyze level sets of $\mathcal{I}$ because interpolation of the discrete grid values $\mathcal{I}(i,j,k)$ provides the only evaluable manifold approximation $\widehat{\mathcal{M}}_c \approx \mu^{-1}(c)$.

Geometric Distortion in Volumetric Projections

Overview

The Geometry Problem in Raw Volumetric Imaging

Anatomical boundaries (e.g. brain cortex, colon wall) are curved 2D Riemannian manifolds $\mathcal{M} \subset \mathbf{R}^3$. Direct 2D planar projections cause severe angular shear, distorting local shapes and intersection angles. Non-conformal projections stretch circular functional regions into elongated ellipses, corrupting diagnostic metrics.

$$\text{Conformal Condition: } \angle(u, v) = \angle(df(u), df(v)), \quad \forall u,v \in T_p \mathcal{M}$$
Planar Projection Shear
Planar Projection Shear (Orthogonal Grid to Sheared 2D Grid)

1. Isometric Mappings (Length Preserving)

Isometric

Definition & Constraints

An isometric mapping preserves all geodesic distances exactly. The metric tensor under mapping remains identical to the identity matrix:

$$g_{ij} = \delta_{ij}$$

Physical Limitation: Gauss's Theorem

By Gauss's Theorema Egregium, a curved surface can be mapped isometrically to a flat plane if and only if its Gaussian curvature is zero ($K = 0$).

Since anatomical structures (e.g. cortex, colon) possess non-zero Gaussian curvature ($K \neq 0$), isometric flattening is physically impossible without introducing cuts or tears.

2. Conformal Mappings (Angle Preserving)

Conformal

Isotropic Metric Scaling

Conformal maps preserve local angles and map infinitesimal circles to circles. The metric tensor scales isotropically by a conformal factor $e^{2\lambda}$:

$$g_{ij} = e^{2\lambda} \delta_{ij}$$

Anatomical Suitability

Unlike isometric maps, conformal mappings accommodate arbitrary curvature without introducing artificial cuts or tears. By preserving local inner angles, they prevent shape distortion in critical areas, making them ideal for mapping cortical sulci or colon folds.

Conformal vs Non-Conformal Comparison
Conformal (Circle to Circle) vs. Non-Conformal (Circle to Ellipse)

3. Authalic Mappings (Area Preserving)

Authalic

Area Element Preservation

Authalic mappings preserve local area elements exactly across the domain. The determinant of the metric tensor remains constant:

$$\det(g_{ij}) = 1$$

Clinical Drawback

While authalic mappings are essential for counting cell densities or computing integrals over surfaces, they introduce severe angular shear.

This shear distorts local shape profiles (e.g. circles to ellipses), rendering them poor for shape-based visual diagnosis.

Isosurface Extraction

Level-Set Surfaces

Level Set Definition

Extract a 2D surface $\mathcal{M}$ from the volume by defining an implicit level set:

$$\mathcal{M} = \left\{(x,y,z) \in \mathbf{R}^3 \;\middle|\; \tilde{\mu}(x,y,z) = c \right\}$$

where $c$ is a designated threshold separating distinct tissue types (e.g. bone vs. soft tissue) and $\tilde{\mu}$ is the continuous interpolant proxy of the discrete voxel intensities.

Isosurface Extraction
2D Manifold Extracted from 3D Voxel Grid

Marching Cubes Algorithm

Triangulation

Discrete Surface Triangulation

  • Evaluates density values at 8 corners of cubic voxel cells.
  • Linearly interpolates grid edge crossings to compute 3D vertex coordinates.
  • Connects vertices with edges to bound flat triangular faces.
  • Forms a closed, watertight 2D manifold $\mathcal{M} \subset \mathbf{R}^3$.
Marching Cubes Configurations
Simplicial Complex $\mathcal{M} \subset \mathbf{R}^3$

Brain Imaging & Cortical Surface Mapping

Genus Zero Manifold
  • Anatomical & Topological Model
    The 2D physical interface separating gray matter from white matter (or CSF) forms a closed, orientable Riemannian manifold of genus zero ($g=0$) with no topological boundary ($\partial \mathcal{M} = \emptyset$).
  • Uniformization Theorem
    Guarantees a conformal diffeomorphism $f \colon \mathcal{M} \to \mathbb{S}^2$.
  • Pull-Back Metric Scaling
    $g_{\mathcal{M}} = e^{2\lambda(z)} g_{\mathbb{S}^2}$, where $e^{2\lambda(z)} > 0$ is the conformal factor dictating local area distortion.
  • Necessity of Conformality (Isotropy)
    Preserves local angles & geometry. Circular functional regions remain circular on the sphere; sulcal intersection angles are preserved.
Human Cerebral Cortex
PULL-BACK METRIC
$g_{\mathcal{M}} = e^{2\lambda(z)} g_{\mathbb{S}^2}$

Virtual Colonoscopy & Flattening

Haustral Fold Flattening
  • Tubular Surface Topology:
    Inner colon surface is a cylinder with two boundaries ($g=0, b=2$), featuring sharp haustral folds.
  • Conformal Unrolling:
    Surface is sliced and conformally flattened into a 2D rectangle or annulus for polyp inspection.
  • Gauss's Theorema Egregium:
    Intrinsic curvature makes metric distortion mathematically unavoidable.
  • Preserving Polyps for Diagnosis:
    Polyps are precursors to tumors and appear as small spherical structures. Area-preserving maps introduce intense shear distortion, making them undetectable. Conformal mapping preserves their local round shape, ensuring they are caught by radiologists.
Virtual Colonoscopy 3D Mesh
TOPOLOGICAL DOMAIN
$g=0, b=2 \implies \text{2D Rectangle / Annulus}$

Harmonic Maps & Discrete Cotangent Laplacian

Finite Element Discretization

Dirichlet Energy & Euler-Lagrange

Holomorphic functions $f = u+iv$ are harmonic ($\Delta u = 0$). Critical points of Dirichlet energy $E_D(f) = \frac{1}{2} \int_{\mathcal{M}} |\nabla f|^2 dA$ yield the Laplace-Beltrami PDE $\Delta_{\mathcal{M}} f = 0$.

Discrete Cotangent Formula

Evaluated over vertex $i$'s 1-ring neighborhood $\mathcal{N}(i)$ on a triangular mesh:

$$(\Delta u)_i = \frac{1}{2 A_i} \sum_{j \in \mathcal{N}(i)} (\cot \alpha_{ij} + \cot \beta_{ij}) (u_j - u_i)$$

Pinning boundary vertices yields sparse linear system $Lu = 0$.

COTANGENT LAPLACIAN
$$Lu = 0$$
Piecewise linear Finite Element Method discretization

FEM Derivation of Cotangent Weights

Variational Geometry (Part 1)

1. Piecewise Linear Hat Functions

For a triangular face $T = (i,j,k)$, any scalar function $u$ is interpolated using hat basis functions $\phi_i$:

$$u(x,y) = \sum_{i} u_i \phi_i(x,y)$$

The gradient $\nabla \phi_i$ on triangle $T$ is constant and perpendicular to the opposite edge $jk$.

FEM Derivation of Cotangent Weights

Variational Geometry (Part 2)

2. Inner Product of Gradients

Integrating the inner product of basis gradients over triangle $T$ yields the cotangent of the interior angle $\gamma_k$:

$$\int_{T} \nabla \phi_i \cdot \nabla \phi_j \, dA = -\frac{1}{2} \cot \gamma_k$$

Summing over the two triangles sharing edge $ij$ yields the cotangent weight $w_{ij} = \frac{1}{2}(\cot \alpha_{ij} + \cot \beta_{ij})$.

Model 2: Cartography

Map Projections

Gauss's Theorema Egregium

The Earth is modeled as a sphere $\mathbb{S}^2$ with constant Gaussian curvature $K=1$. The flat map lies in the Euclidean plane $\mathbf{R}^2$ with constant Gaussian curvature $K=0$.

Isometric (distance-preserving) mappings between surfaces of unequal Gaussian curvature are mathematically impossible.

Therefore, every map projection must distort distance, area, or angle. For maritime navigation, preserving local angles (course headings) is paramount, naturally requiring mappings to be strictly conformal. A different angle will take you nowhere!

Maritime Navigation Angle

Stereographic Projection

Conformal Mapping to the Plane

Maps a point $P = (X, Y, Z)$ on the Riemann sphere $\mathbb{S}^2$ conformally to $z = x + iy$ on the complex plane $\mathbb{C}$. By intersecting the plane exactly at the equator, the lower hemisphere maps to the interior of the unit disk $\mathbf{D} = \{z \in \mathbb{C} : |z| < 1\}$.

$$z = \frac{X + iY}{1 - Z}$$

The projection preserves angles but distorts distances, quantified by the scalar conformal factor: $ds^2_{\mathbb{S}^2} = \lambda(x, y) ds^2_{\mathbf{R}^2}$.

Stereographic Projection

Characterizing Conformal Maps

Mathematical Agenda

The Mathematical Objective

We have seen that conformal mappings naturally arise in both Medical Imaging and Cartography. By the Uniformization Theorem, any simply connected Riemann surface is conformally equivalent to either the Riemann Sphere $\hat{\mathbb{C}}$, the Complex Plane $\mathbb{C}$, or the Upper Half Plane $\mathbb{H}$.

For the remainder of this lecture, we will explicitly characterize the algebraic structure of all conformal automorphisms mapping these canonical domains to themselves.

Conformal Maps as Biholomorphisms

Complex Analysis

In the context of the canonical domains ($\hat{\mathbb{C}}, \mathbb{C}, \mathbb{H}$), the geometric property of conformality is inextricably linked to complex analyticity.

Fundamental Equivalence

Any angle-preserving (conformal) mapping between these spaces is necessarily a holomorphic function with a non-zero derivative.

Therefore, characterizing the conformal automorphisms of these spaces is mathematically identical to finding all their biholomorphisms (bijective holomorphic mappings).

The Algebraic Proof of Conformality

Part 1: Vector Calculus

Let $f: \mathbf{R}^2 \to \mathbf{R}^2$ be a differentiable mapping given by $f(x,y) = (u(x,y), v(x,y))$. Geometrically, preserving angles and orientation requires the Jacobian matrix $J_f$ to act purely as a scaled rotation at every point.

The Jacobian Matrix

$$ J_f = \begin{pmatrix} u_x & u_y \\ v_x & v_y \end{pmatrix} = \lambda \begin{pmatrix} \cos \theta & -\sin \theta \\ \sin \theta & \cos \theta \end{pmatrix} $$
where $\lambda > 0$ strictly preserves orientation and prevents degeneracy.

Reduction to Complex Analyticity

Part 2: Algebraic Equivalence

Substituting a single matrix $\begin{pmatrix} a & -b \\ b & a \end{pmatrix}$ for the scaled rotation yields a direct algebraic equivalence.

The Cauchy-Riemann Equations

Equating the components of $J_f = \begin{pmatrix} a & -b \\ b & a \end{pmatrix}$ instantly recovers:
$$ u_x = v_y \quad \text{and} \quad u_y = -v_x $$

Because $u$ and $v$ satisfy the Cauchy-Riemann equations, the function $f(z) = u + iv$ is strictly holomorphic. The scaling factor $\lambda > 0$ ensures $f'(z) \neq 0$.

The Projective Space $\mathbf{P}^1$

Automorphisms

The conformal automorphisms of the projective space $\mathbf{P}^1 = \mathbf{C} \cup \{\infty\}$ are exactly the fractional linear transformations, also known as Möbius transformations.

Theorem

The group $\operatorname{Aut}(\mathbf{P}^1)$ consists precisely of all rational functions of the form:
$$f(z) = \frac{az + b}{cz + d}$$
where $a,b,c,d \in \mathbf{C}$ and $ad - bc \neq 0$.

This group is isomorphic to the projective special linear group $PSL(2, \mathbf{C})$. The condition $ad - bc \neq 0$ ensures the mapping is invertible and non-degenerate. Möbius transformations are uniquely determined by their action on any three distinct points in $\mathbf{P}^1$.

Proof of Characterization

Riemann Sphere

Any conformal automorphism of $\mathbf{P}^1$ is a meromorphic function globally defined on the extended complex plane. Because the mapping is bijective, it can possess exactly one pole, and this pole must be simple.

Formal Proof

Let $f$ have a simple pole at $z_0 \neq \infty$. Its Laurent expansion is $f(z) = \frac{c}{z-z_0} + g(z)$, where $g(z)$ is entire because $f$ has no other poles. Since $f$ is bijective, $z_0$ is the unique pole, meaning $f$ has no pole at $\infty$. Consequently, $g(z)$ must be bounded as $z \to \infty$. By Liouville's Theorem, $g(z)$ is a constant $k$.

Thus, $f(z) = \frac{c}{z-z_0} + k = \frac{kz + (c - kz_0)}{z-z_0}$, which is exactly a Möbius transformation. If $z_0 = \infty$, the principal part is $cz$, and $f(z) - cz$ is an entire bounded function, thus a constant $d$, yielding $f(z) = cz + d$.

The Complex Plane $\mathbf{C}$

Automorphisms

The conformal automorphisms of the complex plane $\mathbf{C}$ are the affine transformations.

Theorem

The group $\operatorname{Aut}(\mathbf{C})$ consists precisely of all functions of the form:
$$f(z) = az + b$$
where $a, b \in \mathbf{C}$ and $a \neq 0$.

This group is a subgroup of the full Möbius group. It contains all Möbius transformations that strictly fix the point at infinity. Geometrically, these mappings represent translations, rotations, and dilations.

Proof of Characterization

Complex Plane

Let $f \in \operatorname{Aut}(\mathbf{C})$. To preserve bijectivity, $f$ must be an entire function with a well-defined behavior at infinity. If $\infty$ were a removable singularity, $f$ would be bounded and thus constant by Liouville's theorem. If $\infty$ were an essential singularity, the Casorati-Weierstrass theorem dictates the image of $\{z : |z| > R\}$ is dense in $\mathbf{C}$. By the Open Mapping Theorem, the image of $\{z : |z| < R\}$ is a non-empty open set. Their intersection is therefore non-empty, violating injectivity.

Formal Proof

Therefore, the isolated singularity at $\infty$ is strictly a pole. Because $f(\infty) = \infty$, the map extends to an automorphism of the projective space $\mathbf{P}^1$.

By the previous theorem, every such automorphism is a Möbius transformation $f(z) = \frac{az+b}{cz+d}$. For this transformation to map $\infty \mapsto \infty$, we must strictly have $c = 0$, immediately reducing the function to the affine form $f(z) = az + b$.

The Upper Half Plane $\mathbf{H}$

Automorphisms

The conformal automorphisms of the upper half plane $\mathbf{H} = \{z \in \mathbf{C} : \operatorname{Im}(z) > 0\}$ correspond naturally to the real projective transformations.

Theorem

The group $\operatorname{Aut}(\mathbf{H})$ consists precisely of all functions of the form:
$$f(z) = \frac{az + b}{cz + d}$$
where $a, b, c, d \in \mathbf{R}$ and $ad - bc > 0$.

By scaling the coefficients such that $ad - bc = 1$, we obtain an explicit isomorphism to the projective special linear group $PSL(2, \mathbf{R})$.

Proof of Characterization

Upper Half Plane

To rigorously prove that every conformal automorphism takes this form, we must show that any $f \in \operatorname{Aut}(\mathbf{H})$ is necessarily a fractional linear transformation. We establish this by extending the domain of $f$.

Formal Proof

Let $f \colon \mathbf{H} \to \mathbf{H}$ be a biholomorphism. Because $f$ maps the real axis $\mathbf{R}$ bijectively to itself, the Schwarz Reflection Principle allows us to analytically continue $f$ across $\mathbf{R}$ into the lower half plane. This symmetric extension produces a bijective meromorphic function on the entire projective space $\mathbf{P}^1$.

By our previous characterization, any such global automorphism must be a Möbius transformation $f(z) = \frac{az+b}{cz+d}$. For this transformation to preserve the extended real axis $\hat{\mathbf{R}}$, its coefficients must be real (up to a global complex scalar), yielding $a,b,c,d \in \mathbf{R}$.

Finally, evaluating the imaginary part yields $\operatorname{Im}(f(z)) = \frac{(ad-bc)\operatorname{Im}(z)}{|cz+d|^2}$. For $f$ to map the upper half plane onto itself (rather than the lower half plane), we must strictly have $ad - bc > 0$.

Isomorphism to the Unit Disk

Cayley Transform

The upper half plane $\mathbf{H}$ and the open unit disk $\mathbf{D}$ are conformally equivalent domains. This equivalence establishes a strict group isomorphism between their respective automorphism groups.

The Cayley Transform

The explicit biholomorphism mapping $\mathbf{H} \to \mathbf{D}$ is the Cayley transform:
$$W(z) = \frac{z - i}{z + i}$$
Consequently, $\operatorname{Aut}(\mathbf{D})$ is completely characterized by conjugating elements of $PSL(2, \mathbf{R})$ via this specific transformation.

3-Point Rigidity

Properties of Möbius Transformations

A Möbius transformation is uniquely determined by its action on three distinct points.

Theorem

Given three distinct points $(z_1, z_2, z_3)$ in $\mathbf{P}^1$ and three distinct target points $(w_1, w_2, w_3)$ in $\mathbf{P}^1$, there exists a unique Möbius transformation $f$ such that $f(z_k) = w_k$ for $k=1,2,3$.

Proof of 3-Point Rigidity

Proof

First, construct a transformation $T$ mapping $(z_1, z_2, z_3)$ to $(0, 1, \infty)$. The explicit formula is:

$$T(z) = \frac{(z - z_1)(z_2 - z_3)}{(z - z_3)(z_2 - z_1)}$$

Similarly, construct $S$ mapping $(w_1, w_2, w_3)$ to $(0, 1, \infty)$. The composition $f = S^{-1} \circ T$ maps $z_k \mapsto w_k$.

For uniqueness, assume $g$ also maps $z_k \mapsto w_k$. The composition $h = S \circ g \circ T^{-1}$ is a Möbius transformation $h(z) = \frac{az+b}{cz+d}$ fixing $0, 1,$ and $\infty$. Since $h(\infty) = \infty$, $c = 0$, so $h(z) = \frac{a}{d}z + \frac{b}{d}$. Since $h(0) = 0$, $b = 0$, yielding $h(z) = \frac{a}{d}z$. Finally, $h(1) = 1$ implies $a = d$. Thus, $h(z) = z$, the identity map. This strictly proves $g = S^{-1} \circ T = f$.

Application: Anatomical Registration

Medical Imaging

Landmark Pinning on the Cortex

When flattening a 3D cortical surface onto a sphere, the mapping retains the three complex degrees of freedom of $PSL(2, \mathbf{C})$.

By fixing exactly three distinct anatomical landmarks on the cortex, the 3-point rigidity theorem ensures the resulting Möbius transformation is uniquely determined. This eliminates arbitrary rotational freedom.

Cortical Surface

Cross-Ratio Preservation

Properties of Möbius Transformations

The cross-ratio of four distinct points $(z_1, z_2, z_3, z_4)$ in $\mathbf{P}^1$ is defined as:

$$(z_1, z_2; z_3, z_4) = \frac{(z_1 - z_3)(z_2 - z_4)}{(z_1 - z_4)(z_2 - z_3)}$$

Theorem

For any Möbius transformation $f$ and any four distinct points $z_k \in \mathbf{P}^1$, the cross-ratio is invariant: $(f(z_1), f(z_2); f(z_3), f(z_4)) = (z_1, z_2; z_3, z_4)$.

Proof of Invariance

Proof

Define $T(z) = (z, z_2; z_3, z_4)$. By construction, $T$ is the unique Möbius transformation mapping $(z_2, z_3, z_4)$ to $(1, 0, \infty)$.

Let $f$ be an arbitrary Möbius transformation. The composition $T \circ f^{-1}$ is a Möbius transformation mapping $(f(z_2), f(z_3), f(z_4))$ to $(1, 0, \infty)$.

By the 3-point rigidity theorem, $T \circ f^{-1}$ is the unique map sending those three points to $(1, 0, \infty)$. Evaluating this map at $f(z_1)$ computes the cross-ratio of the images: $(T \circ f^{-1})(f(z_1)) = T(z_1) = (z_1, z_2; z_3, z_4)$. Thus, the cross-ratio is preserved.

Application: Fluoroscopic Invariants

Medical Imaging

X-Ray Pose Estimation

In 2D X-ray fluoroscopy, the imaging process acts as a projective transformation on the 2D detector plane.

By placing four collinear radio-opaque markers on a surgical instrument, their cross-ratio remains strictly invariant on the resulting radiograph. This invariant allows the tracking system to identify the instrument independently of the C-arm camera angle.

Medical Scanner

Generalized Circles

Properties of Möbius Transformations

A generalized circle in $\mathbf{P}^1$ is either a standard circle or a straight line (a circle passing through the point at infinity).

Theorem

The action of $PSL(2, \mathbf{C})$ on $\mathbf{P}^1$ maps the set of all generalized circles bijectively onto itself.

Proof of Circle Preservation

Proof

The general equation for a generalized circle in $\mathbf{C}$ is $A|z|^2 + B z + \bar{B} \bar{z} + C = 0$, where $A, C \in \mathbf{R}$ and $|B|^2 > AC$.

Since any Möbius transformation can be decomposed into translations $z \mapsto z + a$, dilations/rotations $z \mapsto k z$, and inversion $z \mapsto 1/z$, it suffices to check invariance under these generators.

Translations and dilations trivially map circles to circles and lines to lines. For inversion $w = 1/z$, substitute $z = 1/w$ into the equation and multiply by $|w|^2$ to obtain $C|w|^2 + \bar{B} w + B \bar{w} + A = 0$. This preserves the form, yielding a generalized circle.

Application: Stereographic Projection

Cartography

Mapping the Sphere to the Plane

The stereographic projection is exactly a Möbius transformation mapping $\mathbf{S}^2$ to the extended complex plane $\mathbf{P}^1$.

Due to the generalized circle-preserving property, geometric circles on the Earth's surface (such as lines of latitude) map exclusively to circles or straight lines on the flat cartographic projection.

Stereographic Projection