Notebook

これは日々の作業を通して学んだことや毎日の生活で気づいたことをを記録しておく備忘録である。

HTML ファイル生成日時: 2025/11/23 20:31:30.065 (台灣標準時)

Astropy の astropy.coordinates.get_body_barycentric () について

Astropy の astropy.coordinates を使うと、太陽や月や惑星などの位置を簡単に取得することができて便利でござる。 astropy.coordinates.get_body () を使うと太陽や月や惑星の赤道座標系での位置が得られるでござる。赤道座標系での位置がわかれば、座標変換して地平座標系での位置を求めることも簡単でござる。

三次元直交座標系 (X, Y, Z) での位置を求めるにはどうすればよいのか調べてみた所、 astropy.coordinates.get_body_barycentric () というものがあると知ったでござる。使ってみたところ、「なんだか結果がおかしい」と思ったでござる。惑星の座標を取得してみたところ、 Z 軸方向の座標がゼロ付近ではないのでござる。

JPL Horizons の出してくる値と比べてみたところ、 X 軸方向の座標はほぼ同じなのでござるが、 Y 軸と Z軸方向の座標の値が大きく違っていたでござる。

ここで、 JPL Horizons の state vector についての設定を見てみたでござる。 "Reference plane" が "ecliptic x-y plane derived from reference frame (standard obliquity, inertial)" となっていたでござる。もしや、と思い、他にどのような選択肢があるのか見てみたでござる。すると、 "x-y axes of reference frame (equatorial or equatorial-aligned, inertial)" というものもあったでござる。もしかしたら get_body_barycentric () が返す値は、黄道座標系での直交座標系の (X, Y, Z) ではなく、赤道座標系での直交座標系の (X, Y, Z) なのかもしれないと思い、試してみたでござる。すると、 get_body_barycentric () の返す値とほぼ同じ座標が得られたでござる。

get_body で赤道座標系での座標を取得して、それを barycentric な黄道座標系に変換し、それを三次元直交座標系に変換したら、 JPL Horizons とだいたい同じ結果が得られたでござる。

太陽系天体については、三次元直交座標系の X 軸と Y 軸を黄道面に置くことが多いと思うのでござるが、どうして get_body_barycentric () は X 軸と Y 軸を赤道面に置いた座標を返すようになっているのでござろうか。

About this article:

author: daisuke
file: 20251015_00.html
category: Misc___Python
title: Astropy の astropy.coordinates.get_body_barycentric () について
mode: public
last modified: 2025/10/15 08:32:17 (UT+08:00)
html generated: 2025/11/23 20:31:30.065 (Taiwan Standard Time)

Frequently accessed files

Misc___Taiwan/20240207_00.html

100 page views
title: 台灣から台灣の外に EMS で荷物を発送する方法

Computer___TeX/20231107_00.html

55 page views
title: LaTeX での counter の利用方法について

Misc___Taiwan/20240819_00.html

42 page views
title: 住所から台灣の郵便番号を調べる方法

Computer___TeX/20230726_01.html

31 page views
title: Beamer の Metropolis テーマで block の色を変える方法

Misc___Taiwan/20240903_01.html

28 page views
title: 台湾の郵便局で EMS を利用して荷物を海外に送る方法

Computer___Python/20250330_00.html

26 page views
title: Jupyter Notebook で Matplotlib を使うときのこと

Misc___Japan/20240718_00.html

25 page views
title: 日本の運転免許証が失効してしまったときの再取得の手続きについて

Computer___TeX/20240411_00.html

25 page views
title: LuaTeX を使って PDF ファイルを作成する方法

Computer___TeX/20240414_01.html

24 page views
title: LuaTeX を使って日本語と中国語の漢字が混在した PDF ファイルを作成する方法

Computer___Python/20220715_0.html

22 page views
title: SciPy による最小二乗法

Computer___Network/20240130_00.html

22 page views
title: Google Colaboratory で Python 3.12 を使う方法

Computer___NetBSD/20230119_00.html

22 page views
title: NetBSD でバイナリーパッケージを利用する方法

Computer___NetBSD/20250307_00.html

21 page views
title: uim と mozc による日本語入力の設定について

Computer___TeX/20240410_00.html

21 page views
title: pdfTeX を使って PDF ファイルを作る方法

Computer___Hardware/20240820_00.html

21 page views
title: Raspberry Pi 5 の演算性能

Computer___NetBSD/20250301_01.html

21 page views
title: yt-dlp で YouTube の動画をダウンロードするときのこと

Misc___Japan/20240610_00.html

21 page views
title: NHK ラジオの「聞き逃し」の変更点 (2024 年 06 月)

Computer___FreeBSD/20220621_0.html

20 page views
title: FreeBSD での X.org の設定の仕方

Computer___NetBSD/20240805_03.html

19 page views
title: NetBSD 10 のインストール

Computer___NetBSD/20241223_00.html

18 page views
title: 2024 年 12 月に購入した新しい玩具

Computer___NetBSD/20240810_00.html

18 page views
title: Raspberry Pi 5 に NetBSD-10 をインストールしてみた

Computer___Python/20220518_0.html

17 page views
title: Matplotlib で作る図の縦横比

Computer___NetBSD/20241102_00.html

17 page views
title: ImageMagick の convert コマンドについて

Computer___NetBSD/20230121_00.html

16 page views
title: pkgin を使ってバイナリーパッケージをインストールする際の注意点

Computer___TeX/20240414_00.html

16 page views
title: LuaTeX での fontspec を使ったフォントの設定の方法

Computer___NetBSD/20240101_02.html

16 page views
title: ffmpeg を使って動画に音声を追加する方法

Computer___Network/20251123_00.html

16 page views
title: Nano Banana について

Computer___NetBSD/20220406_0.html

16 page views
title: pkgsrc-wip に新たなパッケージを追加する方法

Computer___Network/20240416_00.html

15 page views
title: git push としたときの error: RPC failed

Computer___TeX/20230503_00.html

14 page views
title: LaTeX CJK で日本語や中国語を取り扱うための準備について

HTML file generated by Kinoshita Daisuke.