2016-06-13 58 views
1

我有一個12級的分類柵格圖像,並希望用它從另一個柵格提取信息。我想從分類柵格中的每個類中提取信息。我該怎麼做?使用R中的分類柵格提取信息

library(raster) 
R1 <- raster(nrows=108, ncols=21, xmn=-19, xmx=19,ymn=2, xmx=29) 
###Let's assume the classified raster is D1, I'll like to use the classes of D1 within the boundaries of R1 to extract classified information form R1. 
RD1<-extract(R1, D1) ### where D1 is a raster with classes 

分類圖像可在下面的鏈接: Koppen-Geiger climate classes

+0

如果您想從基於另一個柵格的柵格中提取信息(如您的問題所述),您可以使用'mask'而不是'extract'。 – rar

+0

hi @rar。我所面臨的問題是遮掩正在普遍而不是按類別進行。如果有一種方法可以通過課堂掩飾,它應該可以工作。我該怎麼做,請嗎? –

+0

一種方法是你可以爲每個班級製作面具。例如。 'm1 <-D1','m1 [!m1 == 1] < - NA'將爲1級創建蒙版'm1'。因此,您可以爲每個課程製作蒙版。您還可以爲選定的課程製作蒙版。例如。 'm <-D1','m1[m1> 5] < - NA'將爲1至5級創建掩碼。 – rar

回答

0

@Joke○:聽起來你正在尋找的rasterzonal功能。使用此函數,您可以指定要爲每個類提取哪些信息(例如,sum,min,mean,max等)。

參見raster包文檔:

zonal((x, z, fun='mean', digits=0, na.rm=TRUE, ...)) 

其中x是光柵對象,並且y是與12度的類分類的柵格。

+0

謝謝蒂姆。這工作。 –