Notebook

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

HTML ファイル生成日時: 2026/01/10 09:34:17.487 (台灣標準時)

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: 2026/01/10 09:34:17.487 (Taiwan Standard Time)

Frequently accessed files

Misc___Taiwan/20240207_00.html

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

Computer___TeX/20231107_00.html

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

Misc___Taiwan/20240819_00.html

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

Book___Chinese/20240424_00.html

306 page views
title: 曹操の詩

Computer___TeX/20240410_00.html

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

Computer___TeX/20240411_00.html

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

Computer___TeX/20230726_01.html

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

Misc___Taiwan/20240903_01.html

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

Computer___Network/20230516_00.html

197 page views
title: OpenVPN 2.6 を使い VPN Gate に接続するときの注意点

Computer___NetBSD/20250301_01.html

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

Computer___TeX/20240414_01.html

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

Computer___Network/20241214_00.html

184 page views
title: Google Colab で Python 3.12 を使う方法

Computer___NetBSD/20230119_00.html

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

Computer___FreeBSD/20220621_0.html

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

Computer___Network/20240130_00.html

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

Computer___TeX/20240414_00.html

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

Computer___Python/20250330_00.html

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

Misc___Japan/20240718_00.html

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

Misc___Japan/20240610_00.html

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

Computer___NetBSD/20250113_00.html

136 page views
title: TigerVNC での X server の解像度の指定について

Computer___NetBSD/20240805_03.html

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

Computer___Python/20220518_0.html

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

Computer___Python/20220715_0.html

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

Computer___NetBSD/20220818_1.html

134 page views
title: Emacs の markdown-mode について

Computer___Network/20220413_1.html

134 page views
title: HTML 文書の中の一部の文字を点滅させる方法

Computer___Python/20240101_00.html

133 page views
title: Matplotlib の 3D plot においての注意点

Computer___NetBSD/20240810_00.html

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

Computer___Hardware/20240820_00.html

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

Computer___NetBSD/20250409_00.html

126 page views
title: Raspberry Pi 5 に NetBSD-10.1 をインストールする手順

Computer___NetBSD/20241102_00.html

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

HTML file generated by Kinoshita Daisuke.