HAMLOG プログラム開発室


新着表示
タイトル
記事No
投稿日
投稿者

Re: 緯度経度  方位と距離を表示する
966
2021/11/08(Mon) 07:54:17
JA9AOB 銅子
> JA9AOB 銅子さん、ありがとうございます。
>
>
> C# から VB.net への変換ミスを修復し切れていない問題も残っています。
> GlCalculate() にて
>
>
> If isEast = True Then
> longitude += CDbl(180)
> Else
> longitude = CDbl(180) - longitude
> End If
> ....の次が
>
> ''' Input:
> ''' longitude /= 20;
> ...の様に未変換のままです。 ここは必要なのかどうか?
>
> 上記のIf文と同じコードが出て来ますが、そこにも
>
> ''' Input:
> ''' longitude /= 20;
>
> -------------------------
> GlToTransit() にも同じ箇所があります。

複合代入演算子 と言われるものです。
longitude /= 20; は longitude = longitude/20; と同じです
他に
x += y
x = x + y
x -= y
x = x - y
x *= y
x = x * y
x /= y
x = x / y
等もあります

>
>
> GlCalculate() 関数の呼び出しでエラー。 パラメターの与え方が悪い?
> コールの仕方が悪いのか?
>
> 'エラー
> BC30451 'GlCalculate' は宣言されていません。アクセスできない保護レベルになっています。
>
> 下記の様にしているのですが??
> Public Class SurroundingClass
>
>   Public Function GlCalculate(ByVal edeg As Integer, 以下省略
>
> Class から Module に変えたら上記のエラーは出なくなりました。
> デバッグで動かすと
>  
> '// 1文字目 経度情報の整数部 0~17 をAからRに変換する
> gl += characters(CInt(longitude))

gl += ((int)longitude).ToString();
はdoubleをintに変換してstringに変換 すなわち整数部分の文字表現を取り出し gl の後ろに付け足す処理です。 

10進数に直して 整数部分を取り出しその値を文字化する方法もあります。
BigDecimal bigDecimal = new BigDecimal(String.valueOf(doubleNumber));
int intValue = bigDecimal.intValue();

その他に doubleを文字化して 小数点以前を取り出す方法でも可能です
String doubleAsString = String.valueOf(doubleNumber);
int indexOfDecimal = doubleAsString.indexOf(".");
gl += doubleAsString.substring(0, indexOfDecimal);


>
> ...にて「ユーザーが処理していない例外」と言うエラー
>
>
> 長くなりますので、後は落ち着いたらメールっさせて頂きます。
このスレッドは管理者によってロックされていますので、返信並びに編集は出来ません!
pagetop
タイトル
記事No
投稿日
投稿者

Re: 緯度経度  方位と距離を表示する
965
2021/11/07(Sun) 21:10:49
JG1MOU浜田
訂正します。
Ido2Radian(138, 14, 16 or &H80)
でした。

> > 最上位ビットが立っていればとか、下記は難しい。
> > int Houi_Kyori(const double Ido1, Kdo1, Ido2, Kdo2, char *sbuff)
>
> VBだとこんな感じでしょうか。
> Ido2Radian(138, 14, 16 and &H80)
このスレッドは管理者によってロックされていますので、返信並びに編集は出来ません!
pagetop
タイトル
記事No
投稿日
投稿者

Re: 緯度経度  方位と距離を表示する
963
2021/11/07(Sun) 20:00:32
JG1MOU浜田
> 最上位ビットが立っていればとか、下記は難しい。
> int Houi_Kyori(const double Ido1, Kdo1, Ido2, Kdo2, char *sbuff)

VBだとこんな感じでしょうか。
Ido2Radian(138, 14, 16 and &H80)
南緯、西経 の場合ですので、とりあえず気にしないでやってみてください。

> 実際に動かして結果を見るにはどうすればよいでしょうか。
>
> 自分の経度を度、分、秒で入れる。 Ido2Radian(138, 14, 16)
> 緯度は?
>
> 相手の経度は?緯度は?  どこにセットするのでしょうか?

自分の緯度・経度、相手の緯度・経度 で、
Ido2Radian()を4回呼び出します。
このスレッドは管理者によってロックされていますので、返信並びに編集は出来ません!
pagetop
タイトル
記事No
投稿日
投稿者

Re: 緯度経度  方位と距離を表示する
962
2021/11/07(Sun) 19:43:52
JA2BQX 太田
お世話になります。

取り敢えずは VB.net に変換しました。
ソース上ではエラーが出ない状態にはしました。(表面上はかな??)

最上位ビットが立っていればとか、下記は難しい。
int Houi_Kyori(const double Ido1, Kdo1, Ido2, Kdo2, char *sbuff)


実際に動かして結果を見るにはどうすればよいでしょうか。

自分の経度を度、分、秒で入れる。 Ido2Radian(138, 14, 16)
緯度は?

相手の経度は?緯度は?  どこにセットするのでしょうか?
このスレッドは管理者によってロックされていますので、返信並びに編集は出来ません!
pagetop
タイトル
記事No
投稿日
投稿者

Re: 緯度経度  方位と距離を表示する
961
2021/11/07(Sun) 18:13:38
JA2BQX 太田
JG1MOU 浜田さん、こんばんは。

> // Pascal言語をC言語に直しました。

ソースをありがとうございます。 C から VB.net に変換してみます。
上手く動くと良いのですが。
このスレッドは管理者によってロックされていますので、返信並びに編集は出来ません!
pagetop
タイトル
記事No
投稿日
投稿者

Re: 緯度経度  方位と距離を表示する
960
2021/11/07(Sun) 18:11:33
JA2BQX 太田
こんばんは。

> これは気がつきませんでした。
> 今まで誰も気がつかなかった? Hi

出来れば計算して方位を
 117 = 入力ウインドウの方位と距離の文字列を得る。
...様にして頂けると助かります。

ユーザーリスト登録者だけが 「免許Get7s」利用時だけの表示でしょうか?
Hamlog 使用している局が全てが利用出来ると言う分けでは無いのですね。
このスレッドは管理者によってロックされていますので、返信並びに編集は出来ません!
pagetop
タイトル
記事No
投稿日
投稿者

Re: 緯度経度  方位と距離を表示する
959
2021/11/07(Sun) 18:03:13
JA2BQX 太田
JA9AOB 銅子さん、ありがとうございます。


C# から VB.net への変換ミスを修復し切れていない問題も残っています。
GlCalculate() にて


If isEast = True Then
longitude += CDbl(180)
Else
longitude = CDbl(180) - longitude
End If
....の次が

''' Input:
''' longitude /= 20;
...の様に未変換のままです。 ここは必要なのかどうか?

上記のIf文と同じコードが出て来ますが、そこにも

''' Input:
''' longitude /= 20;

-------------------------
GlToTransit() にも同じ箇所があります。


GlCalculate() 関数の呼び出しでエラー。 パラメターの与え方が悪い?
コールの仕方が悪いのか?

'エラー
BC30451 'GlCalculate' は宣言されていません。アクセスできない保護レベルになっています。

下記の様にしているのですが??
Public Class SurroundingClass

  Public Function GlCalculate(ByVal edeg As Integer, 以下省略

Class から Module に変えたら上記のエラーは出なくなりました。
デバッグで動かすと
 
'// 1文字目 経度情報の整数部 0~17 をAからRに変換する
gl += characters(CInt(longitude))

...にて「ユーザーが処理していない例外」と言うエラー


長くなりますので、後は落ち着いたらメールっさせて頂きます。
このスレッドは管理者によってロックされていますので、返信並びに編集は出来ません!
pagetop
タイトル
記事No
投稿日
投稿者

Re: 緯度経度  方位と距離を表示する
958
2021/11/07(Sun) 17:53:03
JG1MOU浜田
これは気がつきませんでした。
今まで誰も気がつかなかった? Hi

> ユザーリストにヒットしないけど「免許状gets」から code が転記される局は方位計算はしていないのでしょうか?
このスレッドは管理者によってロックされていますので、返信並びに編集は出来ません!
pagetop
タイトル
記事No
投稿日
投稿者

Re: 緯度経度  方位と距離を表示する
957
2021/11/07(Sun) 17:45:58
JG1MOU浜田
// Pascal言語をC言語に直しました。
// 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; // 小数点第2位を四捨五入する

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; // 方位を返す}
このスレッドは管理者によってロックされていますので、返信並びに編集は出来ません!
pagetop
タイトル
記事No
投稿日
投稿者

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

太田さん どこでエラーになりますか。
from の緯度経度
to の緯度経度
をラージアンに変換して 次の式に代入するだけです。
ラージアンとは角度180度 をπラージアン と言いますπは約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;
このスレッドは管理者によってロックされていますので、返信並びに編集は出来ません!
pagetop