grdtrend¶
官方文档: | grdtrend |
---|---|
简介: | 拟合网格的趋势面并计算残差 |
该命令会读取一个2D网格文件,并用最小二乘方法拟合一个低阶多项式趋势面。多项式趋势面的定义为:
m1 + m2*x + m3*y + m4*x*y + m5*x*x + m6*y*y + m7*x*x*x + m8*x*x*y + m9*x*y*y + m10*y*y*y.
必选选项¶
<gridfile>
- 2D网格文件名
-N<n_model>[+r]
指定要拟合的模型。
<n_model>
指定要拟合的模型的参数个数。例如-N3
表示bilinear趋势,-N6
表示 quadratic趋势面。加上+r
表示robust拟合,此时,程序会根据 robust scale estimate多次迭代,给数据重新赋予权重,以得到一个对outliers 不敏感的解。
可选选项¶
-D<diff.nc>
- 将残差(输入减去拟合)结果写到网格文件中
-T<trend.nc>
- 将拟合得到的趋势文件写到网格文件
<trend.nc>
中 -W<weight.nc>
- 若
<weight.nc>
存在,则读取该文件,并求解一个有权重的最小二乘问题。默认为常规的最小二乘拟合。若-N
选项中指定了robust拟合,则robust拟合中所使用的权重会写到文件<weight.nc>
中。
示例¶
从网格文件中移除线性趋势,并将结果写到残差文件中:
gmt grdtrend hawaii_topo.nc -N3 -Dhawaii_residual.nc
对网格文件做bicubic面的robust拟合:
gmt grdtrend hawaii_topo.nc -N10r -Thawaii_trend.nc -Whawaii_weight.nc -V