grdtrend
- 官方文档:
- 简介:
拟合网格的趋势面并计算残差
该命令会读取一个2D网格文件，并用最小二乘方法拟合一个低阶多项式趋势面。多项式趋势面的定义为：
\[m_1 + m_2x + m_3y + m_4xy + m_5x^2 + m_6y^2 + m_7x^3 + m_8x^2y + m_9xy^2 + m_{10}y^3.\]
必选选项
<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