2016-04-09 3 views
0

この例では、ラスタとしてプロットすること.binファイル形式での海氷データを使用していますラスターに.binファイルNSIDCの海氷()

「極性の例」の下hereを与えます。私は、National Snow and Ice Data Centerのoriginal ftp serverから入手できる別のファイルを使って同じことを試みています。したがって、問題はないはずです。 。このファイルからRasterLayerオブジェクトを作成できません

:私はプロンプトしようとすると、ただし、R次のスクリプト

# from NSIDC sea ice concentration file 
baseurl <- "ftp://sidads.colorado.edu/pub/DATASETS/" 

f2 <- paste(baseurl, 
     "nsidc0051_gsfc_nasateam_seaice/final-gsfc/north/daily/2013/nt_20130111_f17_v1.1_n.bin", 
     sep='') 

if (!file.exists(basename(f2))) download.file(f2, basename(f2), mode = "wb") 
ice2 <- raster(basename(f2)) 

を= "RasterLayer" .rasterObjectFromFile(X、バンド=バンド、オブジェクトタイプでのエラー、実行私は間違っています

は?破損した.binファイルです?任意の助け感謝!

ありがとう!OK、働くgithubの上の解決策を見つけた

答えて

0

すごくいい。新しい名前の構造はnamly一桁のドットで、元の1若干異なるとしてだけ

hemi <- tolower(substr(bx, 22L, 22L)), 

によってスクリプト

hemi <- tolower(substr(bx, 21L, 21L)) 

の14行目に

https://github.com/cran/raster/blob/master/R/nsidcICE.R

置き換えますv1.1シーケンス!

は比較:

"nt_19781119_f07_v01_s.bin" 

私は

"nt_20130111_f17_v1.1_n.bin" 
に興味があったバージョンに
関連する問題