News
This Month
Introduction
Observing
Science
Gallery
Information
Site Map
|
[Japanese][English]
■ HDSデータのQuick Look Reduction on IRAF
すばる望遠鏡 HDS (High Dispersion Spectrograph)の
IRAF上でのQuick Reductionについての解説です。
基本的には山頂(もしくはリモート観測環境)でのANA環境を前提に説明し
ていますが、ユーザー各自のIRAF環境でも同様にデータの一次解析を進
めることも可能です。(その場合はこちらから
CLスクリプトをダウンロードし、各自の環境にインストールしてください)
スペクトルは一次元化されるため、空間情報が必要な場合はこれ以外の
手段を用いる必要があります。
|
|
| [1] |
ユーザーアカウント(u05XXXX など)でANAのコンソールからログインする。
|
| [2] |
デスクトップの右クリックメニューの上から二番目『Subaru Observatory
Software System (ANA Only) login』を選択。
しばらく待つとANA Menuのウィンドウが現れる。
|
| [3] |
ANA Menuのなかの『HDS』のボタンを押す。
|
| [4] |
IRAFのターミナル(xgterm)と"hdsdir"のターミナルが起動する。
初回起動時はIRAFのターミナルで自動的にmkirafのコマンドが立ち上がって
いるので、"xgterm"と入力する。
"hdsdir"は/data/o0XXXX/の下の生データを見張っているソフトで、
簡易ログとして使える。
|
| [5] |
起動したIRAF上ではoverscanなどのタスクや
set stdimage=imt4096, set imtype="fits"がすでに登録されている他、
これから解説するクイックリダクション用タスクの"hdsql"も登録されている。
|
| [6] |
各タスクをHDS用に最適化したパラメータを/opt/share/hds/uparm/の下に
保存しているので、それをコピーしておく。
cp /opt/share/hds/iparm/* ~/uparm/
|
| [7] |
IRAF上で
cl> cd /work/o0XXXX
このディレクトリの下で作業を行う。
/data/o0XXXX はユーザー書き込み不可のデータディレクトリで、
イメージはここに保存される。
|
|
|
▲"hdsql"を使ったリダクション過程
図の左側がhdsqlの内部処理。
これを行う前に右側1〜3の準備を行う。
●は
一次元化されたスペクトルを得るために最低限必要なプロセス。
ピンク
で書かれているのはそれぞれに使用される(する)IRAFのタスク。
|
|
以下ではANA上のIRAFにて、タスク"hdsql"を使っての
クイックリダクションの過程を解説する。
上図にてこの過程の概要をまとめているが、
hdsqlを使うことにより、生データから最終的には一次元化され、フラット補正および
波長較正までされたデータを作ることができる。
hdsql内の各過程は必要に応じてそれぞれスキップすることも可能である。
hdsqlを使用する前に上図右側の1〜3の準備が必要となるが、
フラット補正や波長較正が必要でない場合は最低限1のみの準備でいい。
これらの準備を完了させ、
各天体フレームについてhdsqlを実行するというのがクイックリダクションの
スタイルになる。
なお、HDSでは一回の露出で赤(CCD1=L: HDSAで奇数番号)と青(CCD2=R: HDSAで偶数番号)の
二枚のフレームが得られるが、これらはそれぞれ独立して解析する必要がある。
|
|
|
HDSの各CCDはそれぞれ二つの読み出し口を持っており、
フレームの中央にオーバースキャン領域を持っている。
BIASフレームの差し引きを行うかわりに
このオーバースキャン領域の平均値をフレームから差し引き、
フレームからオーバースキャン領域を消すタスクが"overscan"である。
また、このタスクは内部でADUからelectron number(約1.7倍)への
変換も行っている。
オーバースキャンを行う前の生データは、
FITSのExtensionテーブルを持っているため、
"hogehoge.fits[0]"というように、
IRAF上では"[0]"をファイル名につけないと取り扱うことができない。
hdsqlの処理でもまず最初に生データからオーバースキャンをした
データを生成するところから処理を始める。
[準備その1]でオーダーのテンプレートを作成する前にも
テンプレートにするフレームの生データに対してこの処理をしなくてはいけないので、
hdsqlでオーバースキャン領域処理のみをする方法を紹介しておく。
|
| [1] |
生データ(HDSA000XXXXX.fits)の保存されているディレクトリを確認する。
通常は/data/o05XXXX/ になる。
|
| [2] |
IRAF上にて
cl> eparam hdsql
タスクhdsqlのパラメータを指定する。
"indirec"でデータ格納ディレクトリおよびファイル名の頭を指定し、
それ以外の番号部分を"inid"で指定するというカタチになる。
なお、オーバースキャンしたデータはタスクを実行したディレクトリに保存
されるので、生データをコピーする必要はなくなる。
|
PACKAGE = echelle
TASK = hdsql
inid = 22045 Input frame ID
indirec = /data/o05129/HDSA000 directory of Input data
overscan= yes Overscan?
cosmicra= no Cosmicray-event rejection?
scatter = no Scattered light subtraction?
flat = no Flat fielding?
apall = no Extract spectra with apall?
wavecal = no Wavelength calibration?
(os_save= yes) Save overscaned data?
......(以下略)
|
|
| [2] |
この状態でhdsqlを実行すると "H22045o.fits"というファイルが作成されている。
これがオーバースキャン済みのデータになる。
|
▲OverScan前のイメージ(Red CCD)
|
→
|
▲OverScan後のイメージ(Red CCD)
|
|
|
[準備その1]
|
オーダートレース用のテンプレートは標準星などの明るい星のデータ、
もしくはフラットのデータを用いて作成することができる。
標準設定(StdYb)を使っている場合は、"/opt/sharey/hds/"
の下にある"StdYbLshsl2x2.fits"
(StdYb, L[Red側,(Blue側はR)], 2x2binning)などのデータを流用し、
以下の作業はスキップすることができる。
このときは database/ディレクトリの下にある apStdYbLshsl2x2 などの
ファイルも作業しているディレクトリの database/ 下ににコピーする必要がある。
|
| [1] |
hdsqlを使ってテンプレートにするデータについてoverscanのみを実行。
できた HXXXXXo.fits をわかりやすい名前にリネーム
(ここではapYbL2x2.fitsとしておく)する。
|
| [2] |
これに関してapallを実行。
パラメータの例としては以下のようになる。
オーダーの抽出時にbad columnが重なってうまくいかないときは、
"line"の値を少し変える(デフォルトはディスパージョン方向の中央の
pixel=1000[2x2bin]なので、たとえば1200とか800とかにしてみる)
ことで対処するとよい。
オーダーの抽出画面では、ライン上"m"キーでマニュアルのオーダー抽出。
1番となるべきオーダー上で"o"キー→Aperture(1)="1"で
オーダーの順番が整列。
その後"q"キーで次の作業に移る。
(yes/no)と聞かれるところは基本的に"yes"(そのままリターンキー押下)
でよい。
フィッティング時にはbad columnの影響を受けているところについては
"s"→"s"で抽出範囲の決定、"t"で抽出範囲のクリア、"f"で再フィット
などの小技を使う以外はそれほど細かい作業はない(sky等を気にしなければ)。
|
PACKAGE = echelle
TASK = apall
input = apYbL2x2 List of input images
(output = ) List of output spectra
(apertur= ) Apertures
(format = echelle) Extracted spectra format
(referen= ) List of aperture reference images
(profile= ) List of aperture profile images
(interac= yes) Run task interactively?
(find = yes) Find apertures?
(recente= yes) Recenter apertures?
(resize = yes) Resize apertures?
(edit = yes) Edit apertures?
(trace = yes) Trace apertures?
(fittrac= yes) Fit the traced points interactively?
(extract= yes) Extract spectra?
(extras = yes) Extract sky, sigma, etc.?
(review = yes) Review extractions?
(line = INDEF) Dispersion line
(nsum = 20) Number of dispersion lines to sum or median
# DEFAULT APERTURE PARAMETERS
(lower = -20) Lower aperture limit relative to center
(upper = 20) Upper aperture limit relative to center
(apidtab= ) Aperture ID table (optional)
# DEFAULT BACKGROUND PARAMETERS
(b_funct= chebyshev) Background function
(b_order= 1) Background function order
(b_sampl= -10:-6,6:10) Background sample regions
(b_naver= -3) Background average or median
(b_niter= 0) Background rejection iterations
(b_low_r= 3.) Background lower rejection sigma
(b_high_= 3.) Background upper rejection sigma
(b_grow = 0.) Background rejection growing radius
# APERTURE CENTERING PARAMETERS
(width = 15.) Profile centering width
(radius = 30.) Profile centering radius
(thresho= 0.) Detection threshold for profile centering
# AUTOMATIC FINDING AND ORDERING PARAMETERS
nfind = 22 Number of apertures to be found automatically
(minsep = 40.) Minimum separation between spectra
(maxsep = 1000.) Maximum separation between spectra
(order = increasing) Order of apertures
# RECENTERING PARAMETERS
(aprecen= ) Apertures for recentering calculation
(npeaks = INDEF) Select brightest peaks
(shift = no) Use average shift instead of recentering?
# RESIZING PARAMETERS
(llimit = -17.) Lower aperture limit relative to center
(ulimit = 17.) Upper aperture limit relative to center
(ylevel = 0.05) Fraction of peak or intensity for automatic width
(peak = yes) Is ylevel a fraction of the peak?
(bkg = no) Subtract background in automatic width?
(r_grow = 0.) Grow limits by this factor
(avglimi= yes) Average limits over all apertures?
# TRACING PARAMETERS
(t_nsum = 10) Number of dispersion lines to sum
(t_step = 3) Tracing step
(t_nlost= 10) Number of consecutive times profile is lost befor
(t_funct= legendre) Trace fitting function
(t_order= 3) Trace fitting function order
(t_sampl= *) Trace sample regions
(t_naver= 1) Trace average or median
(t_niter= 2) Trace rejection iterations
(t_low_r= 3.) Trace lower rejection sigma
(t_high_= 3.) Trace upper rejection sigma
(t_grow = 0.) Trace rejection growing radius
# EXTRACTION PARAMETERS
(backgro= none) Background to subtract
(skybox = 1) Box car smoothing length for sky
(weights= none) Extraction weights (none|variance)
(pfit = fit1d) Profile fitting type (fit1d|fit2d)
(clean = no) Detect and replace bad pixels?
(saturat= INDEF) Saturation level
(readnoi= 0.) Read out noise sigma (photons)
(gain = 1.) Photon gain (photons/data number)
(lsigma = 4.) Lower rejection threshold
(usigma = 4.) Upper rejection threshold
(nsubaps= 1) Number of subapertures per aperture
(mode = ql)
|
|
| [3] |
apYbL2x2.ec.fits という一次元化されたファイルができる。
またトレースされたアパーチャーの情報はdatabase/apapYbL2x2というファイルに保存される。
|
▲apallによるオーダー抽出
|
▲apallによるオーダーフィッティング
|
|
|
[準備その2]
|
上記のオーダートレース用のテンプレートとフラットフレーム
から規格化されたフラットを作成する。
もちろんフラットフィールディングをしない場合はこの作業はスキップできる。
オーダーテンプレートや波長較正データと違って、
標準設定でもフラットは保存していないので、必要な場合は
各自で作業されたい。
|
| [1] |
hdsqlを使って各フラットフレームをoverscanする。
HDSではたいていの場合Blue用とRed用に最適化したフラットのシリーズを
とるため、使える方のものだけを使う。
|
| [2] |
overscan済みのHXXXXXo.fitsデータのテキストリストを作成し、
imcombineで平均化されたフラットイメージを作成する
(combine=averageでよい)。
|
| [3] |
タスク apflattenで規格化されたフラットの作成。
refrenにはオーダートレース用テンプレートのファイル名を入れる。
フラットをフィッティングする関数は比較的高次で13〜7次程度くらい
(":order 11"のようにして変更)。
bad columnの影響部分や端のピクセルなどを取り除きながら行う行程は
apallとほぼ一緒である。
outputのFlatYbL2x2.nmをhdsqlでフラットとして使う。
|
PACKAGE = echelle
TASK = apflatten
input = FlatYbL2x2 List of images to flatten
output = FlatYbL2x2.nm List of output flatten images
(apertur= ) Apertures
(referen= apYbL2x2) List of reference images
(interac= yes) Run task interactively?
(find = yes) Find apertures?
(recente= yes) Recenter apertures?
(resize = yes) Resize apertures?
(edit = yes) Edit apertures?
(trace = yes) Trace apertures?
(fittrac= yes) Fit traced points interactively?
(flatten= yes) Flatten spectra?
(fitspec= yes) Fit normalization spectra interactively?
(line = INDEF) Dispersion line
(nsum = 100) Number of dispersion lines to sum or median
(thresho= 10.) Threshold for flattening spectra
(pfit = fit1d) Profile fitting type (fit1d|fit2d)
(clean = no) Detect and replace bad pixels?
(saturat= INDEF) Saturation level
(readnoi= 0.) Read out noise sigma (photons)
(gain = 1.) Photon gain (photons/data number)
(lsigma = 4.) Lower rejection threshold
(usigma = 4.) Upper rejection threshold
(functio= spline3) Fitting function for normalization spectra
(order = 3) Fitting function order
(sample = *) Sample regions
(naverag= 1) Average or median
(niterat= 5) Number of rejection iterations
(low_rej= 3.) Lower rejection sigma
(high_re= 3.) High upper rejection sigma
(grow = 0.) Rejection growing radius
(mode = q)
|
|
▲apflattenでのオーダー抽出
(テンプレートからのresize、recenterで行う)
|
▲apflatternでのオーダーフィッティング
(spline3のorder=9でフィッティング)
|
|
|
[準備その3]
|
上記のオーダートレース用のテンプレートと波長較正用のTh-Arフレーム
を取得していれば、波長較正用のデータを作成することができる。
もちろん波長情報がいらない場合はこの作業はスキップできる。
またオーダートレースのときと同様に、
標準設定(StdYb)を使っている場合は、"/opt/sharey/hds/"
の下にある"THARStdYbL2x2.fits"
(StdYb, L[Red側,(Blue側はR)], 2x2binning)などのデータを流用し、
以下の作業はスキップすることができる。
このときは 上記のfitsファイルの他に、
database/ディレクトリの下にある ecTHARStdYbL2x2 などの
ファイルも作業しているディレクトリの database/ 下ににコピーする必要がある。
|
| [1] |
hdsqlをThArのフレームに対して使用。
今回はoverscanとapallを実行することになる。
apallのアパーチャーのリファレンスには先ほどテンプレートを作ったときの
元のデータを指定する。
また、apallでのアパーチャーのresize, recenter, editは
ThArのデータなのでやらなくてよい。
apall後のデータをセーブするために"ap_save=yes"としておく。
|
PACKAGE = hds
TASK = hdsql
inid = 22073 Input frame ID
indirec = /data/o05129/HDSA000 directory of Input data
overscan= yes Overscan?
cosmicra= no Cosmicray-event rejection?
scatter = no Scattered light subtraction?
flat = no Flat fielding?
apall = yes Extract spectra with apall?
wavecal = no Wavelength calibration?
(os_save= yes) Save overscaned data?
(...中略...)
(ap_save= yes) Save apalled data?
(ap_in = ) Input frame for apall (if necessary)?
(ap_refe= apYbL2x2) Reference frame for apall
(ap_inte= yes) Run apall interactively?
(ap_rece= no) Recenter apertures?
(ap_resi= no) Resize apertures?
(ap_edit= no) Edit apertures?
(ap_trac= no) Trace apertures?
(ap_fitt= no) Fit the traced points interactively?
(ap_llim= INDEF) Lower aperture limit relative to center
(ap_ulim= INDEF) Upper aperture limit relative to center
(ap_ylev= 0.05) Fraction of peak for automatic width determinatio
(ap_peak= yes) Is ylevel a fraction of the peak?
(...後略...)
|
|
| [2] |
H22073o_ec.fits という一次元化されたファイルができる。
このファイルをわかりやすい名前 "ThArYbL2x2.ec.fits" などにかえておく。
|
| [3] |
リネームしたファイルについて、タスク"ecidentify"で、波長較正。
オーダーの数にもよるが、1オーダーにつき5〜6点、2〜3オーダーとばしで
波長情報を記入して"f"でフィット。
":xorder 4", ":yorder 4"などとして4次程度の関数でフィットするのが
いいようである。
いらない点を"d"でつぶして再フィットし、エラーが小さいとうなら、"l"で
リスト中の他のラインも同定し再フィット→いらいない点をけす。
……というような手順で完了。
|
PACKAGE = echelle
TASK = ecidentify
images = ThArYbL2x2.ec Images containing features to be identified
(databas= database) Database in which to record feature data
(coordli= linelists$thar.dat) User coordinate list
(units = ) Coordinate units
(match = 1.) Coordinate list matching limit in user units
(maxfeat= 1000) Maximum number of features for automatic identifi
(zwidth = 10.) Zoom graph width in user units
(ftype = emission) Feature type
(fwidth = 4.) Feature width in pixels
(cradius= 5.) Centering radius in pixels
(thresho= 10.) Feature threshold for centering
(minsep = 2.) Minimum pixel separation
(functio= chebyshev) Coordinate function
(xorder = 2) Order of coordinate function along dispersion
(yorder = 2) Order of coordinate function across dispersion
(niterat= 0) Rejection iterations
(lowreje= 3.) Lower rejection sigma
(highrej= 3.) Upper rejection sigma
(autowri= no) Automatically write to database?
(graphic= stdgraph) Graphics output device
(cursor = ) Graphics cursor input
(mode = q)
|
|
▲ecidentifyによるTh-Ar lineの同定
(1orderにつき6本くらい)
|
▲ecidentifyでのフィッティング誤差
(2x2 binning, 0.6 arcsec slitで0.01A以内程度)
|
|
|
|
上記の[準備](1は必須、2,3はスキップ可)がすんだら、
いよいよ天体のフレームのクイックルックをhdsqlでおこなう。
波長較正やフラット、Cosmic Rayや散乱光の処理についても
スキップすることが可能であるので、
スキップしたい場合はhdsqlのパラメータで該当する部分に
"no"を記述する。
|
| [1] |
hdsqlが行う作業は(1)overscan (2)Cosmic Ray除去 (3)Scattered Light除去
(4)Flat Fielding (5)Aperture Trace (6)波長較正
となっている。
最低 (1)と(5)を行えば一次元化されたクイックルックとなる。
ここでは全行程を行うパラメータを示すが、スキップしたい項目に関しては
パラメータで"no"を指定すればいいだけである。
各工程のなかの作業は以下のようなものになる。
| (1) |
overscan
各フレームのオーバースキャン領域の平均値を出して、
それを画像全体から差し引く。
いわゆるBIAS Subtructど同等の作業であるが、
同時にADUからelectron (photon) numberへの変換(約1.7倍)も
行っている。
|
| (2) |
cosmic-ray noise の処理
wacosm11.cl をhdsql内部で用いている。
これは、天体のフレームにメディアンフィルターをかけ、
それで天体のフレームを割り、飛び抜けた値を持つピクセルを同定し、
そのピクセルについては周囲からフィットした値を用いるという処理を行うスクリプトである。
パラメータとしては、
メディアンフィルターをかける前にゲタをはかせる際の定数がある
(ゲタをはかせるのは、オーダー間のカウントがゼロに近いと、
あとで割り算を行う時に問題になるため)。
この定数は、 cr_base として与えるが、
天体データのピークカウントと同程度以上の値を与えておくと問題ないようである
(小さい値の方が cosmic-ray hit を同定しやすいが、
リアルなデータまで拾ってしまう恐れがあるので注意。)
この処理を行うには、cosmicra を yes にする。その結果を保存したい場合
には、cr_save を yes にする。
結果は、H22045oc.fits のように、'c' のついたファイル名として保存される。
|
| (3) |
散乱光の除去
apscatter を用いる。
これは、天体の光の当たっていない範囲を指定して曲面フィットを行い、
全体から差し引くタスクである。
この処理のためには、scatter を yes にする。
結果を保存するには sc_save を yes にする。
すると、H22045ocs.fits のように、's'のついたファイルが保存される。
それ以外のパラメータ(sc_*)は、apscat のパラメータである。sc_refe
には、
オーダトレース用の点プレート(ここでは apYbL2x2)を指定する。
sc_inte によって、インタラクティブに処理を進めるかどうか指定できる。
|
| (4) |
フラット補正
フラット補正を行うには、まず apnorm や apflat などを用いて規格化され
たフラットフレームを作成する。ここでは、その結果を FlatYdL1x1sn.fitsとする。
パラメータ 'flat' を yes にする。規格化されたフラットを fl_refe に与
える。結果を保存するには、fl_save を yes にする。結果は、
H22045oscf.fits のように、'f' のついたファイル名で保存される。
|
| (5) |
オーダートレース
apallを使用している。
前述のように事前に用意したオーダートレース用のテンプレートを
ap_refeに与える必要がある。
また作業を対話的に進めるか(ap_inte)、
recenter, resize, edit等を行うかどうかも"ap_*"のパラメータで
指定できる。
結果を保存するにはap_saveをyesにする。
結果はH22045oscf_ec.fitsのように'_ec'のついたファイル名で保存される。
|
| (6) |
波長較正
波長較正を行うには、前述のようにまず ecidentify により
ThArデータの処理を行っておく必要がある。
その結果をここでは ThArYbL2x2.ec.fits とする。
波長較正を行うには、wavecal を yes にし、
ThArデータを wv_refe に与える。
結果を保存するには、wv_save を yes にする。
結果は、H22045oscf_ecw.fits のように、'w' ついたファイル名で保存される。
|
すべての工程を実行する場合のパラメータ指定例はいかのとおり。
オーダートレース用のテンプレート(apYbL2x2)、
規格化フラット(FlatYbL2x2.nm)、
波長較正用Th-Arフレーム(ThArYbL2x2.ec)の各フレームを
パラメータで指定している。
これらの作業を行うと実際にはフルリダクションとかわらない結果が
得られるはずである。
|
PACKAGE = hds
TASK = hdsql
inid = 22081 Input frame ID
indirec = /data/o05129/HDSA000 directory of Input data
overscan= yes Overscan?
cosmicra= yes Cosmicray-event rejection?
scatter = yes Scattered light subtraction?
flat = yes Flat fielding?
apall = yes Extract spectra with apall?
wavecal = yes Wavelength calibration?
(os_save= yes) Save overscaned data?
(cr_in = ) Input frame for wacosm1 (if necessary)
(cr_save= yes) Save cosmicray processed data?
(cr_base= 2000.) Baseline for wacosm1
(sc_in = ) Input frame for scattered light subtraction (if
(sc_refe= apYbL2x2) Reference for aperture finding
(sc_save= yes) Save scattered light subtracted data?
(sc_inte= yes) Run apscatter interactively?
(sc_rece= yes) Recenter apertures for apscatter?
(sc_resi= yes) Resize apertures for apscatter?
(sc_edit= yes) Edit apertures for apscatter?
(sc_trac= no) Trace apertures for apscatter?
(sc_fitt= no) Fit the traced points interactively for apscatte
(fl_in = ) Input frame for flat fielding (if necessary)
(fl_refe= FlatYbL2x2.nm) Flat frame
(fl_save= yes) Save flat-fielded data?
(ap_save= yes) Save apalled data?
(ap_in = ) Input frame for apall (if necessary)?
(ap_refe= apYbL2x2) Reference frame for apall
(ap_inte= yes) Run apall interactively?
(ap_rece= yes) Recenter apertures?
(ap_resi= yes) Resize apertures?
(ap_edit= yes) Edit apertures?
(ap_trac= no) Trace apertures?
(ap_fitt= no) Fit the traced points interactively?
(ap_llim= INDEF) Lower aperture limit relative to center
(ap_ulim= INDEF) Upper aperture limit relative to center
(ap_ylev= 0.05) Fraction of peak for automatic width determinati
(ap_peak= yes) Is ylevel a fraction of the peak?
(wv_save= yes) Save wavelength-calibrated data?
(wv_in = ) Input frame for wavelength calibration (if neces
(wv_refe= ThArYbL2x2.ec) Reference frame for refspectra
(mode = q)
|
|
| [2] |
タスクがうまく機能すれば"H22081ocsf_ecw.fits"のようなファイルが
作成されるはずである。
付加される"ocsf_ecw"の意味は、
"o"=overscan, "c"=cosmicray removed, "s"=scatterd light removed,
"f"=flat fielded, "ec"=aparture extracted, "w"=wavelength caliblated
で、途中のファイルをセーブするようにパラメータ内で指定していたら、
それに応じた途中のファイルも存在しているはずである。
|
▲apscatterでの散乱光のフィッティング
|
|
|
|
| [1] |
HDSは、一回の露出で2枚の独立なファイルをつくる。
このため、キャリブレーションファイルも2セットつくることになる。
これらを参照するパラメータの内容を書きかえる手間を省くため、
ANAではhdsqlと同様のタスクを"hdsql1", "hdsql2"としても
登録してあり、カラー毎にタスクを切り替えて使うことができる。
|
| [2] |
どの処理を行うか選択できるため、
途中まで処理済のデータを入力ファイルとして扱うことも可能ではある。
ただし少しスクリプトの書き換えが必要(だと思われる)。
|
| [3] |
出力ファイルは基本的に上書きされない。
同じファイルに同じ処理を施すときは、前につくられたファイルは別名にするか消去しておく。
|
| [4] |
エラーなどで処理が中断されたときは、中間ファイルが消されずに残って
しまうことがある。tmp*という名前になっているはずなので、消去しておく。
|
|
|
|
hdsql.cl,
wacosm11.cl,
overscan.cl
|
上記の３つのcl scriptファイルをダウンロードして、
以下のような設定をlogin.clに書き加えれば、
各自の環境にてhdsqlを使用することもできる。
"/opt/share/hds/hdsql/"はANAでの設定なので、
各自スクリプトをおいたディレクトリに書き換えが必要。
|
set stdimage = imt4096
task overscan = /opt/share/hds/hdsql/overscan.cl
task hdsql = /opt/share/hds/hdsql/hdsql.cl
task wacosm11 = /opt/share/hds/hdsql/wacosm11.cl
set imtype = "fits"
imred
echelle
|
|
Akito Tajitsu
Last modified: 2007/05/22 14:38:59 HST
|
ANAへのlogin, IRAFの起動