.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だと思うので、データの読み込みに問題があるのでしょうか、それともデータが理解できませんか?