我想创建一个这样的地图:
我插入数据,然后使用filled.contour()
绘制它如您所见,地图未正确绘图。我需要找到一种方法将输出转换回纬度和经度,以便地图正确绘制。我该怎么做?
我查看了这个论坛:Plotting interpolated data on map但它没有回答我的问题。
地图代码:
df <-read.table("Petra_phytoplankton+POM_xydata_minusNAs_noduplicates.txt",header=T)
attach(df)
names(df)
fld <- with(df, interp(x = longitude, y = latitude, z = d13C))
filled.contour.ungeoreferenced<-(filled.contour(x = fld$x,
y = fld$y,
z = fld$z,
color.palette =
colorRampPalette(c("blue", "green","yellow","orange","red")),
xlab = "Longitude",
ylab = "Latitude",
key.title = title(main = "d13C", cex.main = 1)))
map(add=TRUE, fill=TRUE)
这里的FYI是我的数据集的片段:
latitude longitude d13C
-65 -70 -27.7
-61 150 -32.2
-61 150 -28.3
-60 116 -26.8
-60 116 -24.7
-47 38 -24.8
-38 150 -20.5
19 -65.7 -19.9
19 -65.5 -18.5
18 -60.7 -20
18 -58.5 -18.2
18 -57.8 -19
17 -55.4 -18.6
17 -50.8 -18
17 -47.1 -18.3
17 -45.5 -19.4
16 -43.3 -17.9
15 -40.7 -18.5
14 -39.3 -19.9
12 -36.7 -19.9
12 -36.2 -19.9
11 -34.4 -19.2
10 -32 -18.5
9 -30.3 -19.3
8 -29.2 -19.4
7 -26.6 -18.2
7 -25.5 -19.3
6 23.9 -20
3 -21.3 -20.4