Question

官方文档

mgd77list

简介

从 MGD77 文件中提取数据

mgd77list 读取 MGD77 文件并生成 ASCII 或二进制表数据。MGD77 文件包括测线信息，比如 ID，时间，位置，控制代码和地球物理观测，其中地球物理观测数据可能包括：重力，磁，水深等，因此还可能包括对这些观测数据的改正，例如：Eotvos 改正和 diurnal 改正。MGD77+ NetCDF 文件还可能包含其他用户自定义信息，见 mgd77info 的 -C 选项；若想添加自定义列，参考 mgd77manage 。用户可以提取包括上面信息以及下面信息的任意组合，其中包括 8 个计算量（距离、航向、航向变化、速度、Carter 改正、Eotvos 改正以及重力和磁的全球参考场）和时间的子部分（年、月、日、时、分、秒）以及 GEODAS id 以及权重（见 -W ）等。通过指定沿轨的累积时间或距离或者指定地理区域，可以选定特定的数据。提取出来的数据可以使用一定的逻辑判断来进一步筛选。如果需要提取多个测线的信息，则多个测线之间由段信息分隔。

语法

gmt mgd77list GEODAS-ids

NGDC-ids

可以是 5 种说明符种的一种或者几种:

1. 8 个字符的 GEODAS ID，例如：01010083，JA010010

2. 2 个字符的机构 ID，将会返回所有该机构的测线

3. 4 个字符代码，返回对应机构和测量船的测线，机构和测量船的代码见 mgd77info -L

4. =list ，list 文件中包含 GEODAS ID 列表，每行一个记录

5. 什么都不指定，返回所有的测线

如果不指定文件扩展名，则搜索下面 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 ，单位见 -Ac|d|f|m|tcode[+f] 默认情况下，改正后的测深（ depth ），残余磁异常（ mag ），自由空气重力异常（ faa ），以及 Carter 深度改正（ carter ）如果被 从 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 重新计算的异常 当通过 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 为 MGD77 (Marine Geophysical Data Exchange Format) 格式见 http://www.ngdc.noaa.gov/mgg/dat/geodas/docs/mgd77.txt. IGRF 见 https://www.ngdc.noaa.gov/IAGA/vmod/igrf.html 相关模块 mgd77convert, mgd77header, mgd77magref mgd77manage, mgd77path, mgd77track, mgd77sniffer, mgd77track