[GNSSに興味] 電子基準点を使った中国から購入したアンテナの座標の校正

こんにちは

またまたGNSSについての話です。　他にも色々チャレンジしているのですが、思った以上にレベルが高くて、どれも失敗の連続です。ブログとして公開したいことはあるのですが、解決していないので公開まで至りません。

中国から購入したアンテナAとB

でGNSSの話しですが、

自宅のアンテナは「GNSSに興味を持ってしまいました　その4」で紹介したように中国から購入したアンテナが２個と奮発した TallysmanのTW2710の１個を含めた計３個です。お互いのアンテナ間は40cmほど離れています。　

アンテナAとBを使ってほぼ同時にデータを取得し電子基準点を使ってFIX解をなんとか得ることができました。

データ取得時間は、２時間と区切って３回ほどチャレンジしたのですが、３回のうち１回だけFIX解が使えました。　残りはFIX解が得られていても、バラけており収束していないことが多いです。

たまたま解が得られたものをプロットしてみました。

さて得られたデータを使って２つのアンテナ間の距離を計算してみたいと思います。

つまり、「GNSSに興味を持ってしまいました　その4」で示した結果と比較したいということです。

ネットを使って調べてみると多くのサイトで詳しい説明がなされています。

緯度は北極と南極で考えれば良いようで

単純に単純に地球の半径が使えるとのことです。

地球の半径(赤道の半径の実測値) r ＝6,378,136m だそうです。

円周の長さは　2πr で求められるので

円周 = 2 * π * 6,378,136 = 40,075,010(m)

1度は

40075010/360 = 111,319 (m)

ここでは10進法を使っているので

小数点以下7桁目が1.1cm　となります。

ここらへんの話はトランジスタ技術2018.01 73Pに説明があります。

軽度の場合の距離は緯度と違って面倒なようです。

つまり赤道上は緯度と同じ計算でもとまるようですが、

北極は０となります。

軽度によって半径が変わるとのことで

計算式は

r = Rcosθ　( R: 地球の半径、θ = 経度))、r = 経度での切断面の半径)

例えば北緯35度の場合は

6,378,136 * cos(35) = 5224663(m)

円周 = 2 * π * 5224663 = 32810884(m)

1度は 32810884/360 = 91141(m)

アンテナAとB間の距離を求めるにはネット上に多くの解説があります。

その中の一つを使って計算して見ました。　

済みませんが一部伏字にさせていただきます。