HAMLOG プログラム開発室

// Pascal言語をＣ言語に直しました。
// Delphiでは、Double型よりも大きいExtended型を使っています。
double Ido2Radian(const int i1, i2, i3)

{// 緯度・経度をラジアン値に変換する
// i1=時、i2=分、i3=秒
double Ido1;
Ido1 = i2 * 60 + (i3 & $7f); // Max 3599, 3600秒＝1度
Ido1 = Ido1 / 3600.0; // 分・秒を 0.nnnn度に直す
Ido1 = Ido1 + i1;
Ido1 = DegToRad(Ido1); // ラジアンを返す
if (i3 & $80) // 最上位ビットが立っていれば
return -Ido1; // 南緯・西経
else
return Ido1; // 東経・北緯}


int Houi_Kyori(const double Ido1, Kdo1, Ido2, Kdo2, char *sbuff)

{// Ido1～Kdo2は、Ido2Radian()で求めた値を渡す
int FAngle;
double al, al2;
double U, V, W, U2, V2, W2, DI;
// Ido1, Kdo1 から Ido2, Kdo2 への距離、方位角を求める
if (Ido1==Ido2 && Kdo1==Kdo2)
return -1; // 同じ位置です

U = cos(Ido2) * cos(Kdo2);
V = cos(Ido2) * sin(Kdo2);
W = sin(Ido2);

U2 = U * sin(Ido1) * cos(Kdo1) +
V * sin(Ido1) * sin(Kdo1) - W * cos(Ido1);
V2 = -U * sin(Kdo1) + V * cos(Kdo1);
W2 = U * cos(Ido1) * cos(Kdo1) +
V * cos(Ido1) * sin(Kdo1) + W * sin(Ido1);

DI = 63700.0 * ( Pi / 2 - ArcSin( W2 )); // 距離を求める
DI = RoundOff(DI) / 10.0; // 小数点第２位を四捨五入する

Al = ArcTan( V2 / U2 );
if (U2 >= 0.0)
Al2 = Al;
else
{ if (Al <= 0.0)
Al2 = Al + Pi;
else
Al2 = Al - Pi;
}
FAngle = Trunc(RadToDeg(Pi - Al2)); // 方位角を求める
if (DI >= 100.0) then // 三桁以上の場合
wsprintf(sbuff, "方位:%d度%5dkm", FAngle, RoundOff(DI));
else
wsprintf(sbuff, "方位:%d度 %.1fkm", FAngle, DI);
return FAngle; // 方位を返す}

> JA9AOB　銅子さん、ありがとうございます。
>
> > GLを使っていますが、緯度・経度に直して方位と距離を出します。
> > 送って置きました。参照下さい。
>
> C# からVB.net に変換して見ましたが数か所エラーが出ます。
> 調査中ですが対処出来るかどうか？
>
> どなたか、他にも２点間の緯度経度から方位を求めるコードがありませんか？
>
> VB.net版が良いですが、C++ や C# でも変換努力しますが。

太田さん　どこでエラーになりますか。
from の緯度経度
to の緯度経度
をラージアンに変換して　次の式に代入するだけです。
ラージアンとは角度１８０度　をπラージアン　と言いますπは約3.14です　Math.PI　です
緯度経度をラージアンに直すのは　角度*π/180　で得られます　変数は
double で十分です
πが字が変わって読みにくいですがパイ　です

const double r=6378.137; // km----　これは地球の半径
double sa = tolongitude - fromlongitude;
distance = r * Math.Acos(Math.Sin(fromlatitude) * Math.Sin(tolatitude) + Math.Cos(fromlatitude) * Math.Cos(tolatitude) * Math.Cos(sa)); // Km
origin = 180.0 * Math.Atan2(Math.Sin(sa), Math.Cos(fromlatitude) * Math.Tan(tolatitude) - Math.Sin(fromlatitude) * Math.Cos(sa)) / Math.PI; // 度
if (origin < 0)
origin = origin + 360;

JG1MOU 浜田さん、ありがとうございます。

> ちょっとバタバタしてまして、明日の夕方までにはメールで送ります。

お手数をお掛けします。浜田さんが参照したと言う Nifty フォーラムはクローズしていて中身は既に削除されていてソースにたどり着けませんでした。

> VB.netは、浮動小数点が doubleよりも大きい変数は扱えるのでしょうか。

どうでしょうか、通常は Int32 , Long 位しか使わないので。
Ｗｅｂで調べると下記のデータがありました。

Double System.Double 倍精度浮動小数点型 64bit 8バイト

負の値：
-1.79769313486231570×10^308 ～ -4.94065645841246544×10^-324

正の値：
4.94065645841246544×10^-324 ～ 1.79769313486231570×10^308

ちょっとバタバタしてまして、明日の夕方までにはメールで送ります。
VB.netは、浮動小数点が doubleよりも大きい変数は扱えるのでしょうか。
浮動小数点演算の精度が若干落ちるのかな、と思いまして。
C#なら、long doubleとかあると思いますが。

大昔に勉強したような三角関数、ラジアンとか使ってます。
今ではさっぱり理解できません。
教えていただいた計算式を、Delphiのソースに変えただけです。
三角関数系のライブラリがVB.netにあれば使えると思います。

> どなたか、他にも２点間の緯度経度から方位を求めるコードがありませんか？
>
> VB.net版が良いですが、C++ や C# でも変換努力しますが。

JA9AOB　銅子さん、ありがとうございます。

> GLを使っていますが、緯度・経度に直して方位と距離を出します。
> 送って置きました。参照下さい。

C# からVB.net に変換して見ましたが数か所エラーが出ます。
調査中ですが対処出来るかどうか？

どなたか、他にも２点間の緯度経度から方位を求めるコードがありませんか？

VB.net版が良いですが、C++ や C# でも変換努力しますが。

何時もお世話になります。

確認です。ユザーリスト、オープンさせていて、「免許状gets」[ From ]　で使用中に
ユザーリストにヒットすればＣｏｄｅが入り、なおかつ方位も表示されます。

ユザーリストにヒットしないけど「免許状gets」から code が転記される局は方位計算はしていないのでしょうか？

JA9AOB　銅子さん、こんにちは。何時もありがとうございます。

> 役に立かは分かりませんが、C#のもの送りますか。
> GLを使っていますが、緯度・経度に直して方位と距離を出します。Hamlogとほぼ同じの結果が得られます。
> 送って置きました。参照下さい。

いま、受取りました。C#から VB.net への変換をこれからします。

>
> それにしてもWeb上に VB.net で検索で「緯度経度　方位」のサンプルが出てこない。
>
AOB　銅子です
役に立かは分かりませんが、C#のもの送りますか。
GLを使っていますが、緯度・経度に直して方位と距離を出します。Hamlogとほぼ同じの結果が得られます。

送って置きました。参照下さい。

浜田さん、こんばんは。

> WM_COPYDATA で
> 117 = 入力ウインドウの方位と距離の文字列を得る。
> というのがありましたけど、使えますか？

あらら、これは見落としていました。先日までは何回か HamlogMs.txt を
見ていたのですが、注意不足でした。ありがとうございます。

それにしてもWeb上に VB.net で検索で「緯度経度　方位」のサンプルが出てこない。

WM_COPYDATA で
117 = 入力ウインドウの方位と距離の文字列を得る。
というのがありましたけど、使えますか？

入力ウインドウで適当に W6AAA と入れてEnterキーを押し、
117番を送ったら、
「方位:44度 9302km」という文字列が返ってきました。

> 何時もお世話になります。
> Hamlog ユーザーからの依頼なのですが
> メニューの「オプション」「環境設定」「設定３」にての
> 緯度経度 「方位と距離を表示する」の設定をしておいて、
> データ入力ウインドウの右上に表示される「方位」データを VB.net にて取得出来る方法はありますか？
> VB.net にて２点間の緯度経度から方位を取得するサンプルコードは未発見です。
> 有っても直ぐに利用出来るかどうかですが。
>
> 国内局の時にはユーザーリストにマッチする等でｃｏｄｅ欄に記入されないと表示されませんが。
>
> ＤＸ局の時にはコールサインから Entity を確定して Entity と方位を
> 関連つけたファイルから取得は出来ました。