嘗試將NetCDF文件讀取爲柵格，但NA值發生變化

我試圖讀取NetCDF文件作爲R中的柵格。我正在使用的數據是來自NOAA海洋顏色實驗室的SST數據。我已經在程序「SeaDAS」中打開了這個文件，所以我知道這個問題不在原始文件中。嘗試將NetCDF文件讀取爲柵格，但NA值發生變化

我的代碼讀取該文件看起來是這樣的：

library(raster) 
sst.nc<-raster("may2015_sst_monthly.nc") 
proj4string(sst.nc)<-CRS("+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0") 
plot(sst.nc)

出於某種原因，柵格已經賦值到陸地（見圖片），這應該是「不適用」。任何人都知道我需要做什麼來解決這個問題？

來源

2016-12-02 Splash1199

您是否已經檢查過這些土地面積的價值，而不是僅僅在一個情節中觀察它們？ – rosscova

@rosscova是的，我做過 – Splash1199

它有什麼價值？你指定它不是「不適用」，但知道這個價值是什麼，它說明了這個故事的很大一部分。如果NetCDF分配一個數字值，似乎肯定它是一個相當合理的（因爲它應該可以與所有其他數據可靠區分），這使@Jack Ballard的答案看起來很合理。 – rosscova

它可能不是NA，而是標記爲一些非常大的數字，比如67676（或者可能是一個非常大的負數）。然後，你可以只是做：

sst.nc[sst.nc==67676]=NA

如果你不能找出不管是什麼原因的確切數量（有時它也有小數點走出去到16位數字），你可以隨便挑了一大批顯然超出範圍您的數據：

sst.nc[sst.nc>60000]=NA

來源

2016-12-02 21:14:07

你是對的！真的很奇怪，當我在SeaDAS中打開文件時，值是「NA」。但在R中，所有基於陸地的小區都被分配了值'45'。 – Splash1199

嘗試將NetCDF文件讀取爲柵格，但NA值發生變化

回答

相關問題