House Systems: Mathematical Foundations and Cusp Calculation #

Every astrologer uses houses, but few understand why different house systems give different cusps. The reason is geometric: there is no unique way to divide the space around an observer into twelve sectors. Different house systems divide different things (the ecliptic, the equator, the prime vertical, or time itself), and each produces a different mapping from sky geometry to house cusps.

This article classifies house systems by their underlying geometry, gives the exact algorithms for the most common systems, and explains the polar-latitude failures that affect time-based systems.

The Division Problem #

An astrological house system must produce twelve cusps (ecliptic longitudes) from three inputs: RAMC, geographic latitude $\phi$, and obliquity $\varepsilon$. The MC and ASC are fixed (computed as in Chart Angles), and they always form cusps 10 and 1 in quadrant-based systems. The challenge is computing the six intermediate cusps (11, 12, 2, 3 and their opposites 5, 6, 8, 9).

The intermediate cusps differ between systems because each system answers a different question about what “equal division” means.

Ecliptic-Based Systems #

These systems divide the ecliptic circle itself. They are simple, unambiguous, and work at all latitudes.

Whole Sign #

Each house occupies exactly one zodiac sign. The 1st house is the entire sign containing the ASC:

$$\text{Cusp}1 = 30° \times \left\lfloor \frac{\lambda{\text{ASC}}}{30°} \right\rfloor$$

$$\text{Cusp}_N = (\text{Cusp}_1 + (N-1) \times 30°) \bmod 360°$$

The MC is not a cusp in this system — it floats somewhere in the 9th, 10th, or 11th house depending on the chart.

Equal House #

Twelve 30° sectors starting from the exact ASC degree:

$$\text{Cusp}N = (\lambda{\text{ASC}} + (N-1) \times 30°) \bmod 360°$$

Like Whole Sign, the MC is not necessarily on cusp 10.

Porphyry #

Trisects the ecliptic arc between the four angles. For the quadrant from MC to ASC:

$$\Delta = (\lambda_{\text{ASC}} - \lambda_{\text{MC}}) \bmod 360°$$

$$\text{Cusp}{11} = (\lambda{\text{MC}} + \Delta / 3) \bmod 360°$$ $$\text{Cusp}{12} = (\lambda{\text{MC}} + 2\Delta / 3) \bmod 360°$$

And similarly for the other three quadrants. Porphyry is the only system that trisects ecliptic longitude directly. It works at all latitudes and never fails, but it has no physical or temporal justification beyond simplicity.

Space-Based Systems #

These systems divide a great circle in 3D space, then project the division points onto the ecliptic.

Campanus #

Divides the prime vertical (the great circle through the zenith and the east/west points) into twelve 30° arcs. The division points are projected onto the ecliptic along vertical circles (great circles through the zenith and nadir).

For arm angle $A_k = k \times 30°$ measured along the prime vertical from the east point:

$$\tan\lambda_k = \frac{\sin A_k}{\cos A_k \cos\varepsilon - \tan\phi,\sin\varepsilon}$$

Use atan2 for quadrant resolution. The formula has the same structure as the ASC formula because the Ascendant is itself the ecliptic’s intersection with the horizon, which is the $A = 0°$ case of the prime vertical division.

Regiomontanus #

Divides the celestial equator into twelve 30° arcs of right ascension, starting from RAMC. Through each division point, draws a great circle passing through the north and south points of the horizon (a “position circle”). The ecliptic longitude where each position circle intersects the ecliptic is the cusp.

For the $k$-th division point at RAMC + $k \times 30°$:

$$\tan\lambda_k = \frac{\sin(\text{RAMC} + 30k)}{\cos(\text{RAMC} + 30k),\cos\varepsilon - \tan\phi,\sin\varepsilon \cdot \cos(30k)}$$

Regiomontanus cusps are closely related to the mathematics of primary directions, which is why this system was historically preferred by astrologers who used directions extensively.

Time-Based Systems #

These systems trisect the time it takes for an ecliptic degree to travel through a quadrant (from horizon to meridian or vice versa). They produce cusps that depend on the diurnal motion of the sky, which is why they are called “time-based.”

