gmtspatialgmt - 云端在线

程序：

您的姓名

gmtspatial - 对线和多边形进行地理空间操作

概要

空间 [ 表 ] [ [a最小距离][单元]] [ ] [
[+f文件][+a最大][+dDMAX][+c|C最大值][+s事实] ] [ +|- ] [ [l] ] [ -I[e|i] ] [
文件[+a][+p开始][+r][+z] ] [ [[-|+]*单元*][+h][+l][+p] ] [ 地区 [ i|u|s|j [
[多边形] ] [ [水平] ] [ -b] [ -d] [ -f] [ -g] [
-h] [ -i] [ -o] [ -:[i|o]]

请注意： 选项标志和相关参数之间不允许有空格。

商品描述

空间 读取一个或多个数据文件（可能是多段文件），其中包含
封闭的多边形并以指定的方式操作这些多边形。 操作包括
面积计算、手性反转和多边形相交。

所需 争论

没有。

可选 争论

表 一个或多个 ASCII（或二进制，见 -双[恩科斯][类型]) 数据表文件保存一个
数据列数。 如果没有给出表，那么我们从标准输入中读取。

-A[一个最小距离][单元]
执行空间最近邻 (NN) 分析：确定
每个点并报告每对中涉及的 NN 距离和点 ID（ID
是从 0 开始的输入记录号）。 用 -AA 抽取数据集
没有神经网络距离低于阈值 最小距离. 在这种情况下，我们写
输出（可能是平均的）坐标和更新的 NN 距离和点 ID。
负点数表示原点被加权平均值代替
（绝对 ID 值给出要包含的第一个原始点 ID 的 ID
在平均。）

-C 将多边形裁剪到地图区域，根据需要包括多边形的地图边界。
结果是一个封闭的多边形（见 -T 代替截断）。 需要 -R.

-D[+f文件][+一个最大][+dDMAX][+c|C最大值][+s事实]
检查输入线或多边形之间的重复项，或者，如果 文件 是通过
+f, 检查输​​入特征是否已经存在于 文件。 我们
考虑精确（相同的数字和坐标）和近似匹配的情况
（两个特征的最近点之间的平均距离小于阈值）。
我们还认为某些功能可能已被颠倒。 特征被考虑
近似匹配，如果它们的最小距离小于 DMAX [0]（见单位）和
它们的接近度（定义为两个人之间的平均距离之比）
特征除以其平均长度）小于 最大值 [0.01]。 对于每个
发现重复，输出记录以单个字母 Y 开头（完全匹配）或
～（近似匹配）。 如果两个匹配段的长度相差超过一个
因子为 2，那么我们将重复项视为子集 (-) 或超集
(+)。 最后，我们还要注意两条线是否是连续分割的结果
跨越日期变更线 (|) 的线。 对于多边形，我们还考虑分数
地区差异； 重复项的差异必须小于 最大 [0.01]。 默认情况下，
我们计算平均线间隔。 用 +C厘米 改为计算中线
分离，因此具有稳健的接近度值。 此外，默认情况下，我们考虑所有
一条线上的点与另一条线上的点之间的距离。 附加 +p 限制比较
垂直投影到另一条线上的点的点（而不是它的
延期）。

-E+|- ]
重置所有多边形的手性以匹配给定的 + （逆时针）或 -
（顺时针）。 暗示 -Q+.

-F[l] 强制输入数据在输出时变成多边形，即，如果不是，则明确关闭它们
已经关闭。 可选地，附加 l 强制线几何。

-我[e|我]
确定所有多边形对之间的交点位置。 附加 i 至
仅计算内部（即自相交多边形）交叉或 e 只
计算外部（即，在多边形巴黎之间）交叉[默认为两者]。

-N文件[+a][+p开始][+r][+z]
确定一个（或全部，与 +a) 输入数据中每个特征的点是
在任何给定的多边形内 文件. 如果在里面，则报告哪个多边形
这是; 多边形 ID 要么取自分配给 Z 的空间值，要么
段头（第一个 -Z， 然后 -L 被扫描），或者被分配运行
初始化为的数字 开始 [0]。 默认情况下，输入段是
发现位于多边形内的多边形 ID 被写入标准输出，其中多边形 ID 编码为
段头作为 -ZID. 或者，附加 +r 只报告哪个多边形
包含一个功能或 +z 将 ID 添加为输出的额外数据列。
不在多边形内的线段不会被写出。 如果不止一个
多边形包含我们跳过第二个（和更远的）场景的相同部分。

