パッケージの導入方法、地球上の2点間の距離・方位角の求め方の説明。パッケージについて説明していきたいと思う。関する計算を行うパッケージであり、地球上の2点間の距離を計算したり、大圏航路を求めたりすることができる。今回はまずパッケージの導入と、2点間の距離と方位角の計算方法について説明していきたいと思う。パッケージはCRANに入っているので、以下のようにして簡単にインストールすることができる。install.packages("geosphere")実際に使うときには、library()を使って読みこんでおこう。使って表現する必要がある。具体的には、経度を第1の要素とし、緯度を第2の要素とする長さ2の数値ベクトルで表現する。経度が先に来るので要注意。西経・南緯については負数を用いて表す。3地点を定義しておこう。パッケージには、地球上の2点間の距離を求める函数がいくつか用意されている。いう函数を使うのが良さそうである。東京とロンドンの間の距離を出すには以下のようにする。計算結果の単位はメートルであり、東京とロンドンの間の距離は9,585,300メートル、すなわち約9,585キロメートルということが分かる。パッケージに含まれている2点間の距離を求める函数としては以下のものがある。いずれも結果はメートルで出力される。地球を真球と見なして距離を求めるものdistCosine()distVincentySphere()distHaversine()地球を楕円体と見なして距離を求めるものdistMeeus()distGeo()distVincentyEllipsoid()地球は真球というよりも楕円体に近いので、楕円体と見なしたものの方が正確な値が出る。正確であるという話がこのパッケージのドキュメンテーションに載っている。