DPAM: 空中磁気探査計測データ処理プログラム

English

  ・先頭232バイトのヘッダに続いて，データブロックが羅列する．
  ・ヘッダ(232バイト)の内容は，29個の8バイト実数の羅列であり，その内容は，
        Version番号，日時(Windows形式)，Julian日，年，月，日，
        時，分，秒，1秒未満，センサタイプ(0/1)，Flight番号，計測周期(秒)，
        および 16個のADC較正係数 (勾配とオフセット ×8ch)．
  ・データブロックのブロック長は，68バイトまたは98バイトであり，
    第4バイトが0のとき68バイト長，0でないとき98バイト長．
  ・データブロックの第5-68バイトの内容は，
        Fiducial番号，System時刻，t200，全磁力値
    を示す4個の8バイト実数と，
        8チャネルADCのデータ(Fluxgate 3成分，電波高度，気圧高度，ほか3成分)
    を示す8個の4バイト実数からなる．
    （「t200」の項目は，ある種のFiducial番号と見られるが，詳細は不明．） 
  ・データブロックの第69-98バイトの内容は，GPSの
        UTC時刻(hhmmss.ssss)，緯度値(dddmmss.ssss)，経度値(dddmmss.ssss)
    の「整数部のみ10進化擬似60進形式」の値(秒単位)を示す8バイト実数と，
    高度値の4バイト実数，測位法コード(1バイト)，および 使用衛星数(1バイト)
    からなる．

/AMDAQ     ver.2014-09-19  by  T.Nakatsuka (GSJ)
/DateTime:  2014-11-26 14:53:44  (Logging start)
/CyclingInterval:  100 mSec  (60 mSec Gate-time)
/GyromagneticRatioUsed:  3.498577 Hz/nT (Scintrex)
M 53624.35 00027.11 45451.232 -3.657+2.026+1.548+0.010+4.243+0.010+0.005+0.005
M 53624.45 00027.21 45451.061 -3.652+2.036+1.548+0.000+4.238+0.010+0.005+0.005
M 53624.55 00027.31 45451.150 -3.652+2.031+1.548+0.005+4.248-0.005+0.005+0.005
M 53624.65 00027.41 45450.871 -3.657+2.026+1.548+0.000+4.248+0.005+0.005+0.005
S 53624.59 00027.35$GPGGA,055345.00,3515.0110,N,13655.3893,E,1,16,0.8,15.38,M,37.80,M,,*55^M        
M 53624.74 00027.51 45451.003 -3.652+2.036+1.543+0.010+4.248+0.010+0.005+0.005
M 53624.85 00027.61 45450.811 -3.657+2.031+1.543+0.005+4.258+0.005+0.005+0.005
M 53624.95 00027.71 45451.083 -3.657+2.031+1.548+0.005+4.238+0.010+0.005+0.005
M 53625.05 00027.81 45451.211 -3.657+2.031+1.543-0.005+4.243-0.005+0.005+0.005
          :
          :
先頭4行は，ヘッダコメント．
Mレコードは，'M'，受信時刻(秒)，Fiducial番号，全磁力値，8個のADC値．
Sレコードは，'S'，受信時刻(秒)，Fiducial番号，GPS受信機による $GPGGAレコード
  （'$GPGGA'，UTC時刻，緯度値，'N'， 経度値，'E'，測位法コード，使用衛星数，
    水平精度HDOP，アンテナ高，'M'，ジオイド高，'M'，，'*' チェックサム '^M'）．

/AMBDAQ     by  T.Nakatsuka (GSJ)  ver.2007NOV22
/DateTime:  2007-11-23 09:12:04  (Logging start)
          :
          :
36411.16,0,  46741.827,     0.000,10:06:52.45,11/23/07,  0
36411.26,0,  46743.765,     0.000,10:06:52.55,11/23/07,  0
36411.36,0,  46745.239,     0.000,10:06:52.65,11/23/07,  0
36411.46,0,  46746.836,     0.000,10:06:52.75,11/23/07,  0
36411.56,0,  46748.750,     0.000,10:06:52.85,11/23/07,  0
36411.66,0,  46750.397,     0.000,10:06:52.95,11/23/07,  0,36411.54,$GPGGA,010653.00,3526.0033,N,13850.1248,E,1,08,1.2,1206.08,M,38.12,M,,*51
36411.76,0,  46751.727,     0.000,10:06:53.05,11/23/07,  0
36411.86,0,  46753.094,     0.000,10:06:53.15,11/23/07,  0
36411.96,0,  46754.179,     0.000,10:06:53.25,11/23/07,  0
36412.06,0,  46755.981,     0.000,10:06:53.35,11/23/07,  0
36412.16,0,  46757.653,     0.000,10:06:53.45,11/23/07,  0
          :
          :
