Measuring Shape

Measuring Shape

Generative Radial Geometry (GRG): measuring shape to change perspective

This page presents concise formulas and proofs showing how GRG replaces constants like $\pi$ with shape measures (circles, hemispheres, ellipses, polygons) and with unit-sum quadratures. The goal is twofold: measurement precision and operational advantage (stability, variance/error reduction, traceability), even in asymmetric scenarios.

1) Principle

In GRG, many occurrences of $2\pi$, $4\pi$ or $\pi$ are not “mysteries” but the area/volume of a shape in the natural domain (e.g., directions, frequencies). By changing the measure on that domain (circular units: turn, hemi, sphere) the constant is absorbed and the integral becomes a normalized average.

2) Core theorems

Theorem 1 — Circular measure equivalence

Let $S^1$ be the unit circle of directions with arc-length measure $\lambda$ so that $\lambda(S^1)=2\pi$. Define the normalized measure $\mu=\lambda/(2\pi)$ with $\mu(S^1)=1$. Then for all $g\in L^1(S^1)$:

$$ \int_{S^1} g \, d\lambda \;=\; 2\pi \int_{S^1} g \, d\mu . $$

Proof. By definition of $\mu$. For $S^2$ the same holds with $\lambda(S^2)=4\pi$ and for the hemisphere with $2\pi$.

Theorem 2 — Anisotropic measure (no symmetry required)

Let $\rho(\theta)\ge 0$ with $\int \rho(\theta)\,d\theta = 1$ be an angular density (possibly non-uniform). Define the weighted measure $d\nu(\theta):=\rho(\theta)\,d\theta$. For all $g\in L^1$:

$$ \int g(\theta)\,\rho(\theta)\,d\theta \;=\; \mathbb{E}_{\theta\sim \rho}\!\big[g(\theta)\big]. $$

Interpretation. The “shape” of the domain is deformed according to $\rho$. Precision comes from encoding asymmetry in the measure, not from any pre-existing symmetry.

Theorem 3 — Radial quadratures with unit-sum weights

There exist sets of directions $\{\hat u_k\}_{k=1}^{N}\subset S^2$ and weights $w_k>0$ with $\sum_k w_k = 1$ such that, for every $g$ in a target family (e.g., polynomials in $\mu=\cos\theta$ up to degree $t$, or spherical harmonics up to degree $t$):

$$ \int_{S^2} g(\hat u)\, d\Omega \;=\; \sum_{k=1}^{N} w_k\, g(\hat u_k). $$

Note. Such cubatures (e.g., spherical t-designs or refined symmetric polyhedra) give controllable accuracy with $N=O(t^2)$.

3) Operational propositions: “removing $\pi$”

Proposition A — Fourier pair in turns (Hz), no $2\pi$

With frequency in cycles/s (Hz) and phase in turns, the Fourier pair

$$ X(f)=\int_{-\infty}^{\infty} x(t)\,e^{-i\,2 f t}\,dt , $$
$$ x(t)=\int_{-\infty}^{\infty} X(f)\,e^{+i\,2 f t}\,df . $$

is equivalent to the classical pair with $\omega=2\pi f$. Proof. Substitute $\omega=2\pi f$ and $d\omega=2\pi\,df$.

Proposition B — RC integral: rad/s vs Hz

$$ \int_{0}^{\infty} \frac{d\omega}{1+(\omega\tau)^2} \;=\; \frac{\pi}{2\,\tau} $$
$$ \Longleftrightarrow\quad \int_{0}^{\infty} \frac{df}{1+(2\pi f \tau)^2} \;=\; \frac{1}{4\,\tau} . $$

Proof. Use $\omega=2\pi f$. The constant $\pi$ is a unit-conversion factor.

4) Demonstrative examples (non–“pure geometry”, explicit shapes)

4.1 Radiometry (Lambert): the number is the area of a disk

For a perfectly diffuse surface,

$$ E \;=\; L \int_{\text{hemisphere}} \cos\theta \, d\Omega . $$

In spherical coordinates $d\Omega=\sin\theta\, d\theta\, d\varphi$, $\varphi\in[0,2\pi)$, $\theta\in[0,\tfrac{\pi}{2}]$. Then

$$ \int_{\text{hemisphere}} \cos\theta \, d\Omega \;=\; \int_0^{2\pi}\! d\varphi \int_0^{\pi/2}\! \cos\theta \sin\theta\, d\theta \;=\; 2\pi \cdot \tfrac{1}{2} \;=\; \pi . $$

