Каким образом может быть спроецировано CDF CDR в полярное изображение в R?Проектирование блока исследований климата (CRU) CDF растра в стереополюсный вид в R: How-To?
В пакет cruts и его функция cruts2raster() получает данные временных рядов хорошо, и я могу визуализировать данные за один месяц (или больше), используя растровый пакет и участок() или image() функции. После ввода всех данных и проверки RasterStack он отображает и дает справочную систему координат при запросе, но я не могу получить ее на полярную проекцию.
Я использую rgdal (установив последнюю proj.4 в моей Мака/USR/местные/Библиотека/pkg-config), но, кажется, не проецировать в этом полярном зрения:
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
«Ошибка, если (значение [1]! = nrow (x) | значение [2]! = ncol (x)) {: отсутствует значение, в котором требуется TRUE/FALSE Кроме того: Предупреждающее сообщение: В dim<- (*tmp*, value = c (nr, nc)): NA, вводимые путем принуждения в целочисленный диапазон "
Мне не удалось выяснить «отсутствующие значения».
В документации projectRaster() есть это: «projectExtent не очень хорошо работает при преобразовании проецируется циркумполярные данные в (например) долготы/широты С такими данными может потребоваться настроить возвращаемый объект Eg сделать.. ymax (объект) < - 90 ", но я не знаю, как реализовать их советы.