先頭2行は，ヘッダコメント．
以後の各レコードは，コンマ区切りの7項目
（受信時刻(秒)，0，全磁力値1，全磁力値2，時刻hh:mm:ss.tt，日付MM/DD/YY，0）
からなり，さらにGPSデータがあるとき，続いて 受信時刻(秒) と GPS受信機による
$GPGGAレコード が置かれる．($GPGGAレコードの内容は，上の「AMDAQ 収録生データ」
の項を参照のこと．)

6         0.000        0.000 14:50:47.500 11/30/16       7599           1           1   5          0
0   46214.218      0.000 14:50:47.500 11/30/16   0
21 $GPGGA,055047.00,3533.4110,N,14023.9125,E,1,15,0.7,3.28,M,33.30,M,,*6F 14:50:47.400 11/30/16  65 0
0   46214.335      0.000 14:50:47.400 11/30/16   0
0   46214.455      0.000 14:50:47.300 11/30/16   0
0   46214.646      0.000 14:50:47.200 11/30/16   0
0   46214.663      0.000 14:50:47.100 11/30/16   0
0   46214.932      0.000 14:50:47.00 11/30/16   0
0   46215.002      0.000 14:50:46.900 11/30/16   0
0   46215.400      0.000 14:50:46.800 11/30/16   0
0   46215.644      0.000 14:50:46.700 11/30/16   0
0   46215.778      0.000 14:50:46.600 11/30/16   0
0   46216.190      0.000 14:50:46.500 11/30/16   0
21 $GPGGA,055046.00,3533.4110,N,14023.9126,E,1,15,0.7,3.40,M,33.30,M,,*63 14:50:46.400 11/30/16  66 0
0   46216.671      0.000 14:50:46.400 11/30/16   0
0   46217.048      0.000 14:50:46.300 11/30/16   0
0   46217.265      0.000 14:50:46.200 11/30/16   0
          :
          :
0   45980.890      0.000 14:20:37.800 11/30/16   0
0   45980.962      0.000 14:20:37.700 11/30/16   0
0   45980.293      0.000 14:20:37.600 11/30/16   0
0   45979.364      0.000 14:20:37.500 11/30/16   0
21 $GPGGA,052037.00,3533.4097,N,14023.9110,E,1,13,0.8,7.46,M,33.30,M,,*62 14:20:37.500 11/30/16 208 0
0   45978.507      0.000 14:20:37.400 11/30/16   0
0   45977.638      0.000 14:20:37.300 11/30/16   0
0   45974.065      0.000 14:20:37.200 11/30/16   0
6         0.000        0.000 14:20:36.900 11/30/16          0           0           0  36          0