Hence $E=\pi L$. With the “hemi” measure (hemisphere $=1$) the average is $\langle \cos\theta\rangle_{\text{hemi}} = \tfrac{1}{2}$ and, expressed per-hemi, $E_{\text{hemi}} = L \cdot \tfrac{1}{2}$. With the projected “Lambert” measure $d\Sigma = \tfrac{\cos\theta}{\pi}\, d\Omega$ (unit disk area), $\int d\Sigma=1$ and therefore $E = L$. In all cases the constant is a shape measure.

4.2 Signals/Control (RC): angle as shape

With the substitution $\omega\tau = \tan\theta$, the integral

$$ \int_0^\infty \frac{d\omega}{1+(\omega\tau)^2} $$

becomes a measure of an angle on the unit circle (arc from $0$ to $\tfrac{\pi}{2}$), which explains $\tfrac{\pi}{2}$. In “turn/Hz” units that angle is $\tfrac{1}{4}$, giving $\tfrac{1}{4\tau}$. Same content, different shape and measure.

4.3 Statistics (Cauchy as a projection of the circle)

If $\Theta \sim \operatorname{Unif}(-\tfrac{\pi}{2},\tfrac{\pi}{2})$ and $X=\gamma \tan\Theta$, then

$$ f_X(x) \;=\; \frac{1}{\pi \gamma}\,\frac{1}{1+(x/\gamma)^2} . $$

Interpretation: $X$ is the tangent projection of a point uniformly distributed on a circle. In “turn” units, let $\Phi \sim \operatorname{Unif}(-\tfrac{1}{4},\tfrac{1}{4})$ and $\Theta = 2\pi \Phi$; with respect to the turn-normalized angular measure (total length $=1$) the normalization constant is absorbed. Converting back to radians restores $1/(\pi\gamma)$.

5) Asymmetry: elliptical metrics and weighted measures

For anisotropic phenomena introduce an elliptical metric $M\succ 0$ and define $r_M(x)=\sqrt{x^\top M x}$. The “balls” of $r_M$ are ellipses, and many constants are absorbed by switching to $M$-radial coordinates. Any directional density $\rho(\theta)$ defines $d\nu=\rho\,d\theta$: integration becomes a unit-sum average without invoking symmetry.

6) Precise discrete re-determination (quadratures)

Design sets $\{(\hat u_k,w_k)\}$ with $\sum_k w_k=1$ that integrate exactly a useful family (e.g., $\{1,\cos\theta,\sin\theta,\cos2\theta,\sin2\theta\}$ for low-order anisotropic BRDFs, or spherical harmonics up to degree $t$). This yields a precise discrete scheme: constants like “$4\pi$” or “$2\pi$” become simply “sum of weights = 1”.

7) Practical conversions (GRG ↔ classical form)

Domain Classical form GRG form Observation
Radiometry (Lambert) $E=\pi L$ $E=L$ with unit projected measure $\pi$ = area of the projected disk
Signals (RC) $\displaystyle \int \frac{d\omega}{1+(\omega\tau)^2}=\frac{\pi}{2\tau}$ $\displaystyle \int \frac{df}{1+(2\pi f\tau)^2}=\frac{1}{4\tau}$ $\pi$ = rad↔turn conversion
Statistics (Cauchy) $\frac{1}{\pi\gamma}\frac{1}{1+(x/\gamma)^2}$ $\frac{1}{\gamma}\frac{1}{1+(x/\gamma)^2}$ with angle in turns $\pi$ = length of the angular domain
Fields/Fluxes $4\pi$ factors on $S^2$ Average over sphere with $\sum w_k=1$ $4\pi$ = area of the unit sphere

8) Measurable benefits

  • Numerical stability: fewer fixed large/small factors, better conditioning.
  • Variance reduction: with anisotropic $\rho$, MC error shrinks where it matters.
  • Traceability: results come pre-normalized (averages on 1); conversion back to standards is multiplication by $2\pi$ or $4\pi$.

9) Validation protocol (A/B)

  1. Select 3 cases (Lambert, RC, Cauchy/Maxwell).
  2. Define the GRG measure (turn/hemi/sphere or weighted $\rho$) and the quadratures (degree $t$).
  3. Compare to the classical baseline at equal budget (samples/time).
  4. Metrics: relative error, variance, conservation (energy/flux), stability along pipelines.
  5. Success threshold: ≥50% lower cost at equal error or ≤½ the error at equal cost in ≥2 out of 3 cases.

10) Conclusion

“Removing $\pi$” is not rhetoric: it is the outcome of measuring the domain in the right shape. GRG provides a measurement language (circular units and unit-sum quadratures) that yields cleaner formulas and, crucially, quantifiable benefits. This is a change of perspective that delivers progress: measurable improvement, not stubbornness.

Inside-Out GRG

Inside-Out GRG

GRG — Looking from the Inside: a discrete, π-free language for shapes

Generative Radial Geometry (GRG) — Looking from the Inside

