Calculation Trace — Campanus ハウス計算(B-7)

Layer B-7 — Campanus ハウスシステムの完全トレース
入力例:読み込み中…

⚑ このページは占星術の設計上の選択(ハウス分割)を扱います。 前段の trace_angles.html までが天文計算であり、 Campanus 法で分割するという選択には天文的な必然性はありません。

本ビューでは、実装上は関数に吸収される処理も、計算過程の理解と検証のために段階的に表示しています。

Stella-JS を読み込み中…

Campanus ハウスシステムの定義

② 定義・公式

主垂直圏(Prime Vertical)を東点(方位角 90°)を起点に 30° 等分し、 各分割点を地平線の北点・南点を通る大円で黄道に投影する。

主垂直圏は以下の4点を通る大円:

$$\text{東点(E)} \to \text{天頂(Z)} \to \text{西点(W)} \to \text{天底(N)}$$

ローカルフレームの単位ベクトル(赤道直交座標系):

$$\hat{e}_E = (-\sin\text{RAMC},\;\cos\text{RAMC},\;0) \quad\text{(東点方向)}$$ $$\hat{e}_Z = (\cos\varphi\cos\text{RAMC},\;\cos\varphi\sin\text{RAMC},\;\sin\varphi) \quad\text{(天頂方向)}$$

主垂直圏上の分割点(East 起点、反時計回り):

$$\hat{P}_i = \cos(pv_i)\,\hat{e}_E + \sin(pv_i)\,\hat{e}_Z, \qquad pv_i = (12 - i) \times 30°$$

各分割点から黄道カスプへの投影(ベクトル外積法):

$$\hat{N} = (\sin\varphi\cos\text{RAMC},\;\sin\varphi\sin\text{RAMC},\;-\cos\varphi) \quad\text{(地平線北点)}$$ $$\hat{h}_p = \widehat{\hat{N} \times \hat{P}_i} \quad\text{(ハウス円の極)}$$ $$\hat{e}_p = (0,\;-\sin\varepsilon,\;\cos\varepsilon) \quad\text{(黄道北極)}$$ $$\hat{c}_d = \widehat{\hat{e}_p \times \hat{h}_p} \quad\text{(黄道上の交点方向)}$$ $$\lambda_i = \text{atan2}(c_{d,y}\cos\varepsilon + c_{d,z}\sin\varepsilon,\;c_{d,x}) \bmod 360°$$
出典:Meeus "Astronomical Algorithms" 2nd ed. Ch.24 / Holden "The Elements of House Division"(Campanus の幾何学的定義)。
$\hat{\ }$ は正規化(単位ベクトル化)を表す。
Regiomontanus との違い:Regiomontanus は天の赤道を等分するが、Campanus は主垂直圏を等分する。 投影の手順(地平線北点大円への投影)は両者で共通。

入力値(B-1〜B-4 の出力を受ける)

① 入力値
読み込み中…

Step 1 — ローカルフレームベクトルの構成

③ 代入
計算中…
$\hat{e}_E$(東点)と $\hat{e}_Z$(天頂)は RAMC と観測地緯度 φ から定まるローカル座標系の基底ベクトル。 主垂直圏はこの2ベクトルが張る平面上にある。

Step 2 — 主垂直圏の 30° 等分と分割点ベクトル

③④ 代入・中間値
計算中…
$pv_i = (12 - i) \times 30°$:ハウス番号 $i$ に対応する主垂直圏上の角度。
東点(E)を $pv = 0°$ とし、天頂(Z)方向に正で増加する。

Step 3 — 地平線北点大円投影(H11 を例に全展開)

③④ 代入・外積計算
計算中…
外積 $\hat{N} \times \hat{P}$ は「$\hat{N}$ と $\hat{P}$ の両方に垂直な方向」=ハウス円の極を与える。
さらに $\hat{e}_p \times \hat{h}_p$ で黄道北極とハウス極に垂直な方向=黄道上の交点方向が得られる。

Step 4 — 全カスプの黄経変換(8 中間カスプを一括計算)

④ 全カスプ
計算中…
Campanus では H5・H6・H8・H9(下半球)も _horizonCusp() で直接計算する。
$pv$ が 180°〜360° の場合、$\cos(pv) < 0$ になるため $\hat{P}$ が「天底方向」に向く。 外積 $\hat{N} \times \hat{P}$ の向きが反転し、$\hat{c}_d$ も反転するため、 $\lambda$ が自動的に上半球対応点 $+180°$ に近い値となる。
(Placidus は $\lambda_{\text{lower}} = \lambda_{\text{upper}} + 180°$ として明示的に導出するが、 Campanus の外積法では不要)

Stella-JS 照合(housesCampanus())

計算中…

出力まとめ — 全 12 カスプ

計算中…