mapprojectgmt - 云端在线

程序：

您的姓名

mapproject - 进行正向和反向地图转换、基准转换和大地测量

概要

地图项目 [ 表 ] 参数 地区 [ b|B|f|F|o|O[龙0/纬度0] ] [ [dx/dy] ] [ c|i|p [
[基准] ] [ [单元] ] [ [x0/y0/][[+|-]单元][+|-] ] [ ] [ 线.xy[/[+|-]单元][+] ] [
[a|c|g|m] ] [ [d|e ] [ ] [h]止[/至] ] [ [水平] ] [ [w|h] ] [ -b] [
-d] [ -f] [ -g] [ -h] [ -i] [ -o] [
-p] [ -s] [ -:[i|o]]

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

商品描述

地图项目 读取（经度，纬度）位置 表 [或标准输入]和
使用指定的地图投影和比例计算 (x,y) 坐标。 可选地，它
可以读取（x，y）位置并计算（经度，纬度）值进行逆运算
转型。 这可用于转换通过数字化获得的线性 (x,y) 点
已知投影到地理坐标的地图。 也可以计算沿线的距离
跟踪，到一个固定点，或最接近一条线。 最后，可以用来执行
各种数据转换。 前 2 列后允许附加数据字段
必须有（经度，纬度）或（x，y）。 查看选项 -: 关于如何阅读
（纬度，经度）文件。

所需 争论

-J参数 （更多的 ...）
选择地图投影。

-R[单元]最小/XMAX/YMIN/最大值[R] （更多的 ...）
指定感兴趣的区域。 UTM 投影的特殊情况：如果 -C 用来
和 -R 未给出，则该区域设置为与给定的 UTM 区域重合，因此
以保留完整的椭球解（有关更多信息，请参阅 RESTRICTIONS
信息）。

可选 争论

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

-Ab|B|f|F|o|O[龙0/纬度0]
-AF 从固定点计算（向前）方位角 经纬度 到每个数据点。
使用 -抗体 获得从数据点到固定点的反向方位角。 用 -敖 要得到
方向 (-90/90) 而不是方位角 (0/360)。 大写 F, B or O 将
从大地纬度转换为地心纬度并估计测地线的方位角
（假设当前椭球不是球体）。 如果没有给出固定点，那么
我们从前一点计算方位角（或后方位角）。

-C[dx/dy]
将投影坐标的中心设置为地图投影中心 [默认为
左下角]。 或者，在要添加的投影单位中添加偏移量（或
减去时 -I 设置）到（来自）投影坐标，例如false
特定投影区域的东距和北距 [0/0]。 用于的单位
offsets 是有效的绘图距离单位（请参阅 PROJ_LENGTH_UNIT），除非 -F is
使用，在这种情况下，偏移量以米为单位。

-DC|i|p
暂时覆盖 PROJ_LENGTH_UNIT 并使用 c （厘米）， i （英寸），或 p （分）
反而。 不能与 -F.

-E[基准]
从大地测量（经度、纬度、高度）转换为地心地球固定 (ECEF)
(x,y,z) 坐标（添加 -I 用于逆转换）。 附加数据 ID（请参阅 -Qd)
或给 椭圆体:dx,dy,dz 哪里 椭圆体 可能是一个椭球 ID（见 -Qe），或
给出为 a[,*inv_f*], 其中 a 是半长轴和 反演_f 是逆
展平（如果省略则为 0）。 如果 基准 是 - 或未给出，我们假设 WGS-84。

-F[单元]
强制 1:1 缩放，即输出（或输入，参见 -I) 数据是实际预测的
米。 要指定其他单位，请附加所需的单位（请参阅 UNITS）。 没有 -F,
输出（或输入，见 -I) 使用 PROJ_LENGTH_UNIT 指定的单位（但
看到 -D).

-G[x0/y0/][[+|-]单元][+|-]
计算沿轨道的距离 or 到可选点集 -Gx0/y0. 附加
距离单位（见 UNITS），包括 c （笛卡尔距离使用输入
坐标）或 C （使用投影坐标的笛卡尔距离）。 这 C 单元
需要 -R 和 -J 要设置。 在没有给出固定点的情况下，我们计算累积
沿轨道的距离。 附加 - 获得连续之间的增量距离
点。 附加 + 通过输入文件中的两个额外列指定第二个点。

-I 做逆变换，即从(x,y) 数据得到(longitude,latitude)。

-L线.xy[/[+|-]单元][+]
确定从输入数据点到给定的线的最短距离
ASCII 多段文件 线.xy. 距离和坐标
最近的点将作为三个新列附加到输出中。 附加
距离单位（见单位），包括 c （使用输入坐标的笛卡尔距离）
or C （使用投影坐标的笛卡尔距离）。 这 C 单位要求 -R 和
-J 要设置。 最后，追加 + 报告线段 ID 和小数
点号而不是最近点的 lon/lat。

-N[a|c|g|m]
从大地纬度转换（使用当前椭球；见 PROJ_ELLIPSOID）
到四个不同的辅助纬度之一（经度不受影响）。 选择
止 a乌塔利奇， c非正式的， g以自我为中心，和 m经纬度[地心]。 用 -I
将辅助纬度转换为大地纬度。