//daq2asc v.2017-10-16 by T. Nakatsuka
//InputDataFilename: /home/nktk/data/Tfield/20150622_085808.daq
//PC-Time data were Shifted by  +0.47 sec.
/AMDAQ     ver.2014-09-19  by  T.Nakatsuka (GSJ)
/DateTime: 20150622 08:58:08 (Logging start)
/CyclingInterval:  100 mSec  (60 mSec Gate-time)
/GyromagneticRatioUsed:  3.498577 Hz/nT (Scintrex)
/   FID     SYSTIME      t200       MAG     FGx     FGy     FGz    Ralt    Balt     AD6     AD7     AD8    LTsec         LAT          LON     ALT  Q  N
        :
        :
  134.9 08:59:40.85   134.910 44022.322  -3.965  -0.605   2.588   0.010   0.233   0.024   0.005   0.000   *         *            *            *    *  *
  135.0 08:59:40.94   135.010 44021.805  -3.901  -0.630   2.559   0.005   0.238   0.000   0.010   0.005   *         *            *            *    *  *
  135.1 08:59:41.05   135.110 44024.093  -3.901  -0.601   2.549   0.005   0.238  -0.005   0.005   0.005 32381.00  34.5741050  135.8166467  265.24  1 17
  135.2 08:59:41.14   135.210 44020.111  -3.887  -0.645   2.520   0.005   0.238   0.020   0.005   0.005   *         *            *            *    *  *
  135.3 08:59:41.25   135.310 44018.997  -3.955  -0.605   2.573   0.000   0.233  -0.015   0.005   0.005   *         *            *            *    *  *
  135.4 08:59:41.35   135.410 44022.616  -3.970  -0.654   2.578   0.010   0.243  -0.015   0.005   0.005   *         *            *            *    *  *
  135.5 08:59:41.44   135.510 44025.231  -3.960  -0.615   2.593   0.010   0.233  -0.015   0.005   0.005   *         *            *            *    *  *
  135.6 08:59:41.55   135.610 44018.321  -3.950  -0.571   2.603  -0.005   0.238   0.015   0.005   0.000 32381.50  34.5741050  135.8166467  265.23  1 17
  135.7 08:59:41.64   135.710 44020.228  -3.950  -0.576   2.520   0.000   0.233  -0.010   0.005   0.005   *         *            *            *    *  *
  135.8 08:59:41.75   135.810 44019.916  -3.950  -0.649   2.515   0.000   0.233  -0.010   0.005   0.000   *         *            *            *    *  *
        :
        :
  153.5 08:59:59.44   153.510 44394.005  -3.975  -0.591   2.588   0.000   0.224   0.010   0.005   0.005   *         *            *            *    *  *
  153.6 08:59:59.55   153.610 44933.933  -3.960  -0.645   2.510   0.005   0.243  -0.005   0.005   0.000 32399.50  34.5741050  135.8166483  265.30  1 17
  153.7 08:59:59.64   153.710 45287.924  -3.940  -0.625   2.559   0.000   0.194  -0.010   0.010   0.005   *         *            *            *    *  *
  153.8 08:59:59.75   153.810 44750.706  -3.970  -0.610   2.603   0.005   0.258   0.005   0.005   0.000   *         *            *            *    *  *
  153.9 08:59:59.85   153.910 44848.155  -3.901  -0.605   2.515   0.005   0.263  -0.010   0.005   0.005   *         *            *            *    *  *
  154.0 08:59:59.94   154.010 44516.191  -3.916  -0.635   2.529   0.005   0.248   0.024   0.005   0.005   *         *            *            *    *  *
  154.1 09:00:00.05   154.110 44512.921  -3.970  -0.625   2.607   0.000   0.248  -0.010   0.005   0.005 32400.00  34.5741050  135.8166483  265.31  1 17
  154.2 09:00:00.14   154.210 44779.757  -3.901  -0.552   2.612   0.005   0.175   0.034   0.005   0.005   *         *            *            *    *  *
        :
        :
  156.5 09:00:02.44   156.510 43976.703  -3.882  -0.566   2.529   0.005   0.277  -0.015   0.005   0.005   *         *            *            *    *  *
  156.6 09:00:02.55   156.610 43987.158  -3.887  -0.615   2.515   0.010   0.121   0.034   0.005   0.005 32402.50  34.5741017  135.8166417  264.79  1 16
  156.7 09:00:02.64   156.710 43995.208  -3.906  -0.649   2.622   0.000   0.302   0.020   0.005   0.005   *         *            *            *    *  *
  156.8 09:00:02.75   156.810 44004.666  -3.979  -0.562   2.627  -0.005   0.297  -0.010   0.005   0.000   *         *            *            *    *  *
  156.9 09:00:02.85   156.910 44007.559  -3.950  -0.571   2.593   0.000   0.346   0.015   0.005   0.005   *         *            *            *    *  *
  157.0 09:00:02.94   157.010 43997.328  -3.892  -0.664   2.515   0.005   0.341   0.005   0.005   0.005   *         *            *            *    *  *
  157.1 09:00:03.05   157.110 44007.911  -3.916  -0.654   2.578   0.000   0.414   0.010   0.005   0.005 32403.00  34.5741017  135.8166417  264.79  1 16
  157.2 09:00:03.14   157.210 44012.890  -3.945  -0.654   2.520   0.000   0.331   0.010   0.005   0.005   *         *            *            *    *  *
        :
        :

