用R中的级别打印netcdf文件
我最近开始在R中使用netcdf。示例数据如下: >175种作物的收获面积和产量>单个作物>大豆产量2000 在此文件夹中,有一个名为用R中的级别打印netcdf文件,r,raster,netcdf,netcdf4,R,Raster,Netcdf,Netcdf4,我最近开始在R中使用netcdf。示例数据如下: >175种作物的收获面积和产量>单个作物>大豆产量2000 在此文件夹中,有一个名为soybean\u AreaYieldProduction.nc 这就是我打开netcdf的方式 library(ncdf4) dat <- nc_open("soybean_AreaYieldProduction.nc") print(soy) 1 variables (excluding dimension variables): fl
soybean\u AreaYieldProduction.nc
这就是我打开netcdf的方式
library(ncdf4)
dat <- nc_open("soybean_AreaYieldProduction.nc")
print(soy)
1 variables (excluding dimension variables):
float soybeanData[longitude,latitude,level,time]
LayerDescriptions: struct(5).Data(:,:,1/2/3/4/5/6) to access data layer: 1=Harvested Area fraction, 2=Yield 3=Harvested Area data quality, 4=Yield data quality, 5=Harvested Area in hectares, 6= Production
Units: Harvested Area Fraction(1)=percent of gridcell that was harvested, Yield(2)=metric tons per hectare, Harvested Area Hectares(5)=total hectares harvested per gridcell, Production(6)=Metric Tons
DataQuality: In levels 3 and 4, a value of 1 = county; .75 = state; .5 = interpolated from within 2 degrees lat/long; .25 = country; 0 = missing.
4 dimensions:
longitude Size:4320
units: longitude
latitude Size:2160
units: latitude
level Size:6
time Size:1
使用
光栅
软件包可能是最简单的:
library(raster)
r1 <- raster("soybean_AreaYieldProduction.nc", level=1)
r2 <- raster("soybean_AreaYieldProduction.nc", level=2)
plot(r1)
image(r1)
s <- stack(r1, r2)
plot(s)
查看CRAN上的其他映射包是否有自动堆叠ncdf文件中的级别的功能?
level.1 <- ncvar_get(dat, ????)
image(lon,lat, level)
library(raster)
r1 <- raster("soybean_AreaYieldProduction.nc", level=1)
r2 <- raster("soybean_AreaYieldProduction.nc", level=2)
plot(r1)
image(r1)
s <- stack(r1, r2)
plot(s)
spplot(s)
library(rasterVis)
levelplot(r1)
levelplot(s)