img2grdgmt - 云端在线

程序：

您的姓名

img2grd - 以墨卡托或地理格式提取 img 文件的子集

概要

图像2grd 图像文件 文件 地区 类型 [ ] [ [最小拉/最大拉特] ] [ ] [ 分钟 ] [ ] [ 导航
] [ [规模] ] [ [水平] ] [ 万能龙 [ -n]

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

商品描述

图像2grd 读取 img 格式文件，提取子集，并将其写入网格文件。 这 -M
选项指示 img 文件的球面墨卡托投影是否为
保留或是否应通过撤消墨卡托投影来编写地理网格。 如果
选择了地理网格，您还可以请求重新采样到精确的 -R 给定的。

所需 争论

图像文件
墨卡托 img 格式文件，例如海洋重力或海底地形字段
根据 Sandwell 和 Smith 的卫星高度计数据估计。 如果用户设置了
环境变量 $GMT_DATADIR， 然后 图像2grd 会尝试找到 图像文件 in
$GMT_DATADIR; 否则它会尝试打开 图像文件 直。

-G文件
文件 是输出网格文件的名称。

-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 选项。
在后一种情况下，绘制了平面的透视图，没有第三个
尺寸。

可选 争论

-C 设置相对于投影中心的 x 和 y 墨卡托坐标 [默认为
相对于网格的左下角]。 需要 -M.

-D[最小拉/最大拉特]
使用扩展纬度范围 -80.738/+80.738。 或者，附加
最小拉/最大拉特 作为输入 img 文件的纬度范围。 [默认是
-72.006/72.006]。 通常不需要，因为我们可以确定范围从
检查文件大小。

-E 可以在什么时候使用 -M 未设置为强制最终网格具有完全相同的
区域按要求与 -R. 默认情况下，最终区域是直接投影
原始墨卡托区域，通常会稍微超出
请求的纬度范围，而且纬度的网格增量不
匹配经度增量。 然而，额外的重采样引入了小
插值错误，只有在输出网格必须匹配
请求的区域并具有 x_inc = y_inc。 在这种情况下，区域设置为 -R 必须
以增量的倍数给出（例如， -R0/45/45/72）。

-I 表明 分钟 作为输入 img 像素的宽度（以经度分钟为单位）。
[默认为 2.0]。 通常不需要，因为我们可以从
尺寸的检查。

-M 输出球面墨卡托网格 [默认是地理经纬度网格]。 这
img 文件的球面墨卡托投影被保留，因此该区域 -R
用户设置略有修改； 修改区域对应于边缘
像素数[或一组 导航 像素]。 设置网格文件头，以便 x 和
y 轴长度表示距图像西边和南边的距离，
以用户默认单位衡量，其中 -Jm1 和调整 -R. 通过设置
默认 项目_ 椭圆体 = 球体，用户可以使用调整后的进行叠加 -R
使它们匹配。 看 示例 以下。 调整后的 -R 也写在格子里
header 注释，以便稍后找到。 看 -C 设置相对于的坐标
投影中心。

-N导航 将输入 img 像素中的值平均为 导航 by 导航 正方形，并创建
每个这样的正方形有一个输出像素。 如果与 -T3 它会报告一个平均值
0 到 1 之间的约束。如果与 -T2 输出将是平均数据值
或 NaN 根据平均约束是否 > 0.5。 导航 必须平均分配
到 imgfile 的尺寸（以像素为单位）。 [默认 1 不求平均值]。

-S[规模]
将 img 文件值乘以 规模 在存储在网格文件中之前。 [默认是
1.0]。 对于最近的 img 文件：img topo 文件存储在（更正的）米 [-S1]；
mGal*10 中的自由空气重力文件 [-S0.1 得到 mGal]; 垂直偏转文件
微弧度*10 [-S0.1得到微弧度]，垂直重力梯度文件
欧特沃斯*50 [-S0.02 获得 Eotvos，或 -S0.002 得到 mGal/km]）。 如果不 规模 给出
我们尝试通过检查文件名寻找线索来确定比例。

-T类型 类型 处理约束信息的编码。 类型 = 0 表示没有
此类信息编码在 img 文件中（用于 1995 年之前版本的
重力数据）并获取所有数据。 类型 > 0 表示约束信息为
编码（1995 年和更新（当前）版本的 img 文件），以便人们可以
生成网格文件如下： -T1 获取所有点的数据值， -T2 获取数据
约束点的值和插值点的 NaN； -T3 获得 1 在
约束点和 0 插值点 [默认为 1]。

-V[水平] （更多的 ...）
选择详细级别 [c]。 特别推荐这里，因为它有助于看到
如何调整坐标。

-W万能龙
表明 万能龙 作为输入 img 文件的最大经度范围。 版本
自 1995 年以来 万能龙 = 360.0，而一些早期的文件有 万能龙 = 390.0。
[默认为 360.0]。

-^ or 只是 -
打印一条关于命令语法的短消息，然后退出（注意：在 Windows 上
只用 -).

-+ or 只是 +
打印广泛的使用（帮助）消息，包括对任何
模块特定选项（但不是 GMT 通用选项），然后退出。

-? or 没有 参数
打印完整的使用（帮助）消息，包括选项的解释，然后
退出。

- 版
打印 GMT 版本并退出。

--显示数据目录
打印 GMT 共享目录的完整路径并退出。

地理 示例

