Call 名前	※必須
文字色	■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■
Ｅメール	公開または未記入非公開
タイトル	sage
ＵＲＬ
メッセージ	手動改行 強制改行 図表モード メッセージには上記と同じURLを書き込まないで下さい > // 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; // 方位を返す}
画像File	(130kBまで)
暗証キー	(英数字で8文字以内)
投稿キー	（投稿時 を入力してください）
プレビュー