GRG conventions. Angles are measured in turns \(u\in[0,1)\). Radial computations (perimeter and radial area) run on the full turn; the Cartesian area uses the quarter (and then ×4). Trigonometric forms are optional cross-checks, not inputs to the pipeline.

Work-in-progress. GRG is a discrete, π-free way to measure and compare shapes by radial generation. The turning point was this: I looked at the circle from its inside—a point that sends rays outward and measures what they meet. That inside→out switch is small to say and large to practice. The door is open; we will walk through it.

Why GRG matters

  • Discrete by nature: everything runs in finite steps (N)—heights, chords, bands.
  • π-free computation: the pipeline never needs π to run; π appears only when we compare to classical values.
  • Immediate comparability: two shapes at the same resolution are directly comparable with the same grammar.
  • Multi-scale: increase N; the very same procedure refines itself (coarse → fine).
  • GPU-friendly: short loops and sums, ideal for digital geometry, CAD/graphics, and simulation.

Looking from the inside

Not from a grid of squares, not from external coordinates—from a center that emits rays. A shape becomes a radial profile (a finite vector of values) rather than a formula to solve. From there, measurement and comparison happen discretely and immediately.

The joined method (circle)

Perimeter (radial, full turn): split the full turn \(u\in[0,1)\) into \(4N\) equal steps \(\Delta u=1/(4N)\). Join successive GRG border nodes by chords; the quarter view in the figure is for clarity.

Area (Cartesian, quarter): split one radius into \(N\) parts; at each split, drop the perpendicular to the circle (the “height”). Adjacent heights form a trapezoid; summing them tiles the quarter of a disk (the last one degenerates to a triangle). Multiply by 4 for the full disk.

Figure — joined method (circle)

Blue = angular sectors (perimeter, shown on the quarter). Red = single-radius trapezoids (area on the quarter). Dashed red rays start from the height intersections and generally do not coincide with the blue rays.

Figure 1. Quarter circle built from angular sectors (blue) and single-radius trapezoids (red); dashed rays start at height intersections.

Discrete comparability (the practical power)

At the same resolution N, any two shapes share the same grammar (heights, chords). This makes them comparable immediately—π-free inside the computation. We use simple, discrete indices relative to a baseline circle at the same N:

  • Area Index \(a_N\)
  • Perimeter Index \(c_N\)
  • Isoperimetric-like ratio \(q_N=a_N/c_N^2\) (the circle maps to 1 at the same N)
Minimal math — turns-first (π-free, node-based)

Quarter nodes (turns):

\[ \Delta u \;=\; \frac{1}{4N} \] \[ u_k \;=\; k\,\Delta u \quad (k=0,\dots,N) \]

Quarter inscribed perimeter (node distances):

\[ C_{\mathrm{in}}^{(Q)}(N) \;=\; \sum_{k=0}^{N-1} \bigl\|\,P_{k+1}-P_k\,\bigr\| \]

Full inscribed perimeter:

\[ C_{\mathrm{in}}(N) \;=\; 4\,C_{\mathrm{in}}^{(Q)}(N) \]

Radial area (full turn, shoelace on 4N nodes):

\[ A_{\mathrm{in}}(N) \;=\; \frac{1}{2} \sum_{j=0}^{4N-1} \bigl(x_j\,y_{j+1}-x_{j+1}\,y_j\bigr) \]

Cartesian area (quarter, Δx) → ×4:

\[ \Delta x \;=\; \frac{R}{N} \] \[ x_i \;=\; i\,\Delta x \] \[ y_i \;=\; \sqrt{\,R^2 - x_i^{\,2}\,} \] \[ A_{\mathrm{trap}}^{(Q)}(N) \;=\; \sum_{i=1}^{N} \frac{y_{i-1}+y_i}{2}\;\Delta x \] \[ A_{\mathrm{trap}}(N) \;=\; 4\,A_{\mathrm{trap}}^{(Q)}(N) \]

Discrete indices (baseline = circle at the same N):

\[ a_N \;=\; \frac{A_{\text{shape}}(N)}{A_{\text{circle}}(N)} \] \[ c_N \;=\; \frac{C_{\text{shape}}(N)}{C_{\text{circle}}(N)} \] \[ q_N \;=\; \frac{a_N}{c_N^{\,2}} \]
Classical cross-check (optional)

Inscribed perimeter (quarter step in turns):

\[ \Delta u \;=\; \frac{1}{4N} \] \[ s_N \;=\; 2R\,\sin\!\bigl(\pi\,\Delta u\bigr) \] \[ C_{\mathrm{in}}(N) \;=\; 4\,N\,s_N \;=\; 8\,N\,R\,\sin\!\Bigl(\frac{\pi}{4N}\Bigr) \]

