mgd77list
- 官方文档:
- 简介:
从 MGD77 文件中提取数据
mgd77list 读取 MGD77 文件并生成 ASCII 或二进制表数据。MGD77 文件包括测线信息，比如 ID，时间，位置，控制代码和地球物理观测，其中地球物理观测数据可能包括：重力，磁，水深等，因此还可能包括对这些观测数据的改正，例如：Eotvos 改正和 diurnal 改正。MGD77+ NetCDF 文件还可能包含其他用户自定义信息，见 mgd77info 的 -C 选项；若想添加自定义列，参考 mgd77manage 。用户可以提取包括上面信息以及下面信息的任意组合，其中包括 8 个计算量（距离、航向、航向变化、速度、Carter 改正、Eotvos 改正以及重力和磁的全球参考场）和时间的子部分（年、月、日、时、分、秒）以及 GEODAS id 以及权重（见 -W ）等。通过指定沿轨的累积时间或距离或者指定地理区域，可以选定特定的数据。提取出来的数据可以使用一定的逻辑判断来进一步筛选。如果需要提取多个测线的信息，则多个测线之间由段信息分隔。
语法
gmt mgd77list GEODAS-ids -Fcolumns[,logic][:bittests] [ -Ac|d|f|m|tcode[+f] ] [ -DA|astartdate ] [ -DB|bstopdate ] [ -E ] [ -Gastartrec ] [ -Gbstoprec ] [ -Ia|c|m|t ] [ -L[corrtable] ] [ -Nd|sunit ] [ -Qa|c|vmin/max ] [ -Rregion ] [ -Sastartdist ] [ -Sbstopdist ] [ -T[m|e] ] [ -V[level] ] [ -Wweight ] [ -Zn|p ] [ -bobinary ] [ -hheaders ] [ -jflags ] [ --PAR=value ]
必选选项
- NGDC-ids
可以是 5 种说明符种的一种或者几种:
8 个字符的 GEODAS ID，例如：01010083，JA010010
2 个字符的机构 ID，将会返回所有该机构的测线
4 个字符代码，返回对应机构和测量船的测线，机构和测量船的代码见 mgd77info -L
=list ，list 文件中包含 GEODAS ID 列表，每行一个记录
什么都不指定，返回所有的测线
如果不指定文件扩展名，则搜索下面 4 个扩展名。搜素顺序为：
MGD77+ (“.nc”)
MGD77T (“.m77t”)
MGD77 (“.mgd77”)
普通文本文件 (“.dat”)
使用 -I 选项可以略过其中的某种或几种格式。
搜索文件时首先在当前目录查找，如果未设置 $MGD77_HOME, 则在默认的 $GMT_SHAREDIR/mgd77 目录查找，如果设置了 $MGD77_HOME，则在其下的 $MGD77_HOME/mgd77_paths.txt 文件中保存的路径中查找。
- -Fcolumns[,logic][:bittests]
指定输出列 columns ，columns 中的多列必须使用逗号分隔，列名为缩写，按给定的顺序输出。如果列名为大写字母，则表示该列输出中不能存在 NaN。除非经过设置，否则，输出的 ASCII 文件将使用 FORMAT_FLOAT_OUT 指定的格式。所有的可用的列名的缩写如下表：
列名 |
含义 |
|---|---|
recno |
计数器 |
dst |
数字记录类型，通常为 3 或者 5 （Y2K-compliant 测量） |
id |
测量 ID 字符串 |
ngdcid |
8 字符的 NGDC（现在为 NCEI）测线 ID 字符串（通常是文件前缀） |
time |
时间，从绝对时间（默认），相对时间，以及小数年时间中选择，其中三种时间的参
数设置分别为：
- 绝对时间：FORMAT_DATE_OUT 和 FORMAT_CLOCK_OUT
- 相对时间：FORMAT_FLOAT_OUT 和 TIME_SYSTEM 或者（
TIME_EPOCH 和 TIME_UNIT ）
- 小数年时间：FORMAT_FLOAT_OUT
|
lon |
经度，格式由 FORMAT_GEO_OUT 指定 |
lat |
纬度，格式由 FORMAT_GEO_OUT 指定 |
twt |
双程传播时间，单位 s |
depth |
改正后的测深，单位为 m，正值，表示海面下深度 |
mtf1 |
传感器 1 的总磁场强度，单位 nTesla |
mtf2 |
传感器 2 的总磁场强度，单位 nTesla |
mag |
残余磁异常，单位 nTesla |
gobs |
重力，单位 mgal |
faa |
自由空气重力异常，单位 mgal |
ptc |
位置类型代码，1 为固定位置，3 为内插位置，9 为未指定 |
bcc |
测深改正代码。
- 01-55 Mattews zone 改正，
- 59 Mattews 改正，但未指定 zone，
- 60 S，Kuwahara formula for T-S，
_ 61 Wilson formula for T-S，
- 62 Del Grosso formula for T-S，
- 63 Carter’s tables，88 其他，在头部记录中，
- 99 未指定
|
btc |
测深类型代码，1 观测值，3 内插值，9 未指定 |
msens |
评估残余磁场的传感器，1 第一个或者 leading 传感器，2 第二个或 trailing
传感器，9 未指定
|
msd |
磁传感器深度（或高度），单位 m，向下为正 |
diur |
磁 diurnal 改正，单位 nTesla |
eot |
Eotvos 改正，单位 mgal |
sln |
Seismic Line Number 字符串 |
sspn |
Seismic Shot Point Number 字符串 |
nqc |
导航质量代码，5 怀疑（测量机构），6 怀疑（NCEI）, 9 确定没问题 |
year |
年 |
month |
月 |
day |
日 |
hour |
时 |
min |
分 |
sec |
秒 |
dist |
沿轨距离，计算方法见 -j ，单位见 -N [km] |
az |
方位角，从北方向顺时针计算，单位 度 |
cc |
航向变更，单位 度 |
vel |
速度，单位见 -N [m/s] |
weight |
该数据集权重，见 -W |
carter |
Carter 深度改正（必须有 twt 的前提下） |
igrf |
国际地磁参考场（总场），单位 nTesla |
ngrav |
国际重力参考场，即正常重力，单位 mgal。
- 默认使用头部记录中的 Gravity Theoretical Formula Code 参数
- 如果该参数未设置就使用 IGF80。也可以使用
-Af 指定
|
ceot |
从导航信息中计算 Eotvos 改正，
E = 7.5038 * V * cos(lat) * sin(az) + 0.004154 * V^2
|
处上述列名外，下表中还列出了一些其他可识别的参数：
参数 |
含义 |
|---|---|
mgd77 |
输出所有的 27 个 MGD77 字段 |
mgd77t |
输出所有的 26 个 MGD77T 字段 |
geo |
输出 10 个字段，分别为 time ，lon，lat 和 7 个地球物理观测 twt，depth，mtf1，mtf2，mag，gobs，和 faa |
dat |
和 mgd77t 相同，但是使用 plain 表文件顺序 |
all |
返回文件中所有数据列 |
allt |
与 mgd77t 相同，但是时间项输出为 日期-时间 格式 |
在这些集合后面使用 + 追加 dist，azim，cc，vel，和 weight 将可输出对应信息。对上面提取的结果还可进行逻辑测试，其中 logic 由一个或多个逗号分隔的 parOPvalue 指令组成，其中 par 是上面列出的列名中的一个，OP 为逻辑运算符（<，<=，=，!=，>=，>，|），value 为常数。符点数进行数值比较，字符参数进行词法比较。位比较（ | ）中 value 中必须至少有一个位的结果是 1；且最少有一个逻辑通过测试测可以输出结果。但如果使用大写列名，则必须所有的逻辑测试通过才可输出。
注 ：(1) 指定对某列进行逻辑判断并不意味着该列被输出，输出必须经过指定。(2) 某些运算符是特殊的 Unix 字符，建议在 -F 整个参数外添加引号。(3) 逻辑测试仅适用于观测数据，派生数据，例如距离和速度等，必须使用 -D ，-Q ，-S 等选项。
对于 MGD77+ 文件，可以追加 :bittests ，即一个冒号，后面为一个或多个逗号分隔的 +-col
项。+ 表示所选位必须是 1（打开），- 表示必须为 0（关闭）。所有的位必须都通过测试。默认情况下，具有特殊的 MGD77_flags 列的 MGD77+ 文件将使用这些标志，并且将与打开（ON）位相关的观测设置为 NaN，即表示出现错误。: 后面不添加任何信息将关闭此行为。需要注意的是，这些记录标志和用于沿轨的系统改正不同，后者不能使用 -T 选项。
可选选项
- -Ac|d|f|m|tcode[+f]
默认情况下，改正后的测深（ depth ），残余磁异常（ mag ），自由空气重力异常（ faa ），以及 Carter 深度改正（ carter ）如果被 -F 选项选择的话，都会按原样输出；-A 选项可改变该行为。对于上面的每种数据，都有 2-4 种方法来调整数据。通过 c，d 和 f 或者 m 可分别选择上述对应的变量。也可以通过数字代码相加的和表示选择多种处理。例如，-Ac3 将首先尝试使用 -Ac1 表示的方法来估计 Carter 改正，但如果深度为 NaN，将使用 -Ac2 表示的方法。在所有的情况下，如果计算过程中存在 NaN，则会导致结果为 NaN。使用 +f 选项将会重新计算异常，即使文件中异常为 NaN。此外，用户可以使用 -At 为没有时间的测线通过距离和持续的时间来创建伪时间。
- -Ac
设置计算 carter 改正的方法。C(twt) Carter 改正的深度（同时依赖 lon 和 lat ）；U(twt, v) 是未改正的深度（U = twt * v / 2 ），其中的 v 使用 MGD77 头部信息中的 “Assumed Sound Velocity” 参数，如果该速度为无效值，则假定为 1500 m/s，也可以自己设定 v 并追加 m/s 作为单位；TU(depth, v) 是从未改正的深度计算的双程传播时间；TC(depth) 是通过 Carter 改正后得到的深度计算的双程传播时间。下面介绍可选择的选项：
-Ac1[,v] 返回 U(twt, v) 和 depth 的差 [默认]
-Ac2[,v] 返回 U(twt, v) 和 Carter (twt) 的差
-Ac4[,v] 返回未改正的 depth 和 Carter (TU(depth)) 的差
-Ac8[,v] 返回 U(TC(depth), v) 和 depth 的差
- -Ad
确定如何获取 depth 输出：
-Ad1 返回储存在数据集中的 depth [默认]
-Ad2[,v] 返回计算的未改正的深度 U(twt, v)
-Ad4 返回计算的改正后的深度 C(twt)
- -Af
确定如何获取 faa 输出。如果计算中需要 ngrav （即International Gravity reference Field, IGF 或正常重力），则根据 MGD77 标题参数中的 “Theoretical Gravity Formula Code”计算重力异常，如果上述代码不存在或无效，默认设置为 4。可以在下面 4 个代码中选择来设置重力参考场：
1 代表 Heiskanen 1924
2 代表 IGF 1930
3 代表 IGF 1967
4 代表 IGF 1980
可设置的选项如下:
-Af1[,field] 返回储存在数据集中的 faa [默认]。可使用 -F 选项输出 ngrav
-Af2[,field] 返回 gobs 和 ngrav 的差
-Af4[,field] 返回以下结果 gobs + eot - ngrav
-Af8[,field] 返回以下结果 gobs + pred_eot - ngrav
- -Am
确定如果获取 mag 输出。文件中可能存在一个或者两个测量值（mtf1 和 mtf2 ），msens 可表明其中哪个传感器优先。可设置的选项如下：
-Am1 返回储存在数据集中的 mag [默认]
-Am2 返回 mgfx 和 igrf 的差，其中 x 由 msens 指定，可以为 1 或者 2，表示优先使用的传感器，如果未指定则默认为 1
-Am4 返回 mgfx 和 igrf 的差，其中 x 由 msens 指定，可以为 1 或者 2，表示优先使用的传感器，如果未指定则默认为 2
-Amcoffset 使用由于测量位置和船位置不同（导航引起）导致的改正值。通过将总磁场重新差值到船所在的位置得到（可能存在一定的偏移 offset ）。由于插值以及导航质量差可能会导致出现重复点。GMT 已经采取一定措施处理这些情况，但不能保证全部处理干净。内插算法使用 GMT_INTERPOLANT 设置的算法。offset 的默认单位为 m，可以通过 f，k，M 或 u 设置为英尺，km，英里，海里和 survery 英尺。
- -At
使用已知的测线观测时长来计算测线缺失的时间记录，该计算时间记录为伪时间记录
- -Dastartdate
不列出开始日期 startdata 之前收集的数据，startdata 格式为 yyyy-mm-ddT[hh:mm:ss]。使用 -DA 可排除缺失时间记录的数据。[默认输出所有数据]
- -Dbstopdate
不列出截止日期 stopdata 之后收集的数据，stopdata 格式为 yyyy-mm-ddT[hh:mm:ss]。使用 -DB 可排除缺失时间记录的数据。[默认输出所有数据]
- -E
精确匹配，仅输出指定的列 [默认输出至少一个匹配的列]
- -Gastartrec
不列出在 startrec 之前的记录 [默认为 0，即第一个记录]
- -Gbstoprec
不列出在 stoprec 之后的记录 [默认为最后一个记录]
- -Ia|c|m|t
忽略某些数据格式的文件。a|c|m|t 分别表示忽略 MGD77 ASCII、MGD77+ NetCDF、MGD77T ASCII 、tab 分隔的普通文本文件。该选项可以重复使用以忽略多个格式
- -L[corrtable]
使用改正表 corrtable 中的值来做改正[默认使用 $MGD77_HOME 文件夹下的 mgd77_corrections.txt 文件作为改正表]。改正表格式见 改正表
- -Nd|sunit
设置距离和速度的单位。 d 表示距离，s 表示速度，然后可以分别指定距离和速度的单位。
e m 或 m/s，
f feet 或 feet/s
k km 或 kms/hr
m miles 或 miles/hr
n nautical miles 或 knots
u survey feet 或 survey feet/s).
默认地，设置为 -Ndk -Nse （km 和 m/s）
- -Qamin/max
指定可接受的方位角范围 min/max ，只输出该方位角范围内的记录 [0/360]。
- -Qcmin/max
指定航向改变的范围 min/max ，只输出该范围内的记录 [-360/+360]。使用 -QC 可在判断前先取绝对值。
- -Qvmin/max
指定速度范围 min/max ，只输出该范围内的记录 [0/infinity]。
- -Rxmin/xmax/ymin/ymax[+r][+uunit] (more …)
指定数据范围
- -Sastartdist
只列出测线累积距离大于 startdist 的数据，可追加单位。
- -Sbstopdist
只列出测线累积距离小于 stopdist 的数据，可追加单位。
- -T[m|e]
关闭基于 MGD77+ 文件中的改正项的自动调整，用于抵消源机构使用 MGD77+ 导出数据创建 MGD77 时使用的错误的单位等问题，普通 MGD77 文件不受该选项影响。m 表示只用于 MGD77 文件中存在的变量，e 表示只用于扩展的变量，默认两者都包括。注 ：基于记录的 E77 标志（仅限于 MGD77+ 文件）不被视为系统误差改正。相反，这些位标志可以通过 -F 选项中的冒号来控制。
- -Wweight
为数据设置权重。权重的输出必须在 -F 选项中设置。如果使用 blockmean ， blockmedian ，blockmode 模块进行加权平均处理数据，可使用该选项
- -Zn|p
为低于海平面的 depth，carter 和 msd 指定符号，-Zn 指定为负的测深，默认使用 -Zp 表示向下为正
- -V[level] (more …)
设置 verbose 等级 [w]
- -bo[ncols][type][w][+l|b] (more …)
设置二进制输出的数据格式
- -h[i|o][n][+c][+d][+msegheader][+rremark][+ttitle] (more …)
跳过或生成指定数目的头段记录
- -je|f|g (more …)
设置球面距离的计算方式
- -^ 或 -
显示简短的帮助信息，包括模块简介和基本语法信息（Windows下只能使用 -）
- -+ 或 +
显示帮助信息，包括模块简介、基本语法以及模块特有选项的说明
- -? 或无参数
显示完整的帮助信息，包括模块简介、基本语法以及所有选项的说明
- --PAR=value
临时修改GMT参数的值，可重复多次使用。参数列表见 配置参数
示例
从 01010047.mgd77 文件中提取 distance，heading，gravity 以及 bathymetry，将数据限定为 June 3 1971 20:45 以后，且累积观测距离小于 5000 km 的范围内
gmt mgd77list 01010047 -Da1971-06-03T20:45 -Sb5000 -Fdist,azim,faa,depth > myfile.txt
使用 Cruises.lis 文件中列出的所有测线中的自由空气重力异常以便后续作为 blockmean
和 surface 模块的输入，数据必须在限定的区域内，且输出为二进制
gmt mgd77list `cat cruises.lis` -Flon,lat,faa -R-40/-30/25/35 -bo > allgrav.b
从 Cruises.lis 中列出的测线数据中提取深度超过 9000 m 且未插值（ btc != 1 ）的数据的位置
gmt mgd77list `cat cruises.lis` -F"depth,DEPTH>9000,BTC!=1" > really_deep.txt
从所有的 MGD77+ NetCDF 文件中提取 dist，faa，和 grav12_2 ，要求这些字段的值不能为 NaN，测线文件的 ID 位于 Cruises.lis
gmt mgd77list `cat cruises.lis` -E -Ia -F"dist,faa,grav12_2,depth<3000" > \
shallow_grav.txt
与上例使用相同的数据，但只提取 faa 的位标志为 BAD 的数据
gmt mgd77list `cat cruises.lis` -E -Ia -F"dist,faa,grav12_2:+faa" > bad_grav.txt
输出所有测线中的 lon，lat，mag 和 faa，根据最新的参考字段计算 mag 和 faa 的残余量
gmt mgd77list `cat cruises.lis` -Flon,lat,mag,faa -Af2,4 -Am2 > data.txt
重新计算的异常
当通过 -A 选项显式地或通过 MGD77+ 文件中 E77 元数据隐式地指定重新计算异常时，GMT 只计算原始异常中不为 NaN 的异常。这样的做的原因是，许多异常列包括手动编辑的修改，这些结果不能被重新计算。
IGRF (International Geomagnetic Reference Field)
IGRF 基于英国地质调查局的 Susan Macmilan 编写的 Fortran 程序，由 Algrave 大学的 Joaquim 通过 f2c 转换到 C 语言，最终 Paul Wessel 改成 GMT 风格。
IGF (International Gravity reference Field)
这里给出代码中使用的公式
(1) g = 978052.0 * [1 + 0.005285 * sin^2(lat) - 7e-6 * sin^2(2*lat) + 27e-6 * cos^2(lat) * cos^2(lon-18)]
(2) g = 978049.0 * [1 + 0.0052884 * sin^2(lat) - 0.0000059 * sin^2(2*lat)]
(3) g = 978031.846 * [1 + 0.0053024 * sin^2(lat) - 0.0000058 * sin^2(2*lat)]
(4) g = 978032.67714 * [(1 + 0.00193185138639 * sin^2(lat)) / sqrt (1 - 0.00669437999013 * sin^2(lat))]
如果您发现上述公式存在问题，请向 GMT 官方 报告。
改正表
改正表是一个 ASCII 文件，其中包含实施改正所需的系数和参数。文件头部可以存在以 # 为开头的注释，改正记录的格式如下
cruiseID observation correction
其中 cruiseID 为测线 ID，NCEI 定义的前缀，observation 为 -F 选项里列名表中的某个列名， correction 由一项或者多项组成，所有的项在输出观测值前都被扣除，每项的格式为
factor[*[function]([scale](abbrev[-origin]))[^power]]
其中方括号中的值是可选的，除不同项之间，其他地方不允许有空格。factor 是基函数的幅值，可选的 function 可以是 sin，cos，或者 exp。可选的 scale 和 origin 用于对参数进行一定的转换。 abbrev 参数是上面列出的列名之一。如果 origin 为 T ，则将其替换为文件中第一个 abbrev 值（通常同于时间）。如果第一条记录为 NaN，则将 origin 设置为 0。在乘以幅度前，可也对整个表达式设置一定的幂指数 power 。
下面是对测线 99999999 的虚构数据的改正示例，三行分别表示：对深度 depth 进行 Carter 改正，对自由空气重力异常 faa 移除线性趋势，对磁异常 mag 进行依赖航向和纬度的改正，对重力 gobs 需要添加 10 mGal 的增量
99999999 depth 1.0*((carter))
99999999 faa 14.1 1e-5*((time-T))
99999999 mag 0.5*cos(0.5*(azim-19))^2 1.0*exp(-1e-3(lat))^1.5
99999999 gobs -10
参考
MGD77 (Marine Geophysical Data Exchange Format) 格式见 http://www.ngdc.noaa.gov/mgg/dat/geodas/docs/mgd77.txt.
相关模块
mgd77convert, mgd77header, mgd77magref mgd77manage, mgd77path, mgd77track, mgd77sniffer, mgd77track