“高分辨率marmap” - 地图范围< 1分钟

时间:2017-10-23 14:35:08

标签: r gis noaa

是否可以将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)

Resolution is insufficient to show the marine study area...

由于我无法通过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)

A little bit better but not enough to show the sample properly...

这稍微好一些。但是,这还不够。这是与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正确地说明我的学习领域,适合出版吗?

谢谢。

0 个答案:

没有答案