1

是否可以将marmap用于高度缩放的地图？这是我挣扎的地方——来自 NOAA 网格的分辨率似乎在高分辨率方面不够……

# Draw a map with getNOAA.bathy

library(marmap)

# sample location
s <- data.frame(y = -44.968276, x = 167.415300)

# map with NOAA 1m data
NOAA.bathy <- getNOAA.bathy(167, 168, -44.5, -45.5, res = 1, keep = TRUE)
plot(NOAA.bathy, image = T, land = F)
points(s$x,s$y, pch = 21, col = "black", bg = "red", cex = 1.3)

分辨率不足以显示海洋研究区域...

由于我无法使用 getNOAA.bathy() 获得足够的分辨率，因此我尝试了 GEPCO 数据库（30 分钟）并下载了相应的网格：

# Same location with GEBCO data =
GEPCO.bathy <- readGEBCO.bathy("~/GEBCO_2014_2D_167.0_-45.5_168.0_-44.5.nc")
plot(GEPCO.bathy, image = T, land = F)
points(s$x,s$y, pch = 21, col = "black", bg = "red", cex = 1.3)

稍微好一点，但不足以正确显示样本......

这只是好一点。然而，这还不够。这是与ggmap的比较。

#comparison to gmaps 
library(ggmap)
gmap = get_map(location = c(lon = s$x, lat = s$y), zoom = 11)
ggmap(gmap) + geom_point(data = s, aes(x = s$x, s$y), color = "red", size = 4)

你们能帮我找到一种方法来用marmap正确说明我的研究领域，适合出版吗？

谢谢你。

4

0 回答 0