我已经在观察到的气候数据(全球表面温度)的netcdf文件中读到了R,并希望在此顶部绘制一些物种位置数据(纬度/经度坐标点),但我很难这样做。 / p>
我使用以下方法生成了我的温度图:
obsula.pdclann <- open.ncdf("obsula.pdclann.nc")
temp_var <- "temp_mm_1_5m"
temp <- get.var.ncdf(obsula.pdclann, varid=temp_var)
image.plot(temp)
我的物种出现数据通过以下方式读入:
speciesall <- read.table("C:/Users/Waterson/Documents/UOB/Data/speciesData/P.subrufa.txt", header=T, sep="\t")
coordinates <- cbind(speciesall$longitude, speciesall$latitude)
然后我创建了一个坐标为
的SpatialPointsDataFramespeciesall_df <- SpatialPointsDataFrame(coordinates, speciesall, proj4string=CRS("+proj=longlat +ellps=WGS84"))
现在我不确定如何将点数据与netcdf地图相结合?