オープンソースプラネタリウムソフトの Stellarium で、選択した天体の（左上に表示されている）情報をクリップボードにコピーしたいときは、その状態で Ctrl+C でできるようです。撮影した天体のメタ情報として記録したり、ブログにアップするときに便利そうです。

太陽
種類: 恒星
等級: -26.77 (大気減光時: -26.39)
絶対等級: 4.83
赤経/赤緯 (J2000): 15h21m20.5s/-18°27'24.7"
赤経/赤緯 (瞬時): 15h22m8s/-18°30'21"
銀経/銀緯: -14°20'58.1"/+31°42'34.5"
時角/赤緯: 1h01m43s/-18°30'21" (幾何学上の)
時角/赤緯: 1h01m42s/-18°29'09" (見かけ上の)
方位角/高度: +197°15'22"/+31°44'22" (幾何学上の)
方位角/高度: +197°15'22"/+31°45'37" (見かけ上の)
Ecliptic Topocentric (of date): +232°55'53.7"/-0°00'07.4"
Obliquity (of date, for Earth): +23°26'15.0"
距離: 0.98910856天文単位
見かけの直径: +0°32'20.4"

太陽高度角と方位角を、与えられた日時と経緯度から計算する関数を作成しました。 直接の元にした計算プログラムは、こちらのページのC言語版です。 計算方法については、こちらのページをご覧ください。

• solang.inc.php (2.8KB) / ソースコード

使い方

list( $elevation, $azimuth ) = solar_angle( $lat, $lon, $dt, $timelon );

与えられた引数における太陽高度角と方位角を返します。 lat、lon はそれぞれ、緯度、経度です。 dt にタイムスタンプ、timelon は標準時における子午線経度を指定してください。 すべての引数は省略でき、その場合、lat、lon は日本経緯度原点、dt は現在日時、timelon は 135 が使われます。 経緯度を六十進で与えたい場合は、139+44/60+28.8869/3600（東経139度44分28秒8869）のような式で書くことができます。 戻り値は、高度、方位の順の配列になりますので、list( $e, $a ) のように受け取って下さい（元プログラムではラジアンですが、度数に変換しています）。

うるう年については特別考慮していません（元のプログラムでは2月は28日固定として通算日を計算していますが、PHP版では date('z') によりうるう年の2月は29日になる点だけ異なります）。 これについては、太陽年との端数 0.2422×n 日を「$A = 2 * pi() * $yday / 365」このあたりの式で調整してあげればいいのかなと思いましたが、ちょっと自信が無いのでこのままとしました。 よりベタ移植がいい人は、通算日の計算を元プログラムに合わせて変更してください。

また、dec2sex、dir16 の関数を同梱しています。 dec2sex は、十進経緯度を六十進表記に変換します。 dir16 は、方位角（度）を16方位（南南東とか）に変換します。

こちらの関数は、LHCAM で使用するために作成しました。