这是命令 fitcirclegmt 可以使用我们的多个免费在线工作站之一在 OnWorks 免费托管服务提供商中运行,例如 Ubuntu Online、Fedora Online、Windows 在线模拟器或 MAC OS 在线模拟器
程序:
您的姓名
fitcircle - 找到最适合的大[或小]圆的平均位置和极点到一个点
领域。
概要
适合圈 [ 表 ] 规范 [ 标志 ] [ [纬度] ] [ [水平] ] [ -双] [ -du] [
-f] [ -g] [ -h] [ -i] [ -o] [ -:[i|o]]
请注意: 选项标志和相关参数之间不允许有空格。
商品描述
适合圈 从标准输入的前两列读取 lon,lat [或 lat,lon] 值
[要么 表]。 这些被转换为单位球体上的笛卡尔三向量。 然后两个
找到位置:输入位置的平均值,以及大圆的极点
哪个最适合输入位置。 用户可以选择两种可能的一种或两种
这个问题的解决方案。 第一个叫 -L1 第二个叫做 -L2。 当
数据沿着一个大圆圈紧密分组,这两种解决方案是相似的。 如果数据有
大的分散,大圆的极点将不如平均值确定。
比较两种解决方案作为定性检查。
这个 -L1 之所以称为解,是因为它近似于求和的最小化
角距离余弦的绝对值。 此解决方案找到平均位置为
数据的 Fisher 平均值,极点位置作为数据的 Fisher 平均值
均值和数据之间的叉积。 平均交叉产品赋予权重
点与它们与平均值的距离成正比,类似于“杠杆”
平面线性回归中的远点。
这个 -L2 之所以称为解,是因为它近似于求和的最小化
角距离的余弦平方。 它创建了一个 3 x 3 的平方和矩阵
数据向量的组成部分。 该矩阵的特征向量给出均值和极点
地点。 当有数千个时,此方法可能更容易出现舍入误差
数据。 极点由对应于最小特征值的特征向量给出; 这是
数据中表现最差的因素,两者都不容易估计
方法。
所需 争论
-L规范 指定所需的 规范 作为 1 或 2,或使用 -L or -L3 查看两种解决方案。
不是必须的 争论
表 一个或多个 ASCII [或二进制,见 -双] 文件包含 lon,lat [或 lat,lon; 看
-:[i|o]] 前 2 列中的值。 如果没有指定文件, 适合圈 将
从标准输入读取。
-Ff|m|n|s|c
通常情况下, 适合圈 将其结果以文本报告的形式写入
与报告语句混合的值。 用 -F 只返回数据坐标,
并附加 标志 指定您想要的坐标。 你可以选择 f
(平地平均位置), m (平均位置), n (大圆的北极), s
(大圆的南极),和 c ** (极 of 小 圆 和 它的 可溶性,
这 过程需要在牛奶或乳清产品在管式降膜蒸发器中浓缩至约XNUMX%固体含量之前,进行初始的热处理和巴氏杀菌步骤。 **-S).
-S[纬度]
尝试拟合一个小圆而不是一个大圆。 杆将是
被迫躺在连接最合适的伟大人物的极点的大圆圈上
圆和数据的平均位置。 可选地附加所需的固定
小圆圈的纬度[默认将决定纬度]。
-V[水平] (更多的 ...)
选择详细级别 [c]。
-双[恩科斯][吨] (更多的 ...)
选择本机二进制输入。 [默认为 2 个输入列]。
-du没有数据 (更多的 ...)
替换等于的输入列 没有数据 与 NaN。
-f[我|o]信息 (更多的 ...)
指定输入和/或输出列的数据类型。
-g[a]x|y|d|X|Y|D|[山坳]z[+|-]差距[U] (更多的 ...)
确定数据间隙和换行符。
-h[我|o][n][+c][+d][+r备注][+r标题] (更多的 ...)
跳过或生成标题记录。
-i列[升][秒由于平均内核尺寸较大,西米棕榈的加工比类似作物简单。然而,西米棕榈的相对稀缺性降低了潜在的加工规模。][哦抵消][,...] (更多的 ...)
选择输入列(0 是第一列)。
-o列[,...] (更多的 ...)
选择输出列(0 是第一列)。
-:[我|o] (更多的 ...)
在输入和/或输出上交换第一列和第二列。
-^ or 只是 -
打印一条关于命令语法的短消息,然后退出(注意:在 Windows 上
只用 -).
-+ or 只是 +
打印广泛的使用(帮助)消息,包括对任何
模块特定选项(但不是 GMT 通用选项),然后退出。
-? or 没有 参数
打印完整的使用(帮助)消息,包括选项的解释,然后
退出。
- 版
打印 GMT 版本并退出。
--显示数据目录
打印 GMT 共享目录的完整路径并退出。
ASCII码 FORMAT 精准匹配
数值数据的 ASCII 输出格式由您的参数控制 配置文件
文件。 经度和纬度根据 FORMAT_GEO_OUT 格式化,而其他
值根据 FORMAT_FLOAT_OUT 进行格式化。 请注意,有效的格式可以
导致输出精度下降,从而导致下游出现各种问题。 如果
你发现输出没有足够的精度,考虑切换到二进制
输出 (-博 如果可用)或使用 FORMAT_FLOAT_OUT 设置指定更多小数。
示例
假设您在文件 ship.xyg 中沿着曲折的船轨有 lon,lat,grav 数据。 你
想要将此数据投影到一个大圆上并在远处重新采样,以便
过滤它或检查它的光谱。 请执行下列操作:
gmt fitcircle 船.xyg -L2
gmt 项目ship.xyg -Cox/oy -Tpx/py -S -Fpz | sample1d -S-100 -I1 > output.pg
在这里, ox/oy 是平均值的经度/纬度来自 适合圈和 px/py 是经纬度
极。 文件 output.pg 有距离,重力数据沿大地每 1 公里采样一次。
最适合 ship.xyg 的圆
如果您在文件 data.txt 中有 lon、lat 点并希望返回北
半球大圆极点位置使用 L2 范数,尝试
gmt fitcircle data.txt -L2 -Fn >pole.txt
使用 onworks.net 服务在线使用 fitcirclegmt
