GRS80

測地学
基本
  • 測地学
  • 地理力学(英語版)
  • ジオマティクス(英語版)
  • 地図学
  • 測地学の歴史(英語版)
概念
技術
基準(歴史(英語版)
NGVD29(英語版) 海面測地系1929年
OSGB36(英語版) イギリス陸上測量1936年
SK-42(英語版) Systema Koordinat 1942 goda
ED50(英語版) 欧州座標系1950年
SAD69(英語版) 南米測地系1969年
GRS80GRS80地球楕円体1980年
NAD83 北米測地系1983年
WGS84 世界測地系1984年
NAVD88(英語版) 北米垂直測地系1988年
ETRS89(英語版) 欧州地球基準座標系システム1989年
GCJ-02 中国の暗号化された座標系2002年
Geo URI(英語版) 地点へのインターネットリンク 2010年

Geodetic Reference System 1980GRS80)は、測地学における地球重力ポテンシャル・地球楕円体のモデルの一つである。GRS80は一つの幾何学定数と三つの物理定数を基本的に定義している。

GRS80準拠楕円体は現在世界の測地系で最もよく使われている。その長半径は 6 378 137 m、扁平率は 1/298.257 222 101 である。

なおGPSで用いられる測地系であるWGS84の準拠楕円体はこのGRS80の長半径に等しく、扁平率はほぼ等しい。

概要

重力ポテンシャルの観点からの地球の正確な形は、平均海水面を大陸まで延長した仮想面によって表現される。これをジオイドと呼ぶ。しかし、ジオイドを直接に測地測量に用いると計算が極めて複雑になるため、通常はジオイドに最適に近似した回転楕円体を用いて地球の形を表現する。このような回転楕円体地球楕円体と呼ばれるが、そのうちで個々の測地系が準拠すべき地球楕円体準拠楕円体 (Reference ellipsoid) と呼ばれる。GRS80は、全地球的測地系を構成する準拠楕円体の一つであり、現在では世界的に最も広く使われているものである。

このシステムは、第17回国際測地学・地球物理学連合 (IUGG) において採択された。これは、本質的にGPSによる測地位置の基礎であるばかりでなく、測地学コミュニティ以外においても広範に使用される[1]

地図海図の制作のために様々な国々で使われてきた数多くの他の準拠楕円体は、次第に世界的に共通のGRS80やWGS84を用いる全地球的測地系に移行しつつある。

日本

日本における測地系も、測量法と関係法令を改正することにより、2002年4月1日から、GRS80をその一部とする全地球的測地系に移行した。具体的には、測量法施行令第3条[2]の改正によって法的に位置づけられた。

ただし、水路業務法の場合は、国際的な取り決めが異なるために、GRS80とはわずかに異なるWGS84を基にして、水路業務法施行令第2条第2号[3]を改正している。

GRS80の定義数値

基準楕円体面は、長半径(赤道半径) a {\displaystyle a} と、短半径(極半径) b {\displaystyle b} アスペクト比 ( b / a ) {\displaystyle (b/a)} 扁平率 f {\displaystyle f} かのいずれかとで通常は定義される。

しかしGRS80は、 a {\displaystyle a} G M {\displaystyle GM} J 2 {\displaystyle J_{2}} および ω {\displaystyle \omega } を基本的な定義数値として選んでいる。

定義する幾何学定数
長半径 = 赤道半径 = a {\displaystyle a} = 6 378 137 m
定義する物理定数
空気塊を含む地心重力定数 (Geocentric gravitational constant, including mass of the atmosphere) G M {\displaystyle GM} = 3 986 005・108 m3/s2
力学的形状係数 (Dynamical form factor) J 2 {\displaystyle J_{2}} = 108 263・ 10-8
回転角速度 (Angular velocity of rotation) ω {\displaystyle \omega } = 7 292 115・10-11 s-1

上記の定義値から各種の量が導出される。ただし扁平率 f {\displaystyle f} については、導出値の逆数を小数以下9桁に丸めて定義値の扱いとしてよいことになっている。

f {\displaystyle f} = 1/298.257 222 101

GRS80の導出値

扁平率 f {\displaystyle f} の導出

GRS80の離心率を与える式は、

e 2 = a 2 b 2 a 2 = 3 J 2 + 4 15 ω 2 a 3 G M e 3 2 q 0 {\displaystyle e^{2}={\frac {a^{2}-b^{2}}{a^{2}}}=3J_{2}+{\frac {4}{15}}{\frac {\omega ^{2}a^{3}}{GM}}{\frac {e^{3}}{2q_{0}}}} [4]

で与えられる。ここで、

