Question

我正在尝试从包marmap执行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)。

我可以成功执行以下内容（来自here）：

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)

但是，当我执行以下内容时：

data<-data.frame(code= c(345, 346), name= "Yemen", v1= c("", "text1"), v2= c("text2", ""),v3= c("text3", ""),v4= c("", "text4"))
code  name    v1    v2    v3    v4
345   Yemen         text2 text3      
346   Yemen   text1             text4

aggregate(x=data[c("v1","v2","v3","v4")], by=list(name=data$name), paste)
name v1.1  v1.2  v2.1 v2.2  v3.1 v3.2 v4.1  v4.2
1 Yemen      text1 text2      text3           text4

我收到了错误：

getNOAA.bathy出错（lon1 = -80，lon2 = -79.833333，lat1 = 32.7，lat2 = 32.833333，：无法联系到NOAA服务器

问题：

对LAT / LON值有特殊限制吗？我错误估算了这里的东西？
是否有更好的＆＃34;可以支持我的LAT / LON值的包吗？

Answer 1

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

getNOAA.bathy()

这绝对不是高分辨率（你得到80个值的网格：10个纬度值乘以8个经度值），但这是getNOAA.bathy(lon1=-80, lon2= 79.833333, lat1=32.7, lat2=32.833333, res=1)得到的最大值。

如果您需要更高的分辨率，则应尝试其他来源，例如getNOAA.bathy()的{{1}}和readGEBCO()功能。您还应该看一下列出其他来源的marmap插图的第2和第3部分。

marmap-ImportExport

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

问题：

1 个答案: