[In the later half of this, there is the article written in English.] ========================================== 選択した星の方位・高度をテキストで吐き出すソフトです。 二番目の星を選択し、一番目の星を基準とした相対座標も吐き出せます。 あまり考えないでプログラムを書いたので遅いです。 ======================================== ◎内容物 strdrt.exe:プログラム本体 ◎初期設定 hayami.exeと同じディレクトリに放り込みます。 基本的にhayami.exeのファイルを流用します。 メニューの File->Data Files で現れるウィンドウで、 第一星表と太陽、月、惑星のパラメータファイルを設定します。 歳差、章動、平均軌道要素、力学時換算のファイルはファイル名固定です。 ======================================== ◎使用方法 1)まず観測地の経緯度・標高、タイムゾーンを設定します。 コメント欄込で、hayami.exeの観測地ファイルを流用できます (File->Load,Save から)。 2)出力する日時を、開始日時、間隔(分)、出力行数を設定します。 実際に出力するのは(設定した行数+1)行です。 したがって最初の行と最後の行の時間差は間隔×指定行数(分)です。 3)出力するファイル名を指定します。同じ名前のファイルがある状態で Outputを押すと、警告を出して止まります。 4a)一番目の星のための座標系を選びます。 「Shift of Azimuth」は、例えば天頂座標系で真東(標準では+270degree)を 基準に観測する際、+90と入力すれば真東が0degreeになります。 4b)二番目の星を選択するかしないかを選びます。 「Object1 and 2」か「Object2 only」を選択すれば、 相対座標のためのドロップとObject2の入力枠が有効になります。 4c)二番目の星を有効にした場合、 座標系は原則として一番目の星が基準となります。 [1]Polar Coord. ひとつめの星を極とする極座標系(方位角と離角)で表示します。 方位角 0degree は「the Upper」の方向で、右回りが正です。 [2]P.Rect.Coord. 天球でいえば春分点にあたる方向をひとつめの星の位置とします。 直交座標に近いので「擬似直交座標」としておきます。 極方向を「the Upper」の方向とします。 [3]Diff. of Coord. 4a)で選択した座標系において、 (2番目の星の座標値)−(1番目の星の座標値) を計算します。「the Upper」は影響しません。 5)一番目、二番目の星をそれぞれ選択します。 選択方法は [1]座標値直接入力 座標系を選んで数値を入力、必要に応じて単位等を選択しますが、 距離を0(または数値として無効な入力)にすると無限遠点とみなします。 [2]第一星表の星表コード 5.5等星までの星に関しては、第一星表の星表コードで指定出来ます。 hayami.exeのクリック検索機能オプションの「コード等表示」を選択し、 目的の星をクリック、またはキーワード検索すればコードが分かります。 固有運動記述のある星表であれば固有運動を補正します。 [3]太陽、月、惑星を直接選択 太陽、月、惑星を選択した場合は視半径も計算・出力します (視半径を非表示にも出来ます[11)])。 [4]地球の影(本影、半影) 影までの距離を設定してください(File->Data Files)。 デフォルトは「月までの距離」です。 地球の影は、大気の影響により、固体地球で計算した場合よりも 大きくなります。 その係数はデフォルトで1.018ですが、変更する事も出来ます。 次のサイトなどを参考にして適当に決めてください。 http://eclipse.gsfc.nasa.gov/LEcat5/shadow.html 6)Output を押せば計算を始めます。 あまり沢山計算すると時間がかかります。 特に太陽、月、惑星を2つ選んで計算すると、 いちいちパラメータファイルを読み込んで計算するので時間がかかります。 そもそも一晩の観測の計画用なので、ほどほどにしておいてください。 7)終了する際に「strdrt.ini」に観測位置と時刻の情報とウィンドウ位置を記録し、 次回起動時に読み込みます。 ============================================ ◎オプション 8)メニューの Units で単位を選択します。デフォルトはdegree(10進)です。 9)一番目の星から二番目の星までの距離を、一番目の星の半径を単位として 表示出来ます(100倍まで)。「内縁」「中心」「外縁」までの距離の3種類です。 二番目の星の中心がある方向が正で、「内縁」はマイナスになる場合があります。 日食の食分は「(1−(内縁までの距離))÷2」です。 二番目の星が半径を持たない星の場合は、中心距離を表示します。 10)出力の際の日付に月名を付けるかどうかを選択できます。 デフォルトは無しです。 なお、最初の行と最後の行の時間差が30日を越える場合、 年月表示のために20行ごとにコメント行 (計算値表示の直前にある年月が表示された部分と同じもの) を挿入します。 11)2つの星両方の情報を表示する際に、2行に分割するかどうかを 選択できます。デフォルトは1行表示です。 12)太陽、月、惑星においては視半径を計算して表示しますが、 非表示にする事も出来ます(二つ選んだ場合は両方とも)。 ============================================ ◎技術情報 13)計算の手順は 一番目、二番目の星の元々の座標系(星表コード指定ではICRF、 太陽、月、惑星ではJ2000.0黄経座標系)から出発し、 一旦指定地表からの天頂座標を計算し、 指定に応じて年周光行差、日周光行差、地心視差、大気差を補正し、 その後に指定した出力座標に変換する という形です。 日周光行差+大気差、地心視差の指定は他の指定から独立しています。 「地心からの値」と「日周光行差+大気差」の両方を指定して J2000.0黄経座標を選んでも、大気差を計算します、今のところ。 14)精度はhayami.exeとほぼ同じです。 恒星は年周視差を計算してないので、明るい星で1秒程度と見ておきましょう。 太陽、月、木星より内側の惑星は3秒以内程度です。 星食計算に使えないこともないけど参考値程度に、というレベルです。 15)地球の形状はGRS80です。ジオイド補正は各自で行うように。 よほど変な使い方をしない限りジオイドの凸凹の影響は無いと思いますが。 ============================================= ◎使用言語 使用言語:F-BASIC V6.3L10 U0003(富士通) ============================================= ◎履歴 H21/05/30:超暫定版 H21/06/05:出力形式の調整,2行分割の選択肢,惑星等の視半径計算関連, レイアウト微調整,ウィンドウ位置保存, 数値表示を端数切捨てから四捨五入にして桁を1桁減 H21/06/12:計算順序による誤差蓄積が大きかったので一部入れ替え (特にassigned Drct.関連),「地球の影」関連 ============================================= http://www2s.biglobe.ne.jp/~ftceter ftcenter@mth.biglobe.ne.jp June 12th, 2009 ========================================================= ========================================================= This program outputs a text file of Azimuth and Elavation data. When you select the second star, this program outputs the relative coordinate data from the first star. ============================================ **Contents strdrt.exe:the program **Setting Put strdrt.exe into the same directory for hayami.exe. This program uses files same as hayami.exe. You select "File->Data files", then a window apears. Then you set the first star catalogue and parameter files for Sun, Moon and planets. Parameter files for procession, nutation, mean orbital elements and Dynamical Time are used with fixed names. ============================================ **How to use 1)You set the longitude, the latitude, the altitude and Time Zone of your obsevational point. you can read the obsevational point file for hayami.exe. (File->Load,Save.) 2)You set the date and time for outputing. You need to set: Year, month, day, hour, minute for starting time, interval (the unit is minutes), the number of lines to output. The program outputs (the assigned lines number +1) lines. Thus the difference from the starting-time to the ending-time is equal to (the assigned lines number)*(interval time). 3)You set the file name to output. If there is the file with the same name and you push "Output", the program alerts. 4a)You select a coordinate system for the first star. If you set Zenithal coordinate system and "Shift of Azimuth"=90degree, then the Azimuth of the true east (the normal azimuth=270 degree) becomes 0 degree. 4b)You choose which you output data for the second star. When you set "Object1 and 2" or "Object2 only", Then the program enable "the relative coordinate" and "Object 2". 4c)If you enable the second star, the coordinate system for second star is the relative coordinate system from the first star. [1]Polar Coord. The polar coordinate system that the polar is the first star. Azimuth=0 degree is the direction of "the Upper" and the clockwise is plus. [2]P.Rect.Coord. Comparing to the equatorial coordinate system of Celestial Grobe, the first star come to the equinox point. Around the origin, this coordinate is nearly equal to the Rectangular coordinate system. Thus we call "Pseudo-Rectangular Coordinate system" for convenience sake. [3]Diff. of Coord. The program outputs (the coordinate of the second star)-(the coordinate of the first star) in the coordinate system for the first star. "the Upper" do not affect it. 5)You set the first star and the second star. You can set stars by the following ways: [1]directly assigning You select a coordinate system and enter the coordinates and the distance of the star directly. If you enter "0" (or meaningless data as number) in "distance", The program calculate as a infinity point. [2]Code in the first star catalogue You can appoint a fixed atar brighter than Magnitude 5.5 by the star code in the first star catalogue. In "Click option" of Hayami.exe, you set "Star-code". Then you can see the star code, after you click on the star or you search by keyword. If you use the star catalogue with proper motion data, the program adjusts the proper motion. [3]Sun, Moon, Planets If you choose Sun Moon or Planets, the program calculates the radius too. [4]The Shadow of Earth(Umbra, Penumbra) You need to set the distance from the center of Earth to the shadow. (The default is "the distance to Moon".) The shadow of Earth is slightly enlarged by the atmosphere. The defalt factor is 1.018, but you can chage it. Please decide it with the following reference. http://eclipse.gsfc.nasa.gov/LEcat5/shadow.html 6)You push "Output", then the program output the file. If you order many lines, the program take time. Especially you assign two planets, you need a very long time. 7)When you close the window, the program save the data of the observational point, the time and the location of the main window (into "strdrt.ini"). When you boot the program the next time, the program loads the data. ============================================ **Option 8)You can choose units of data(in "Unit" menu). The default is "degree(decimal)". 9)You can output the distance from the first star to the second star by the radius of the first star (if less than 100). You can output "inside", "center", "outside" respectively. The direction to the center of the second star is "positive". The "inside" may become negative. The phase of eclipse is equal to (1-"inside")/2. If the second star do not have the radius, the program outputs "center". 10)You can choose to output the month in each line. The default is "no display". On the other hand, if the difference from the time of the first line to the time of the last line is more than 30day, the comment to display the year and month is inserted at every 20 lines. 11)When you choose "Object 1 and 2", you can divede data into 2 lines. The default is "don't devide". 12)When you select Sun, Moon or a planet, the program calculates the radius. But you can choose to not display the radius data. (If you choose two planets, you can set to display or to not display the both.) ============================================= **Technical information 13)The sequence of the calculation is: The coordinate in the original coordinate (ICRF equatorial coordinate system for fix stars, J2000.0 ecliptic coordinate system for Sun, Moon, and planets) ->translating to Zenithal coordinate system (from the assigned surface point) ->adjusting Annual aberration, Diurnal aberration, Geocentric parallax and Refraction if you assign ->translating to the assigned coordinate system. "Diurnal aberration + Refraction" and "From Geocenter" are independent from each other. Even if you choose "Diurnal abberation + Refraction", "From Geocenter" and "J2000.0 ecliptic Coordinate system", the program adjusts Refraction (for now). 14)The accuracy is same as hayami.exe. The program do not adjust annual parallax, thus fix stars may be 1 arcsecond accracy. Sun, Moon planets inside of Jupitor is 3 arcsecond accuracy. 15)The shape of the Earth is GRS80. Please adjust geoid-shape youself, if you need. =========================================== **Language F-BASIC V6.3L10 U0003(Fujitsu) =========================================== **History May 30th,2009:the temporary version Jun 05th,2009:Adjusting Output format, Option of "Two lines", calculating the radius of planets, Adjusting layout, Saving the location of the main window, Changing the number-displaying from "rounding toward 0" to "Round Half Up", taking off 1 digit from Azimuth and elavation. Jun 12th,2009:Errors occured by ordering calculations is large, so I exchange some calculations. (especially on "assigned Drct."), On the shadow of Earth. ========================================== http://www2s.biglobe.ne.jp/~ftceter ftcenter@mth.biglobe.ne.jp June 12th, 2009