/Base:  46490               ← 基準磁力値(一般に静穏日の夜間値) (*)
/Date: 20030215             ← 日付（第7カラム以降に yyyymmdd 形式の8桁）(*)
130008 464795
130023   46480
130038  464795
130053 4648030              ← 定点磁力計収録データそのままを想定．
130108 46480.5              ← HHMMSS の6桁, 1桁以上のスペース, 
130123 46480.71             ←   磁力値(1nT／0.1nT／0.01nT 単位のいずれか)
130138 464806.1
130153 4648052.
   .......

        (*) 基準磁力値, 日付 の指定は，ファイル途中に指定されると
            その位置から後ろに向って有効となる.

   12:51:11.000 35.0960951 137.6353737 983.83
   12:51:12.000 35.0963192 137.6353688 985.57
   12:51:13.000 35.0965443 137.6353639 987.19
        .......
  14:14:46.00  36.4404884  138.5358380  1408.51
  14:14:47.00  36.4404885  138.5358380  1408.50
  14:14:48.00  36.4404886  138.5358380  1408.47
        .......

     すべての行は，フリーフォーマットのGPS位置データであり，時刻(Localtime)，
       緯度(度単位)，経度(度単位)，高度(m単位) で構成される．
       時刻は，"HH:MM:SS.tt" の形式でなければならない．

=030217f1.asc  030217f1.pnav
220      095250 100100    0
210      100330 101000  180
200      101210 101950    0
     .......
140      111210 112050    0
130      112340 113120  180
=030217f2.asc  030217f2.pnav
120      121300 122010    0
110      122235 122840  180
100      123105 123850    0
     .......

        先頭桁が '=' の行は，データファイルの指定を示し，共通 ASCII測定data と
          GPS測位データ（通常はPostNav処理結果）の両ファイル名を与える．
          （xldam に対しては，測位データファイル名はなくてよい．）
        それ以外は，個別の測線情報データであり，行頭から，
          測線名，開始／終了時刻（hhmmss.tt の形式，通常 .tt は省略）と，
          測線飛行方位角（度単位） を，空白区切りで並べる．磁力センサーの位置
          ズレ補正を行わない場合は「測線飛行方位角」の指定がなくても構わない．

　測線データファイル形式（AM Line data type）には，
　　1) DPAM測線データファイル（DPAM Line data）
　　2) GSJ測線データファイル（AMDB-GSJ line data）
　　3) NEDO統合測線データファイル（AMDB-NEDO line data）
　　4) 標準測線データファイル（StdLIN data）
　　0) 一般ランダム点データ（generic random point data）
のバリエーションがある．
　2 および 3 のデータ形式では，原則的に，緯経度値は世界測地系での値を 高度値は
ジオイドからの標高を用いる．しかし，その他の形式では 測地系・高度基準とも規定
されておらず，そのデータの利用にあたっては測地系等の種別を指定する必要がある．
　また，0 は一般的な非固定フォ－マットのデータファイルを扱うための形式であり，
緯経度値の単位や各データの配列順をも指定する必要がある．そのときデータ行は数値
データのみで構成されている必要があるが，特例として ":" は空白と同等のデータの欄
区切りとみなされる．
　なお，いずれの場合も，(測線を構成する)個別点データは 1点/1行 のデータであり，
その先頭欄は 空白(もしくは数値先頭桁) である．また，測線の区切りや属性情報を
示す測線開始指示行は，先頭欄が "&", "#" または "%" となっており，その他のコメ
ント情報として，先頭欄が "#" の行が ファイル先頭に置かれることがある．
（従って，先頭欄が "#", "&", "%" の行を無視すると，点データの羅列となる．）

# Areaname: Kobe-Kyoto
# Survey Date: 1995.12.07-12.27
&A-01
  2079.02221N  8116.27649E   277.87m   -45.15nT
  2079.04052N  8116.31640E   278.58m   -44.66nT
  2079.05883N  8116.35641E   279.19m   -44.47nT
        .......
  2087.39584N  8134.25592E   275.40m   -48.38nT
  2087.41585N  8134.29643E   275.41m   -53.69nT
