Форумы GIS-Lab.info

Код: Выделить всё

library("raster")
#генерируем "подопытный" растр r
r <- raster(nrow=seq(from =  -90,  to =  90,  by = 1),
                ncol=seq(from = -180,   to = 180, by = 2))

#0 или 1 присваиваем значениям ячееек "случайным образом": либо "0", либо "1"
values(r) <- sample(c(0,1), replace=TRUE, size=ncell(r))

#"подопытный" растр r готов
str(r)

#значения ячеек с "0" "для удобоваримости" в NA "обращаем"
values(r)[values(r)==0] <- NA
str(r)

#получаем количество ячеек со значением "1"
countCells1 <- cellStats(r, stat='sum', na.rm=TRUE)
str(countCells1)

#получаем растр a с площадями ячеек растра r (NA-значения - они же "0" - " игнорируем")
a <- area(r, na.rm = TRUE)
str(a)

#вот так, получаем собственно, общую площадь ячеек растра r со значением "1"
AreaCellsWith1 <- cellStats(a, stat='sum', na.rm=TRUE)
str(AreaCellsWith1)

Код: Выделить всё

library(dplyr)
df <- as.data.frame(freq(r))
df <- mutate(df, area=count*30/1e4) # площадь в га, если пространственное разрешение растра 30 м