SLA lib(有名な位置天文ライブラリ)

1. SLA libのインストール
2. Astro::SLAのインストール
3. 黄経黄緯<=>赤経赤緯 on PDL


#######################
1. SLA libのインストール
#######################
C vertionは売り物！(TPOINT http://www.tpsoft.demon.co.uk/)、
fortran versionはstarlink projectから得られる。

http://www.starlink.rl.ac.uk/
http://www.starlink.rl.ac.uk/store/store.html
http://www.starlink.rl.ac.uk/store/storelibs.html
http://www.starlink.rl.ac.uk/cgi-store/storeform1?SLA

<TIPS>
 install情報は $HOME/.star_config
 に記録されるようだ。

いろいろ選べるようだが、デフォルトの
$HOME/star
にインストールすることにしてみよう。

GUIにて、"source only"を選び、"download"すると、
自己解凍shell script "storeform1hdlr"が得られた。

自己解凍とインストール
$HOME/star にインストールする場合は、
sh storeform1hdlr
に続いて、returnを押していけばOK。
~/star 以下に、binだのlibだの、いろいろ生成した。

使い方は、
http://www.starlink.rl.ac.uk/star/docs/sun67.htx/sun67.html
参照。


#######################
2. Astro::SLAのインストール
#######################
CPAN www.cpan.org とかからmodule "Astro::SLA"のソース Astro-SLA-0.99.tar.gz
を取り寄せ、インストールする。
# SLA C versionは売り物！(TPOINT http://www.tpsoft.demon.co.uk/)
# SLA starlink/fortran versionを使う
# そのためには、starlinkの場所を環境変数"STARLINK"で指示しておく必要がある。
tar zxvf Astro-SLA-0.99.tar.gz
cd Astro-SLA-0.99
setenv STARLINK ~/star        #<=これがポイント!
perl Makefile.PL
make
make test (passしない。ほんのすこし精度に問題があるようだ。が、十分!)
make install


#######################
3. 黄経黄緯<=>赤経赤緯 on perldl/PDL
#######################
Astro::SLAをつかってのlamdabeta=>j2000への変換を行なってみる。

>perldl
# 後で、epochをModified Juliandayで与える必要があるので、、、
perldl> use Time::JulianDay
perldl> $jd= julian_day(2000, 1, 1);
perldl> $mjd=$jd-2400000.5 
perldl> print $mjd
51544.5


perldl> use Astro::SLA
perldl> print DPI                # man Astro::SLA
3.14159265358979

(lamda,beta=>alpha,delta)
perldl> $lamda0=0;$beta0=90
perldl> slaEcleq($lamda0/180*DPI,$beta0/180*DPI,$mjd,$ra2000,$dec2000);
perldl> print $ra2000/DPI*180,$dec2000/DPI*180
270 66.5607088888889

(lamda,beta<=alpha,delta)
perldl> slaEqecl($ra2000,$dec2000,$mjd,$lamda,$beta)
perldl> print $lamda/DPI*180, $beta/DPI*180
214.877529091438 90
うーんまずまず。

see also
man Time::JulianDay
man Astro::SLA
file:///home/wada/star/docs/sun67.htx/sun67.html
http://www.starlink.rl.ac.uk/star/docs/sun67.htx/sun67.html