Placidus #

The Placidus system is the most widely used quadrant system. Its principle: an intermediate cusp is the ecliptic degree whose semi-arc has been trisected in proportion to the cusp number.

Every point on the celestial sphere traces a diurnal arc from rising to upper culmination (the semi-diurnal arc, SA). The key quantity is the oblique ascension (OA) or oblique descension (OD) of an ecliptic degree — the RAMC at which that degree rises or sets at latitude $\phi$.

The oblique ascension of an ecliptic degree with right ascension $\alpha$ and declination $\delta$ is:

$$\text{OA} = \alpha - \text{AD}$$

where the ascensional difference is:

$$\text{AD} = \arcsin(\tan\phi,\tan\delta)$$

The semi-diurnal arc is:

$$\text{SA} = 90° + \text{AD}$$

The Placidus condition for cusp 11: the ecliptic degree $\lambda_{11}$ is the one for which:

$$\text{RAMC} = \text{OA}(\lambda_{11}) + \frac{1}{3},\text{SA}(\lambda_{11})$$

For cusp 12:

$$\text{RAMC} = \text{OA}(\lambda_{12}) + \frac{2}{3},\text{SA}(\lambda_{12})$$

The problem is that OA and SA depend on $\delta$, which depends on $\lambda$, which is what we’re solving for. This makes the equation implicit — it cannot be solved in closed form.

The Iterative Algorithm #

The standard approach is fixed-point iteration:

1. Initial guess: set $\text{RA}_{\text{cusp}} = \text{RAMC} + f \times 90°$ where $f$ is the fractional position (1/3 for cusp 11, 2/3 for cusp 12).
2. From RA, find $\lambda$: assuming the cusp lies on the ecliptic ($\beta = 0$), use the equatorial-to-ecliptic transformation:

$$\tan\lambda = \frac{\tan\text{RA}}{\cos\varepsilon}$$

