最近、R で netcdf を使い始めました。サンプル データは次のとおりです。
http://www.earthstat.org/data-download/ > 175 作物の収穫面積と収量 > 個々の作物 > 大豆_HarvAreaYield2000_NetCDF
このフォルダには、と呼ばれる netcdf ファイルがあります。soybean_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
各レベルをプロットしたいのですが、各レベルのデータを抽出する方法がわかりません。
これは、経度と緯度のデータを抽出する方法です。
lon <- ncvar_get(dat,"longitude") # extract long
lat <- ncvar_get(dat,"latitude") # extract lat
しかし、個々のレベルを抽出するにはどうすればよいでしょうか?
level.1 <- ncvar_get(dat, ????)
最終的な目的は、次のコマンドを使用して、視覚化する各レベルを視覚化することです。
image(lon,lat, level)