&    C-2r
  2088.27126N  8134.37994E   279.12m   -44.90nT
  2088.25637N  8134.33845E   279.33m   -40.81nT
  2088.24078N  8134.29646E   279.64m   -40.12nT
        .......

     先頭桁が，'#' の行は，単なるコメント行であり，通常，ファイルの先頭にのみ
       置かれる．（一連の個別点データ(測線データ)の途中に，置いてはならない．）
     先頭桁が，'&' または '%' の行は測線の開始を示す．
       第2～9桁に，測線名を記述する．それ以外の記述内容とその形式に制約は
       設けない．（開始時刻や含まれる測点数などを記述する場合が多い．）
       ファイルの末尾に，含まれる測点数が 0 の測線を表示して，データの末尾を
       明示する場合がある．
     それ以外は，個別の点データであり，緯度(分単位)・経度(分単位)・高度(m単位)
       ・残差磁力値(nT単位)の値が，
          format(1x, f11.5, 1hN,1x, f11.5, 1hE,1x, f8.2, 1hm,1x, f8.2, 2hnT)
       の形式（計47桁）で入っている．
     北緯を意味する 'N'，東経を意味する 'E' の文字を含むが，南半球・西半球の
       場合には 緯度値・経度値をマイナス値で表現する．
     入力データの場合は数値部分の桁数は任意である．

# Areaname: Ootoge
# Survey Date: 2003.02.17
&220       20030217  95250.00 100100.00
  418860 20030217  95250.02  3  35.0885765  137.7122326 1033.28 46445.27   -50.13  -3.535   2.783   1.099  35570.02
  418870 20030217  95250.09  3  35.0885932  137.7122356 1033.31 46445.02   -50.39  -3.525   2.783   1.108  35570.09
  418880 20030217  95250.17  3  35.0886091  137.7122389 1033.34 46445.90   -49.51  -3.545   2.781   1.084  35570.17
        .......
  494660 20030217 100059.89  3  35.2059824  137.7122510 1258.37 46440.41  -115.48  -3.560   2.891   0.713  36059.89
  494670 20030217 100059.95  3  35.2059986  137.7122531 1258.30 46439.93  -115.95  -3.555   2.874   0.708  36059.95
&210       20030217 100330.00 101000.00
  517780 20030217 100330.09  3  35.2047662  137.7067634 1247.53 46418.68  -138.34  -3.281  -3.225  -0.249  36210.09
  517790 20030217 100330.19  3  35.2047375  137.7067660 1247.46 46418.48  -138.52  -3.296  -3.230  -0.273  36210.19
  517800 20030217 100330.29  3  35.2047093  137.7067705 1247.39 46418.52  -138.47  -3.286  -3.210  -0.288  36210.29
        .......

     先頭桁が，'#' の行は，単なるコメント行であり，通常ファイルの先頭にのみ
       置かれる．（一連の個別点データ(測線データ)の途中に，置いてはならない．）
     先頭桁が，'&' または '%' の行は測線の開始を示し，測線名・日付yyyymmdd
       ・開始／終了時刻が記述される．時刻は，HHMMSS.tt 形式の 倍精度実数値．
     それ以外は，個別の点データであり，
           1) Fiducial number, 2) Date (yyyymmdd), 3) Time (HHMMSS.tt), 
           4) データ種別(*), 5) 緯度(Deg.), 6) 経度(Deg.), 7) 高度(m), 
           8) 磁力値(nT), 9) IGRF残差値(nT), 10-12) Fluxgate 3成分値(V), 
           13) Localtime秒(00:00からの秒数) が，
             format(i8, 1x,i8, 1x,f9.2, 1x,i2, 1x,f11.7, 1x,f12.7, 1x,f7.2,
            *       2(1x,f8.2), 3(1x,f7.3), 1x,f9.2)
       の形式（計115桁＋改行コード）で入っている．
       ecomp, fcomp の出力では，さらに続いて '(4(1x,f8.2))' 形式で，
           14) 補正前残差値(tres), 15) 機体磁気補正量(corr), 
           16) ランダム成分値(rand), 17) リニアトレンド分(trend) 
       の各データ [tres ＝ corr ＋ rand ＋ trend] が入る．
       (*) ここに，データ種別(k) の意味は次のとおりである．
            k = 4～7 :   5) 緯度, 6) 経度, 7) 高度 が RealTimeGPS による
            k = 0～3 :   5) 緯度, 6) 経度, 7) 高度 が PostNav処理 による
            k = 2,3,6,7 :   8) 磁力値, 9) IGRF残差値 とも 日変化が未補正
            k = 0,1,4,5 :   8) 磁力値, 9) IGRF残差値 とも 日変化補正済
            k = 1,3,5,7 :   8) 磁力値, 9) IGRF残差値 とも 機体磁気が未補償
            k = 0,2,4,6 :   8) 磁力値, 9) IGRF残差値 とも 機体磁気補償済

        .......
 20051012 09:14:51.10  99999.99  46949.90
 20051012 09:14:51.20  99999.99  46949.89
 20051012 09:14:51.30  99999.99  46949.85
 20051012 09:14:51.40  46781.33  46949.85
 20051012 09:14:51.50  99999.99  46949.79
 20051012 09:14:51.60  99999.99  46949.77
 20051012 09:14:51.70  46795.16  46949.70
 20051012 09:14:51.80  46850.40  46949.78
 20051012 09:14:51.90  46851.73  46949.73
 20051012 09:14:52.00  46851.76  46949.70
        .......

     各データは，Date (yyyymmdd), Time (HH:MM:SS.tt),
       第1センサ全磁力値(nT), 第2センサ全磁力値(nT) が
          format(1x, i8, 1x, i2,1h:,i2,1h:,f5.2, 1x, f9.2, 1x, f9.2)
       の形式（計41桁＋改行コード）で入っている．

