Calculation Trace — Regiomontanus ハウス計算(B-8)

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

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

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

Stella-JS を読み込み中…

Campanus との比較(設計選択の違い)

② 定義の対比
項目Campanus(B-7)Regiomontanus(B-8)
等分する大円主垂直圏(東点 → 天頂)天の赤道(赤緯 0° の大円)
起点東点(方位角 90°)RAMC(子午線と赤道の交点)
分割点ベクトル $\hat{P}$$\cos(pv)\hat{e}_E + \sin(pv)\hat{e}_Z$(z ≠ 0)$(\cos\theta,\;\sin\theta,\;0)$(赤道面上・z = 0)
下半球の導出外積の向きで自動処理明示的に $+180°$
投影手順共通 — 地平線北点大円への投影(ベクトル外積法)
どちらも「ある大円を 30° 等分し、地平線北/南点を通る大円で黄道に投影する」という同じ手順を使う。
唯一の設計上の違い:等分する大円が主垂直圏か赤道か。
出典:Holden, J.H. "The Elements of House Division"(各システムの幾何学的定義); Meeus "Astronomical Algorithms" 2nd ed. Ch.24

Regiomontanus ハウスシステムの定義

② 定義・公式

天の赤道(Celestial Equator)を RAMC を起点に 30° ずつ等分し、 各分割点を地平線の北点・南点を通る大円で黄道に投影する。

$n$ の定義:RAMC(子午線と赤道の交点 = MC の赤経)を起点に、 赤道を東方向へ 30° 刻みで数えた通し番号。

nθ = RAMC + n×30°対応するハウス導出方法
0RAMC + 0°H10(MC)軸点 — 子午線×黄道交点(球面三角法)
1RAMC + 30°H11赤道等分 → 黄道投影
2RAMC + 60°H12赤道等分 → 黄道投影
3RAMC + 90°(H1 には使わない)RAMC+90° を投影しても ASC とは一致しない。
H1 = ASC は地平線×黄道交点として別途計算(球面三角法)
4RAMC + 120°H2赤道等分 → 黄道投影
5RAMC + 150°H3赤道等分 → 黄道投影
6RAMC + 180°H4(IC)軸点 — IC = MC + 180°
n=3 が使われない理由:赤道を純粋に等分すると n=3 は「RAMC + 90°」の赤道点であり、 これを地平線北点大円で投影した黄経は ASC とは一般に一致しない(緯度・ε によって差が生じる)。 Regiomontanus は「H1 = ASC」という占星術の取り決めを優先し、n=3 の投影値を捨てて 球面三角法で定義された ASC を H1 に採用する。
これは「等分」の原則を部分的に破る設計選択であり、天文的な必然性はない。

赤道上の分割点(n=1,2,4,5 のみ使用):

$$\theta_n = \text{RAMC} + n \times 30°, \qquad \hat{P}_n = (\cos\theta_n,\;\sin\theta_n,\;0)$$

地平線北点大円への投影(Campanus と共通):

$$\hat{N} = (\sin\varphi\cos\text{RAMC},\;\sin\varphi\sin\text{RAMC},\;-\cos\varphi)$$ $$\hat{h}_p = \widehat{\hat{N} \times \hat{P}_n},\quad \hat{e}_p = (0,\;-\sin\varepsilon,\;\cos\varepsilon),\quad \hat{c}_d = \widehat{\hat{e}_p \times \hat{h}_p}$$ $$\lambda_n = \text{atan2}(c_{d,y}\cos\varepsilon + c_{d,z}\sin\varepsilon,\;c_{d,x}) \bmod 360°$$

下半球(H5/H6/H8/H9)は対称的に導出:

$$\lambda_{H5} = \lambda_{H11} + 180°,\quad \lambda_{H6} = \lambda_{H12} + 180°,\quad \lambda_{H8} = \lambda_{H2} + 180°,\quad \lambda_{H9} = \lambda_{H3} + 180°$$
$n=3$ は H1(ASC)、$n=0$ は H10(MC)に相当する。これらは赤道等分ではなく天文的に定義された軸点のため、 Regiomontanus の等分式ではなく前段の MC・ASC 計算から取得する。

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

① 入力値
読み込み中…

Step 1 — 赤道上の分割点ベクトル(θ = RAMC + n×30°)

③④ 代入・中間値

① $\theta_n = \text{RAMC} + n \times 30°$(一般式)に RAMC の数値を代入(②)し、 cos・sin を評価(③)して $\hat{P}_n$ を構成する。 Campanus と異なり、$\hat{P}_n$ の z 成分は常に 0(赤道面上)。

計算中…
$|\hat{P}_n| = 1$(単位ベクトル確認)。z = 0 は赤道面上にあることの定義的帰結であり、補正ではない。
Campanus の $\hat{P}_i$ は z ≠ 0(主垂直圏は赤道面の外に出る)のに対し、 Regiomontanus の $\hat{P}_n$ は赤道面内に留まる。

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

③④ 代入・外積計算

① 一般式(記号のまま) → ② 数値代入(φ・RAMC・ε の値を代入した式) → ③ 評価結果(数値計算後)

計算中…
この投影手順($\hat{N}$・外積・正規化・atan2)は Campanus と完全に同一。 入力となる $\hat{P}$ の構成方法だけが異なる。

Step 3 — 全カスプの黄経変換(上半球 + 下半球の明示的導出)

④ 全カスプ
計算中…
Regiomontanus では下半球(H5/H6/H8/H9)を「上半球 + 180°」として明示的に導出する。
Campanus では $\hat{P}$ の z 成分の符号が反転することで外積が自動的に逆向きになるが、 Regiomontanus の $\hat{P}$ は z = 0 固定のため、この自動処理が成立しない。
そのためコード上も cusps[4] = normAngle(cusps[10] + 180) のように明示される。

Stella-JS 照合(housesRegiomontanus())

計算中…

出力まとめ — 全 12 カスプ

計算中…