柵格堆棧正在繪製經緯度座標

我已經從TRMM（熱帶降雨）衛星下載了降水數據作爲netCDF文件，並且一直試圖在R中繪製數據作爲rasterstack。然而，R堅持不正確地繪製緯度和經度座標軸，使得經度被繪製在x軸上（因爲它應該是），但是使用緯度座標，而緯度在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度經度的熱帶降雨量。

來源

2017-06-12 cbur

如果座標到任何處理之前切換很奇怪。也許你想評估下載數據的來源，如果還有更好的數據。

反正，（在此期間）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

來源

2017-06-13 07:29:11 Val

感謝您的回覆。看起來奇怪的是，座標出現了問題，我嘗試下載一組新的數據，併發生同樣的問題。但是，由於您的輸入，我可以通過使用transpose（）和flip（）函數來糾正問題。我不得不調換數據，然後沿着x和y方向翻轉它，因爲圖像是「鏡像」的。這裏是我的情況下，其他人使用的代碼遇到此問題與TRMM數據集：

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)

來源

2017-06-13 10:31:33 cbur

這些數據集旋轉「設計」。另見這裏：https://gis.stackexchange.com/questions/170137/converting-binary-trmm-to-geotiff-using-r/170173#170173 – lbusett

謝謝@LoBu，這是超級有用的，我想我會用未來的Rsenal軟件包 – cbur

柵格堆棧正在繪製經緯度座標

回答

相關問題