3. From $\lambda$, find $\delta$: using $\sin\delta = \sin\varepsilon,\sin\lambda$.
4. From $\delta$, compute AD: $\text{AD} = \arcsin(\tan\phi,\tan\delta)$.
5. Update RA: $\text{RA}_{\text{new}} = \text{RAMC} + f \times (90° + \text{AD})$.
6. Check convergence: if $|\text{RA}{\text{new}} - \text{RA}{\text{cusp}}| < 10^{-7}$ radians (~0.02"), stop. Otherwise, set $\text{RA}{\text{cusp}} = \text{RA}{\text{new}}$ and go to step 2.

Convergence is typically achieved in 3–6 iterations for moderate latitudes. The nocturnal cusps (2 and 3) use the semi-nocturnal arc ($90° - \text{AD}$) instead.

Opposite cusps are obtained by adding 180°: cusp 5 = cusp 11 + 180°, cusp 6 = cusp 12 + 180°, etc.

Koch #

Koch divides the diurnal arc of the MC degree (not the cusp degree) into three equal time segments. At each division point, it computes what RAMC would be, then finds the ASC for that RAMC. This produces cusps that are really “the Ascendant at a fraction of the MC’s semi-arc ago.”

Koch cusps can be computed without iteration (they use the MC’s declination, which is known), but the system has stronger polar-latitude failures than Placidus because it requires the MC degree itself to have a finite diurnal arc.

Topocentric (Polich-Page) #

The Topocentric system, introduced by Polich and Page in 1961, produces cusps very close to Placidus but uses a closed-form algorithm based on “artificial latitudes.” The intermediate cusps are computed by evaluating the standard ASC formula with modified latitudes:

$$\tan\phi’{11} = \frac{\tan\phi}{3}, \qquad \tan\phi’{12} = \frac{2\tan\phi}{3}$$

$$\lambda_{11} = \text{ASC}(\text{RAMC} + 30°,; \phi’{11},; \varepsilon)$$ $$\lambda{12} = \text{ASC}(\text{RAMC} + 60°,; \phi’_{12},; \varepsilon)$$

where ASC() is the standard Ascendant formula from Chart Angles. This is computationally efficient (no iteration) and numerically stable. The resulting cusps differ from Placidus by at most a few arcminutes at moderate latitudes.

Gauquelin Sectors #

Gauquelin sectors divide the diurnal arc into 36 equal time segments (18 above the horizon, 18 below). They are not a house system in the traditional sense but use the same semi-arc mathematics as Placidus. A planet’s Gauquelin sector is computed from its mundane position — its fractional progress through its own semi-diurnal arc.

Polar Latitude Failure #

At geographic latitudes above $90° - \varepsilon \approx 66.56°$, portions of the ecliptic become circumpolar (never setting) or permanently invisible (never rising). When this happens:

• Placidus fails because $\tan\phi,\tan\delta > 1$ for some ecliptic degrees, making $\text{AD} = \arcsin(\tan\phi,\tan\delta)$ undefined. The iteration diverges or returns a domain error.
• Koch fails even more aggressively because it depends on the MC degree’s diurnal arc, which may not exist.
• Campanus and Regiomontanus also produce degenerate cusps at extreme latitudes, though the failure is less severe (cusps cluster rather than becoming undefined).
• Whole Sign, Equal, and Porphyry never fail, because they divide the ecliptic without reference to the horizon or diurnal motion.

Robust software must detect the failure condition and fall back gracefully. Common strategies:

• Fall back to Porphyry or Equal at $|\phi| > 66°$
• Compute as many cusps as possible and flag the degenerate ones
• Use the Topocentric system, which degrades more gracefully than Placidus

Worked Example #

Given: RAMC = 75.996°, $\phi = 41.9028°$, $\varepsilon = 23.4393°$.

Cusp 11 (Placidus, first iteration):

1. Initial RA = 75.996° + 30° = 105.996°
2. $\lambda = \text{atan2}(\sin 105.996°,; \cos 105.996° \times \cos 23.4393°)$ $= \text{atan2}(0.9613,; -0.2756 \times 0.9175) = \text{atan2}(0.9613,; -0.2529)$ $= 180° - 75.26° = 104.74°$
3. $\delta = \arcsin(\sin 23.4393° \times \sin 104.74°) = \arcsin(0.3977 \times 0.9662) = \arcsin(0.3842) = 22.62°$
4. $\text{AD} = \arcsin(\tan 41.9028° \times \tan 22.62°) = \arcsin(0.8973 \times 0.4163) = \arcsin(0.3735) = 21.92°$
5. RA$_{\text{new}}$ = 75.996° + (1/3)(90° + 21.92°) = 75.996° + 37.31° = 113.30°

Second iteration: repeat from step 2 with RA = 113.30°. The value converges within 4–5 iterations to the final cusp 11 longitude.

Cusp 11 (Topocentric, closed form):

$$\phi’_{11} = \arctan!\left(\frac{\tan 41.9028°}{3}\right) = \arctan(0.2991) = 16.64°$$

$$\lambda_{11} = \text{ASC}(75.996° + 30°,; 16.64°,; 23.4393°)$$ $$= \text{atan2}(\cos 105.996°,; -\sin 105.996° \times \cos 23.4393° - \tan 16.64° \times \sin 23.4393°)$$ $$= \text{atan2}(-0.2756,; -0.9613 \times 0.9175 - 0.2988 \times 0.3977)$$ $$= \text{atan2}(-0.2756,; -0.8820 - 0.1188) = \text{atan2}(-0.2756,; -1.0008)$$ $$= 180° + 15.40° = 195.40°$$

This is 15°24’ Libra. The Placidus iteration would converge to a value within a few arcminutes of this.

References #

• Holden, R. W. (1977). The Elements of House Division. L. N. Fowler.
• Polich, W., & Page, A. P. N. (1961). “The Topocentric System of Houses.” Astrologischer Auskunftsbogen.
• Dalton, J. (1893). The Spherical Basis of Astrology.
• Hand, R. (1981). Essays on Astrology: The Equal House System. Para Research.
• Meeus, J. (1998). Astronomical Algorithms, 2nd ed. Willmann-Bell. Chapter 13.