`shift_x_y_to_zk1_2_wavediff()` function has dimensional inconsistencies
まだ誰も着手していません。
評価
- 難易度
- 4/5
- 見積もり時間
- 3〜5日
- 初心者へのやさしさ
- 35/100
- issue の種類
- バグ
- 明瞭さ
- 説明が足りない
- 活発さ
- 停滞
- 技術スタック
- python
調査の方向性
Read src/wf_psf/utils/preprocessing.py around shift_x_y_to_zk1_2_wavediff() and src/wf_psf/utils/centroids.py around compute_zernike_tip_tilt. Trace the displacement units through both functions and clarify the intended physics and Zernike wavediff convention. Done means the arctan input and conversion path are dimensionally consistent, with the affected behavior verified.
索引モデルが issue の本文から書いたものです。
説明
Hi @tobias-liaudat,
There seems to be a dimensional analysis issue in the code you originally wrote. Since you're most familiar with the intended physics here, could you take a look?
In the code snippet below, the function takes arctan() of a quantity with units [1/m], which is dimensionally incorrect:
Moreover, this function is called incompute_zernike_tip_tilt. There, the arg passed in is the displacement pixel shift multiplied by pixel_sampling=12e-6 m, which will do the conversion from the [pixel] units to [m].
But, then this is undone by division by reference_pix_sampling = 12e-6 in shift_x_y_to_zk1_2_wavediff in the arctan function above. Both values are the same, making this conversion roundtrip unnecessary.
Just to add: this pattern also appears in the original code from last year:
current_shifts = centroid_calc.return_shifts() # In [pixel]
dx = reference_shifts[1] - current_shifts[1] # In [pixel]
dy = reference_shifts[0] - current_shifts[0] # In [pixel]
z1_wavediff = shift_x_y_to_zk1_2_wavediff(dx * pixel_sampling) # Input in [m]
z2_wavediff = shift_x_y_to_zk1_2_wavediff(dy * pixel_sampling) # Input in [m]
return (
z1_wavediff,
z2_wavediff,
) # Output in Zernike coefficients in wavediff convention
Could you help clarify the intended physics and suggest how to make this dimensionally consistent? I can think of a couple of approaches, but you're the expert and so I'd value your input. Thanks!
- 主要言語
- Python
- スター
- 24
- フォーク
- 10
- 平均マージ
- 2日 13分
- マージ済み PR（30日）
- 2
環境構築
- Dockerfile・Docker Compose ファイルなし
- プルリクエストのテンプレートあり
- コントリビューションガイドを読む
はじめの一歩
- issue を最後まで読み、次にプロジェクトのコントリビューションガイドを読みます。
- 着手することを issue にコメントします — 二人が同じ作業をするのを防げます。
- リポジトリをフォークし、ブランチを切って変更します。
- issue 番号を参照したプルリクエストを送ります。
CosmoStat/wf-psf のほかの issue
-
bug documentation
難易度 1/5 1〜3時間 初心者へのやさしさ 72/100
-
[DOCUMENTATION] Installation instructions list different environment name than in `environment.yml`オープンbug documentation
難易度 1/5 1時間未満 初心者へのやさしさ 78/100
-
bug
難易度 1/5 1時間未満 初心者へのやさしさ 72/100
-
[TASK] Implement a configurable ThresholdRejectionPolicy対応中かも @jeipollack が 9 日前に担当しました。 オープンneeded
-
[BUG] Update `QualityControlPipeline` to consume the dictionary-based `QualityMetric.compute()` API.対応中かも @jeipollack が 22 日前に担当しました。 オープンbug
CosmoStat/wf-psf の issue をすべて見る
似ている issue
-
難易度 2/5 1〜3時間 初心者へのやさしさ 77/100
lidofinance/diffyscan#239 ·
メンテナーはふだん 1 日以内に返信
-
docs: the dev-setup table and AGENTS.md name the reviewed_by guard that CI's Changelog job runsオープンarea: tools documentation good first issue ready for a contributor
難易度 2/5 1〜3時間 初心者へのやさしさ 82/100
eyalgolan/statecraft-israel#873 · コメント 1 件 ·
メンテナーはふだん 1 日以内に返信
-
GET rooms API returns 500 due to nonexistent Room.is_active attribute対応中かも @Mars-60 が今日担当しました。 オープンbug
難易度 2/5 1〜3時間 初心者へのやさしさ 84/100
メンテナーはふだん 1 日以内に返信
-
enhancement
難易度 2/5 1〜3時間 初心者へのやさしさ 72/100
メンテナーはふだん 5 日以内に返信
-
Harmony OPeNDAP SubSetter (HOSS) Geographic production rssmif16d 7 XYZ_PROV
難易度 2/5 1〜3時間 初心者へのやさしさ 62/100
nasa/harmony-autotester#341 ·