在R中投影氣候研究單位（CRU）CDF柵格爲立體極地視圖：方法？

如何將氣候研究組（CDU）CDF投影到R的極地視角？在R中投影氣候研究單位（CRU）CDF柵格爲立體極地視圖：方法？

包cruts其功能cruts2raster（）很好地獲取時間序列數據，我可以使用光柵包和圖（）顯現一個月的數據（或更多）或圖像（）函數。引入所有數據並檢查RasterStack後，繪圖並在查詢時給出座標參照系，但我無法將其轉換爲極座標投影。

我用rgdal（已安裝了最新的PROJ.4到我的Mac的在/ usr/local/lib目錄/ pkgconfig），但似乎無法投射到這個極面觀：

library(cruts) 
library(raster) 
## Bring the downloaded file into R 
r <- cruts2raster("cru_ts3.24.2011.2015.pre.dat.nc", 
        timeRange=c("2011-01-01","2015-01-01")) 
## get the first layer 
r1 <- r[[1]] 

## admire its colors (beautiful!) 
plot(r1) 

# check its coordinate reference system and projection 
r1 
# gives this 
# coord. ref. : +init=epsg:4326 
# .    +proj=longlat 
# .    +datum=WGS84 
# .    +no_defs 
# .    +ellps=WGS84 +towgs84=0,0,0 

## define the stereo polar projection to make it really gorgeous 
newproj <- "+proj=stere 
      +lat_0=90 
      +lat_ts=70 
      +lon_0=-45 +datum=WGS84" 
## try to project 
x <- projectRaster(r1,crs=newproj) 
# and FAIL

「錯誤的，如果（value [1]！= nrow（x）| value [2]！= ncol（x））{：缺少值，其中TRUE/FALSE需要另外：警告消息： dim<-（*tmp*，value = c nr，nc））：通過強制整數範圍引入的NAs「

我一直無法弄清楚「缺失值」。

在爲projectRaster（）的文檔有這樣的：「projectExtent不工作時轉換預計極地數據（例如）經度/緯度非常好，這樣的數據，你可能需要調整返回的對象如做。 ymax（物體）< - 90「但我不知道如何實施他們的建議。

來源

2017-01-04 Dial

這解決了這個問題對我來說：

R1E < - 程度（-180,180,60,90）

R，C < - 作物（R1，R1E）

圖（R，C）

X < - projectRaster（R1C，CRS =取名爲NewProj）

積（X）

來源

2017-01-07 06:18:48

在R中投影氣候研究單位（CRU）CDF柵格爲立體極地視圖：方法？

回答

相關問題