问题描述
使用R中的功能定位器，我可以按netcdf4文件映射中的一个点，它会给我它们的坐标。但是，我也希望显示该点是否在陆地或水上。该信息在变量mask_rho上，但是在按下坐标的情况下，我不知道如何提取它。
library(ncdf4)
library(fields)
nc_file <- 'C:/Users/dell/Desktop/Proyecto/ich/input/newperush_avg.Y2011.M12.newperush.nc'
lon <- 280.81572993
lat <- -8.148526246
nc <- nc_open(nc_file)
lon_rho <- ncvar_get(nc,'lon_rho')
lat_rho <- ncvar_get(nc,'lat_rho')
mask_rho <- ncvar_get(nc,'mask_rho')
xlim <- c(lon-1.5,lon+1.5)
ylim <- c(lat-1.5,lat+1.5)
X11()
image.plot(lon_rho,lat_rho,mask_rho,xlim = xlim,ylim = ylim,xlab =
'LON',ylab = 'LAT')
points(lon,lat)
locator(n = 1,type = 'p')
解决方法
暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！
如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@）