Skip to content

Conversation

fukudakazuya
Copy link
Contributor

Related issues

#600

Description

  • 疑似距離、搬送波位相観測量に電離層遅延量(クロブッチャーモデル)を追加した。
  • 疑似距離も搬送波位相と同様に周波数帯別に計算、ログに保存するようにした。
  • (搬送波位相にすでに実装済みのクロックバイアス項を追加した。)

Test results

L1, L2, L5それぞれの電離層遅延量
スクリーンショット 2025-06-18 185657
スクリーンショット 2025-06-18 185706
スクリーンショット 2025-06-18 185716

L1の疑似距離、搬送波位相
スクリーンショット 2025-06-18 185953
スクリーンショット 2025-06-18 190033

MATLABですでに実装、公開されているコードを実行し、同じinputで以下の計算値が同じになることを確認した。

input
klobuchar_alpha_[0] = 3.82E-8; klobuchar_alpha_[1] = 1.49E-8; klobuchar_alpha_[2] = -1.79E-7; klobuchar_alpha_[3] = 0.0; klobuchar_beta_[0] = 1.43E5; klobuchar_beta_[1] = 0.0; klobuchar_beta_[2] = -3.28E5; klobuchar_beta_[3] = 1.13E5; elevation_rad = 20.0 * math::deg_to_rad; azimuth_rad = 210.0 * math::deg_to_rad; receiver_latitude_rad = 40.0 * math::deg_to_rad; receiver_longitude_rad = 260.0 * math::deg_to_rad; gps_time_s_ = 593100.0;
S2Eの計算結果
earth_centered_angle_semii: 0.039960 latitude_ipp_semi: 0.187616 longitude_ipp_semi: 1.420409 geomagnetic_latitude_ipp_semi: 0.239793 local_time_ipp_s: 49661.667431 slant_factor: 2.176025 ionospheric_delay_gps_l1_m/c_m_s: 0.0000000793

MATLABの計算結果
psi = 0.03996 lat_i = 0.18762 long_i = 1.4204 lat_m = 0.23979 t = 49661.6674 sF = 2.176 dIon1 = 7.9335e-08

Impact

Describe the scope of influence of the changes, e.g., The behavior of feature ** changes.

Supplementary information

Provide any supplementary information.

@fukudakazuya fukudakazuya self-assigned this Jun 18, 2025
@fukudakazuya fukudakazuya requested review from sksat and a team as code owners June 18, 2025 10:12
@fukudakazuya fukudakazuya added the priority::medium priority medium label Jun 18, 2025
@fukudakazuya fukudakazuya added this to S2E Jun 18, 2025
@fukudakazuya fukudakazuya added component component emulation major update incompatible API changes labels Jun 18, 2025
@fukudakazuya fukudakazuya requested review from 200km, ogoogo, Hiro-0110, seki-hiro, suzuki-toshihir0 and conjikidow and removed request for a team June 18, 2025 10:12
Copy link
Member

@200km 200km left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

コメントつけました。

double azimuth_rad = elevation_azimuth_rad[0];
double elevation_rad = elevation_azimuth_rad[1];

double earth_centered_angle_semi = 0.0137 / (elevation_rad * rad2semi + 0.11) - 0.022;
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

[WANT] これは、こちらに定義するのが良いかなと思います。

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

この値はelevation_radに応じて変化するので定数として定義できないと思うのですがいかがでしょうか。

Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@fukudakazuya こちらのコメント確認をお願いします。

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@200km すみません、pendingになっていて返信できていませんでした。
earth_centered_angle_semiはelevation_radに応じて変化するので定数として定義できないと思うのですがいかがでしょうか。

Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

確かにそうですね。
でしたら、この式の中に入っているマジックナンバーを説明するコメントや参照などを書いてもらえると嬉しいです。

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

298行目に
// Ref: Klobuchar, J.A., (1996) "Ionosphercic Effects on GPS", in Parkinson, Spilker (ed), "Global Positioning System Theory and Applications,
// pp.513-514.
と一応参照を入れていますが、不十分でしょうか?

Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

可読性を考えると、マジックワードの近くにより細かな数値の意味を書いてもらえると嬉しいです。
例えば、「参考資料の式***参照」などあると便利だと思います。(英語で書いてもらいたいですが)

@fukudakazuya
Copy link
Contributor Author

コメントありがとうございました。修正したので確認よろしくお願いいたします。

Copy link
Member

@200km 200km left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

一つだけ、前回指摘分が残っているので確認をお願いします。それ以外は問題ないです。

@github-project-automation github-project-automation bot moved this to 🏗 In Progress in S2E Jul 7, 2025
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
component component emulation major update incompatible API changes priority::medium priority medium
Projects
Status: 🏗 In Progress
Development

Successfully merging this pull request may close these issues.

2 participants