-M 如果您需要输出网格位于地理区域，则应排除选项
坐标。 提取区域内的数据 -R-40/40/-70/-30 从 world_grav.img.7.2 和
重新投影以产生地理坐标，您可以尝试

img2grd world_grav.img.16.1 -Gmerc_grav.nc -R-40/40/-70/-30 -V

因为 img 文件中的纬度间距是以墨卡托为单位的等距，所以
结果网格将与指定的不匹配 -R 准确地说，纬度间距不会
等于经度间距。 如果您需要与您的完全匹配 -R 和相同的间距
在经度和纬度中，使用 -E 选项​​：

img2grd world_grav.img.16.1 -Gmerc_grav.nc -R-40/40/-70/-30 -E -V

机甲 示例

由于 img 文件位于墨卡托投影中，因此您不应提取地理网格
如果您的计划是制作墨卡托地图。 如果你这样做了，你就结束了投影和
重新投影网格，丢失短波长细节。 更好用 -M 并绘制网格
使用与所需墨卡托投影具有相同比例的线性投影（参见 GMT
例 29)。 提取区域内的数据 -R-40/40/-70/-30 从 world_grav.img.7.2， 跑

gmt img2grd -M world_grav.img.7.2 -Gmerc_grav.nc -R-40/40/-70/-30 -V

请注意 -V 选项告诉我们范围已调整为
-R-40/40/-70.0004681551/-29.9945810754。 我们也可以使用 数据信息 找到网格文件
标题显示其区域为 -R0/80/0/67.9666667 这是我们将得到的 x,y 的范围
从球面墨卡托投影使用 -R-40/40/-70.0004681551/-29.9945810754 和
-Jm1. 因此，要获取 ship.lonlatgrav 并使用它来对 merc_grav.nc 进行采样，我们可以这样做：

gmt 设置 PROJ_ELLIPSOID 球体

gmt 地图项目-R-40/40/-70.0004681551/-29.9945810754 -Jm1i ship.lonlatgrav | \
gmt grdtrack -Gmerc_grav.nc | gmt 地图项目\
-R-40/40/-70.0004681551/-29.9945810754 -Jm1i -I > ship.lonlatgravsat

建议使用上述投影和反投影数据的方法。
一个应用程序，因为只有一个插值步骤（在 轨道）。 如果一个
先尝试把grid文件转换成lon,lat然后采样，有两个
插值步骤（在转换和采样中）。

要从上面的网格制作一个 lon,lat 网格，我们可以使用

gmt grdproject Merc_grav.nc -R-40/40/-70.0004681551/-29.9945810754 -Jm1i -I -D2m -Ggrav.nc

在某些情况下，这并不容易，因为 -R 在两个坐标系中可能不会对齐
好。 出现这种情况时，我们也可以使用（其实用起来可能总是更好）

gmt grd2xyz merc_grav.nc | gmt 地图项目 \
-R-40/40/-70.0004681551/-29.994581075 -Jm1i -I | \
gmt 表面 -R-40/40/-70/70 -I2m -Ggrav.nc

要制作上述区域的墨卡托地图，假设我们的 gmt.conf 值为
PROJ_LENGTH_UNIT 是英寸。 那么由于上面的 merc_grav.nc 文件是用 -Jm1i
它有 80 英寸宽。 我们可以使用以下方法制作 8 英寸宽的地图 -Jx0.1i 在任何地图程序上
应用于这个网格（例如， 等高线, 格图, 视图)，然后对于覆盖
在 lon,lat 工作（例如， 精神力, 海岸线) 我们可以使用上面的调整 -R 和 -Jm0.1 得到
两个系统相匹配。

然而，我们可以比这更聪明。 意识到输入的 img 文件有像素 2.0
分钟宽（或使用 grdinfo merc_grav.nc 检查 nx 和 ny）我们意识到
merc_grav.nc 使用了 img 文件的全分辨率，它有 2400 x 2039 像素，并且
在 8 英寸宽时，这是每英寸 300 像素。 我们决定不需要那么多，我们
将满足每英寸 100 像素，因此我们希望将数据平均为 3 x 3
正方形。 （如果我们想要一个等高线图，我们可能会选择更多地平均数据
（例如，6 x 6）以获得平滑的轮廓。）由于 2039 不能被 3 整除，我们将得到一个
不同的调整 -R 这次：

gmt img2grd -M world_grav.img.7.2 -Gmerc_grav_2.nc -R-40/40/-70/-30 -N3 -V

这次我们发现调整后的区域是 -R-40/40/-70.023256525/-29.9368261101 和
输出为 800 x 601 像素，对我们来说是更好的尺寸。 现在我们可以创建一个人工
为此使用照明文件 梯度:

gmt grdgradient Merc_grav_2.nc -Gillum.nc -A0/270 -Ne0.6

如果我们还有一个名为“grav.cpt”的 CPT 文件，我们可以创建一个彩色阴影浮雕图
喜欢这个：

gmt grdimage merc_grav_2.nc -Iillum.nc -Cgrav.cpt -Jx0.1i -K > map.ps
gmt psbasemap -R-40/40/-70.023256525/-29.9368261101 -Jm0.1i -Ba10 -O >> map.ps

假设您只想从 img 文件中获取受约束的数据值，以纬度/经度表示
坐标。 然后运行 图像2grd 与 -T2 选项，使用 坐标系 转储值，管道
通过 grep -v NaN 消除 NaN，并通过管道 地图项目 与逆
投影如上。

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