-Q[d|e 列出所有投影参数。 要仅列出基准，请使用 -Qd. 只列出
椭球体，使用 -Qe.

-S 抑制区域外的点。

-T[小时]止[/至]
基准之间的坐标转换 止 和 至 使用标准的莫洛登斯基
转型。 用 -Th 如果第三个输入列的高度高于椭球 [默认
假设高度 = 0，即在椭圆体上]。 使用基准 ID 指定基准（请参阅
-Qd) 或给 椭圆体:dx,dy,dz 哪里 椭圆体 可能是一个椭球 ID（见 -Qe），或
给出为 a[,*inv_f*], 其中 a 是半长轴和 反演_f 是逆
展平（如果省略则为 0）。 如果 基准 是 - 或未给出，我们假设 WGS-84。 -T 也许
结合使用 -R -J 在坐标投影之前更改基准
（加 -I 在逆投影后应用基准转换）。 确保
PROJ_ELLIPSOID 设置适合您的情况。

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

-W[w|h]
在标准输出上打印地图宽度和高度。 不读取输入文件。 只到
输出宽度或高度，追加 w or h， 分别。 的单位
尺寸可以通过改变 -D.

-双[恩科斯][吨] （更多的 ...）
选择本机二进制输入。 [默认为 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 是第一列）。

-p[x|y|z]阿齐姆/提升[/水平][+w龙0/纬度0[/z0]][+vx0/y0] （更多的 ...）
选择透视图。

-s[列][一|r] （更多的 ...）
设置对 NaN 记录的处理。

-:[我|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 设置指定更多小数。

示例

要将 UTM 坐标（以米为单位）转换为地理位置，给定文件 utm.txt 和
知道 UTM 区域（以及区域或半球），尝试

gmt 地图项目 utm.txt -Ju+11/1:1 -C -I -F

将具有 (longitude,latitude) 的文件转换为 (x,y) 在墨卡托上以 cm 为单位的位置
网格为每度 0.5 厘米的给定比例，运行

gmt 地图项目 lonlatfile -R20/50/12/25 -Jm0.5c > xyfile

使用（纬度，经度）转换几个 2 列、二进制、双精度文件
进入横轴墨卡托网格（中心经度 75W）上以英寸为单位的 (x,y) 位置，用于
scale = 1:500000 并抑制那些落在地图区域之外的点，运行

gmt 地图项目轨迹。* -R-80/-70/20/40 -Jt-75/1:500000 -: -S -Di -bo -bi2 > tmfile.b

转换 NAD27 文件 old.dat 中的大地坐标（lon、lat、height）
CONUS 基准（使用 Clarke-131 椭球的基准 ID 1866）到 WGS 84，运行

gmt mapproject old.dat -Th131 > new.dat

计算输入文件 quakes.dat 中每个点之间的最近距离（以公里为单位）
和多段ASCII文件coastline.xy中给出的线段，运行

gmt mapproject quakes.dat -Lcoastline.xy/k > quake_dist.dat

限制

矩形输入区域设置为 -R 通常会被映射到非矩形
网格。 除非 -C 设置后，此网格上最左边的点的 xvalue = 0.0，并且
最低点的 yvalue = 0.0。 因此，在将地图数字化之前，请运行极端
地图坐标通过 地图项目 使用适当的比例并查看 (x,y) 值
它们被映射到。 在设置数字化时使用这些值，以便
逆变换正常工作，或者，使用 AWK 缩放和移动
(x,y) 转换前的值。

对于某些投影，尽管用户选择了一个球面解
椭圆体。 这发生在用户 -R 设置意味着超出域的区域
其中椭球级数展开是有效的。 这些是条件：（1）
兰伯特等角圆锥曲线 (-JL) 和阿尔伯斯等面积 (-JB) 将使用球形解
当地图比例超过 1.0E7 时。 (2) 横轴墨卡托 (-JT) 和 UTM (-菊) 会
当西边界或东边界在 -R 超过
距中央子午线 10 度，与 (3) 卡西尼号 (-JC) 但有一个限制
只有4度。

椭圆体 AND 球体

GMT 将使用椭圆体公式，如果它们被实现并且用户选择了一个
椭球体作为参考形状（参见 PROJ_ELLIPSOID）。 用户需要了解一些
潜在的陷阱：(1) 对于某些预测，例如横轴墨卡托、阿尔伯斯和
Lambert 的共形圆锥我们使用椭球表达式，当映射的区域是
小，并切换到球面表达式（并替换适当的辅助
纬度）用于更大的地图。 椭球公式的使用如下： (a) 横向
墨卡托：当所有点都在中央子午线 10 度以内时，(b) 圆锥投影
当纵向范围小于 90 度时，(c) Cassini 投影，当所有点
位于中央子午线 4 度以内。 (2) 当你试图匹配一些历史
数据（例如，通过某个投影和某个参考获得的坐标
椭圆体）您可能会发现 GMT 给出的结果略有不同。 一个可能
这种不匹配的根源在于较旧的计算通常使用较少的有效数字。 为了
例如，斯奈德的例子经常使用克拉克 1866 椭球（他定义为具有
扁平化 f = 1/294.98）。 从 f 我们得到偏心率的平方为 0.00676862818
（这是 GMT 使用的），而 Snyder 舍入并使用 0.00676866。 这种差异可以
给出几十厘米的差异。 如果您需要再现坐标投影
有了这种稍微不同的离心率，您应该使用
与 Clarke 1866 相同的参数，但 f = 1/294.97861076。 另外，请注意，年龄较大的
数据可能会引用不同的数据，除非您知道使用了哪个数据并且
将所有数据转换为通用数据，您可能会遇到数十到数百个不匹配的情况
米。 (3) 最后，请注意 PROJ_SCALE_FACTOR 有某些默认值
投影，因此您可能必须覆盖设置以匹配生成的结果
与其他设置。

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