С вашей помощью в другой ветке мне удалось построить несколько глобальных карт. Сначала я конвертирую метеорологические данные 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