&F13       20051012 142618.00 152530.00
 1893030 20051012 142618.00 80  36.4324162  138.4260987 2285.58 46934.03   -58.27 46936.09   -56.09    -2.06
 1893040 20051012 142618.10 80  36.4323834  138.4261012 2285.41 46934.69   -57.58 46936.12   -56.03    -1.43
 1893050 20051012 142618.20 80  36.4323596  138.4261236 2285.44 46935.28   -56.97 46936.75   -55.38    -1.47
        .......
 2248210 20051012 152529.80 80  36.4200127  138.5840531 2700.97 47043.23   112.45 47048.76   118.09    -5.53
 2248220 20051012 152529.90 80  36.4200489  138.5840760 2700.40 47043.27   112.47 47048.77   118.09    -5.50
 2248230 20051012 152530.00 80  36.4200661  138.5840904 2699.93 47042.84   112.02 47048.52   117.81    -5.67
&F23a      20051013 125215.00 130130.00
 1388080 20051013 125215.00 80  36.4079470  138.5286123 2771.76 47531.58   593.48 47536.52   598.54    -4.94
 1388090 20051013 125215.10 80  36.4079393  138.5285881 2772.09 47531.73   593.64 47536.70   598.73    -4.98
 1388100 20051013 125215.20 80  36.4079166  138.5285668 2772.32 47531.83   593.75 47536.73   598.77    -4.90
 1388110 20051013 125215.30 80  36.4079055  138.5285373 2772.65 47531.95   593.87 47536.86   598.90    -4.91
        .......

     先頭桁が，'&' の行は測線の開始を示し，測線名・日付yyyymmdd・開始／終了
       時刻が記述される．時刻は，HHMMSS.tt 形式の 倍精度実数値．
     それ以外は，個別の点データであり，
           1) Fiducial number, 2) Date (yyyymmdd), 3) Time (HHMMSS.tt), 
           4) データ種別(*) [傾度測定HGAM測線データ]
           5) 緯度(Deg.), 6) 経度(Deg.), 7) 高度(m), 
           8) 磁力値1(nT), 9) IGRF残差値1(nT), 10) 磁力値2(nT), 
           11) IGRF残差値2(nT), 12) 磁力値差(磁力値1－磁力値2)(nT) が，
             format(i8, 1x,i8, 1x,f9.2, 1x,i2, 1x,f11.7, 1x,f12.7, 1x,f7.2,
            *       5(1x,f8.2))
       の形式（計108桁＋改行コード）で入っている．
       (*) データ種別(k) は，緯経度値が PostNAV処理によるものであれば k＝80，
          RealTimeGPS による場合は k＝84 である．
          (いずれも日変化補正済・機体磁気補正不要)
     このデータは，DPAM測線データファイルの 1変種の形態をとっており，磁力値1
       を扱う場合は，DPAM測線データと同様に取り扱うことができる．磁力値2 を扱
       う場合には，一般ランダム点データとして処理すればよい．