arrays - R - 等高线图

Question

我已经绘制了等高线图，但我需要做一些改进。这是使用的数据的结构:

str(lon_sst)
# num [1:360(1d)] -179.5 -178.5 -177.5 -176.5 -175.5 ...

str(lat_sst)
# num [1:180(1d)] -89.5 -88.5 -87.5 -86.5 -85.5 -84.5 -83.5 -82.5 -81.5 -80.5 ...

dim(cor_Houlgrave_SF_SST_JJA_try)
# [1] 360 180

require(maps)
maps::map(database="world", fill=TRUE, col="light blue")
maps::map.axes()
contour(x=lon_sst, y=lat_sst, z=cor_Houlgrave_SF_SST_JJA_try[c(181:360, 1:180),],
        zlim=c(-1,1), add=TRUE)
par(ask=TRUE)
filled.contour(x = lon_sst, y=lat_sst,
               z=cor_Houlgrave_SF_SST_JJA_try[c(181:360, 1:180),],
               zlim=c(-1,1), color.palette=heat.colors)

因为大多数相关性都接近于 0，所以很难看到大的。

我可以让它更容易看，或者我可以改变分辨率以便我可以放大吗？目前轮廓间距太小，所以我看不到轮廓水平是多少。

我在哪里可以看到增量，我将范围设置为(-1,1)，我不知道如何手动设置间隔。

有人能告诉我如何绘制 map 的特定区域，例如从 100 到 160 的经度和从 -50 到 -80 的纬度吗？我试图更换 lon_sst和 lat_sst ，但它有一个尺寸错误。谢谢。

Answer 1

要回答似乎是相同请求的 1 和 3，请尝试:

maps::map(database="world", fill=TRUE, col="light blue", 
                            ylim=c(-80, -50), xlim=c(100,160) )

要解决 2:您的范围比 [-1,1] 小得多。这些等高线上的标签是像 .06、-.02 和 .02 这样的数字。 contour函数将接受“nlevels”或“levels”参数。一旦你有一个放大的部分，你可以用它来调整轮廓的 z 分辨率。