ラスタスタックのプロットの緯度と経度が正しくない

Question

TRMM（熱帯雨量）衛星の降水量データをnetCDFファイルとしてダウンロードし、Rでデータをラスタスタックとしてプロットしようとしています。しかし、Rは緯度と経度の軸を間違ってプロットすることを主張しています。そのため、緯度はY軸上にありますが、経度は座標を使用しています。私はplot()とlevelplot()の両方の関数を使ってみましたが、いずれもうまくいかないようです。誰も私がこれを修正するのを助けることができる？あなたが現在の出力を見ることができ、次の画像で

class  : RasterStack 
dimensions : 1440, 186, 267840, 12 (nrow, ncol, ncell, nlayers) 
resolution : 0.25, 0.25 (x, y) 
extent  : -23.25, 23.25, -180, 180 (xmin, xmax, ymin, ymax) 
coord. ref. : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
names  : X2016.01.16, X2016.02.15, X2016.03.16, X2016.04.15, X2016.05.16, X2016.06.15, X2016.07.16, X2016.08.16, X2016.09.15, X2016.10.16, X2016.11.15, X2016.12.16 
Date  : 2016-01-16, 2016-02-15, 2016-03-16, 2016-04-15, 2016-05-16, 2016-06-15, 2016-07-16, 2016-08-16, 2016-09-15, 2016-10-16, 2016-11-15, 2016-12-16 

：

これらは、スタックの特徴です。それは、緯度-23〜23度、経度-180〜180度の熱帯雨量を示すはずです。

Answer 1

座標は、任意の処理の前に切り替えられる場合、それは奇妙です。データをダウンロードしたソースと、おそらくより良いソースがあるかどうかを調べたいと思うかもしれません。

とにかく、（その間に）rasterパッケージはあなたのために助けになることができます。具体的には転位t()機能です。ここでは例です：

# data before transpose 
x <- getData('worldclim',var='tmean',res=10) 
plot(x) 

# data after transpose 
y <- t(x) 
plot(y) 

はあなたのために興味がある可能性がありraster内の他の機能のいくつもあります。flipとrotate

HTH

あなたのお返事ありがとうございます。3210

Answer 2

座標が正しくボックスの外にぶち壊されてしまって、データの新しいセットをダウンロードしようとしたが、同じ問題が発生したことは奇妙なことです。しかし、あなたのおかげで、私はtranspose（）とflip（）関数を使って問題を解決することができました。データを転置してから、イメージが「ミラー化」されたときに、xとyの両方の次元に沿ってそれを反転しなければなりませんでした。

a.t = t(test.rasterstack) a.flipy = flip(a.t, direction = 2) a.t.flipxy = flip(a.t.flipy, direction = 1) levelplot(a.t.flipxy)