ncdf文件的尺寸长度不同

Question

我是R的新手，并且尝试使用R绘制网格数天，但未成功。 我从我所在地区的ecmwf ERA临时2013年获得了降水数据，并且由于ncdf文件的尺寸长度不同，因此无法绘制网格。 如果有人可以帮助我，那将真的很棒。

我在Google云端硬盘中的.nc数据： https ://drive.google.com/file/d/0B0X5ssvd281gUnIyUkZHRGs1a0k/edit?usp=sharing

我使用代码：

> library(ncdf)
> precipitation = open.ncdf("precipitation_2013.nc")
> lonmat  = get.var.ncdf(nc= precipitation,varid="longitude")
> latmat  = get.var.ncdf(nc= precipitation,varid="latitude")
> plot(lonmat, latmat, main= "Precipitation grid 2013")

谢谢你的帮助：）

Answer 1

我从未处理过geo-NetCDF文件，但是当您查看该文件时，您会看到实际情况：

print(precipitation)
[1] "file ~/Downloads/precipitation_2013.nc has 3 dimensions:"
[1] "longitude   Size: 8"
[1] "latitude   Size: 5"
[1] "time   Size: 730"
[1] "------------------------"
[1] "file ~/Downloads/precipitation_2013.nc has 1 variables:"
[1] "short tp[longitude,latitude,time]  Longname:Total precipitation Missval:-32767"

您的实际数据具有longitude ， latitude和time三个维度的short类型的变量。 尺寸值存储为单独的变量，而收集值的变量仅引用尺寸。 因此，该问题与NetCDF或地理信息系统问题无关：它是什么，当您拥有值时实际上要绘制的图

> latmat
[1] 60.0 58.5 57.0 55.5 54.0
> lonmat
[1] 19.5 21.0 22.5 24.0 25.5 27.0 28.5 30.0

您可以通过提取实际变量来构成数据，将其采用长格式，然后将维度值分配​​给结果data.frame ：

library(reshape2)
tp.molten <- melt(get.var.ncdf(nc= precipitation,varid="tp"), varnames=c("longitude", "latitude", "time"))
tp.molten$longitude <- lonmat[tp.molten$longitude]
tp.molten$latitude <- lonmat[tp.molten$latitude ]
tp.molten$time<- get.var.ncdf(nc= precipitation,varid="time")[tp.molten$time ]

close.ncdf( precipitation ) # don't forget

head(tp.molten)
longitude latitude   time       value
1      19.5     19.5 990564 0.002813700
2      21.0     19.5 990564 0.004143370
3      22.5     19.5 990564 0.006473583
4      24.0     19.5 990564 0.008655690
5      25.5     19.5 990564 0.007778240
6      27.0     19.5 990564 0.005223891

现在，您可以例如按一个时间点对这些数据进行子集并将其绘制在lon-lat-g​​rid上

library(ggplot2)
ggplot( data = tp.molten[ tp.molten$time == 992052, ], aes( x = latitude, y = longitude, fill = value ) ) +
  geom_raster()