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