NSIDC 海冰 .bin 文件到光栅（）

Question

我正在尝试按照 "A polar example"

下的 here 示例

此示例使用 .bin 格式的海冰数据绘制为栅格。我正在尝试使用国家冰雪数据中心 original ftp server 提供的不同文件进行同样的操作。因此，我认为应该没有问题。但是，当我尝试提示 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))

Error in .rasterObjectFromFile(x, band = band, objecttype = "RasterLayer", : Cannot create a RasterLayer object from this file.

我哪里错了？ .bin 文件是否已损坏？任何帮助表示赞赏！

谢谢！

Answer 1

好的，在 github 上找到了一个非常有效的解决方案。

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

只需替换脚本的第 14 行

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

来自

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

由于新的名称结构与原来的名称结构略有不同，在 v1.1 序列中只有一个点！

比较：

"nt_19781119_f07_v01_s.bin"

到我感兴趣的版本

"nt_20130111_f17_v1.1_n.bin"

NSIDC 海冰 .bin 文件到光栅（）

NSIDC sea ice .bin file to raster()

gis

geospatial

bin

r-raster