我的地图（R）倒过来。我如何翻转它？

我有一个.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。这是我做过什么：我的地图（R）倒过来。我如何翻转它？

文件1 < -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'）

时间= ncvar_get（pr_file， '时间'）

雨= ncvar_get（pr_file， '雨'）

image.plot（LON，经纬度，雨[,, 1]，主要= “沉淀”）

错误image.default（...，休息休息=，=添加添加，列= COL）：增加 'X' 和 'Y' 值预计

LAT1 < -sort （lat）

image.plot（LON，LAT1，雨[,, 1]，主要= “沉淀”，zlim = C（0,20））

然后，地图出来，但上下颠倒。即，Antartica位居榜首。我如何正确转动我的地图？

来源

2017-04-17 Kim Heey Jin

请提供[重复的例子（http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example）当你问一个题。在这种情况下，预期的结果也会有很大的帮助。 –

快速入侵是翻转数据本身。

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))

更新： 我，不知何故，切换纬度和经度。请现在试试这个。而且您需要将轴标签更改为有意义的内容。

来源

2017-04-17 04:17:30 kitsune

感谢您的回答。你的地图正是我想要的，但对我而言这并不是太好。所以我做： > LAT = SEQ（1360，length.out = nrow（雨）） > LON = SEQ（1,720，length.out = NcoI位（雨）） > image.plot（LON，LAT ，降水[，ncol（precipitation）：1,1]，main =“降水”）这就是我得到的： image.default错误（...（x）（ - 1）乘以长度（y）（ - 1）我得到的图像在不同位置有lon和lat与你的相比，也没有地图。标题为“降水”的空白轴。怎么了？ –

没关系，我明白了。刚刚混合了经纬度的序列。非常感谢！ –

我的地图（R）倒过来。我如何翻转它？

回答

相关问题