2016-04-07 11 views
2

RにNetCDFファイルをラスタとして読み込もうとしています。 netCDFファイルには、経度、緯度、深度の関数としての海洋の年間平均気温が記述されています。私は表面の海(netCDFファイルの最初のレベル)に関心があり、t_anはnetCDFファイルの変数の名前です。したがって、私は次のコードを使用しR:netCDFファイルをラスタとして読む

Warning message: 
In .getCRSfromGridMap4(atts) : cannot process these parts of the CRS: epsg_code=EPSG:4326 

あなたが見ることができるように、netCDFファイルがCRS EPSG 4326(またはWGS 84)が、まだラスタ:

MyRast <- raster("Temperature.nc", level = 1, varname = "t_an") 

をこれは私に次の警告を与えますそれが作成され、次のCRSがあります

+proj=longlat +lon_0=0 +a=6378137 +rf=298.257232666016 

私は正しいCRSとのファイルのNetCDFを読み込む方法任意のアイデア?

+0

あなたはNetCDFファイルへのリンクを与えることができますか? –

+0

https://www.nodc.noaa.gov/cgi-bin/OC5/woa13/woa13.plフォーマットとしてnetCDFを選択し、使用可能なグリッドとして1/4インチを選択し、このファイルの名前はt00_04v2.ncとなります。年次見出し。 – Lyngbakr

+0

あなたはnetCDFを使って死んでいるのですか? – MikeJewski

答えて

3

あなたがデータを読み込み後の投影を設定するだけです:

r <- raster("Temperature.nc", varname = "t_an") 
proj4string(r)=CRS("+init=EPSG:4326") 
+1

または 'crs(r)< - " + init = EPSG:4326 "' – RobertH

関連する問題