2 q 0 = ( 1 + 3 ( e ) 2 ) arctan ( e ) 3 e {\displaystyle 2q_{0}=\left(1+{\frac {3}{\left(e'\right)^{2}}}\right)\arctan \left(e'\right)-{\frac {3}{e'}}}

なお、 e = e / 1 e 2 {\displaystyle e'=e/{\sqrt {1-e^{2}}}} は第二離心率である。すると、

e 2 = 0.006 694 380 022 903 415 749 574 948 586 289 306 212 443 890 {\displaystyle e^{2}=0.006\,694\,380\,022\,903\,415\,749\,574\,948\,586\,289\,306\,212\,443\,890\ldots }

e 2 = f ( 2 f ) {\displaystyle e^{2}=f(2-f)} であるから

1 / f = 298.257 222 100 882 711 {\displaystyle 1/f=298.257\,222\,100\,882\,711\,\ldots }

この値を小数点第9位まで丸めて

f = 1 / 298.257 222 101 {\displaystyle f=1/298.257\,222\,101}

を地球楕円体GRS80の定義値として使用している。

なお、グローバル・ポジショニング・システムで使用される測地系は、WGS84 (World Geodetic System 1984) と呼ばれる。WGS84楕円体では、計算方法の違いにより極めてわずかだが値が異なる。WGS84楕円体の扁平率は、

f = 1 / 298.257 223 563 {\displaystyle f=1/298.257\,223\,563}

と定義されている。詳細は扁平率#地球の扁平率地球楕円体#WGS84楕円体測地系#WGS84を参照。

その他の導出値

楕円または回転楕円体長半径a、短半径を b とすると、扁平率

f = a b a = 1 b a {\displaystyle f={\frac {a-b}{a}}=1-{\frac {b}{a}}}

であるので、上記の定義から、

極半径 b = 6 356 752.314 140 356mとなる。

また扁平率は、0.003 352 810 681 182 319である。

離心率の2乗 e 2 {\displaystyle e^{2}} = f(2-f) であるから、 e 2 {\displaystyle e^{2}} = 0.006 694 380 022 900 788、

離心率 e は、0.081 819 191 042 815 791となる。

第3扁平率 n は(子午線弧長の計算に用いられる)、

n = a b a + b = f 2 f {\displaystyle \quad n={\frac {a-b}{a+b}}={\frac {f}{2-f}}}

であり、n = 0.001 679 220 394 628 744 689 667である。

導出された幾何学定数
扁平率 = f {\displaystyle f} = 0.003 352 810 681 225
扁平率の逆数 = 1 / f {\displaystyle 1/f} = 298.257 222 101
短半径 = 極半径 = b {\displaystyle b} = 6 356 752.314 14 m
アスペクト比 = b / a {\displaystyle b/a} = 0.996 647 189 318 816
国際測地学・地球物理学連合が定義する平均半径 (Mean radius) R1= (2a+b)/3 = 6 371 008.7714 m
正積の平均半径 (Authalic mean radius) = 6 371 007.1810 m
同じボリュームの球の半径 (Radius of a sphere of the same volume) = ( a 2 b ) 1 / 3 {\displaystyle (a^{2}b)^{1/3}} = 6 371 000.7900 m
リニア偏心 (Linear eccentricity) = ( a 2 b 2 ) .5 {\displaystyle (a^{2}-b^{2})^{.5}} = 521 854.0097 m
両極を通る楕円断面の離心率 (Eccentricity of elliptical section through poles) = a 2 b 2 / a {\displaystyle {\sqrt {a^{2}-b^{2}}}/a} = 0.081 819 191 0435;
極の曲率半径 (Polar radius of curvature) = a 2 / b {\displaystyle a^{2}/b} = 6 399 593.6259 m
赤道上の子午線曲率半径 (Equatorial radius of curvature for a meridian) = b 2 / a {\displaystyle b^{2}/a} = 6 335 439.3271 m
子午線象限 (Meridian quadrant) = 赤道から極までの子午線弧長 = 10 001 965.7293 m;

脚注

[脚注の使い方]
  1. ^ “GRS80楕円体とは何か。”. 世界測地系移行に関する質問集(Q&A). 国土地理院. 2012年6月3日閲覧。
  2. ^ “測量法施行令”. e-Gov法令検索. 総務省行政管理局. 2012年6月3日閲覧。
  3. ^ “水路業務法施行令”. e-Gov法令検索. 総務省行政管理局. 2012年6月3日閲覧。
  4. ^ Moritz 1980, pp. 395, 398.

参考文献

  • その他の派生物理定数および測地学の数式は、次を参照されたい。Moritz, Helmut (September 1980). “Geodetic Reference System 1980”. Bulletin Géodésique 54 (3): 395–405. doi:10.1007/BF02521480.  Republished (with corrections) in Moritz, Helmut (March 2000). “Geodetic Reference System 1980”. Journal of Geodesy 74 (1): 128–162. doi:10.1007/S001900050278. 

関連項目

外部リンク

  • GRS 80 Specification (英語)
スタブアイコン

この項目は、建築土木に関連した書きかけの項目です。この項目を加筆・訂正などしてくださる協力者を求めています(プロジェクト:建築/Portal:建築)。

  • 表示
  • 編集