Вопрос по dictionary, r, raster – R растровые данные и заданные пределы оси

9

С вашей помощью в другой ветке мне удалось построить несколько глобальных карт. Сначала я конвертирую метеорологические данные GRIB2 в Netcdf, а затем строю глобальные карты.

Теперь я хочу построить только субрегион карты. Я попытался обрезать команду и успешно извлек субрегион глобального файла nc. Но при заговоре я могуне найти, как контролировать пределы оси. Он строит карту больше, чем область данных, поэтому с обеих сторон появляются большие пробелы.

Это сценарий, который яиспользую для построения карт

library("ncdf")
library("raster")
library("maptools")

DIA=format(Sys.time(), "%Y%m%d00") # Data d'avui
url=sprintf("ftp://ftp.ncep.noaa.gov/pub/data/nccf/com/gfs/prod/gfs.%s/gfs.t00z.pgrb2f00", DIA) # Ruta del ftp
loc=file.path(sprintf("%s",url))
download.file(loc,"gfs.grb",mode="wb")

system("/usr/bin/grib2/wgrib2/wgrib2 -s gfs.grb | grep :TMP: | /usr/bin/grib2/wgrib2/wgrib2 -i gfs.grb -netcdf temp.nc",intern=T)

t2m 
Просто добавьтеxlim иylim на вашplot команды, напримерplot(.... , xlim = c(-10,30) , ylim = c(30 , 80 ) ) И хороший сюжет, кстати, +1 Simon O'Hanlon
Вы должны указать свой размер холста png соответственно. Например.png(width=200,height=600, ... ) Robert Hijmans
Может быть, посмотрите на пакет googleVis. Не уверен, что это поможет, но это довольно аккуратно. Он содержит функции IntensityMap, GeoMap и Map, которые могут помочь. Geoffrey Absalom
Привет @ SimonO101 Это была моя первая попытка, прежде чем посмотреть на урожай, не уверен, пробовал ли и то и другое. Не на работе сейчас, попробую. Большое спасибо. pacomet
@ GeoffreyAbsalom Привет и спасибо. Я'У меня есть несколько скриптов, использующих растр, поэтому я предпочитаю не пробовать другой пакет. Но слышал много googleVis, должно быть в моем списке желаний. Еще раз спасибо. pacomet

Ваш Ответ

1   ответ
9

После загрузки файла данных, я могу читать напрямую с растром. Я выбираю группу 221, которая (если я не ошибаюсь) это то, что вам нужно в соответствии сэтот стол:

library("raster")
t2mc  t2mc
class       : RasterLayer 
band        : 221  (of  315  bands)
dimensions  : 361, 720, 259920  (nrow, ncol, ncell)
resolution  : 0.5, 0.5  (x, y)
extent      : -0.25, 359.75, -90.25, 90.25  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +a=6371229 +b=6371229 +no_defs 
data source : /home/oscar/gfs.grb 
names       : gfs 

Ты нет нужна вся степень, поэтому вы используетеcrop чтобы получить желаемую степень:

e 
Я все еще буду искать решение с сюжетом. pacomet
Gracias @ oscar-perpinan Ваше предложение работает нормально. Я'я буду искать дополнительные параметры spplot для оси, заголовка и так далее ... Спасибо pacomet

Похожие вопросы