2

最近、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)
4

1 に答える 1