解説用の星図をさっさかつくるプロジェクト　メモ9　Javascript+fireworks+ＶＢ２００５編

前のメモに引き続き、用紙を決めて、円形の星図を作ってみる。

座標変換

とりあえず北半球用とし、星座早見などで使う、正距方位図法を採用する。

まず、天の北極を中心にする（RA：0時，DEC：90度）とする。これは、

DEC　中心からの距離Rに対応：　R=（90度−DEC）＊A　　　Aは適当な倍率

RA　中心からの方位θに対応：　θ＝360／24＊ＲＡ（角度単位）　＝2＊π／24＊ＲＡ（ラジアン）＝π／12＊ＲＡ

そして、Ｘ，Ｙ，との関係は、

Ｘ＝Ｒ＊ｃｏｓθ＋Ｘ０

Ｙ＝Ｒ＊ｓｉｎθ＋Ｙ０

となる。なお、Ｘ０，Ｙ０ならびに、Ａは、紙のサイズで決まる

X0=maxpaperX/2

Y0=maxpaperY/2

Aについては

ＤＥＣの値域を（90度〜−９０度＋緯度）　となるので緯度をφとすると、

ＤＥＣ：φ−９０〜９０度　→　Ｒ：Ａ＊（０〜180-φ）であり、

Ｒの最大値は maxpaperX または maxpaperYの大きい方 maxpaper の半分でなければいけないので

Ｒmax＝maxpaper/2=A*（180-φ）　なので、整理すると

Ａ＝maxpaper / 2 *(180-φ）

となる。以上をふまえて、北半球、北緯φのばあい

θ＝　＝π／12＊ＲＡ

Ｘ＝　（90−DEC）＊　maxpaper / 2 *(180-φ）＊ｃｏｓθ＋maxpaperX/2

Ｙ＝　（90−DEC）＊　maxpaper / 2 *(180-φ）＊sinθ＋maxpaperY/2

ついてに

starR = maxstarR / 10^(0.2(Vmag))

とする。

これをVB2005 のルーチンにする のはおあずけ

2008年04月12日 記す

続く・・・のかな