0

「InterruptedGoode'sHomolosineProjection BinaryFormat」のGISデータをRにロードしようとしていますが、その方法がわかりません。

データの例は、次のリンケージから取得できます。ftp: //ftp.glcf.umd.edu/glcf/Continuous_Fields_Tree_Cover/Global/gl-goodes-treecover/gl-goodes-deciduous.bin.gz

Rパッケージ「caTools」の関数「read.ENVI」を使おうとしましたが、動作しません。

> library(caTools)
Loading required package: bitops
> r <- read.ENVI("gl-goodes-deciduous")
Error in read.ENVI("gl-goodes-deciduous") : 
  read.ENVI: Could not open input header file: gl-goodes-deciduous.hdr

助言がありますか？ありがとう〜

4

1 に答える 1

4

600M以上のバイナリファイル全体を窒息させることなく読み取ることができるマシンをお持ちの場合は、すばらしいでしょう。それ以外の場合は、ディスク上でGDALイメージを使用するラスターパッケージのスキルを利用することをお勧めします。

gl-goodes-deciduous.vrtまず、次のようなファイルを作成します。

<VRTDataset rasterXSize="40031" rasterYSize="17347">
<VRTRasterBand dataType="Byte" band="1" subClass="VRTRawRasterBand">
    <SourceFilename relativetoVRT="1">gl-goodes-deciduous.bin</SourceFilename>
    <ImageOffset>0</ImageOffset>
    <PixelOffset>1</PixelOffset>
    <LineOffset>40031</LineOffset>
    <ByteOrder>MSB</ByteOrder>
  </VRTRasterBand>
</VRTDataset>

これにより、GDALライブラリにデータを読み取るのに十分な情報が提供されます。また、VRTのドキュメントを読んで、座標を正しく取得することもできます。次に、Rで、VRTを開き、その低解像度サンプルをプロットします。

> require(raster)
> r = raster("gl-goodes-deciduous.vrt")
> plot(r,maxpixels=100*100)

ProTip：.tifファイルが利用できる場合は、それを使用してください。これは、必要なすべてのメタデータが組み込まれたGeoTIFFである可能性があるためです。

Nrows x Ncolumnsがファイルより163バイト短いのではないかと少し心配していますが、フッターがあるか、これを相殺するためにImageOffsetを163にする必要があります。

于 2013-01-02T09:07:51.937 に答える