我在一组点上有压力值(给定纬度和经度值)。我正在使用以下代码在R中绘制它们:
data("wrld_simpl")
lon = dat$X - 180
lat = dat$Y
phi = dat$phi[,,1]
#quartz(width=11, height=4)
plot.new()
image(lon, rev(lat), phi, xlab="Longitude", ylab="Latitude")
title(main="Phi")
plot(wrld_simpl, add=TRUE)
结果如下:
我希望图像在完整的地图上呈现,而不是截断的版本。请让我知道如何做到这一点。
由于