我的地图(在 R 中)颠倒了。我如何将其翻转回来?
My map (in R) is upside down. How do I flip it back up?
我有一个 .nc 文件,我正试图从中绘制地图。我在这里得到了我的数据:https://www.esrl.noaa.gov/psd/data/gridded/data.UDel_AirT_Precip.html (precip.mon.ltm.v301.nc)。纬度从89.75
开始到-89.75
。经度从 0.25
到 359.75
。这就是我所做的:
File1<-Sys.glob("precip.mon.ltm*.nc")
tmfile<-File1
pr_file<-nc_open(tmfile)
lat=ncvar_get(pr_file,'lat')
lon=ncvar_get(pr_file,'lon')
time=ncvar_get(pr_file,'time')
precip=ncvar_get(pr_file,'precip')
image.plot(lon,lat,precip[,,1],main="Precipitation")
错误 image.default(..., breaks = breaks, add = add, col = col) :
增加 'x' 和 'y' 预期值
lat1<-sort(lat)
image.plot(lon,lat1,precip[,,1],main="Precipitation",zlim=c(0,20))
然后地图出来了,但是颠倒了。即,南极洲在最上面。如何正确转动我的地图?
一个快速的破解方法是翻转数据本身。
image.plot(lon,lat,precip[,ncol(precip):1,1],main="Precipitation")
我认为 lon 和 lat 变量也需要根据错误提示进行不同的定义,例如
(参见 the related question):
lon = seq(1, 360, length.out = nrow(precip))
lat = seq(1, 360, length.out = ncol(precip))
更新:
不知何故,我切换了纬度和经度。请,现在试试这个。您需要将轴标签更改为有意义的内容。
我有一个 .nc 文件,我正试图从中绘制地图。我在这里得到了我的数据:https://www.esrl.noaa.gov/psd/data/gridded/data.UDel_AirT_Precip.html (precip.mon.ltm.v301.nc)。纬度从89.75
开始到-89.75
。经度从 0.25
到 359.75
。这就是我所做的:
File1<-Sys.glob("precip.mon.ltm*.nc")
tmfile<-File1
pr_file<-nc_open(tmfile)
lat=ncvar_get(pr_file,'lat')
lon=ncvar_get(pr_file,'lon')
time=ncvar_get(pr_file,'time')
precip=ncvar_get(pr_file,'precip')
image.plot(lon,lat,precip[,,1],main="Precipitation")
错误 image.default(..., breaks = breaks, add = add, col = col) : 增加 'x' 和 'y' 预期值
lat1<-sort(lat)
image.plot(lon,lat1,precip[,,1],main="Precipitation",zlim=c(0,20))
然后地图出来了,但是颠倒了。即,南极洲在最上面。如何正确转动我的地图?
一个快速的破解方法是翻转数据本身。
image.plot(lon,lat,precip[,ncol(precip):1,1],main="Precipitation")
我认为 lon 和 lat 变量也需要根据错误提示进行不同的定义,例如 (参见 the related question):
lon = seq(1, 360, length.out = nrow(precip))
lat = seq(1, 360, length.out = ncol(precip))
更新: 不知何故,我切换了纬度和经度。请,现在试试这个。您需要将轴标签更改为有意义的内容。