grdredpol
- 官方文档:
- 简介:
计算连续化极（Continuous Reduction To the Pole，RTP）
grdredpol 读取一个磁异常文件并计算化极后的异常。 This anomaly is the one that would have been produce if the bodies were magnetized vertically and the anomalies were observed at the geomagnetic pole. 标准的 RTP 程序假定磁化方向和地磁场在整个区域是一致的。这些假设对于小区域 来说是合理的，但是大区域不成立。
本模块中的方法是同时在频域和空间域中计算的。这里，将一个大区域分成多个小 区域，在小区域中，磁化和地磁场的变化非常小，可以认为是恒定的。对于每个小 区域，计算其滤波参数并对单独的点使用一阶泰勒级数来重构。
语法
gmt grdredpol
ingrid
-Grtp_grd
[ -Cdec/dip ]
[ -Eiinc_grd ]
[ -Eddec_grd ]
[ -Fm/n ]
[ -Mm|r ]
[ -N ]
[ -Wwin_width ]
[ -V[level] ]
[ -Tyear ]
[ -Zfiltergrd ]
[ -V[level] ]
[ -nflags ]
[ --PAR=value ]
输入数据
ingrid[=ID|?varname][+bband][+ddivisor][+ninvalid][+ooffset][+sscale]
输入网格名。通过追加 =ID 可指定 网格格式 [默认为 =nf]。 追加 ?varname 可指定 NetCDF 变量 [默认为 GMT 找到的第一个 2-D 网格]。 参数详细介绍请参考 读 netCDF 文件。
必须选项
- -G
- -Grtp_grd
RTP 结果网格文件名
可选选项
- -C
- -Cdec/dip
使用恒定的偏角和倾角
- -E
-Eiinc_grd
- -Eddec_grd
从 inc_grd 和 dec_grd 中获取磁化倾角和偏角，默认使用 IGRF 计算。 注 ：这两个网格不需要与输入的磁异常一致，范围和分辨率都可以不同。
- -F
- -Fm/n
滤波窗口大小，m 和 n 分别为行和列，默认为 25/25
- -M
- -Mm|r
设置边界条件，m|r 分别表示镜像对称或者以边缘对称
- -N
- -N
不使用泰勒展开
- -R
- -Rwest/east/south/north
设置输入点的区域，默认与输入相同
- -T
- -Tyear
IGRF 程序用于计算每个点的偏角和请教的年份，默认为 2000 年
- -W
- -Wwidth
滑动窗口的大小，单位为度，默认为 5 度
- -Z
- -Zfilter_grd
输出滤波文件
- -V
- -V[level]
设置 verbose 等级 [w]。 (参数详细介绍)
- -n
- -n[b|c|l|n][+a][+bBC][+c][+tthreshold]
设置网格文件的插值方式。 (参数详细介绍)
网格重采样的后果
在对网格进行重采样或重新取样时，会使用多种算法（参见 -n）， 这些算法可能导致重采样后的数值出现一定的失真或意外结果。 使用样条函数进行重采样的一个常见现象是：新的采样值可能会略微超出原始网格的全局最小值或最大值。 如果这种超出范围的现象不可接受，可以在 -n 选项中添加 +c 参数， 以对重采样结果进行裁剪，从而确保输出的数值不会超过输入网格的最小值和最大值。
示例
假定磁异常文件 anom.grd 化极到 2010 历元，dec.grd 和 dip.grd
分别为磁偏角和磁倾角文件，使用 2 度划分整个区域，并用 45 个系数的滤波器计算 RTP
gmt grdredpol anom.grd -Grtp.grd -W2 -F45/45 -T2010 -Edec.grd/dip.grd -V
计算相同的 RTP，但磁场和磁化向量是平行的，并使用 IGRF 计算
gmt grdredpol anom.grd -Grtp.grd -W2 -F45/45 -T2010 -V
参考文献
Luis, J.L. and Miranda, J.M., 2008, Reevaluation of magnetic chrons in the North Atlantic between 35N and 47N: Implications for the formation of the Azores Triple Junction and associated plateau. JGR, VOL. 113, B10105, doi:10.1029/2007JB005573