* .adfファイルをRに読み込む方法は?


9

.adfファイルをRにロードしたい。データは次のページからのものである:http ://www.fao.org/geonetwork/srv/en/metadata.show?id=14057

インターネットでの調査の結果、次のコードを見つけました。問題は、RasterLayerクラスで、そこにあるべきではないネガティブな値を取得することです。なぜこれが起こるのかわかりませんので、誰かが私を助けてくれるといいのですが!?

コード:

library(rgdal)
library(RColorBrewer)
dpath<- path...

x <- new("GDALReadOnlyDataset", dpath)
getDriver(x)
getDriverLongName(getDriver(x))
xx<-asSGDF_GROD(x)
r <- raster(xx)

「r」の出力は次のとおりです。

rクラス:RasterLayer寸法:2160、4320、9331200(nrow、ncol、ncell)解像度:0.08333333、0.08333333(x、y)範囲:-180、180、-90、90(xmin、xmax、ymin、ymax)座標。ref。:+ proj = longlat + ellps = WGS84 + towgs84 = 0,0,0,0,0,0,0 + no_defsデータソース:メモリ名内:band1値:-997、16(最小、最大)

値の「16」は、成長期間の長さの16クラスを指します。しかし、それらの「-997」はどこから来たのでしょうか。たぶん、コーディネートに問題があります。ref?

ここにも「xx」のデータの概要があります。

データの概要:最小。第1四半期 中央値第3四半期 マックス。NAの-997 3 5 -9 8 16 7123158

そして、xxのデータをさらに詳しく見てみると、

テーブル(xx $ band1)

-997 1 2 3 4 5 6 7 8 9 10 11 12 31711 429643 83011 166674 207228 270161 240958 183342 118608 98795 88473 73743 56022 13 14 15 16 30104 45521 52216 31832

実際には、まさにこの '-997'のものが含まれています。NAはozeansだと思うので、データの読み込みに問題があるのでしょうか、それともデータが理解できませんか?


1
タイトルをインテリジェントで正確なものに変更してください。
mdsumner 2015年

fwiw、rgdalを直接使用する必要はなく、raster(dpath)を使用するだけです
mdsumner '29 / 07/29

回答:


3

あなたはほぼ正しいです:

海洋の場合、NODATAは-32768に設定されます。さらに、海岸線から除外されていない大湖には-997が設定されています。

ピクセルコンテンツ(成長期間)は湖では意味がないため、-997もNODATAとして安全に扱うことができます。

弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.