0

NetCDFファイルからバイナリファイルを書きたい

library(ncdf)
download.file("http://gswp/Fixed/SoilDepth.nc", destfile="SoilDepth.nc")
soil <- open.ncdf("SoilDepth.nc")
soil$var[[3]] -> var3 
get.var.ncdf(soil, var3) -> SoilDepth
download.file("http://gswp/Fixed/landmask_gswp.nc", destfile="landmask.nc")
landmask <- open.ncdf("landmask.nc")
landmask$var[[3]] -> varland
get.var.ncdf(landmask, varland) -> land
land = t(land)
land[land==1] <- SoilDepth
land[land==0] <- NA
land = t(land)
image(land)

このコードの結果は次のようになります：！[土壌図] [1]

今、私はそれをバイナリファイルに書きたいです：

結果は逆さまの画像になります。

4

1 に答える 1

1

を指定して、自分で逆向きをリクエストしましたylim=c(1,0)。その最後の行を次のように変更するだけです

image(y)

二重転置はリソースを食べる以外に何もしないので、大丈夫です。

詳細情報：元のデータと再読み込みされたデータの唯一の違いは、前者が（利用不可）NaNであった場所に後者が（つまり番号ではない）あるという事実です。NAこの変更を元に戻すと、完全に同一のデータが生成されます。

y[is.nan(y)] <- NA

この後、yと見分けがつかなくなりlandます。

于 2012-09-07T11:50:56.873 に答える