1

我正在尝试使用“levelplot”制作点光栅图。我希望将我的 x 和 y 刻度标签打印为地理坐标(例如 101°W、100°W、...34°N、35°N 等)以下是我迄今为止尝试的测试数据和代码:

library(rasterVis)
library(raster)
library(colorRamps)
x <- c(-98.02,-96.66,-98.29,-100.53,-95.86,-102.49,-99.64,-99.27,-97.00,-96.33,-96.80,-98.36,-99.72,-96.32,-98.03,-99.80,-98.49,-96.42,-101.60,-95.64,-96.00,-98.48,-99.83,-97.76,-98.11,-95.99,-97.21,-99.01,-94.84,-97.60,-97.46,-96.76,-97.04,-98.96,-99.04,-96.03,-97.09,-96.06,-99.13,-98.52,-97.98,-94.68,-99.41)

y <- c(34.80,34.79,34.91,36.80,35.96,36.69,36.83,35.59,34.87,34.60,35.65,36.74,35.54,33.92,35.54,35.20,36.26,36.84,36.60,35.78,35.84,35.48,34.65,34.52,36.38,34.30,36.06,36.98,36.88,36.15,35.23,36.36,35.95,35.04,36.19,36.41,36.12,34.87,34.43,35.84,34.16,34.98,36.42)

z <- c(0.36,0.78,0.32,0.75,0.91,0.11,0.52,0.08,0.6,0.84,0.53,0.58,0.12,0.45,0.44,0.26,0.3,0.78,0.15,0.83,0.71,0.23,0.13,0.73,0.23,0.81,0.61,0.18,0.68,0.34,0.72,0.84,0.79,0.55,0.38,0.92,0.94,0.8,0.19,0.5,0.96,0.71,0.41)

col=colorRampPalette(matlab.like2(255))
levelplot(z ~ x + y, xlim=c(-103.5,-94), ylim=c(33.2,37.5), panel = panel.levelplot.points, cex = 2, contour = FALSE, margin = FALSE, col.regions = col,xlab = list("Longitude", cex = 2), ylab = list("Latitude", cex=2),scales = list(tck = c(1,0), cex=1.5)) 
4

0 回答 0