地球楕円体上の点
地表の位置は緯度、経度で示されます。地表はデコボコしているので高さも必要です。地球楕円体上の点を楕円体高ゼロで考えます。
常用されている緯度は地理的緯度φです。これは観測点で「おもり」を垂らした垂線と赤道面のなす角です。地球を回転楕円体と見なすので地表の観測点の垂線は地心を通りません。
経度をλとして、(φ、λ) から点A (x,y,z) を計算することを考えます。
地理的緯度をφ、地心緯度をθ とすると両者は図のような関係です。
図は地球楕円体を極軸を含む平面で切断した断面です。また、点Aを通る経線です。地球楕円体は球が極軸を軸に回転して膨らんだ回転楕円体なので赤道は真円です。したがって、どの経度の経線も同じ楕円を描きます。z座標は緯度φで決まり、経度λに依存しません。
断面をXZ平面と見なして考えます。赤道半径に当たる長半径aは常にa=1で計算します。
|
楕円をRで描くには、上のようなスクリプトが簡便です。0から2πの角度tに対して、
点
を、計算して結びます。
上の図で、中心角がβのとき、円周上の点です。点Pを通る垂線と、楕円の交点Aは、 です。
楕円は、 です。したがって、図の範囲では、 です。
とき、
(1) |
です。
この作図は、角度βに対して、点Aを描いています。この媒介変数βは、更成緯度です。
更成緯度βと、地理的緯度φの関係を考えます。
長半径が1の楕円の点における接線は、 です。
なので、接線の傾きは です。したがって、これに直交する法線の傾きは、 です。法線の傾きは なので、 において、
更成緯度βと、地心緯度θの関係を考えます。
したがって、地理的緯度φと地心緯度θの関係は、
(2) |
計算はすべて弧度で計算します。しかし、説明には適宜「度」を使います。
緯度は、赤道をゼロとして、北緯90゜、南緯90゜まで測ります。計算上は南緯を負で表します。
経度は、ゼロとする経線から、東経180°、西経180°まで測ります。計算上は0から360°で表し、東周りで測ります。計算上の179°は東経179°です。181°は、西経179°です。
図はRGLで描画するので右手系で考えます。
図は、長半径1、短半径0.6の地球楕円体です。 更成緯度で北緯45°、東経30°の地点Aを描いています。 経度はX軸から東に測ります。北極方向から見て反時計回りで、(0,0,0)を通る経線を0°として360°まで測ります。 緯度は、赤道が0°で、北極が90°、南極が-90°です。 更成緯度の45°は、地理的緯度で59°、地心緯度で31°と計算されます。 |
扁平率で準拠楕円体を示すと、 です。
GRS80 では、赤道半径 a = 6,378,137 m、扁平率 f = 1/298.257222101 なので、
m
bは極半径となります。
ここでは、地球楕円体の長半径をa=1としたときの短半径bで準拠する楕円体を示します。
この値は、ほぼ1であって、作図では真球と変わらないものになります。
小さなbの値を使って強調して描きます。
前述の式(2)から、
なので、
です。
Z軸と、点Aで決まる平面による断面、X’Z平面で考えます。
です。点AのZ座標は、前述の式(1)のように、楕円の式に を代入して求めました。
したがって、地心緯度θは、
北極から赤道面を見ると図のようになります。X’Z平面の は、XYZ座標で、
(3) |