南卓铜 ([email protected])

一些有名的Land cover数据集并没有使用通用的GIS格式，比如在以下链接，

http://www.geog.umd.edu/landcover/tree-cover-poster/download.html

下载broadleaf Layer，gzip解压缩后为 694,417,920 bytes，扩展名为.img，但不同于Erdas Imagine .img格式。

据介绍Goode binary file是plain binary，即相当于BSQ （Band sequential），其metadata有以下信息：

40031 pixels by 17347 lines，且是8bit unsigned integer，应当总大小为 694,417,757bytes，而解压缩后的大小大了163bytes，目前我还不知道原因。

根据metadata信息，可以添加头文件，broadleaf.hdr，内容如下，

nrows      17347
ncols     40031
nbands     1
nbits     8
layout     bsq
skipbytes     0
ulxmap     -20015500.000
ulymap     8673500.000
xdim     1000
ydim     1000

将broadleaf.img改名为 broadleaf.bsq，在ArcGIS中才可以正常访问，也可以导出为其它格式，比如Grid。不过注意其空间参考系仍为unknown。在ArcCatalog中打开其属性页，在spatial reference中选projected coordinates/World/Goode Homolosine (Land).prj。可以将bsq设置为 Goode投影。

如果希望转换成其它常用格式，

上图将bsq转成ESRI binary grid格式。图形化结果如下图。

在Workstation下好像没有找到Goode投影，但在Desktop下可以找到.prj。因为如果bsq是unknown时，导成esri grid，在ArcCatalog里对grid进行定义投影（属性->define spatial reference），没法子定义到Goode投影，很奇怪。所以这里，我是先对bsq定义投影，再转成grid。此外，通过ArcToolbox的Data Management下的projection define来做，也是可以。