国土地理院の作業規定の準則の計算式集にある
「1.2 楕円体の諸公式」
にある子午線曲率半径および、卯酉(ぼうゆう)線曲率半径を導出する記事である。
これから、ガウス・クリューゲル等角投影法などを理解しようと試みる人への手助けである。
子午線曲率半径および卯酉線曲率半径の意味
微分幾何学の空間曲線論の話になる。
数学では、曲率を以下のように定義している。
空間内のある曲線のある点での曲率とは、
その点において、曲線を円で近似したときの半径の逆数である。
円で近似したとき、その円の半径を曲率半径といい、
その円を、曲率円という。
つまり、
曲率=1/(曲率半径)
である。
例としては、
直線は、曲率0である。
曲率半径は、無限大。
つまり、直線は、半径無限大の円と考えるわけである。
真円は、
曲率円はその円であり、曲率半径はその円の半径である。
道路などにあるR=400は、
次のカーブは半径400mの円の円弧ですよ、
という意味である。
逆数の1/400が曲率になる。
子午線曲率半径と卯酉(ぼうゆう)線曲率半径の導出
それでは、(地理学的、測地学的)緯度における子午線曲率半径を求めていく。
(地理学的)緯度は、地球楕円体上の点における接平面に直行する線(法線)と赤道面とのなす角である。
楕円をパラメトライズするときの、パラメータとは異なるので注意。
これが話をややこしくしている。
(なんで緯度と楕円のパラメタをわざわざ別のものにするのだろう?)
地球を回転楕円体と見立てて、
下図のように座標を入れて考える。
しかし、子午線曲率半径を考えるには、子午線だけを考えればよく、
結局のところ、楕円だけを考えればいいことになる。
というわけで、楕円の曲率を求める。
楕円の曲率を導出する
弧長パラメタとは限らない任意のパラメタに関して、
空間曲線の曲率は以下のように書ける。
この公式を利用する。
楕円を、
とパラメトライズする。
1階と2階の微分を計算して、
空間曲線の曲率の公式に代入して、
となる。
ここで、
離心率と、
を用いて、ルートの中を整理すると、
となる。
ここで、を消去してで整理したが、これは後で楕円のパラメタから緯度に戻すときに、
作業既定の準則の付録の計算式集の、
とつじつまが合うので、どちらで整理してもよい。
これで楕円の曲率(したがってその逆数の曲率半径)が求められた。
楕円のパラメタと緯度の換算式
この楕円のパラメタに関する曲率の式を、
緯度に関する式に換算するために、との関係式を求める。(さえわかればいい)
緯度と回転楕円体の接平面が直行するので、
接ベクトルと、
緯度の方向ベクトルが直行することを使って、
となる。
だから2乗して、をに直して、
より、
だから、
となる。
結局、
を得る。
(厳密なことを書くと、のとき、となるので、としておく。ゼロ割防止のため)
子午線曲率半径を求めよう
以上により得られた楕円の曲率
とパラメタと緯度の換算式
より、緯度における子午線曲率は、
となる。
ここで、
離心率より、
を使って、
を消去した。
卯酉線曲率半径を求める
卯酉(ぼうゆう)線が何かということと、卯酉円曲率半径が以下の図の赤線部になることは前提として、この赤線部の長さを求めていく。
経度と楕円のパラメタの換算式が求まった今、それほどややこしいことはない。
流れとしては、赤線を表示するベクトル方程式を求めて、
軸との交点を求めて、線長を計算するだけである。
赤線のベクトル方程式は、
楕円上の点を通り、
傾きは、接ベクトルを反時計回りに90度回転させたものになるので、
となる。(は任意の実数)
を代入して、
を得る。
よって、赤線部の線長は、
(数式の=をそろえられなくてすいません。)
まとめ
子午線曲率が、
だから、
子午線曲率半径はその逆数
卯酉線曲率半径は、
となる。