如何在netcdf地图上绘制纬度/经度点

时间:2014-04-03 09:40:35

标签: r plot coordinates netcdf

我已经在观察到的气候数据(全球表面温度)的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)

然后我创建了一个坐标为

的SpatialPointsDataFrame
speciesall_df <- SpatialPointsDataFrame(coordinates, speciesall, proj4string=CRS("+proj=longlat +ellps=WGS84"))

现在我不确定如何将点数据与netcdf地图相结合?

0 个答案:

没有答案