3

USDA Crop Data Layer (データはこちらから入手できます。大きなファイルに注意してください: http://www.nass.usda.gov/research/Cropland/SARS1a.htm ) をダウンロードし、分析のために R にインポートしたいと考えています。他の大きな GTIFF ファイル (LANDSAT データなど) をインポートしましたが、関数を呼び出すと次のエラーが発生します。

r <- raster(file)

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

ここで、「ファイル」は USDA が提供する .tif ファイルのいずれかです (中西部地域のコンパイルを使用しています)。

問題は、GTiff が真の序数/スペクトル データではなく、カテゴリ データ (土地利用カテゴリ) であることだと思いますが、画像エディター (プレビューなど) で正常にレンダリングされ、R / GDAL でボークする理由がわかりません。 .

助けてくれてありがとう。

編集: これが GIS フォーラムにある場合は、私に知らせてください。一般的に、空間問題に対するプログラミングのアプローチはここにあるようです。

4

2 に答える 2

1

同じファイルでエラー メッセージを再現できました。少なくとも Mac では、GDAL にはホーム ディレクトリよりも深いパス名が必要なようです。つまり、パス名は、USDA トリミング ファイルの例の 1 つとして、~/cdl_tm_r_nv_2010_utm11.tif ではなく、/Users/myname/cdl_tm_r_nv_2010_utm11.tif に戻る必要がありました。でファイルを開くことができました

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

そして、GeoTiff に移動するだけです。

于 2011-06-24T10:12:30.617 に答える
1

うーん、ファイルの名前を「temp.tif」などに変更すると、問題は解決しました。ファイル名に az、0-9、および _ 文字だけが含まれていたため、非常に奇妙です。

GDAL の問題を経験している他の人へ: これで遊ぶことをお勧めします。誰かが本当の根本的な問題を知っているなら、私に知らせてください。元のファイル名は、、cdl_tim_r_ia_2010_utm15.tifなどcdl_tm_r_il_2010_utm16.tifでした。

ありがとう。

于 2011-06-24T16:24:11.897 に答える