使用 getNOAA.bathy 时出错，坐标值有限制吗？

【问题标题】：Error using getNOAA.bathy, are there restrictions on coordinates values?使用 getNOAA.bathy 时出错，坐标值有限制吗？
【发布时间】：2016-01-22 00:08:04
【问题描述】：

我正在尝试从包marmap 执行getNOAA.bathy。

我可以成功执行以下操作（来自here）：

library(marmap)
getNOAA.bathy(lon1=-20,lon2=-90,lat1=50,lat2=20, resolution=10) -> a
plot(a, image=TRUE, deep=-6000, shallow=0, step=1000)

但是，当我执行以下操作时：

getNOAA.bathy(lon1=-80,lon2=-79.833333,lat1=32.7,lat2=32.833333, resolution=10) -> a
plot(a, image=TRUE, deep=-6000, shallow=0, step=1000)

我收到错误：

getNOAA.bathy(lon1 = -80, lon2 = -79.833333, lat1 = 32.7, lat2 = 32.833333, 错误：无法访问 NOAA 服务器

问题：

对 LAT/LON 值有特殊限制吗？我是不是这里算错了什么？
是否有“更好”的软件包可以支持我的 LAT/LON 值？

【问题讨论】：

你说的是哪个包？
@Pascal：我相信是marmap
我猜你的区域对于数据集分辨率来说太小了。
@Pascal：当分辨率为
您无法更改数据集的分辨率。扩展您的区域。

标签： r

【解决方案1】：

如getNOAA.bathy() 的帮助文件中所述，resolution 参数以分钟表示。所以resolution=10 意味着您的网格单元格的维度为经度 10 分钟乘纬度 10 分钟。数字越大，分辨率越差。因此，考虑到您所在的地区，您需要为 ETOPO1 数据集（即getNOAA.bathy() 获取的数据库）使用尽可能高的分辨率：

getNOAA.bathy(lon1=-80, lon2= 79.833333, lat1=32.7, lat2=32.833333, res=1)

这绝对不是高分辨率（您会得到一个包含 80 个值的网格：10 个纬度值乘 8 个经度值），但这是使用 getNOAA.bathy() 可以获得的最大值。

如果您需要更高的分辨率，您应该尝试其他来源，例如marmap 的GEBCO database 和readGEBCO() 功能。您还应该查看marmap-ImportExport 小插图的第 2 和第 3 部分，其中列出了其他来源。

vignette("marmap-ImportExport")

【讨论】：