Circumscribed perimeter (outer bound):

\[ C_{\mathrm{out}}(N) \;=\; 8\,N\,R\,\tan\!\Bigl(\frac{\pi}{4N}\Bigr) \]

Radial area via triangles (full turn):

\[ A_{\mathrm{tri}}(N) \;=\; 2\,N\,R^{2}\, \sin\!\Bigl(\frac{\pi}{2N}\Bigr) \]

Cartesian area (quarter, Δx) — same as above, then ×4.

Figure — 3D comparator (hemisphere bands)

Clean comparator for surfaces: baseline hemisphere on the left, target surface on the right. Blue = GRG parallels on each latitude; red = bands (average of neighboring circumferences). Below, π-free indices (schematic bars).

Figure 2. Hemisphere bands comparator: left = baseline, right = target. Blue parallels indicate perimeter along latitude; red band shows area via average neighboring circumferences times slant. Bars below sketch π-free indices.

From the circle to the sphere (concept)

We stack latitude bands along the vertical axis and apply the same grammar on each parallel (GRG circumference on the band’s ends, averaged across the slant). Errors shrink fast with balanced vertical and angular splits. The comparator above is designed to show differences band-by-band without clutter.

Applications

  • Graphics/CAD: π-free arc lengths and areas that run fast and parallel.
  • Sampling/tiling on spheres: controllable bands, equal-area approximations.
  • Gear-like and elliptical profiles via directional radial generation.
  • Education: a concrete bridge from discrete geometry to classical values.

Status & invitation

GRG is young and grows with this site. The aim is simple: a practical, comparable, discrete language for shapes and surfaces. If the idea helps, let’s refine proofs, diagrams, and applications together.

ERG – Trapezoidal

ERG – Trapezoidal

Document Status and Current Reference

This page preserves an earlier formulation or research method. Its notation, formulas, and embedded demonstrations belong to that version and should not be treated as the specification of the current binary radial kernels.

Read the current formulation and its analytical verification. Historical performance statements and implementation behaviour require separate validation. Earlier references to circular constants remain in this historical document; they are not inputs to the current geometric kernel definitions.

Δx Trapezoidal (Cartesian)

Area from the inside on the quarter of a circle. Sample uniformly on the x-axis with $\Delta x=r/N$, build trapezoids strictly under the arc $y(x)=\sqrt{r^2-x^2}$, then multiply by 4 for the full circle. No angles and no external constants appear in the primary operators.

Setup (quarter)

$$ \Delta x=\frac{r}{N},\qquad x_k=k\,\Delta x\ \ (k=0,\dots,N),\qquad y_k=\sqrt{\,r^2-x_k^2\,}. $$

Quarter area (→ ×4 for the circle)

$$ A_{\text{quarter}}(N)\;=\;\sum_{k=0}^{N-1}\frac{\Delta x}{2}\,\bigl(y_k+y_{k+1}\bigr), \qquad A_{\text{circle}}(N)\;=\;4\,A_{\text{quarter}}(N). $$

Because $y(x)$ is concave on $[0,r]$, these trapezoids give a monotone under-approximation. The error typically decays like $O(N^{-2})$. A quick check is step-doubling: compare $A_{\text{circle}}(N)$ vs $A_{\text{circle}}(2N)$.

Quarter with Δx trapezoids (Cartesian x–y) Strict Cartesian axes (equal units). r=100, N=5 ⇒ Δx=20. Trapezoids are inside the quarter arc. x y
Fig. 1 — Δx trapezoids strictly inside the quarter arc (example: $r=100$, $N=5$).

Single strip geometry

Between $x_k$ and $x_{k+1}$ the trapezoid has bases $y_k$ and $y_{k+1}$ and width $\Delta x$; the arc segment lies strictly above it.

One Δx trapezoid under the arc Strip between x=60 and x=80 for r=100. Heights are y_k=80 and y_{k+1}=60. Dashed arc segment shown above the trapezoid. x y
Fig. 2 — One Δx strip (here $k=3$): the arc is strictly above the trapezoid.

Pseudocode

Input: r, N
Δx = r / N
Aq = 0
for k = 0..N-1:
    yk   = sqrt(r*r - (k*Δx)*(k*Δx))
    yk1  = sqrt(r*r - ((k+1)*Δx)*((k+1)*Δx))
    Aq  += 0.5 * Δx * (yk + yk1)
A_circle = 4 * Aq
  

Notes

  • Monotone: $A_{\text{circle}}(N)$ increases with $N$ and approaches the true area from below.
  • Accuracy: smooth $O(N^{-2})$ decay; step-doubling is a simple diagnostic.
  • Implementation: works cleanly with fixed-point or vectorized code (just sums and square roots).