1

次のGISデータファイルをRに読み込もうとしています。

いくつかのデータ仕様はここにあります：

画像タイプ：汎用フラットバイナリ、行によるバイトインターリーブ（BIL）

投影：中断されたグード図法

Rパッケージ「ラスター」を使おうとしましたが失敗しました。

library(raster)
r <- raster(file.choose())

Error in .local(.Object, ...) : 
  `C:\global_forest_cover.img' not recognised as a supported file format.

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

では、このGISデータをRにどのようにロードする必要がありますか？また、投影法「Interrupted GoodeHomolosine」をLongLatに変換しますか？

4

1 に答える 1

5

というファイルを作成し、httpglobal_forest_cover.hdr ：//edc2.usgs.gov/glcc/fao/header_file.phpのヘッダー情報をそのファイルに貼り付けます。

BYTEORDER   M
LAYOUT  BIL
NROWS   15059
NCOLS   36543
NBANDS  1
NBITS   8
BANDROWBYTES    36543
TOTALROWBYTES       36543
ULXMAP  -17359000
ULYMAP  8673000
XDIM    1000
YDIM    1000

次に、ファイルを読み取り.imgます。

> forest = raster("global_forest_cover.img")
Geospatial Data Abstraction Library extensions to R successfully loaded
Loaded GDAL runtime: GDAL 1.6.3, released 2009/11/19
Path to GDAL shared files: /usr/share/gdal16
Loaded PROJ.4 runtime: Rel. 4.7.1, 23 September 2009, [PJ_VERSION: 470]
Path to PROJ.4 shared files: (autodetected)
> forest
class       : RasterLayer 
dimensions  : 15059, 36543, 550301037  (nrow, ncol, ncell)
resolution  : 1000, 1000  (x, y)
extent      : -17359500, 19183500, -6385500, 8673500  (xmin, xmax, ymin, ymax)
coord. ref. : NA 
data source : /home/nobackup/rowlings/Downloads/global_forest_cover.img 
names       : global_forest_cover 
values      : 0, 255  (min, max)

GDALは.hdrファイルを確認し、それを使用してファイル内の15059*36543バイトの構造を計算し.imgます。

ただし、投影情報はありません。ただし、それが行われたとしても、IGHプロジェクションを緯度経度にワープすることは多少問題になる可能性があります。他のIGH投影マップを見ると、地球を平らにするために海を破壊することが含まれていることがわかります。緯度経度に戻るには、これらすべての平坦化と涙を元に戻す必要があります。それは可能であるはずであり、最新のPROJ4ライブラリはIGHプロジェクションをサポートしていると思いますが、おそらく一方向にしかサポートされていません。ただし、データの作成に使用したのとまったく同じコードを使用しない限り、逆に正しい答えが得られない可能性があります。

ファイルがzipに含まれていなかったこと、およびプロジェクションがさらに分析するのに実際には適切ではないことを考えると、.hdr別のデータソースを探しに行きます。このデータは明らかに表示目的のみを目的としています。また、IGHの投影によって地球が分割された場所でデータ値が欠落しているようには見えないため、これも貧弱です。http：//edc2.usgs.gov/glcc/fao/index.phpページの写真を他のIGHと比較してください。予測http://en.wikipedia.org/wiki/Goode_homolosine_projection

ここのデータはもっと適切かもしれないと思いました：http：//www.fao.org/forestry/32203/en/ 

> r=raster("./fceurope/europe/w001001.adf")
> r
class       : RasterLayer 
dimensions  : 20000, 40000, 8e+08  (nrow, ncol, ncell)
resolution  : 0.009, 0.009  (x, y)
extent      : -180, 180, -90, 90  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +ellps=clrk66 +no_defs 
data source : /home/nobackup/rowlings/Downloads/fceurope/europe/w001001.adf 
names       : w001001 
values      : 1, 6  (min, max)
Raster Attribute Table
     fields : ID   COUNT
        min :  1  129626
        max :  6 4963775

しかし、私は名前の「ヨーロッパ」と世界的な広がりの間で少し混乱しています。また、非常に大きなラスターです。ああ、それはどこの真ん中でヨーロッパを手に入れました。投影が緯度経度になっていることに注意してください。グリッドをまとめる前にグリッドを細くして、そこに指定されているすべてのファイルからグローバルラスターを作成することをお勧めします。

于 2013-01-14T09:42:12.530 に答える