-Q[[-|+]*单位*][+h][+l][+p]
测量所有多边形的面积或线段的长度。 用 -Q+h 附加
区域到每个多边形段标题 [默认只是将区域写入标准输出]。
对于多边形，我们还计算质心位置，而对于线数据，我们计算
中点（半长）位置。 附加距离单位以选择单位
使用（见单位）。 请注意，该区域将取决于当前的设置
PROJ_椭圆形； 这应该是最近的椭球以获得准确的结果。 这
质心是使用组成的 3-D 笛卡尔向量的平均值计算的
多边形顶点，而面积是通过等积投影获得的。 为了
您可以预先添加的行长度 -|+ 到单位，计算将使用 Flat
地球或测地线算法，分别[默认为大圆距离]。 经过
默认情况下，我们将开放多边形视为线。 附加 +p 关闭打开的多边形和
因此将所有输入视为多边形，或附加 +l 将所有输入视为行，
即使关闭。

-R[单元]西/东/南/北[/最小/最大][r]
西, 东, 南及 北 指定感兴趣的区域，您可以指定
它们以十进制度数或 [+-]dd:mm[:ss.xxx][W|E|S|N] 格式。 附加 r 如果更低
给出了左右地图坐标而不是 w/e/s/n。 他们俩
速记 -Rg 和 -Rd 代表全球域（经度为 0/360 和 -180/+180
分别为 -90/+90 纬度）。 或者对于网格创建，给
R代码/纬度/nx/ny，其中 码 是 L、C、R 的 2 个字符组合（对于左、
中心或右侧）和 T、M、B 表示顶部、中间或底部。 例如，BL 表示左下角。
这表示矩形区域上的哪个点 LON/纬度 坐标是指
和网格尺寸 nx 和 ny 通过网格间距 -I 用于创建
对应的区域。 或者，指定现有网格文件的名称
和 -R 设置（和网格间距，如果适用）是从网格复制的。
运用 -R单元 预计投影（笛卡尔）坐标与所选的兼容 -J
并且我们反向投影以确定实际的矩形地理区域。 为了
透视图（-p), 可选附加 /最小/最大. 在透视图的情况下
(-p), z 范围 (最小, 最大) 可以附加以指示第三维。 这个
只需要在使用 -Jz 选项，而不是仅使用 -p 选项。
在后一种情况下，绘制了平面的透视图，没有第三个
尺寸。 将多边形裁剪到地图区域，包括多边形的地图边界
如所须。 结果是一个封闭的多边形。

-Si|j|s|u
多边形的空间处理。 从中选择 如果 它返回的交集
多边形（封闭）， -苏 返回多边形的并集（闭合）， -SS 会的
分割跨越日期变更线的多边形，以及 -Sj 这将连接多边形
被日期变更线分割。 注意：仅 -SS 已经实施。

-T[多边形]
针对给定的指定多边形截断多边形，可能导致打开
多边形。 如果没有给出参数 -T 我们从 -R 这
那么是必需的。 请注意，当 -R 剪裁生效，我们还将寻找
长度为 4 或 5 的多边形完全匹配 -R 剪裁多边形。

-V[水平] （更多的 ...）
选择详细级别 [c]。

-双[恩科斯][吨] （更多的 ...）
选择本机二进制输入。 [默认为 2 个输入列]。

-博[恩科斯][类型] （更多的 ...）
选择本机二进制输出。 [默认与输入相同]。

-d[我|o]没有数据 （更多的 ...）
替换等于的输入列 没有数据 与 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 共享目录的完整路径并退出。

单位

对于地图距离单位，附加 单元 d 对于弧度， m 对于弧分，和 s 对于弧
第二，或 e 对于仪表 [默认]， f 对于脚， k 公里， M 对于法定英里， n 航海用
英里，和 u 为美国测量脚。 默认情况下，我们使用球面计算这些距离
与大圆近似。 前置 - 到一个距离（或单位是没有距离是
给出）执行“平面地球”计算（更快但不太准确）或预先 + 至
执行精确的测地线计算（更慢但更准确）。

ASCII码 FORMAT 精确

数值数据的 ASCII 输出格式由您的参数控制 配置文件
文件。 经度和纬度根据 FORMAT_GEO_OUT 格式化，而其他
值根据 FORMAT_FLOAT_OUT 进行格式化。 请注意，有效的格式可以
导致输出精度下降，从而导致下游出现各种问题。 如果
你发现输出没有足够的精度，考虑切换到二进制
输出 （-博 如果可用）或使用 FORMAT_FLOAT_OUT 设置指定更多小数。

例

要将多段文件lines.txt 中的所有行转换为闭合多边形，请运行

gmt 空间线.txt -F > 多边形.txt

要计算多段文件polygons.txt 中所有地理多边形的面积，请运行

gmt 空间多边形.txt -Q > area.txt

相同的数据，但现在将所有多边形定向为逆时针方向并将它们的区域写入
段标题，运行

gmt 空间多边形.txt -Q+h -E+ > area.txt

要确定多边形 A.txt 和 B.txt 之间的交点，请运行

gmt 空间 A.txt B.txt -Ce > crossovers.txt

使用 onworks.net 服务在线使用 gmtspatialgmt