IP-VLBI 相関処理ソフトウェア群

最近の更新履歴
2021-02-10 新アーカイブ公開（gcc9.3.0対応）

ソフトウェア一覧

Version 2017-10-09以降

予測値計算関係

apri_calc	予測値計算（ＶＥＸファイルもサポート）
skdchk	スケジュールファイルのチェック(SKEDおよびVEXファイルをサポート）

相関処理関係

cor 相関処理（高速版:１ビットサンプリングデータ専用）（１観測のみの処理
cor_all 相関処理（高速版:１ビットサンプリングデータ専用）（複数観測処理）
fx_cor 相関処理（１観測のみの処理）
fx_cor_all 相関処理（複数観測処理）
cor_new 相関処理（高速版:１ビットサンプリングデータ専用）（１観測のみの処理）各種フォーマットサポート
cor_all_new 相関処理（高速版:１ビットサンプリングデータ専用）（複数観測処理）各種フォーマットサポート
fx_cor_new 相関処理（１観測のみの処理）各種フォーマットサポート
fx_cor_all_new 相関処理（複数観測処理）各種フォーマットサポート
cor_mon 相関データ（cout形式およびKSP形式）の相関関数のダイナミック表示（PGPLOTでコンパイルした時のみ有効）
sdelay フリンジサーチ（粗決定サーチ）（PP毎の結果出力機能、２次サーチ機能、フリンジ位相プロット、PCALプロット機能）

cor	相関処理（高速版:１ビットサンプリングデータ専用）（１観測のみの処理
cor_all	相関処理（高速版:１ビットサンプリングデータ専用）（複数観測処理）
fx_cor	相関処理（１観測のみの処理）
fx_cor_all	相関処理（複数観測処理）
cor_new	相関処理（高速版:１ビットサンプリングデータ専用）（１観測のみの処理）各種フォーマットサポート
cor_all_new	相関処理（高速版:１ビットサンプリングデータ専用）（複数観測処理）各種フォーマットサポート
fx_cor_new	相関処理（１観測のみの処理）各種フォーマットサポート
fx_cor_all_new	相関処理（複数観測処理）各種フォーマットサポート
cor_mon	相関データ（cout形式およびKSP形式）の相関関数のダイナミック表示（PGPLOTでコンパイルした時のみ有効）
sdelay	フリンジサーチ（粗決定サーチ）（PP毎の結果出力機能、２次サーチ機能、フリンジ位相プロット、PCALプロット機能）

データチェック関係

oscillo サンプリングデータの時系列ダイナミック表示 (PGPLOTでコンパイル時のみ有効）
speana オフラインスペクトル表示 (VSSPフォーマット専用）
speana2 オフラインスペクトル表示(高機能版）（VSSPフォーマット専用）
speana_n オフラインスペクトルのダイナミック表示 (PGPLOTでコンパイル時のみ有効）（VSSPフォーマット専用）
g_speana 各種フォーマットデータのスペクトル表示
datachk データチェック（ヘッダー部を頼りにヘッダー間のデータビット数をカウント）。オプションでサンプリング統計表示機能 (VSSPフォーマット専用）
pcalcheck K5/VSSPおよびVSSP32データファイル中のPCAL位相と強度をグラフ表示する
vdifcheck VDIFのファイル形式をチェックする
speana_vdif VDIFデータのスペクトル表示
m5check Mark5のファイル形式チェック（8|16|32|64トラックモードの自動判定）およびシンクブロックチェックなしでのヘッダー部表示
m5time Mark5データの時刻表示
m5btime Mark5Bデータの時刻表示
m5vex_ana VEXファイル解析
g_pulsar パルサープロファイル（パルサー周期に同期したアンサンブル平均）を表示する（全データフォーマットをサポート）

oscillo	サンプリングデータの時系列ダイナミック表示 (PGPLOTでコンパイル時のみ有効）
speana	オフラインスペクトル表示 (VSSPフォーマット専用）
speana2	オフラインスペクトル表示(高機能版）（VSSPフォーマット専用）
speana_n	オフラインスペクトルのダイナミック表示 (PGPLOTでコンパイル時のみ有効）（VSSPフォーマット専用）
g_speana	各種フォーマットデータのスペクトル表示
datachk	データチェック（ヘッダー部を頼りにヘッダー間のデータビット数をカウント）。オプションでサンプリング統計表示機能 (VSSPフォーマット専用）
pcalcheck	K5/VSSPおよびVSSP32データファイル中のPCAL位相と強度をグラフ表示する
vdifcheck	VDIFのファイル形式をチェックする
speana_vdif	VDIFデータのスペクトル表示
m5check	Mark5のファイル形式チェック（8\|16\|32\|64トラックモードの自動判定）およびシンクブロックチェックなしでのヘッダー部表示
m5time	Mark5データの時刻表示
m5btime	Mark5Bデータの時刻表示
m5vex_ana	VEXファイル解析
g_pulsar	パルサープロファイル（パルサー周期に同期したアンサンブル平均）を表示する（全データフォーマットをサポート）

データ編集関係（VSSPフォーマット専用）

datacut サンプリングデータファイルからある時間範囲のデータを抜き出す
datatime_edit サンプリングデータファイルの時刻を修正（変更）する
extdata サンプリングデータファイルから特定のｃｈのデータを抽出しテキストファイルに書き出す
four2one ４chモードで収集したデータから任意のチャンネルのデータを切り出し、１chモードのデータに変換する
one2four ４台のＩＰ－ＶＬＢＩボードから１ｃｈモードで取得したデータファイル４つを結合し、４ｃｈモードのデータファイルを作成する
data_half サンプリングデータを間引きすることにより見かけのサンプリング周波数を半分にする
data_double サンプリングデータを繰り返すことにより見かけのサンプリング周波数を倍にする
adbitconv サンプリングデータファイルのADビット数を任意のADビット数（ただし、1,2,4,8）に変換する
data_recov K5/VSSPおよびK5VSSP32データヘッダー部の修復を行う
aux_recov K5/VSSP32データでデータチェック時にAUX MISALIGN (misaligned AUX field)と表示されるデータの修復を行う
vssplogana sampling（autoobsも可）のログファイルおよびdatachkのサマリー出力ファイルを解析し、発生エラー状況の統計結果を表示する

datacut	サンプリングデータファイルからある時間範囲のデータを抜き出す
datatime_edit	サンプリングデータファイルの時刻を修正（変更）する
extdata	サンプリングデータファイルから特定のｃｈのデータを抽出しテキストファイルに書き出す
four2one	４chモードで収集したデータから任意のチャンネルのデータを切り出し、１chモードのデータに変換する
one2four	４台のＩＰ－ＶＬＢＩボードから１ｃｈモードで取得したデータファイル４つを結合し、４ｃｈモードのデータファイルを作成する
data_half	サンプリングデータを間引きすることにより見かけのサンプリング周波数を半分にする
data_double	サンプリングデータを繰り返すことにより見かけのサンプリング周波数を倍にする
adbitconv	サンプリングデータファイルのADビット数を任意のADビット数（ただし、1,2,4,8）に変換する
data_recov	K5/VSSPおよびK5VSSP32データヘッダー部の修復を行う
aux_recov	K5/VSSP32データでデータチェック時にAUX MISALIGN (misaligned AUX field)と表示されるデータの修復を行う
vssplogana	sampling（autoobsも可）のログファイルおよびdatachkのサマリー出力ファイルを解析し、発生エラー状況の統計結果を表示する

データフォーマット変換関係

k5tom5b K5からMark5Bへの変換(K5/VSSP32フォーマットもサポート）
k5tom5 K5からMark5への変換(K5/VSSP32フォーマットもサポート）
k5tovdif K5からVDIFへのフォーマット変換を行う
ads2k5 ADS3000+(DBBCモード)からK5/VSSPへの変換
m5btok5 Mark5BからK5への変換
m5tok5 Mark5からK5への変換（Mark5全モード対応高速版：K5/VSSP32への変換は未対応）
vdif2k5 VDIFからK5/VSSP32フォーマットへの変換を行う
vdif2m5b VDIFからMark5Bへの変換
octad2vdif 生サンプリングデータからVDIFへの変換

k5tom5b	K5からMark5Bへの変換(K5/VSSP32フォーマットもサポート）
k5tom5	K5からMark5への変換(K5/VSSP32フォーマットもサポート）
k5tovdif	K5からVDIFへのフォーマット変換を行う
ads2k5	ADS3000+(DBBCモード)からK5/VSSPへの変換
m5btok5	Mark5BからK5への変換
m5tok5	Mark5からK5への変換（Mark5全モード対応高速版：K5/VSSP32への変換は未対応）
vdif2k5	VDIFからK5/VSSP32フォーマットへの変換を行う
vdif2m5b	VDIFからMark5Bへの変換
octad2vdif	生サンプリングデータからVDIFへの変換

IP-VLBI 相関処理ソフトウェア群のインストール

IP-VLBIボードがなくてもコンパイル可能です。
【注：Ver.2016-10-04以降インストール法が大幅変更になりました】
Ver.2016-08-08およびそれ以前のアーカイブのインストール法は install_cor.20161003.htmlを参照してください。

K5ソフトウェアインストール環境のセットアップ(Ubuntu, Debian向け)

UbuntuおよびDebianでは以下のパッケージをインストールするとコンパイルの環境が整います。

     sudo apt install build-essential  --  gcc, make などがインストールされる
     sudo apt install libfftw3-dev     --  FFTW3パッケージのインストール
     sudo apt install pgplot5          --  PGPLOTのインストール
     sudo apt install libx11-dev       --  PGPLOTで必要なパッケージ
     sudo apt install gfortran         --  f77のサポート

インストール法

1.	適当なディレクトリにipvlbi_cor20211104.tar.gz（Ver.2021-11-04) を置く
2.	tar xvzf ipvlbi_corXXXXXXXX.tar.gz で解凍する
3.	cd ipvlbiXXXXXXXX で作業ディレクトリを移動し make T=COR を実行する。以下で説明するFFTWパッケージがすでにインストールされている場合は make T=COR F=FFTW を実行する。グラフィック表示はデフォルトではPGPLOTですが、GNUPLOTを使用したい場合はオプションG=GNUPLOTをつけてmakeする。 (追記 2018/09/27：最近のバージョンはGNUPLOTを使ったグラフィック表示のチェックを行なっていません。グラフィック表示が不具合を起こす場合はPGPLOTを使うようにしてください。) コンパイルしなおす場合は make clean を実行後、makeし直す。
5.	【環境の設定を行う】インストールしたソフト群のディレクトリが例えば$HOME/K5/ipvlbiXXXXXXXXである場合シンボリックリンクで ln -s $HOME/K5/ipvlbiXXXXXXXX $HOME/ipvlbi すでにリンクがある場合新しく貼り直す場合は"f"オプションをつけて、 ln -sf $HOME/K5/ipvlbiXXXXXXXX $HOME/ipvlbi でリンクを作成し、.bashrcなどのリソースファイル中で export PATH=$PATH:$HOME/ipvlbi/bin (.bashrcの場合） setenv PATH $PATH:$HOME/ipvlbi/bin (.cshrcの場合）を記述しておけば、K5ソフト群のバージョンアップを行った際もシンボリックリンクを張り直すだけで良くなる。 Ver.2010-02-08以降のアーカイブではmanページファイルも含まれているので、.bashrcなどに export MANPATH=$HOME/ipvlbi/man:$MANPATH (.bashrcの場合） setenv MANPATH $HOME/ipvlbi/man:$MANPATH (.cshrcの場合）を記述しておけば、manコマンド(例えば"man cor"等)でK5ソフト群の使用方法を得ることができる。 (追記2018/07/11: 最近のlinuxのディストリビューションではMANPATHに追加の記述をしなくても.basrcなどでPATHの追加を行なっておけば K5ソフトのmanページへのパスが自動的に追加される。実際のMANPATHはmanpathコマンドで確認できる。) なお、"man k5cor"を実行することにより、K5相関処理ソフト群全体の説明を得ることができる。

ＦＦＴＷパッケージサポート

ＦＦＴＷパッケージがすでにインストールされている場合は make F=FFTW を実行します。
ＦＦＴＷのインストールについてはここまたはアーカイブ解凍後に作成されるfftw3_install.txtをご覧下さい。

相関処理プログラムのチェック法（チェック用生データのダウンロード）

相関処理プログラムcorおよびfx_corチェック用のK5/VSSPデータファイル（10秒間データ）および予測値ファイルのダウンロードが可能です。

オリジナルデータ諸元:

Ｘ局	鹿島１１ｍアンテナ
Ｙ局	苫小牧１１ｍアンテナ
観測日	２００３年７月１６日　（通日：１９７日）
観測時間	02:40:20 - 02:42:00 (UT)（内10秒間のデータ）
電波星	３Ｃ２７３Ｂ
サンプリングモード	８ＭＨｚ、１ビットサンプリング、４チャンネル
ＲＦ周波数	CH1: 8209.99 MHz USB
	CH2: 8219.99 MHz USB
	CH3: 8249.99 MHz USB
	CH4: 8309.99 MHz USB
ＰＣＡＬ周波数	１０ｋＨｚ

１０秒間データ
- 予測値ファイル(apeXY10.txt)
- Ｘ局K5/VSSPデータファイル(Xk5data.10.dat 約３９ＭＢ)
- Ｙ局K5/VSSPデータファイル(Yk5data.10.dat 約３９ＭＢ)
- cor処理結果(cor_10.gif)
- fx_cor処理結果(fx_cor_10.gif)

チェック法
（１）apeXY10.txt、Xk5data.10.dat、Yk5data.10.datをダウンロードしcorrディレクトリに置く。
（２）corrディレクトリにて
　cor ./apeXY10.txt
　や
　fx_cor ./apeXY10.txt
　を実行する。
（３）実行後、相関処理結果図（corrディレクトリに作成されるポストスクリプトファイルpgplot.ps）と"cor_10.gif"や"fx_cor_10.gif"と比較する。

実際の相関処理

実際の相関処理はここ「相関処理の実際」をクリック

ドキュメント

ダウンロード

更新履歴

2021-11-04
アーカイブ(2021-11-04)公開
2021-02-10
新アーカイブ公開（gcc9.3.0対応）
2020-10-20
apri_calcの予測値計算にVer.2020-09-29の改修時にバグが混入していた（星位置が赤経方向にずれる）。そのバグを修正。なお相関処理結果（予測値＋残差）は正しく求まるので観測量には影響はない。
2020-09-29
apri_calcで計算する予測値が0.1psecレベルでの桁落ちが発生する現象を防ぐための改修を行った。モデル遅延とTaylor係数で求めた遅延の差をプロットするプログラムapri_error_plotを作成した。
2020-09-24
apri_calcにオプション'-4th'を追加。予測値のTaylor係数を4次まで求めることができるようにした。モデル値と予測値値のスキャナ内の最大誤差も表示するようにした。 fx_cor_newでTaylor係数4次まで対応できるよう改修した。スペースVLBI以外や国内基線および大陸間基線でもスキャン長が3分以内の場合はこのオプションを使う必要はない。またKSP相関ファイルフォーマットおよびKOM3出力フォーマットは未対応。
2020-09-11
sdelayのグラフ表示でPP#軸をTIME軸に変更した。
2020-09-09
rdf2m5bを開発。RadioAstron観測生データファイルをMark5Bフォーマットに変換する
2020-09-02
g_pulsarを開発。パルサー周期に同期したパワーのアンサンブル平均表示。全データフォーマットをサポート。
2020-08-18
fx_cor_new, cor_new, fx_cor, corへの'-frmid'オプション（フリンジストップを帯域の中央で行う）の実装。今まで実装されていたが機能していなかったが、デバグを行い機能するようにした（相関強度が約10％上昇する）。
重要：Ver.2020-08-18以前の相関処理では例え'-frmid'を指定してもフリンジ回転はベースバンド周波数で行われている。過去バージョンとの互換性からデフォルトはベースバンドでフリンジストップとするが、デフォルトの変更は cor_head.h（fx_corとcor用）およびcor_head_new.h(fx_cor_newとcor_new用）中の#define DEFLT_FRINGE_STOP_MODEで行うことができる（0: ベースバンド、 1: バンド中央)。
2020-01-17
sdelayで64chのcoutデータを処理できなかった不具合を修正。
2020-01-16
アーカイブ2018-09-27以降のfx_cor_newでIQ VDIFデータ処理が出来なくなっていたバグを改修した。
apri_calcでSKEDファイルを入力とした場合 -threadY 0 オプションがエラーを起こすバグを修正した。

更新日２０２１年１１月４日